Comparative analysis of mitochondrial genomes of Rhizophagus irregularis - syn. Glomus irregulare - reveals a polymorphism induced by variability generating elements

Characterization of the complete mitochondrial genome of Mucor indicus Lendn. 1930 (Mucorales: Mucoraceae), isolated from the wine fermentation system

Mucor indicus Lendn. 1930 has been widely used in food fermentation; however, its mitochondrial genome characteristics are not well understood. In this study, the complete mitochondrial genome of M. indicus was obtained, which was 61,400 bp in length with a GC content of 33%. The M. indicus mitochondrial genome was found to contain 14 core protein-coding genes, four free-standing ORFs, 18 intronic ORFs, 26 tRNAs, and two rRNA genes. Phylogenetic trees were generated for 25 early-differentiated fungi using the Bayesian inference (BI) method, which demonstrated that M. indicus is closely related to Mucor piriformis. This study provides useful information for the classification and evolution of Mucor species or other early-differentiated fungi.

Mitochondrial genomic characteristics and phylogenetic analysis of a brewing fungus, Rhizopus microsporus Tiegh. 1875 (Mucorales: Rhizopodaceae)

Rhizopus microsporus Tiegh. 1875 is widely used in a variety of industries, such as brewing, wine making, baking, and medicine production, as it has the capability to break down proteins and generate surface-active agents. To date, the mitochondrial genome features of early evolved fungi from the Rhizopus genus have not been extensively studied. Our research obtained a full mitochondrial genome of R. microsporus species, which was 43,837 bp in size and had a GC content of 24.93%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 7 intronic ORFs, 24 tRNAs, and 2 rRNA genes. Through the use of the BI phylogenetic inference method, we were able to create phylogenetic trees for 25 early differentiation fungi which strongly supported the major clades; this indicated that R. microsporus is most closely related to Rhizopus oryzae.

Diversity of Arbuscular Mycorrhizal Fungi in Distinct Ecosystems of the North Caucasus, a Temperate Biodiversity Hotspot

BACKGROUND

Investigations that are focused on arbuscular mycorrhizal fungus (AMF) biodiversity is still limited. The analysis of the AMF taxa in the North Caucasus, a temperate biodiversity hotspot, used to be limited to the genus level. This study aimed to define the AMF biodiversity at the species level in the North Caucasus biotopes.

METHODS

The molecular genetic identification of fungi was carried out with ITS1 and ITS2 regions as barcodes via sequencing using Illumina MiSeq, the analysis of phylogenetic trees for individual genera, and searches for operational taxonomic units (OTUs) with identification at the species level. Sequences from MaarjAM and NCBI GenBank were used as references.

RESULTS

We analyzed >10 million reads in soil samples for three biotopes to estimate fungal biodiversity. Briefly, 50 AMF species belonging to 20 genera were registered. The total number of the AM fungus OTUs for the "Subalpine Meadow" biotope was 171/131, that for "Forest" was 117/60, and that for "River Valley" was 296/221 based on ITS1/ITS2 data. The total number of the AM fungus species (except for virtual taxa) for the "Subalpine Meadow" biotope was 24/19, that for "Forest" was 22/13, and that for "River Valley" was 28/24 based on ITS1/ITS2 data. Greater AMF diversity, as well as number of OTUs and species, in comparison with that of forest biotopes, characterized valley biotopes (disturbed ecosystems; grasslands). The correlation coefficient between "Percentage of annual plants" and "Glomeromycota total reads" r = 0.76 and 0.81 for ITS1 and ITS2, respectively, and the correlation coefficient between "Percentage of annual plants" and "OTUs number (for total species)" was r = 0.67 and 0.77 for ITS1 and ITS2, respectively.

CONCLUSION

High AMF biodiversity for the river valley can be associated with a higher percentage of annual plants in these biotopes and the active development of restorative successional processes.

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex

Introduction

Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species.

Methods

In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers.

Results and Discussion

Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

Functional annotation and comparative analysis of four Botrytis cinerea mitogenomes reported from Punjab, Pakistan

Botrytis cinerea is one of the top phytopathogenic fungus which ubiquitously cause grey mold on a variety of horticultural plants. The mechanism of respiration in the fungus occurs within the mitochondria. Mitogenomes serve as a key molecular marker for the investigation of fungal evolutionary patterns. This study aimed at the complete assembly, characterization, and comparative relationship of four mitogenomes of Botrytis cinerea strains including Kst5C, Kst14A, Kst32B, Kst33A, respectively. High throughput sequencing of four mitogenomes allowed the full assembly and annotation of these sequences. The total genome length of these 4 isolates Kst5C Kst14A, Kst32B, Kst33A was 69,986 bp, 77,303 bp, 76,204 bp and 55, 226 bp respectively. The distribution of features represented 2 ribosomal RNA genes,14 respiration encoding proteins, 1 mitochondrial ribosomal protein-encoding gene, along with varying numbers of transfer RNA genes, protein-coding genes, mobile intronic regions and homing endonuclease genes including LAGLIDADG and GIY-YIG domains were found in all four mitogenomes. The comparative analyses performed also decipher significant results for four mitogenomes among fungal isolates included in the study. This is the first report on the detailed annotation of mitogenomes as a proof for investigation of variation patterns present with in the B. cinerea causing grey mold on strawberries in Pakistan. This study will also contribute to the rapid evolutionary analysis and population patterns present among Botrytis cinerea.

The Potential Applications of Commercial Arbuscular Mycorrhizal Fungal Inoculants and Their Ecological Consequences

Arbuscular mycorrhizal fungal (AMF) inoculants are sustainable biological materials that can provide several benefits to plants, especially in disturbed agroecosystems and in the context of phytomanagement interventions. However, it is difficult to predict the effectiveness of AMF inoculants and their impacts on indigenous AMF communities under field conditions. In this review, we examined the literature on the possible outcomes following the introduction of AMF-based inoculants in the field, including their establishment in soil and plant roots, persistence, and effects on the indigenous AMF community. Most studies indicate that introduced AMF can persist in the target field from a few months to several years but with declining abundance (60%) or complete exclusion (30%). Further analysis shows that AMF inoculation exerts both positive and negative impacts on native AMF species, including suppression (33%), stimulation (38%), exclusion (19%), and neutral impacts (10% of examined cases). The factors influencing the ecological fates of AMF inoculants, such as the inherent properties of the inoculum, dosage and frequency of inoculation, and soil physical and biological factors, are further discussed. While it is important to monitor the success and downstream impacts of commercial inoculants in the field, the sampling method and the molecular tools employed to resolve and quantify AMF taxa need to be improved and standardized to eliminate bias towards certain AMF strains and reduce discrepancies among studies. Lastly, inoculant producers must focus on selecting strains with a higher chance of success in the field, and having little or negligible downstream impacts.

Arbuscular Mycorrhizal Fungi and the Need for a Meaningful Regulatory Plant Protection Product Testing Strategy

Arbuscular mycorrhizal fungi (AMF) perform key soil ecosystem services and, because of their symbiotic relationship with plant roots, may be exposed to the plant protection products (PPPs) applied to soils and crops. In 2017, the European Food Safety Authority (EFSA) released a scientific opinion addressing the state of the science on risk assessment of PPPs for in-soil organisms, recommending the inclusion of AMF ecotoxicological testing in the PPP regulatory process. However, it is not clear how this can be implemented in a tiered, robust, and ecologically relevant manner. Through a critical review of current literature, we examine the recommendations made within the EFSA report and the methodologies available to integrate AMF into the PPP risk assessment and provide perspective and commentary on their agronomic and ecological relevance. We conclude that considerable research questions remain to be addressed prior to the inclusion of AMF into the in-soil organism risk assessment, many of which stem from the unique challenges associated with including an obligate symbiont within the PPP risk assessment. Finally, we highlight critical knowledge gaps and the further research required to enable development of relevant, reliable, and robust scientific tests alongside pragmatic and scientifically sound guidance to ensure that any future risk-assessment paradigm is adequately protective of the ecosystem services it aims to preserve. Environ Toxicol Chem 2022;41:1808-1823. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Comparative analysis of Malassezia furfur mitogenomes and the development of a mitochondria-based typing approach

Malassezia furfur is a yeast species belonging to Malasseziomycetes, Ustilaginomycotina and Basidiomycota that is found on healthy warm-blooded animal skin, but also involved in various skin disorders like seborrheic dermatitis/dandruff and pityriasis versicolor. Moreover, Malassezia are associated with bloodstream infections, Crohn's disease and pancreatic carcinoma. Recent advances in Malassezia genomics and genetics have focused on the nuclear genome. In this work, we present the M. furfur mitochondrial (mt) genetic heterogenicity with full analysis of 14 novel and six available M. furfur mt genomes. The mitogenome analysis reveals a mt gene content typical for fungi, including identification of variable mt regions suitable for intra-species discrimination. Three of them, namely the trnK–atp6 and cox3–nad3 intergenic regions and intron 2 of the cob gene, were selected for primer design to identify strain differences. Malassezia furfur strains belonging to known genetic variable clusters, based on AFLP and nuclear loci, were assessed for their mt variation using PCR amplification and sequencing. The results suggest that these mt regions are excellent molecular markers for the typing of M. furfur strains and may provide added value to nuclear regions when assessing evolutionary relationships at the intraspecies level.

The Coevolution of Fungal Mitochondrial Introns and Their Homing Endonucleases (GIY-YIG and LAGLIDADG)

Fungal mitochondrial (mt) genomes exhibit great diversity in size which is partially attributed to their variable intergenic regions and most importantly to the inclusion of introns within their genes. These introns belong to group I or II, and both of them are self-splicing. The majority of them carry genes encoding homing endonucleases, either LAGLIDADG or GIY-YIG. In this study, it was found that these intronic homing endonucleases genes (HEGs) may originate from mt free-standing open reading frames which can be found nowadays in species belonging to Early Diverging Fungi as “living fossils.” A total of 487 introns carrying HEGs which were located in the publicly available mt genomes of representative species belonging to orders from all fungal phyla was analyzed. Their distribution in the mt genes, their insertion target sequence, and the phylogenetic analyses of the HEGs showed that these introns along with their HEGs form a composite structure in which both selfish elements coevolved. The invasion of the ancestral free-standing HEGs in the introns occurred through a perpetual mechanism, called in this study as “aenaon” hypothesis. It is based on recombination, transpositions, and horizontal gene transfer events throughout evolution. HEGs phylogenetically clustered primarily according to their intron hosts and secondarily to the mt genes carrying the introns and their HEGs. The evolutionary models created revealed an “intron-early” evolution which was enriched by “intron-late” events through many different independent recombinational events which resulted from both vertical and horizontal gene transfers.

Comparison of mitochondrial genomes provides insights into intron dynamics and evolution in Botryosphaeria dothidea and B. kuwatsukai

Botryosphaeria dothidea is one of the most common fungal pathogens on a large number of hosts worldwide. Botryosphaeria dothidea and B. kuwatsukai are also the main causal agents of apple ring rot. In this study, we sequenced, assembled and annotated the circular mitogenomes of 12 diverse B. dothidea isolates (105.7-114.8 kb) infecting various plants including apple, and five diverse B. kuwatsukai isolates (118.0-124.6 kb) from apple. B. dothidea mitogenomes harboured a set of 29-31 introns and 48-52 ORFs. In contrast, B. kuwatsukai mitogenomes harboured more introns (32-34) and ORFs (51-54). The variation in mitogenome sizes was associated mainly with different numbers of introns and insertions of mobile genetic elements. Interestingly, B. dothidea and B. kuwatsukai displayed distinct intron distribution patterns, with three intron loci showing presence/absence dynamics in each species. Large numbers of introns (57% in B. dothidea and 49% in B. kuwatsukai) were most likely obtained through horizontal transfer from non-Dothideomycetes. The mitochondrial gene phylogeny supported the differentiation of the two species. Overall, this study sheds light into the mitochondrial evolution of the plant pathogens B. dothidea and B. kuwatsukai, and intron distribution patterns could be useful markers for studies on population diversity.

Commercial arbuscular mycorrhizal fungal inoculant failed to establish in a vineyard despite priority advantage

Background

Arbuscular mycorrhizal (AM) fungi associate with most plants and can increase nutrient uptake. As a result, commercial inoculants called “biofertilizers” containing AM fungi have been developed and marketed to increase plant performance. However, successful establishment of these inoculants remains a challenge, and may be negatively impacted by competition with fungi already present (priority effects). Perennial agriculture may be more amenable if inoculants can be successfully established on crops prior to field planting.

Methods

Here, we inoculate grapevine (Vitis vinifera) with a commercial inoculant in three treatments designed to manipulate the strength and direction of priority effects and quantified the abundance of the fungal strain before and after introduction using droplet digital PCR (ddPCR).

Results

We found that the introduced strain did not establish in any treatment, even with priority advantage, and inoculated vines did not differ in performance from non-inoculated vines. Fungal abundance was not greater than in pre-inoculation soil samples during any of the five years sampled and may have been impaired by high available phosphorus levels in the soil. This study highlights the need to understand and evaluate how the management of the agricultural system will affect establishment before introduction of “biofertilizers”, which is often unpredictable.

Arbuscular Mycorrhizal Symbiosis Enhances Photosynthesis in the Medicinal Herb Salvia fruticosa by Improving Photosystem II Photochemistry

We investigated the influence of Salvia fruticosa colonization by the arbuscular mycorrhizal fungi (AMF) Rhizophagus irregularis on photosynthetic function by using chlorophyll fluorescence imaging analysis to evaluate the light energy use in photosystem II (PSII) of inoculated and non-inoculated plants. We observed that inoculated plants used significantly higher absorbed energy in photochemistry (Φ_PSII) than non-inoculated and exhibited significant lower excess excitation energy (EXC). However, the increased Φ_PSII in inoculated plants did not result in a reduced non-regulated energy loss in PSII (Φ_NO), suggesting the same singlet oxygen (¹O₂) formation between inoculated and non-inoculated plants. The increased Φ_PSII in inoculated plants was due to an increased efficiency of open PSII centers to utilize the absorbed light (Fv'/Fm') due to a decreased non-photochemical quenching (NPQ) since there was no difference in the fraction of open reaction centers (q_p). The decreased NPQ in inoculated plants resulted in an increased electron-transport rate (ETR) compared to non-inoculated. Yet, inoculated plants exhibited a higher efficiency of the water-splitting complex on the donor side of PSII as revealed by the increased Fv/Fo ratio. A spatial heterogeneity between the leaf tip and the leaf base for the parameters Φ_PSII and Φ_NPQ was observed in both inoculated and non-inoculated plants, reflecting different developmental zones. Overall, our findings suggest that the increased ETR of inoculated S. fruticosa contributes to increased photosynthetic performance, providing growth advantages to inoculated plants by increasing their aboveground biomass, mainly by increasing leaf biomass.

Fungal Mitogenomes: Relevant Features to Planning Plant Disease Management

Mitochondrial genomes (mt-genomes) are characterized by a distinct codon usage and their autonomous replication. Mt-genomes encode highly conserved genes (mt-genes), like proteins involved in electron transport and oxidative phosphorylation but they also carry highly variable regions that are in part responsible for their high plasticity. The degree of conservation of their genes is such that they allow the establishment of phylogenetic relationships even across distantly related species. Here, we describe the mechanisms that generate changes along mt-genomes, which play key roles at enlarging the ability of fungi to adapt to changing environments. Within mt-genomes of fungal pathogens, there are dispensable as well as indispensable genes for survival, virulence and/or pathogenicity. We also describe the different complexes or mechanisms targeted by fungicides, thus addressing a relevant issue regarding disease management. Despite the controversial origin and evolution of fungal mt-genomes, the intrinsic mechanisms and molecular biology involved in their evolution will help to understand, at the molecular level, the strategies for fungal disease management.

Exploring the Relationship Among Divergence Time and Coding and Non-coding Elements in the Shaping of Fungal Mitochondrial Genomes

The order Hypocreales (Ascomycota) is composed of ubiquitous and ecologically diverse fungi such as saprobes, biotrophs, and pathogens. Despite their phylogenetic relationship, these species exhibit high variability in biomolecules production, lifestyle, and fitness. The mitochondria play an important role in the fungal biology, providing energy to the cells and regulating diverse processes, such as immune response. In spite of its importance, the mechanisms that shape fungal mitogenomes are still poorly understood. Herein, we investigated the variability and evolution of mitogenomes and its relationship with the divergence time using the order Hypocreales as a study model. We sequenced and annotated for the first time Trichoderma harzianum mitochondrial genome (mtDNA), which was compared to other 34 mtDNAs species that were publicly available. Comparative analysis revealed a substantial structural and size variation on non-coding mtDNA regions, despite the conservation of copy number, length, and structure of protein-coding elements. Interestingly, we observed a highly significant correlation between mitogenome length, and the number and size of non-coding sequences in mitochondrial genome. Among the non-coding elements, group I and II introns and homing endonucleases genes (HEGs) were the main contributors to discrepancies in mitogenomes structure and length. Several intronic sequences displayed sequence similarity among species, and some of them are conserved even at gene position, and were present in the majority of mitogenomes, indicating its origin in a common ancestor. On the other hand, we also identified species-specific introns that advocate for the origin by different mechanisms. Investigation of mitochondrial gene transfer to the nuclear genome revealed that nuclear copies of the nad5 are the most frequent while atp8, atp9, and cox3 could not be identified in any of the nuclear genomes analyzed. Moreover, we also estimated the divergence time of each species and investigated its relationship with coding and non-coding elements as well as with the length of mitogenomes. Altogether, our results demonstrated that introns and HEGs are key elements on mitogenome shaping and its presence on fast-evolving mtDNAs could be mostly explained by its divergence time, although the intron sharing profile suggests the involvement of other mechanisms on the mitochondrial genome evolution, such as horizontal transference.

Single-Spore Extraction for Genetic Analyses of Arbuscular Mycorrhizal Fungi

Biomass of arbuscular mycorrhizal fungi (AMF, Glomeromycota) is often only available in small quantities as these fungi are obligate biotrophs and many species are difficult to cultivate under controlled conditions. Here, I describe a simple, efficient approach to produce crude extracts from single or a small number of spores that can be used for genotyping AMF.

The complete mitochondrial genomes of two model ectomycorrhizal fungi (Laccaria): features, intron dynamics and phylogenetic implications

Laccaria amethystine and L. bicolor have served as model species for studying the life history and genetics of ectomycorrhizal fungi. However, the characterizations and variations of their mitogenomes are still unknown. In the present study, the mitogenomes of the two Laccaria species were assembled, annotated, and compared. The two mitogenomes of L. amethystine and L. bicolor comprised circular DNA molecules, with the sizes of 65,156 bp and 95,304 bp, respectively. Genome collinearity analysis revealed large-scale gene rearrangements between the two Laccaria species. Comparative mitogenome analysis indicated the introns of cox1 genes in Agaricales experienced frequent lost/gain eveants, which promoted the organization and size variations in Agaricales mitogenomes. Evolutionary analysis indicated the core protein-coding genes in the two mitogenomes were subject to strong pressure of purifying selection. Phylogenetic analysis using the Bayesian inference (BI) and Maximum likelihood (ML) methods based on a combined mitochondrial gene set resulted in identical and well-supported tree topologies, wherein the two Laccaria species were most closely related to Coprinopsis cinerea. This study severed as the first study on the mitogenomes of Laccaria species, which promoted a comprehensive understanding of the genetics and evolution of the model ectomycorrhizal fungi.

Mobile genetic elements explain size variation in the mitochondrial genomes of four closely-related Armillaria species

Background

Species in the genus Armillaria (fungi, basidiomycota) are well-known as saprophytes and pathogens on plants. Many of them cause white-rot root disease in diverse woody plants worldwide. Mitochondrial genomes (mitogenomes) are widely used in evolutionary and population studies, but despite the importance and wide distribution of Armillaria, the complete mitogenomes have not previously been reported for this genus. Meanwhile, the well-supported phylogeny of Armillaria species provides an excellent framework in which to study variation in mitogenomes and how they have evolved over time.

Results

Here we completely sequenced, assembled, and annotated the circular mitogenomes of four species: A. borealis, A. gallica, A. sinapina, and A. solidipes (116,443, 98,896, 103,563, and 122,167 bp, respectively). The variation in mitogenome size can be explained by variable numbers of mobile genetic elements, introns, and plasmid-related sequences. Most Armillaria introns contained open reading frames (ORFs) that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Insertions of mobile elements were also evident as fragments of plasmid-related sequences in Armillaria mitogenomes. We also found several truncated gene duplications in all four mitogenomes.

Conclusions

Our study showed that fungal mitogenomes have a high degree of variation in size, gene content, and genomic organization even among closely related species of Armillara. We suggest that mobile genetic elements invading introns and intergenic sequences in the Armillaria mitogenomes have played a significant role in shaping their genome structure. The mitogenome changes we describe here are consistent with widely accepted phylogenetic relationships among the four species.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5732-z) contains supplementary material, which is available to authorized users.

Fungal mitochondrial genomes and genetic polymorphisms

Mitochondria are the powerhouses of eukaryotic cells, responsible for ATP generation and playing a role in a diversity of cellular and organismal functions. Different from the majority of other intracellular membrane structures, mitochondria contain their own genetic materials that are capable of independent replication and inheritance. In this mini-review, we provide brief summaries of fungal mitochondrial genome structure, size, gene content, inheritance, and genetic variation. We pay special attention to the relative genetic polymorphisms of the mitochondrial vs nuclear genomes at the population level within individual fungal species. Among the 20 species/groups of species reviewed here, there is a range of variation among genes and species in the relative nuclear and mitochondrial genetic polymorphisms. Interestingly, most (15/20) showed a greater genetic diversity for nuclear genes and genomes than for mitochondrial genes and genomes, with the remaining five showing similar or slower nuclear genome genetic variations. This fungal pattern is different from the dominant pattern in animals, but more similar to that in plants. At present, the mechanisms for the variations among fungal species and the overall low level of mitochondrial sequence polymorphisms are not known. The increasing availability of population genomic data should help us reveal the potential genetic and ecological factors responsible for the observed variations.

Comparison of the Mitochondrial Genome Sequences of Six Annulohypoxylon stygium Isolates Suggests Short Fragment Insertions as a Potential Factor Leading to Larger Genomic Size

Mitochondrial DNA (mtDNA) is a core non-nuclear genetic material found in all eukaryotic organisms, the size of which varies extensively in the eumycota, even within species. In this study, mitochondrial genomes of six isolates of Annulohypoxylon stygium (Lév.) were assembled from raw reads from PacBio and Illumina sequencing. The diversity of genomic structures, conserved genes, intergenic regions and introns were analyzed and compared. Genome sizes ranged from 132 to 147 kb and contained the same sets of conserved protein-coding, tRNA and rRNA genes and shared the same gene arrangements and orientation. In addition, most intergenic regions were homogeneous and had similar sizes except for the region between cytochrome b (cob) and cytochrome c oxidase I (cox1) genes which ranged from 2,998 to 8,039 bp among the six isolates. Sixty-five intron insertion sites and 99 different introns were detected in these genomes. Each genome contained 45 or more introns, which varied in distribution and content. Introns from homologous insertion sites also showed high diversity in size, type and content. Comparison of introns at the same loci showed some complex introns, such as twintrons and ORF-less introns. There were 44 short fragment insertions detected within introns, intergenic regions, or as introns, some of them located at conserved domain regions of homing endonuclease genes. Insertions of short fragments such as small inverted repeats might affect or hinder the movement of introns, and these allowed for intron accumulation in the mitochondrial genomes analyzed, and enlarged their size. This study showed that the evolution of fungal mitochondrial introns is complex, and the results suggest short fragment insertions as a potential factor leading to larger mitochondrial genomes in A. stygium.

Genomic analyses reveal low mitochondrial and high nuclear diversity in the cyclosporin-producing fungus Tolypocladium inflatum

Mitochondrial DNA is generally regarded to evolve faster than nuclear DNA in animals, whereas if this is also true in fungi remains unclear. Herein, we annotate the first complete mitochondrial genome (mitogenome) of the cyclosporin-producing fungus Tolypocladium inflatum and report the genome-wide sequence variations among five isolates originating from distantly separated localities. We found that T. inflatum has among the most compact of fungal mitogenomes; its 25 kb DNA molecule encodes all standard fungal mitochondrial genes and harbors only one intron. Transcriptional analyses validated the expression of most conserved genes. We found several uncommon repetitive elements and evidence of gene transfer from the mitochondrion to the nucleus. Phylogenetic analyses confirmed the placement of T. inflatum in the fungal order Hypocreales although there was uncertainty on its family-level affiliation. Comparative genomic analyses among the five isolates identified an overall lower level of intraspecific variation in mitogenomes than in nuclear genomes; however, both the nuclear and mitochondrial genomes revealed similar isolate relationships, not correlating with geographic sources of these isolates. Our study shed new insights into the evolution of the medicinally important ascomycete T. inflatum.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

Comparative analysis of the mitochondrial genomes of Colletotrichum gloeosporioides sensu lato: insights into the evolution of a fungal species complex interacting with diverse plants

BACKGROUND

The fungal species complex Colletotrichum gloeosporioides sensu lato contains over 20 plant-interacting species. These species exhibit different life styles (e.g., endophytes, foliar and fruit pathogens) and show considerable variation in host and tissue adaptation strategies. Accurate species delimitation in C. gloeosporioides s.l. is very challenging due to nascent lineage boundaries and phenotypic plasticity, which strongly impedes studies of the complex's host-interaction biology. In this study, we first sequenced and compared nine mitogenomes belonging to four C. gloeosporioides s.l. species lineages (C. gloeosporioides, C. fructicola, C. aenigma, and C. siamense s.l.), and evaluated the usefulness of mitogenome sequence in complementing prevailing nuclear markers for species delimitation.

RESULTS

The C. gloeosporioides s.l. mitogenomes ranged between 52,671 and 58,666 bp in size, and each contained an identical set of genes transcribed in the same direction. Compared with previously reported Colletotrichum mitogenomes, these mitogenomes were uniquely featured by: (1) significantly larger genome size due to richer intron content and longer intergenic region; (2) striking GC content elevation at the intergenic region; and (3) considerable intron content variation among different species lineages. Compared with nuclear DNA markers commonly used in phylogeny, the mitogenome nucleotide diversity was extremely low, yet the mitogenome alignment contained the highest number of parsimony informative sites, which allowed the generation of a high-resolution phylogeny recognizing all taxonomic lineages, including ones belonging to the very nascent C. siamense s.l. complex. The tree topology was highly congruent with the phylogeny based on nuclear marker concatenation except for lineages within C. siamense s.l. Further comparative phylogenetic analysis indicated that lineage-specific rapid divergence of GS and SOD2 markers confounded concatenation-based species relationship inference.

CONCLUSIONS

This study sheds light on the evolution of C. gloeosporioides s.l. mitogenomes and demonstrates that mitogenome sequence can complement prevailing nuclear markers in improving species delimitation accuracy. The mitogenome sequences reported will be valuable resources for further genetic studies with C. gloeosporioides s.l. and other Colletotrichum species.

Comparative mitochondrial genomics toward exploring molecular markers in the medicinal fungus Cordyceps militaris

Cordyceps militaris is a fungus used for developing health food, but knowledge about its intraspecific differentiation is limited due to lack of efficient markers. Herein, we assembled the mitochondrial genomes of eight C. militaris strains and performed a comparative mitochondrial genomic analysis together with three previously reported mitochondrial genomes of the fungus. Sizes of the 11 mitochondrial genomes varied from 26.5 to 33.9 kb mainly due to variable intron contents (from two to eight introns per strain). Nucleotide variability varied according to different regions with non-coding regions showing higher variation frequency than coding regions. Recombination events were identified between some locus pairs but seemed not to contribute greatly to genetic variations of the fungus. Based on nucleotide diversity fluctuations across the alignment of all mitochondrial genomes, molecular markers with the potential to be used for future typing studies were determined.

Molecular diagnostic toolkit for Rhizophagus irregularis isolate DAOM-197198 using quantitative PCR assay targeting the mitochondrial genome

Rhizophagus irregularis (previously named Glomus irregulare) is one of the most widespread and common arbuscular mycorrhizal fungal (AMF) species. It has been recovered worldwide in agricultural and natural soils, and the isolate DAOM-197198 has been utilized as a commercial inoculant for two decades. Despite the ecological and economical importance of this taxon, specific markers for quantification of propagules by quantitative real-time PCR (qPCR) are extremely limited and none have been rigorously validated for quality control of manufactured products such as biofertilizers. From the sequencing of 14 complete AMF mitochondrial (mt) genomes, a qPCR assay using a hydrolysis probe designed in the single copy cox3-rnl intergenic region was tested and validated to specifically and accurately quantify the spores of R. irregularis isolate DAOM-197198. Specificity tests were performed using standard PCR and qPCR, and results clearly showed that the primers specifically amplified the isolate DAOM-197198, yielding a PCR product of 106 bp. According to the qPCR analyses on spores produced in vitro, the average copy number of mt genomes per spore was 3172 ± 304 SE (n = 6). Quantification assays were successfully undertaken on known and unknown samples in liquid suspensions and commercial dry formulations to show the accuracy, precision, robustness, and reproducibility of the qPCR assay. This study provides a powerful molecular toolkit specifically designed to quantify spores of the model AMF isolate DAOM-197198. The approach of molecular toolkit used in our study could be applied to other AMF taxa and will be useful to research institutions and governmental and industrial laboratories running routine quality control of AMF-based products.

The large (134.9 kb) mitochondrial genome of the glomeromycete Funneliformis mosseae

Funneliformis mosseae is among the most ecologically and economically important glomeromycete species and occurs both in natural and disturbed areas in a wide range of habitats and climates. In this study, we report the sequencing of the complete mitochondrial (mt) genome of F. mosseae isolate FL299 using 454 pyrosequencing and Illumina HiSeq technologies. This mt genome is a full-length circular chromosome of 134,925 bp, placing it among the largest mitochondrial DNAs (mtDNAs) in the fungal kingdom. A comparative analysis with publically available arbuscular mycorrhizal fungal mtDNAs revealed that the mtDNA of F. mosseae FL299 contained a very large number of insertions contributing to its expansion. The gene synteny was completely reshuffled compared to previously published glomeromycotan mtDNAs and several genes were oriented in an anti-sense direction. Furthermore, the presence of different types of introns and insertions in rnl (14 introns) made this gene very distinctive in Glomeromycota. The presence of alternative genetic codes in both initiation (GUG) and termination (UGA) codons was another new feature in this mtDNA compared to previously published glomeromycotan mt genomes. The phylogenetic analysis inferred from the analysis of 14 protein mt genes confirmed the position of the Glomeromycota clade as a sister group of Mortierellomycotina. This mt genome is the largest observed so far in Glomeromycota and the first mt genome within the Funneliformis clade, providing new opportunities to better understand their evolution and to develop molecular markers.

Independent mitochondrial and nuclear exchanges arising in Rhizophagus irregularis crossed-isolates support the presence of a mitochondrial segregation mechanism

BACKGROUND

Arbuscular mycorrhizal fungi (AMF) are members of the phylum Glomeromycota, an early divergent fungal lineage that forms symbiotic associations with the large majority of land plants. These organisms are asexual obligate biotrophs, meaning that they cannot complete their life cycle in the absence of a suitable host. These fungi can exchange genetic information through hyphal fusions (i.e. anastomosis) with genetically compatible isolates belonging to the same species. The occurrence of transient mitochondrial length-heteroplasmy through anastomosis between geographically distant Rhizophagus irregularis isolates was previously demonstrated in single spores resulting from crossing experiments. However, (1) the persistence of this phenomenon in monosporal culture lines from crossed parental isolates, (2) its correlation with nuclear exchanges and (3) the potential mechanisms responsible for mitochondrial inheritance are still unknown. Using the AMF model organism R. irregularis, we tested whether the presence of a heteroplasmic state in progeny spores was linked to the occurrence of nuclear exchanges and whether the previously observed heteroplasmic state persisted in monosporal in vitro crossed-culture lines. We also investigated the presence of a putative mitochondrial segregation apparatus in Glomeromycota by identifying proteins similar to those found in other fungal groups.

RESULTS

We observed the occurrence of biparental inheritance both for mitochondrial and nuclear markers tested in single spores obtained from crossed-isolates. However, only one parental mitochondrial DNA and nuclear genotype were recovered in each monosporal crossed-cultures, with an overrepresentation of certain mitochondrial haplotypes. These results strongly support the presence of a nuclear-independent mitochondrial segregation mechanism in R. irregularis. Furthermore, a nearly complete set of genes was identified with putative orthology to those found in other fungi and known to be associated with the mitochondrial segregation in Saccharomyces cerevisiae and filamentous fungi.

CONCLUSIONS

Our findings suggest that mitochondrial segregation might take place either during spore formation or colony development and that it might be independent of the nuclear segregation machinery. We present the basic building blocks for a better understanding of the mitochondrial inheritance process and segregation in these important symbiotic fungi. The comprehension of these processes is of great importance since it has been shown that different segregated lines of the same isolate can have variable effects on the host plant.

Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity

How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.

Genome diversity in arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi (Glomeromycota) are the most widespread and important symbionts of plants. They cannot be cultured without plants, are apparently asexual, and have multiple nuclei in a common cytoplasm. There is evidence for genetic variation among nuclei, and for segregation of this variation during growth, but these findings remain contentious. Recently, two papers have reported whole genome sequences for a strain of Rhizophagus irregularis; both suggest that genetic variation among nuclei is low. Genome assembly is very incomplete, though, so significant nuclear diversity cannot be excluded. While the diversity of nuclear genomes remains unresolved, multiple complete mitochondrial genomes are now available; there is virtually no variation within isolates, but significant variation between them.

Comparison of mitochondrial genomes provides insights into intron dynamics and evolution in the caterpillar fungus Cordyceps militaris

Intra-specific comparison of mitochondrial genomes can help elucidate the evolution of a species, however it has not been performed for hypocrealean fungi that form diverse symbiotic associations with other organisms. In this study, comparative analyses of three completely sequenced mitochondrial genomes of a hypocrealean fungus, Cordyceps militaris, the type species of Cordyceps genus, revealed that the introns were the main contributors to mitochondrial genome size variations among strains. Mitochondrial genes in C. militaris have been invaded by group I introns in at least eight positions. PCR assays of various C. militaris isolates showed abundant variations of intron presence/absence among strains at seven of the eight intronic loci. Although the ancestral intron pattern was inferred to contain all eight introns, loss and/or gain events occurred for seven of the eight introns. These introns invaded the C. militaris mitochondrial genome probably by horizontal transfer from other fungi, and intron insertions into intronless genes in C. militaris were accompanied by co-conversions of upstream exon sequences especially for those introns targeting protein-coding genes. We also detected phylogenetic congruence between the intron and exon trees at each individual locus, consistent with the ancestral mitochondria of C. militaris as having all eight introns. This study helps to explain the evolution of C. militaris mitochondrial genomes and will facilitate population genetic studies of this medicinally important fungus.

The mitochondrial genome of the glomeromycete Rhizophagus sp. DAOM 213198 reveals an unusual organization consisting of two circular chromosomes

Mitochondrial (mt) genomes are intensively studied in Ascomycota and Basidiomycota, but they are poorly documented in basal fungal lineages. In this study, we sequenced the complete mtDNA of Rhizophagus sp. DAOM 213198, a close relative to Rhizophagus irregularis, a widespread, ecologically and economical relevant species belonging to Glomeromycota. Unlike all other known taxonomically close relatives harboring a full-length circular chromosome, mtDNA of Rhizophagus sp. reveals an unusual organization with two circular chromosomes of 61,964 and 29,078 bp. The large chromosome contained nine protein-coding genes (atp9, nad5, cob, nad4, nad1, nad4L, cox1, cox2, and atp8), small subunit rRNA gene (rns), and harbored 20 tRNA-coding genes and 10 orfs, while the small chromosome contained five protein-coding genes (atp6, nad2, nad3, nad6, and cox3), large subunit rRNA gene (rnl) in addition to 5 tRNA-coding genes, and 8 plasmid-related DNA polymerases (dpo). Although structural variation of plant mt genomes is well documented, this study is the first report of the presence of two circular mt genomes in arbuscular mycorrhizal fungi. Interestingly, the presence of dpo at the breakage point in intergenes cox1-cox2 and rnl-atp6 for large and small mtDNAs, respectively, could be responsible for the conversion of Rhizophagus sp. mtDNA into two chromosomes. Using quantitative real-time polymerase chain reaction, we found that both mtDNAs have an equal abundance. This study reports a novel mtDNA organization in Glomeromycota and highlights the importance of studying early divergent fungal lineages to describe novel evolutionary pathways in the fungal kingdom.

Intraisolate mitochondrial genetic polymorphism and gene variants coexpression in arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi (AMF) are multinucleated and coenocytic organisms, in which the extent of the intraisolate nuclear genetic variation has been a source of debate. Conversely, their mitochondrial genomes (mtDNAs) have appeared to be homogeneous within isolates in all next generation sequencing (NGS)-based studies. Although several lines of evidence have challenged mtDNA homogeneity in AMF, extensive survey to investigate intraisolate allelic diversity has not previously been undertaken. In this study, we used a conventional polymerase chain reaction -based approach on selected mitochondrial regions with a high-fidelity DNA polymerase, followed by cloning and Sanger sequencing. Two isolates of Rhizophagus irregularis were used, one cultivated in vitro for several generations (DAOM-197198) and the other recently isolated from the field (DAOM-242422). At different loci in both isolates, we found intraisolate allelic variation within the mtDNA and in a single copy nuclear marker, which highlighted the presence of several nonsynonymous mutations in protein coding genes. We confirmed that some of this variation persisted in the transcriptome, giving rise to at least four distinct nad4 transcripts in DAOM-197198. We also detected the presence of numerous mitochondrial DNA copies within nuclear genomes (numts), providing insights to understand this important evolutionary process in AMF. Our study reveals that genetic variation in Glomeromycota is higher than what had been previously assumed and also suggests that it could have been grossly underestimated in most NGS-based AMF studies, both in mitochondrial and nuclear genomes, due to the presence of low-level mutations.

I-OmiI and I-OmiII: two intron-encoded homing endonucleases within the Ophiostoma minus rns gene

The mitochondrial small subunit ribosomal RNA (rns) gene of the ascomycetous fungus Ophiostoma minus [strain WIN(M)371] was found to contain a group IC2 and a group IIB1 intron at positions mS569 and mS952 respectively. Both introns have open reading frames (ORFs) embedded that encode double motif LAGLIDADG homing endonucleases (I-OmiI and I-OmiII respectively). Codon-optimized versions of I-OmiI and I-OmiII were synthesized for overexpression in Escherichia coli. The in vitro characterization of I-OmiII showed that it is a functional homing endonuclease that cleaves the rns target site two nucleotides upstream (sense strand) of the intron insertion site generating 4 nucleotide 3' overhangs. The endonuclease activity of I-OmiII was tested using linear and circular substrates and cleavage activity was evaluated at various temperatures. The I-OmiI protein was expressed in E. coli, but purification was difficult, thus the endonuclease activity of this protein was tested via in vivo assays. Overall this study showed that there are many native forms of functional homing endonucleases yet to be discovered among fungal mtDNA genomes.

Sequencing and comparison of the mitochondrial COI gene from isolates of Arbuscular Mycorrhizal Fungi belonging to Gigasporaceae and Glomeraceae families

Arbuscular Mycorrhizal Fungi (AMF) are well known for their ecological importance and their positive influence on plants. The genetics and phylogeny of this group of fungi have long been debated. Nuclear markers are the main tools used for phylogenetic analyses, but they have sometimes proved difficult to use because of their extreme variability. Therefore, the attention of researchers has been moving towards other genomic markers, in particular those from the mitochondrial DNA. In this study, 46 sequences of different AMF isolates belonging to two main clades Gigasporaceae and Glomeraceae have been obtained from the mitochondrial gene coding for the Cytochrome c Oxidase I (COI), representing the largest dataset to date of AMF COI sequences. A very low level of divergence was recorded in the COI sequences from the Gigasporaceae, which could reflect either a slow rate of evolution or a more recent evolutionary divergence of this group. On the other hand, the COI sequence divergence between Gigasporaceae and Glomeraceae was high, with synonymous divergence reaching saturated levels. This work also showed the difficulty in developing valuable mitochondrial markers able to effectively distinguish all Glomeromycota species, especially those belonging to Gigasporaceae, yet it represents a first step towards the development of a full mtDNA-based dataset which can be used for further phylogenetic investigations of this fungal phylum.

Detection of a transient mitochondrial DNA heteroplasmy in the progeny of crossed genetically divergent isolates of arbuscular mycorrhizal fungi

Nonself fusion and nuclear genetic exchange have been documented in arbuscular mycorrhizal fungi (AMF), particularly in Rhizophagus irregularis. However, mitochondrial transmission accompanying nonself fusion of genetically divergent isolates remains unknown. Here, we tested the hypothesis that mitochondrial DNA (mtDNA) heteroplasmy occurs in the progeny of spores, obtained by crossing genetically divergent mtDNAs in R. irregularis isolates. Three isolates of geographically distant locations were used to investigate nonself fusions and mtDNA transmission to the progeny. We sequenced two additional mtDNAs of two R. irregularis isolates and developed isolate-specific size-variable markers in intergenic regions of these isolates and those of DAOM-197198. We achieved three crossing combinations in pre-symbiotic and symbiotic phases. Progeny spores per crossing combination were genotyped using isolate-specific markers. We found evidence that nonself recognition occurs between isolates originating from different continents both in pre-symbiotic and symbiotic phases. Genotyping patterns of individual spores from the progeny clearly showed the presence of markers of the two parental mtDNA haplotypes. Our results demonstrate that mtDNA heteroplasmy occurs in the progeny of the crossed isolates. However, this heteroplasmy appears to be a transient stage because all the live progeny spores that were able to germinate showed only one mtDNA haplotype.

Rapid mitochondrial genome evolution through invasion of mobile elements in two closely related species of arbuscular mycorrhizal fungi

Arbuscular mycorrhizal fungi (AMF) are common and important plant symbionts. They have coenocytic hyphae and form multinucleated spores. The nuclear genome of AMF is polymorphic and its organization is not well understood, which makes the development of reliable molecular markers challenging. In stark contrast, their mitochondrial genome (mtDNA) is homogeneous. To assess the intra- and inter-specific mitochondrial variability in closely related Glomus species, we performed 454 sequencing on total genomic DNA of Glomus sp. isolate DAOM-229456 and we compared its mtDNA with two G. irregulare isolates. We found that the mtDNA of Glomus sp. is homogeneous, identical in gene order and, with respect to the sequences of coding regions, almost identical to G. irregulare. However, certain genomic regions vary substantially, due to insertions/deletions of elements such as introns, mitochondrial plasmid-like DNA polymerase genes and mobile open reading frames. We found no evidence of mitochondrial or cytoplasmic plasmids in Glomus species, and mobile ORFs in Glomus are responsible for the formation of four gene hybrids in atp6, atp9, cox2, and nad3, which are most probably the result of horizontal gene transfer and are expressed at the mRNA level. We found evidence for substantial sequence variation in defined regions of mtDNA, even among closely related isolates with otherwise identical coding gene sequences. This variation makes it possible to design reliable intra- and inter-specific markers.

Mitochondrial genome rearrangements in glomus species triggered by homologous recombination between distinct mtDNA haplotypes

Comparative mitochondrial genomics of arbuscular mycorrhizal fungi (AMF) provide new avenues to overcome long-lasting obstacles that have hampered studies aimed at understanding the community structure, diversity, and evolution of these multinucleated and genetically polymorphic organisms.AMF mitochondrial (mt) genomes are homogeneous within isolates, and their intergenic regions harbor numerous mobile elements that have rapidly diverged, including homing endonuclease genes, small inverted repeats, and plasmid-related DNA polymerase genes (dpo), making them suitable targets for the development of reliable strain-specific markers. However, these elements may also lead to genome rearrangements through homologous recombination, although this has never previously been reported in this group of obligate symbiotic fungi. To investigate whether such rearrangements are present and caused by mobile elements in AMF, the mitochondrial genomes from two Glomeraceae members (i.e., Glomus cerebriforme and Glomus sp.) with substantial mtDNA synteny divergence,were sequenced and compared with available glomeromycotan mitochondrial genomes. We used an extensive nucleotide/protein similarity network-based approach to investigated podiversity in AMF as well as in other organisms for which sequences are publicly available. We provide strong evidence of dpo-induced inter-haplotype recombination, leading to a reshuffled mitochondrial genome in Glomus sp. These findings raise questions as to whether AMF single spore cultivations artificially underestimate mtDNA genetic diversity.We assessed potential dpo dispersal mechanisms in AMF and inferred a robust phylogenetic relationship with plant mitochondrial plasmids. Along with other indirect evidence, our analyses indicate that members of the Glomeromycota phylum are potential donors of mitochondrial plasmids to plants.