Complete mitochondrial genome of the endophytic fungus Pestalotiopsis fici: features and evolution

Complete mitochondrial genome sequence of the white root rot pathogen Dematophora necatrix (Xylariaceae: Xylariales)

The mitochondrial genome of Dematophora necatrix is 121,350 base pairs in length with a G + C content of 30.19%. Phylogenetic analysis showed that D. necatrix grouped with other members of the Xylariaceae, with which its mitogenome also shares a broadly similar architecture and gene content. The D. necatrix mitogenome contains 14 protein-coding and 26 tRNA-encoding genes, as well as one copy each of the rnl, rns, rps3 and nat1 genes. However, as much as 80% of this genome is intronic or non-coding. This is likely due to expansions and rearrangements caused by the large number of group I introns and the homing endonucleases and reverse-transcriptases they encode. Our study thus provides a valuable foundation from which to explore the mitochondrion's role in the biology of D. necatrix, and also serves as a resource for investigating the pathogen's population biology and general ecology.

Insights into Fungal Mitochondrial Genomes and Inheritance Based on Current Findings from Yeast-like Fungi

The primary functions of mitochondria are to produce energy and participate in the apoptosis of cells, with them being highly conserved among eukaryotes. However, the composition of mitochondrial genomes, mitochondrial DNA (mtDNA) replication, and mitochondrial inheritance varies significantly among animals, plants, and fungi. Especially in fungi, there exists a rich diversity of mitochondrial genomes, as well as various replication and inheritance mechanisms. Therefore, a comprehensive understanding of fungal mitochondria is crucial for unraveling the evolutionary history of mitochondria in eukaryotes. In this review, we have organized existing reports to systematically describe and summarize the composition of yeast-like fungal mitochondrial genomes from three perspectives: mitochondrial genome structure, encoded genes, and mobile elements. We have also provided a systematic overview of the mechanisms in mtDNA replication and mitochondrial inheritance during bisexual mating. Additionally, we have discussed and proposed open questions that require further investigation for clarification.

Metabolome and Genome Analysis of a Novel Endophytic Fungus Aureobasidium pullulans KB3: Discovery of Polyketones and Polyketone Biosynthesis Pathway

Endophytic fungi associated with plants may contain undiscovered bioactive compounds. Under standard laboratory conditions, most undiscovered compounds are inactive, whereas their production could be stimulated under different cultivation conditions. In this study, six endophytic fungi were isolated from the bark of Koelreuteria paniculata in Quancheng Park, Jinan City, Shandong Province, one of which was identified as a new subspecies of Aureobasidium pullulans, named A. pullulans KB3. Additionally, metabolomic tools were used to screen suitable media for A. pullulans KB3 fermentation, and the results showed that peptone dextrose medium (PDM) was more beneficial to culture A. pullulans KB3 for isolation of novel compounds. Sphaerolone, a polyketone compound, was initially isolated from A. pullulans KB3 via scaled up fermentation utilizing PDM. Additionally, the whole-genome DNA of A. pullulans KB3 was sequenced to facilitate compound isolation and identify the biosynthesis gene clusters (BGCs). This study reports the multi-omics (metabolome and genome) analysis of A. pullulans KB3, laying the foundation for discovering novel compounds of silent BGCs and identifying their biosynthesis pathway.

Comparative Mitogenomics Analysis Revealed Evolutionary Divergence among Neopestalotiopsis Species Complex (Fungi: Xylariales)

The genus Neopestalotiopsis consists of obligate parasites that cause ring spot, scab, and leaf blight diseases in higher plant species. We assembled the three complete mitogenomes for the guava fruit ring spot pathogen, Neopestalotiopsis cubana. The mitogenomes are circular, with sizes of 38,666 bp, 33,846 bp, and 32,593 bp. The comparative analyses with Pestalotiopsis fici showed that N. cubana differs greatly from it in the length of the mitogenomes and the number of introns. Moreover, they showed significant differences in the gene content and tRNAs. The two genera showed little difference in gene skewness and codon preference for core protein-coding genes (PCGs). We compared gene sequencing in the mitogenomes of the order Xylariales and found large-scale gene rearrangement events, such as gene translocations and the duplication of tRNAs. N. cubana shows a unique evolutionary position in the phylum Ascomycota constructed in phylogenetic analyses. We also found a more concentrated distribution of evolutionary pressures on the PCGs of Neopestalotiopsis in the phylum Ascomycota and that they are under little selective pressure compared to other species and are subjected to purifying selection. This study explores the evolutionary dynamics of the mitogenomes of Neopestalotiopsis and provides important support for genetic and taxonomic studies.

The complete mitochondrial genome of the rice blast fungus Pyricularia oryzae Cavara 1892 strain Guy11 and phylogenetic analysis

The complete mitochondrial genome of Pyricularia oryzae Cavara 1892 strain Guy11 is 34,865 bp in length (GenBank accession number OP095391), containing 29 tRNA genes, 2 rRNA genes, and 15 protein-coding genes (PCGs). The gene order and orientation are novel compared to other Sordariomycetes species with sequenced mitogenomes in the GenBank database. Phylogenetic analysis suggests that P. oryzae Guy11 and 19 other Sordariomycetes species form a monophyletic group. The complete mitochondrial sequence of P. oryzae Guy11 will be a valuable resource for species identification, population genetics, phylogenetics, and comparative genomics studies in Sordariomycetes and Magnaporthales.

Mitochondrial genome of Cordyceps blackwelliae: organization, transcription, and evolutionary insights into Cordyceps

Cordyceps is a diverse genus of insect pathogenic fungi, with about 180 accepted species, including some well-known ones used as ethnic medicine and/or functional food. Nevertheless, mitogenomes are only available for four members of the genus. The current study reports the mitogenome of Cordyceps blackwelliae, a newly described entomopathogenic fungus. The 42,257-bp mitogenome of the fungus encoded genes typically found in fungal mitogenomes, and a total of 14 introns inserted into seven genes, including cob (1 intron), cox1 (4), cox3 (3), nad1 (1), nad4 (1), nad5 (1), and rnl (3). RNA-Seq analysis revealed differential expression of mitochondrial genes and supported annotations resulting from in silico analysis. There was clear evidence for polycistronic transcription and alternative splicing of mitochondrial genes. Comparison among mitogenomes of five different Cordyceps species (i.e., C. blackwelliae, C. chanhua, C. militaris, C. pruinosa, and C. tenuipes) revealed a high synteny, with mitogenome size expansion correlating with intron insertions. Different mitochondrial protein-coding genes showed variable degrees of genetic differentiation among these species, but they were all under purifying selection. Mitochondrial phylogeny based on either nucleotide or amino acid sequences confirmed the taxonomic position of C. blackwelliae in Cordycipitaceae, clustering together with C. chanhua. This study promotes our understanding of fungal evolution in Cordyceps.

Comparative structural analysis on the mitochondrial DNAs from various strains of Lentinula edodes

The evolution of mitochondria through variations in mitochondrial DNA (mtDNA) is one of the intriguing questions in eukaryotic cells. In order to assess the causes of the variations in mitochondria, the mtDNAs of the 21 strains of Lentinula edodes were assembled for this study, and analyzed together with four published mtDNA sequences. The mtDNAs were within the sizes of 117 kb ~ 122 kb. The gene number was observed consistent except for two mtDNAs, which carry a duplicated trnG1-trnG2 unit or a putative gene deletion. The size variation was largely attributed to the number of introns, repeated sequences, transposable elements (TEs), and plasmid-related sequences. Intron loss and gain were found from cox1, rnl, and rns of three mtDNAs. Loss of two introns in cox1 of KY217797.1 reduced its size by 2.7 kb, making it the smallest cox1 gene (8.4 kb) among the cox1s of the 25 mtDNAs, whereas gain of a Group II intron (2.65 kb) and loss of a Group I intron (1.7 kb) in cox1 of MF774813.1 resulted in the longest cox1 (12 kb). In rnl of L. edodes, we discovered four intron insertion consensus sequences which were unique to basidiomycetes but not ascomycetes. Differential incorporation of introns was the primary cause of the rnl size polymorphism. Homing endonucleases (HEGs) were suggestively involved in the mobilization of the introns because all of the introns have HEG genes of the LAGRIDADG or GIY-YIG families with the conserved HEG cleavage sites. TEs contributed to 11.04% of the mtDNA size in average, of which 7.08% was LTR-retrotransposon and 3.96% was DNA transposon, whereas the repeated sequences covered 4.6% of the mtDNA. The repeat numbers were variable in a strain-dependent manner. Both the TEs and repeated sequences were mostly found in the intronic and intergenic regions. Lastly, two major deletions were found in the plasmid-related sequence regions (pol2-pol3 and pol1-atp8) in the five mtDNAs. Particularly, the 6.8 kb-long deletion at pol2-pol3 region made MF774813.1 the shortest mtDNA of all. Our results demonstrate that mtDNA is a dynamic molecule that persistently evolves over a short period of time by insertion/deletion and repetition of DNA segments at the strain level.

Comparative Mitogenomic Analysis Reveals Intraspecific, Interspecific Variations and Genetic Diversity of Medical Fungus Ganoderma

Ganoderma species are widely distributed in the world with high diversity. Some species are considered to be pathogenic fungi while others are used as traditional medicine in Asia. In this study, we sequenced and assembled four Ganoderma complete mitogenomes, including G. subamboinense s118, G. lucidum s37, G. lingzhi s62, and G. lingzhi s74. The sizes of the four mitogenomes ranged from 50,603 to 73,416 bp. All Ganoderma specimens had a full set of core protein-coding genes (PCGs), and the rps3 gene of Ganoderma species was detected to be under positive or relaxed selection. We found that the non-conserved PCGs, which encode RNA polymerases, DNA polymerases, homing endonucleases, and unknown functional proteins, are dynamic within and between Ganoderma species. Introns were thought to be the main contributing factor in Ganoderma mitogenome size variation (p < 0.01). Frequent intron loss/gain events were detected within and between Ganoderma species. The mitogenome of G. lucidum s26 gained intron P637 in the cox3 gene compared with the other two G. lucidum mitogenomes. In addition, some rare introns in Ganoderma were detected in distinct Basidiomycetes, indicating potential gene transfer events. Comparative mitogenomic analysis revealed that gene arrangements also varied within and between Ganoderma mitogenomes. Using maximum likelihood and Bayesian inference methods with a combined mitochondrial gene dataset, phylogenetic analyses generated identical, well-supported tree topologies for 71 Agaricomycetes species. This study reveals intraspecific and interspecific variations of the Ganoderma mitogenomes, which promotes the understanding of the origin, evolution, and genetic diversity of Ganoderma species.

Characterization of the mitogenome of Gongronella sp. w5 reveals substantial variation in Mucoromycota

Gongronella is a genus of fungi in Mucorales (Mucoromycota). Some of its members have important biotechnological applications, but until now, not a single mitogenome has been characterized in Gongronella. Here, we present the complete mitogenome assembly of Gongronella sp. w5, a soil isolate known to interact with plants and several fungi. Its 36,593-bp circular mitogenome encodes the large and small subunit rRNAs, 14 standard mitochondrial proteins, 24 tRNAs, three free-standing ORF proteins, and the RNA subunit of RNase P (rnpB). These genes arrange in an order novel to known fungal mitogenomes. Three group I introns are present in the cob, cox1, and nad5 genes, respectively, and they are probably acquired by horizontal gene transfer. Phylogenetic analysis based on mitochondrion-encoded proteins supports the grouping of Gongronella sp. w5 with Absidia glauca, forming the Cunninghamellaceae clade within Mucoromycota. Gongronella and most other Mucoromycota species are predicted to use the standard genetic code in mitochondrial translation, rather than code 4 assigned by GenBank. A comparison among seven publicly available mitogenomes in Mucoromycota reveals the presence of the same 14 typical protein-coding genes plus rnpB, yet substantial variation in mitogenome size, intron number, gene order, and orientation. In this comparison, the uniqueness of Gongronella is evident from similarly large differences to its closest phylogenetic neighbor, A. glauca. This study promotes our understanding of fungal evolution in Mucoromycota. KEY POINTS: • This study reports the first mitogenome in Gongronella, which presents a novel gene order. • Different Mucoromycota mitogenomes show substantial variation of gene organizations. • Most Mucoromycota species use the standard genetic code to translate mitochondrial genes.

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.

Revisiting the Neurospora crassa mitochondrial genome

The mitochondrial genome of Neurospora crassa has been less studied than its nuclear counterpart, yet it holds great potential for understanding the diversity and evolution of this important fungus. Here we describe a new mitochondrial DNA (mtDNA) complete sequence of a N. crassa wild type strain. The genome with 64 839 bp revealed 21 protein-coding genes and several hypothetical open reading frames with no significant homology to any described gene. Five large repetitive regions were identified across the genome, including partial or complete genes. The largest repeated region holds a partial nd2 section that was also detected in Neurospora intermedia, suggesting a rearrangement that occurred before the N. crassa speciation. Interestingly, N. crassa has a palindrome adjacent to the partial nd2 repeated region possibly related to the genomic rearrangement, which is absent in N. intermedia. Finally, we compared the sequences of the three available N. crassa complete mtDNAs and found low levels of intraspecific variability. Most differences among strains were due to small indels in noncoding regions. The revisiting of the N. crassa mtDNA forms the basis for future studies on mitochondrial genome organization and variability.

A Comparative Analyses of the Complete Mitochondrial Genomes of Fungal Endosymbionts in Sogatella furcifera, White-Backed Planthoppers

Sogatella furcifera Horvath, commonly known as the white-backed planthoppers (WBPH), is an important pest in East Asian rice fields. Fungal endosymbiosis is widespread among planthoppers in the infraorder Fulgoromorpha and suborder Auchenorrhyncha. We successfully obtained complete mitogenome of five WBPH fungal endosymbionts, belonging to the Ophiocordycipitaceae family, from next-generation sequencing (NGS) reads obtained from S. furcifera samples. These five mitogenomes range in length from 55,390 bp to 55,406 bp, which is shorter than the mitogenome of the fungal endosymbiont found in Ricania speculum, black planthoppers. Twenty-eight protein-coding genes (PCGs), 12 tRNAs, and 2 rRNAs were found in the mitogenomes. Two single-nucleotide polymorphisms, two insertions, and three deletions were identified among the five mitogenomes, which were fewer in number than those of four species of Ophiocordycipitaceae, Ophiocordyceps sinensis, Hirsutella thompsonii, Hirsutella rhossiliensis, and Tolypocladium inflatum. Noticeably short lengths (up to 18 bp) of simple sequence repeats were identified in the five WBPH fungal endosymbiont mitogenomes. Phylogenetic analysis based on conserved PCGs across 25 Ophiocordycipitaceae mitogenomes revealed that the five mitogenomes were clustered with that of R. speculum, forming an independent clade. In addition to providing the full mitogenome sequences, obtaining complete mitogenomes of WBPH endosymbionts can provide insights into their phylogenetic positions without needing to isolate the mtDNA from the host. This advantage is of value to future studies involving fungal endosymbiont mitogenomes.

Characterization of the Complete Mitochondrial Genome of Basidiomycete Yeast Hannaella oryzae: Intron Evolution, Gene Rearrangement, and Its Phylogeny

In this study, the mitogenome of Hannaella oryzae was sequenced by next-generation sequencing (NGS) and successfully assembled. The H. oryzae mitogenome comprised circular DNA molecules with a total size of 26,444 bp. We found that the mitogenome of H. oryzae partially deleted the tRNA gene transferring cysteine. Comparative mitogenomic analyses showed that intronic regions were the main factors contributing to the size variations of mitogenomes in Tremellales. Introns of the cox1 gene in Tremellales species were found to have undergone intron loss/gain events, and introns of the H. oryzae cox1 gene may have different origins. Gene arrangement analysis revealed that H. oryzae contained a unique gene order different from other Tremellales species. Phylogenetic analysis based on a combined mitochondrial gene set resulted in identical and well-supported topologies, wherein H. oryzae was closely related to Tremella fuciformis. This study represents the first report of mitogenome for the Hannaella genus, which will allow further study of the population genetics, taxonomy, and evolutionary biology of this important phylloplane yeast and other related species.

Revealing the high variability on nonconserved core and mobile elements of Austropuccinia psidii and other rust mitochondrial genomes

Mitochondrial genomes are highly conserved in many fungal groups, and they can help characterize the phylogenetic relationships and evolutionary biology of plant pathogenic fungi. Rust fungi are among the most devastating diseases for economically important crops around the world. Here, we report the complete sequence and annotation of the mitochondrial genome of Austropuccinia psidii (syn. Puccinia psidii), the causal agent of myrtle rust. We performed a phylogenomic analysis including the complete mitochondrial sequences from other rust fungi. The genome composed of 93.299 bp has 73 predicted genes, 33 of which encoded nonconserved proteins (ncORFs), representing almost 45% of all predicted genes. A. psidii mtDNA is one of the largest rust mtDNA sequenced to date, most likely due to the abundance of ncORFs. Among them, 33% were within intronic regions of diverse intron groups. Mobile genetic elements invading intron sequences may have played significant roles in size but not shaping of the rust mitochondrial genome structure. The mtDNAs from rust fungi are highly syntenic. Phylogenetic inferences with 14 concatenated mitochondrial proteins encoded by the core genes placed A. psidii according to phylogenetic analysis based on 18S rDNA. Interestingly, cox1, the gene with the greatest number of introns, provided phylogenies not congruent with the core set. For the first time, we identified the proteins encoded by three A. psidii ncORFs using proteomics analyses. Also, the orf208 encoded a transmembrane protein repressed during in vitro morphogenesis. To the best of our knowledge, we presented the first report of a complete mtDNA sequence of a member of the family Sphaerophragmiacea.

Characterization of Two Mitochondrial Genomes and Gene Expression Analysis Reveal Clues for Variations, Evolution, and Large-Sclerotium Formation in Medical Fungus Wolfiporia cocos

Wolfiporia cocos, a precious mushroom with a long history as an edible food and Asian traditional medicine, remains unclear in the genetic mechanism underlying the formation of large sclerotia. Here, two complete circular mitogenomes (BL16, 135,686 bp and MD-104 SS10, 124,842 bp, respectively) were presented in detail first. The salient features in the mitogenomes of W. cocos include an intron in the tRNA (trnQ-UUG²), and an obvious gene rearrangement identified between the two mitogenomes from the widely geographically separated W. cocos strains. Genome comparison and phylogenetic analyses reveal some variations and evolutional characteristics in W. cocos. Whether the mitochondrion is functional in W. cocos sclerotium development was investigated by analyzing the mitogenome synteny of 10 sclerotium-forming fungi and mitochondrial gene expression patterns in different W. cocos sclerotium-developmental stages. Three common homologous genes identified across ten sclerotium-forming fungi were also found to exhibit significant differential expression levels during W. cocos sclerotium development. Most of the mitogenomic genes are not expressed in the mycelial stage but highly expressed in the sclerotium initial or developmental stage. These results indicate that some of mitochondrial genes may play a role in the development of sclerotium in W. cocos, which needs to be further elucidated in future studies. This study will stimulate new ideas on cytoplasmic inheritance of W. cocos and facilitate the research on the role of mitochondria in large sclerotium formation.

Complete Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma atroviride: Genomic Features, Comparative Analysis and Insight Into the Mitochondrial Evolution in Trichoderma

The improvement of biopesticides for use in the agriculture industry requires an understanding of the biological- and ecological principles underlying their behavior in natural environments. The nuclear genomes of members of the genus Trichoderma, which are representative fungal biocontrol agents, have been actively studied in relation to the unique characteristics of these species as effective producers of CAZymes/secondary metabolites and biopesticides, but their mitochondrial genomes have received much less attention. In this study, the mitochondrial genome of Trichoderma atroviride (Hypocreales, Sordariomycetes), which targets wood-decaying fungal pathogens and has the ability to degrade chemical fungicides, was assembled de novo. A 32,758 bp circular DNA molecule was revealed with specific features, such as a few more protein CDS and trn genes, two homing endonucleases (LAGLIDADG-/GIY-YIG-type), and even a putative overlapping tRNA gene, on a closer phylogenetic relationship with T. gamsii among hypocrealean fungi. Particularly, introns were observed with several footprints likely to be evolutionarily associated with the intron dynamics of the Trichoderma mitochondrial genomes. This study is the first to report the complete de novo mitochondrial genome of T. atroviride, while comparative analyses of Trichoderma mitochondrial genomes were also conducted from the perspective of mitochondrial evolution for the first time.

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.

Mitogenome types of two Lentinula edodes sensu lato populations in China

China has two populations of Lentinula edodes sensu lato as follows: L. edodes sensu stricto and an unexcavated morphological species respectively designated as A and B. In a previous study, we found that the nuclear types of the two populations are distinct and that both have two branches (A1, A2, B1 and B2) based on the internal transcribed spacer 2 (ITS2) sequence. In this paper, their mitogenome types were studied by resequencing 20 of the strains. The results show that the mitogenome type (mt) of ITS2-A1 was mt-A1, that of ITS2-A2 was mt-A2, and those of ITS2-B1 and ITS2-B2 were mt-B. The strains with heterozygous ITS2 types had one mitogenome type, and some strains possessed a recombinant mitogenome. This indicated that there may be frequent genetic exchanges between the two populations and both nuclear and mitochondrial markers were necessary to identify the strains of L. edodes sensu lato. In addition, by screening SNP diversity and comparing four complete mitogenomes among mt-A1, mt-A2 and mt-B, the cob, cox3, nad2, nad3, nad4, nad5, rps3 and rrnS genes could be used to identify mt-A and mt-B and that the cox1, nad1 and rrnL genes could be used to identify mt-A1, mt-A2 and mt-B.

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.

The complete mitochondrial genome of the Chan-hua fungus Isaria cicadae: a tale of intron evolution in Cordycipitaceae

Isaria cicadae is an entomogenous fungus of great medicinal value. Its nuclear genome has been reported, while its mitogenome remains unknown. Herein, we first described its mitogenome and then inferred intron evolution from both intraspecific and interspecific perspectives. The fungus represented the largest mitogenome (56.6 kb in strain CCAD02) known in Cordycipitaceae due to the presence of 25 introns interrupting nine genes. Comparison of three I. cicadae strains revealed intron presence/absence dynamics at six intron loci plus a few indels and single nucleotide polymorphisms. Phylogenetic analyses confirmed the placement of I. cicadae in Cordycipitaceae. Comparison of 10 Cordycipitaceae species revealed a high degree of synteny and conserved genetic content. They, however, varied in intron numbers (1-25 per species) with overall 34 intron loci identified, which resulted in more than twofold variations in mitogenome sizes (24.5-56.6 kb). An rnl intron encoding ribosomal protein S3 was present in all species, suggesting its early invasion in Cordycipitaceae, while further divergence occurred for this intron. The other introns identified in this study were present in some, but not all of the species and have undergone multiple gains and losses in Cordycipitaceae. This study greatly enhanced our understanding of intron evolution in Cordycipitaceae.

Mitochondrial genome of the entomophthoroid fungus Conidiobolus heterosporus provides insights into evolution of basal fungi

Entomophthoroid fungi represent an ecologically important group of fungal pathogens on insects. Here, the whole mitogenome of Conidiobolus heterosporus, one of the entomophthoroid fungi, was described and compared to those early branching fungi with available mitogenomes. The 53,364-bp circular mitogenome of C. heterosporus contained two rRNA genes, 14 standard protein-coding genes, 26 tRNA genes, and three free-standing ORFs. Thirty introns interrupted nine mitochondrial genes. Phylogenetic analysis based on mitochondrion-encoded proteins revealed that C. heterosporus was most close to Zancudomyces culisetae in the Zoopagomycota of basal fungi. Comparison on mitogenomes of 23 basal fungi revealed great variabilities in terms of mitogenome conformation (circular or linear), genetic code (codes 1, 4, or 16), AT contents (53.3-85.5%), etc. These mitogenomes varied from 12.0 to 97.3 kb in sizes, mainly due to different numbers of genes and introns. They showed frequent DNA rearrangement events and a high variability of gene order, although high synteny and conserved gene order were also present between closely related species. By reporting the first mitogenome in Entomophthoromycotina and the second in Zoopagomycota, this study greatly enhanced our understanding on evolution of basal fungi.

The first complete mitochondrial genome from the family Hygrophoraceae (Hygrophorus russula) by next-generation sequencing and phylogenetic implications

Hygrophorus russula (Schaeff.) Kauffman is an edible ectomycorrhizal fungus that is widely distributed in the world. In this study, the mitogenome of H. russula was sequenced and assembled. The mitogenome of H. russula is composed of circular DNA molecules, with a total size of 55,769 bp. Further analysis indicated that the frequent use of A and T in codons contributes to the high AT content (80.87%) in the H. russula mitogenome. Comparative analysis indicated that the length and base composition of the core protein-encoding genes, and the number of tRNA genes in the H. russula mitogenome varied from that of other Agaricales mitogenomes. Gene arrangement analysis revealed a novel gene order in the H. russula mitogenome. In addition, the expansion of the mitogenome in Agaricales was found to be closely related to the increase in the number of introns. Phylogenetic analysis of the combined mitochondrial gene set showed strong support for tree topologies, and H. russula was determined to be relatively distant from other Agaricales species. This study is the first report on the mitogenome of a member of genus Hygrophorus as well as family Hygrophoraceae, which improves our understanding of mitochondrial differentiation and evolution in the important ectomycorrhizal fungi Hygrophorus species.

Characterization and comparative mitogenomic analysis of six newly sequenced mitochondrial genomes from ectomycorrhizal fungi (Russula) and phylogenetic analysis of the Agaricomycetes

In this study, the mitochondrial genomes of six Russula species were sequenced using next generation sequencing. The six mitogenomes were all composed of circular DNA molecules, with lengths ranging from 40,961 bp to 69,423 bp. The length and number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among the six mitogenomes. The increased number and total size of introns likely contributed to the size expansion of mitogenomes in some Russula species. Gene synteny analysis revealed some gene rearrangements among the six mitochondrial genomes. The nad4L gene had the lowest K2P genetic distance of the 15 core PCGs among the six Russula species, indicating that this gene was highly conserved. The Ka/Ks values for all 15 core PCGs were <1, suggesting that they were all subject to purifying selection. Phylogenetic analyses based on two gene datasets (15 core PCGs, and 15 core PCGs + rnl + rns) recovered identical and well-supported trees. In addition, cox1 was identified as a potential single-gene molecular marker for the phylogenetic analysis of relationships among Agaricomycetes species. This study provides the first report of mitogenomes from the Russulaceae family and facilitates the investigation of population genetics and evolution of other ectomycorrhizal fungi.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

Comparative analysis of the mitochondrial genome of the fungus Colletotrichum lindemuthianum, the causal agent of anthracnose in common beans

Fungi of the genus Colletotrichum are economically important and are used as models in plant-pathogen interaction studies. In this study, the complete mitochondrial genomes of two Colletotrichum lindemuthianum isolates were sequenced and compared with the mitochondrial genomes of seven species of Colletotrichum. The mitochondrial genome of C. lindemuthianum is a typical circular molecule 37,446 bp (isolate 89 A₂ 2-3) and 37,440 bp (isolate 83.501) in length. The difference of six nucleotides between the two genomes is the result of a deletion in the ribosomal protein S3 (rps3) gene in the 83.501 isolate. In addition, substitution of adenine for guanine within the rps3 gene in the mitochondrial genome of the 83.501 isolate was observed. Compared to the previously sequenced C. lindemuthianum mitochondrial genome, an exon no annotated in the cytochrome c oxidase I (cox1) gene and a non-conserved open reading frame (ncORF) were observed. The size of the mitochondrial genomes of the seven species of Colletotrichum was highly variable, being attributed mainly to the ncORF, ranging from one to 10 and also from introns ranging from one to 11 and which encode a total of up to nine homing endonucleases. This paper reports for the first time by means of transcriptome that then ncORFs are transcribed in Colletotrichum spp. Phylogeny data revealed that core mitochondrial genes could be used as an alternative in phylogenetic relationship studies in Colletotrichum spp. This work contributes to the genetic and biological knowledge of Colletotrichum spp., which is of great economic and scientific importance.

Mitochondrial genome of the nematode endoparasitic fungus Hirsutella vermicola reveals a high level of synteny in the family Ophiocordycipitaceae

Ophiocordycipitaceae is a diverse fungal family comprising multiple ecologically, economically, medicinally, and culturally important fungal species; however, only four species of the family have available mitochondrial genomes (mitogenomes). In this study, the complete mitogenome of the nematode endoparasitic fungus Hirsutella vermicola in Ophiocordycipitaceae was sequenced, and a comparative mitogenomic analysis of Ophiocordycipitaceae was performed. We found that the 53,793-bp circular mitogenome of H. vermicola, except for standard fungal mitochondrial genes, harbors seven introns acquired possibly through lateral transfer from other fungi and three free-standing open reading frames (ORFs) coding for hypothetical proteins. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed its placement in Ophiocordycipitaceae. Comparison on five mitogenomes of Ophiocordycipitaceae revealed great variation on their sizes, from 35.2 kb in Tolypocladium ophioglossoides to 157.5 kb in Ophiocordyceps sinensis, mainly due to variable numbers of introns (from 7 to 54) as well as variable lengths of intergenic regions. The five mitogenomes, however, are highly syntenic to each other in terms of gene order, the presence of an intronic ORF encoding ribosomal protein S3 within rnl, and the nad2/nad3 joining pattern. Our study is the first report of the mitogenome of H. vermicola and has facilitated the understanding of mitogenome evolution of Ophiocordycipitaceae.