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

Comparative Mitogenomics of Fungal Species in Stachybotryaceae Provides Evolutionary Insights into Hypocreales

Stachybotrys chartarum is one of the world’s ten most feared fungi within the family Stachybotryaceae, although to date, not a single mitogenome has been documented for Stachybotryaceae. Herein, six mitogenomes of four different species in Stachybotryaceae are newly reported. The S. chartarum mitogenome was 30.7 kb in length and contained two introns (one each in rnl and cox1). A comparison of the mitogenomes of three different individuals of S. chartarum showed few nucleotide variations and conservation of gene content/order and intron insertion. A comparison of the mitogenomes of four different Stachybotryaceae species (Memnoniella echinata, Myrothecium inundatum, S. chartarum, and S. chlorohalonata), however, revealed variations in intron insertion, gene order/content, and nad2/nad3 joining pattern. Further investigations on all Hypocreales species with available mitogenomes showed greater variabilities in gene order (six patterns) and nad2/nad3 joining pattern (five patterns) although a dominant pattern always existed in each case. Ancestral state estimation showed that in each case the dominant pattern was always more ancestral than those rare patterns. Phylogenetic analyses based on mitochondrion-encoded genes supported the placement of Stachybotryaceae in Hypocreales. The crown age of Stachybotryaceae was estimated to be approximately the Early Cretaceous (141–142 Mya). This study greatly promotes our understanding of the evolution of fungal species in Hypocreales.

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.

Comparative Mitogenomic Analysis and the Evolution of Rhizoctonia solani Anastomosis Groups

Mitochondria are the major energy source for cell functions. However, for the plant fungal pathogens, mitogenome variations and their roles during the host infection processes remain largely unknown. Rhizoctonia solani, an important soil-borne pathogen, forms different anastomosis groups (AGs) and adapts to a broad range of hosts in nature. Here, we reported three complete mitogenomes of AG1-IA RSIA1, AG1-IB RSIB1, and AG1-IC, and performed a comparative analysis with nine published Rhizoctonia mitogenomes (AG1-IA XN, AG1-IB 7/3/14, AG3, AG4, and five Rhizoctonia sp. mitogenomes). These mitogenomes encoded 15 typical proteins (cox1-3, cob, atp6, atp8-9, nad1-6, nad4L, and rps3) and several LAGLIDADG/GIY-YIG endonucleases with sizes ranging from 109,017 bp (Rhizoctonia sp. SM) to 235,849 bp (AG3). We found that their large sizes were mainly contributed by repeat sequences and genes encoding endonucleases. We identified the complete sequence of the rps3 gene in 10 Rhizoctonia mitogenomes, which contained 14 positively selected sites. Moreover, we inferred a robust maximum-likelihood phylogeny of 32 Basidiomycota mitogenomes, representing that seven R. solani and other five Rhizoctonia sp. lineages formed two parallel branches in Agaricomycotina. The comparative analysis showed that mitogenomes of Basidiomycota pathogens had high GC content and mitogenomes of R. solani had high repeat content. Compared to other strains, the AG1-IC strain had low substitution rates, which may affect its mitochondrial phylogenetic placement in the R. solani clade. Additionally, with the published RNA-seq data, we investigated gene expression patterns from different AGs during host infection stages. The expressed genes from AG1-IA (host: rice) and AG3 (host: potato) mainly formed four groups by k-mean partitioning analysis. However, conserved genes represented varied expression patterns, and only the patterns of rps3-nad2 and nad1-m3g18/mag28 (an LAGLIDADG endonuclease) were conserved in AG1-IA and AG3 as shown by the correlation coefficient analysis, suggesting regulation of gene repertoires adapting to infect varied hosts. The results of variations in mitogenome characteristics and the gene substitution rates and expression patterns may provide insights into the evolution of R. solani mitogenomes.

Dynamic evolution of eukaryotic mitochondrial and nuclear genomes: a case study in the gourmet pine mushroom Tricholoma matsutake

Fungi, as eukaryotic organisms, contain two genomes, the mitochondrial genome and the nuclear genome, in their cells. How the two genomes evolve and correlate to each other is debated. Herein, taking the gourmet pine mushroom Tricholoma matsutake as an example, we performed comparative mitogenomic analysis using samples collected from diverse locations and compared the evolution of the two genomes. The T. matsutake mitogenome encodes 49 genes and is rich of repetitive and non-coding DNAs. Six genes were invaded by up to 11 group I introns, with one cox1 intron cox1P372 showing presence/absence dynamics among different samples. Bioinformatic analyses suggested limited or no evidence of mitochondrial heteroplasmy. Interestingly, hundreds of mitochondrial DNA fragments were found in the nuclear genome, with several larger than 500 nt confirmed by PCR assays and read count comparisons, indicating clear evidence of transfer of mitochondrial DNA into the nuclear genome. Nuclear DNA of T. matsutake showed a higher mutation rate than mitochondrial DNA. Furthermore, we found evidence of incongruence between phylogenetic trees derived from mitogenome and nuclear DNA sequences. Together, our results reveal the dynamic genome evolution of the gourmet pine mushroom.

Characterization of the mitochondrial genome of the pathogenic fungus Scytalidium auriculariicola (Leotiomycetes) and insights into its phylogenetics

Scytalidium auriculariicola is the causative pathogen of slippery scar disease in the cultivated cloud ear fungus, Auricularia polytricha. In the present study, the mitogenome of S. auriculariicola was sequenced and assembled by next-generation sequencing technology. The circular mitogenome is 96,857 bp long and contains 56 protein-coding genes, 2 ribosomal RNA genes, and 30 transfer RNA genes (tRNAs). The high frequency of A and T used in codons contributed to the high AT content (73.70%) of the S. auriculariicola mitogenome. Comparative analysis indicated that the base composition and the number of introns and protein-coding genes in the S. auriculariicola mitogenome varied from that of other Leotiomycetes mitogenomes, including a uniquely positive AT skew. Five distinct groups were found in the gene arrangements of Leotiomycetes. Phylogenetic analyses based on combined gene datasets (15 protein-coding genes) yielded well-supported (BPP = 1) topologies. A single-gene phylogenetic tree indicated that the nad4 gene may be useful as a molecular marker to analyze the phylogenetic relationships of Leotiomycetes species. This study is the first report on the mitochondrial genome of the genus Scytalidium, and it will contribute to our understanding of the population genetics and evolution of S. auriculariicola and related species.

Mitochondrial genome in Hypsizygus marmoreus and its evolution in Dikarya

BACKGROUND

Hypsizygus marmoreus, a high value commercialized edible mushroom is widely cultivated in East Asia, and has become one of the most popular edible mushrooms because of its rich nutritional and medicinal value. Mitochondria are vital organelles, and play various essential roles in eukaryotic cells.

RESULTS

In this study, we provide the Hypsizygus marmoreus mitochondrial (mt) genome assembly: the circular sequence is 102,752 bp in size and contains 15 putative protein-coding genes, 2 ribosomal RNAs subunits and 28 tRNAs. We compared the mt genomes of the 27 fungal species in the Pezizomycotina and Basidiomycotina subphyla, with the results revealing that H. marmoreus is a sister to Tricholoma matsutake and the phylogenetic distribution of this fungus based on the mt genome. Phylogenetic analysis shows that Ascomycetes mitochondria started to diverge earlier than that of Basidiomycetes and supported the robustness of the hyper metric tree. The fungal sequences are highly polymorphic and gene order varies significantly in the dikarya data set, suggesting a correlation between the gene order and divergence time in the fungi mt genome. To detect the mt genome variations in H. marmoreus, we analyzed the mtDNA sequences of 48 strains. The phylogeny and variation sited type statistics of H. marmoreus provide clear-cut evidence for the existence of four well-defined cultivations isolated lineages, suggesting female ancestor origin of H. marmoreus. Furthermore, variations on two loci were further identified to be molecular markers for distinguishing the subgroup containing 32 strains of other strains. Fifteen conserved protein-coding genes of mtDNAs were analyzed, with fourteen revealed to be under purifying selection in the examined fungal species, suggesting the rapid evolution was caused by positive selection of this gene.

CONCLUSIONS

Our studies have provided new reference mt genomes and comparisons between species and intraspecies with other strains, and provided future perspectives for assessing diversity and origin of H. marmoreus.

Proposal of a new nomenclature for introns in protein-coding genes in fungal mitogenomes

Fungal mitochondrial genes are often invaded by group I or II introns, which represent an ideal marker for understanding fungal evolution. A standard nomenclature of mitochondrial introns is needed to avoid confusion when comparing different fungal mitogenomes. Currently, there has been a standard nomenclature for introns present in rRNA genes, but there is a lack of a standard nomenclature for introns present in protein-coding genes. In this study, we propose a new nomenclature system for introns in fungal mitochondrial protein-coding genes based on (1) three-letter abbreviation of host scientific name, (2) host gene name, (3), one capital letter P (for group I introns), S (for group II introns), or U (for introns with unknown types), and (4) intron insertion site in the host gene according to the cyclosporin-producing fungus Tolypocladium inflatum. The suggested nomenclature was proved feasible by naming introns present in mitogenomes of 16 fungi of different phyla, including both basal and higher fungal lineages although minor adjustment of the nomenclature is needed to fit certain special conditions. The nomenclature also had the potential to name plant/protist/animal mitochondrial introns. We hope future studies follow the proposed nomenclature to ensure direct comparison across different studies.

Complete mitogenome of the entomopathogenic fungus Sporothrix insectorum RCEF 264 and comparative mitogenomics in Ophiostomatales

The fungal order Ophiostomatales contains numerous species important in medical fields, agriculture, and forestry, and several species have had available mitogenome information. The nuclear genome of the entomopathogenic fungus Sporothrix insectorum has been reported, while its mitogenome remains unknown. Herein, we firstly described the mitogenome of S. insectorum RCEF 264 and then compared Ophiostomatales mitogenomes from both interspecific and intraspecific perspectives. The mitogenome of S. insectorum RCEF 264 was 31,454 bp in length, containing typical fungal mitochondrial genes plus rnpB. Four group I introns interrupted rnl and cox1. Phylogenetic analyses confirmed the placement of S. insectorum RCEF 264 in Ophiostomatales. Comparison of mitogenomes among seven Ophiostomatales species revealed conserved gene contents and a high synteny, although there were also some differences among them. Their mitogenomes showed more than two-fold variations (26.6-65.1 kb) in size, with a total of 37 intron insertional loci from 11 genes (1-25 introns per species). The sole intron shared by all species was an rps3-encoding intron in rnl (mL2450), and this intron-based phylogeny was highly consistent with those constructed using mitochondrial/nuclear genes, suggesting convergent evolution of this intron with Ophiostomatales species. The dendrogram based on presence/absence patterns at all intron loci was quite different from those based on mitochondrial/nuclear genes. Comparison of mitogenomes among two to three intraspecific individuals in Ophiostoma novo-ulmi subsp. novo-ulmi and Sporothrix schenckii revealed mitogenome size variations due to single-nucleotide polymorphisms (SNPs) and indels but without fluctuation of intron numbers for each species. This study greatly enhanced our understanding of mitogenome evolution in Ophiostomatales.

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.