The 256 kb mitochondrial genome of Clavaria fumosa is the largest among phylum Basidiomycota and is rich in introns and intronic ORFs

Trichoderma koningiopsis (Hypocreaceae) has the smallest mitogenome of the genus Trichoderma

Introduction

Fungal mitogenomes exhibit remarkable variation in conformation, size, gene content, arrangement and expression, including their intergenic spacers and introns.

Methods

The complete mitochondrial genome sequence of the mycoparasitic fungus Trichoderma koningiopsis was determined using the Illumina next-generation sequencing technology. We used data from our recent Illumina NGS-based project of T. koningiopsis genome sequencing to study its mitochondrial genome. The mitogenome was assembled, annotated, and compared with other fungal mitogenomes.

Results

T. koningiopsis strain POS7 mitogenome is a circular molecule of 27,560 bp long with a GC content of 27.80%. It harbors the whole complement of the 14 conserved mitochondrial protein-coding genes (PCG) such as atp6, atp8, atp9, cox1, cox2, cox3, cob, nad1, nad2, nad3, nad4, nad4L, nad5, and nad6, also found in the same gene order to other Hypocreales. The mitogenome also contains 26 transfer RNA genes (tRNAs), 5 of them with more than one copy. Other genes also present in the assembled mitochondrial genome are a small rRNA subunit and a large rRNA subunit containing ribosomal protein S3 gene. Despite the small genome size, two introns were detected in the T. koningiopsis POS7 mitogenome, one of them in cox3 gene and the other in rnl gene, accounting 7.34% of this mitogenome with a total size of 2,024 bp. A phylogenetic analysis was done using the 14 PCGs genes of T. koningiopsis strain POS7 mitogenome to compare them with those from other fungi of the Subphyla Pezizomycotina and Saccharomycotina. T. koningiopsis strain POS7 was clustered together with other representatives of Trichoderma lineage, within the Hypocreales group, which is also supported by previous phylogenetic studies based on nuclear markers.

Discussion

The mitochondrial genome of T. koningiopsis POS7 will allow further investigations into the taxonomy, phylogenetics, conservation genetics, and evolutionary biology of this important genus as well as other closely related species.

The First Whole Genome Sequencing of Agaricus bitorquis and Its Metabolite Profiling

Agaricus bitorquis, an emerging wild mushroom with remarkable biological activities and a distinctive oversized mushroom shape, has gained increasing attention in recent years. Despite its status as an important resource of wild edible fungi, knowledge about this mushroom is still limited. In this study, we used the Illumina NovaSeq and Nanopore PromethION platforms to sequence, de novo assemble, and annotate the whole genome and mitochondrial genome (mitogenome) of the A. bitorquis strain BH01 isolated from Bosten Lake, Xinjiang Province, China. Using the genome-based biological information, we identified candidate genes associated with mating type and carbohydrate-active enzymes in A. bitorquis. Cluster analysis based on P450 of basidiomycetes revealed the types of P450 members of A. bitorquis. Comparative genomic, mitogenomic, and phylogenetic analyses were also performed, revealing interspecific differences and evolutionary features of A. bitorquis and A. bisporus. In addition, the molecular network of metabolites was investigated, highlighting differences in the chemical composition and content of the fruiting bodies of A. bitorquis and A. bisporus. The genome sequencing provides a comprehensive understanding and knowledge of A. bitorquis and the genus Agaricus mushrooms. This work provides valuable insights into the potential for artificial cultivation and molecular breeding of A. bitorquis, which will facilitate the development of A. bitorquis in the field of edible mushrooms and functional food manufacture.

The first two mitochondrial genomes from Apiotrichum reveal mitochondrial evolution and different taxonomic assignment of Trichosporonales

Apiotrichum is a diverse anamorphic basidiomycetous yeast genus, and its mitogenome characterization has not been revealed. In this study, we assembled two Apiotrichum mitogenomes and compared them with mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. The mitogenomes of Apiotrichum gracile and A. gamsii comprised circular DNA molecules, with sizes of 34,648 bp and 38,096 bp, respectively. Intronic regions were found contributed the most to the size expansion of A. gamsii mitogenome. Comparative mitogenomic analysis revealed that 6.85-38.89% of nucleotides varied between tRNAs shared by the two Apiotrichum mitogenomes. The GC content of all core PCGs in A. gamsii was lower than that of A. gracile, with an average low value of 4.97%. The rps3 gene differentiated the most among Agaricomycotina, Pucciniomycotina and Ustilaginomycotina species, while nad4L gene was the most conserved in evolution. The Ka/Ks values for cob and rps3 genes were > 1, indicating the two genes may be subjected to positive selection in Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. Frequent intron loss/gain events and potential intron transfer events have been detected in evolution of Agaricomycotina, Pucciniomycotina and Ustilaginomycotina. We further detected large-scale gene rearrangements between the 19 mitogenomes from Agaricomycotina, Pucciniomycotina and Ustilaginomycotina, and fifteen of the 17 mitochondrial genes shared by Apiotrichum varied in gene arrangements. Phylogenetic analyses based on maximum likelihood and Bayesian inference methods using a combined mitochondrial gene dataset revealed different taxonomic assignment of two Apiotrichum species, wherein A. gamsii had a more closely relationship with Trichosporon asahii. This study served as the first report on mitogenomes from the genus Apiotrichum, which promotes the understanding of evolution, genomics, and phylogeny of Apiotrichum.

Evolution and codon usage bias of mitochondrial and nuclear genomes in Aspergillus section Flavi

The fungal genus Aspergillus contains a diversity of species divided into taxonomic sections of closely related species. Section Flavi contains 33 species, many of industrial, agricultural, or medical relevance. Here, we analyze the mitochondrial genomes (mitogenomes) of 20 Flavi species-including 18 newly assembled mitogenomes-and compare their evolutionary history and codon usage bias patterns to their nuclear counterparts. Codon usage bias refers to variable frequencies of synonymous codons in coding DNA and is shaped by a balance of neutral processes and natural selection. All mitogenomes were circular DNA molecules with highly conserved gene content and order. As expected, genomic content, including GC content, and genome size differed greatly between mitochondrial and nuclear genomes. Phylogenetic analysis based on 14 concatenated mitochondrial genes predicted evolutionary relationships largely consistent with those predicted by a phylogeny constructed from 2,422 nuclear genes. Comparing similarities in interspecies patterns of codon usage bias between mitochondrial and nuclear genomes showed that species grouped differently by patterns of codon usage bias depending on whether analyses were performed using mitochondrial or nuclear relative synonymous usage values. We found that patterns of codon usage bias at gene level are more similar between mitogenomes of different species than the mitogenome and nuclear genome of the same species. Finally, we inferred that, although most genes-both nuclear and mitochondrial-deviated from the neutral expectation for codon usage, mitogenomes were not under translational selection while nuclear genomes were under moderate translational selection. These results contribute to the study of mitochondrial genome evolution in filamentous fungi.

Comparative analyses of Flammulina filiformis mitochondrial genomes reveal high length polymorphism in intergenic regions and multiple intron gain/loss in cox1

The golden-needle mushroom Flammulina filiformis is one of the bulk mushroom products in the world. This study obtained complete mitogenomes of 44 wild isolates collected from nine provinces and two artificially bred cultivars of F. filiformis, together with three Flammulina rossica isolates and one Flammulina fennae isolate for comparison. The mitogenome of F. filiformis ranged from 83,540 bp to 90,938 bp, consisting of 14 conserved protein-coding genes (PCGs), two rRNA genes, and 25 tRNA genes. To the best of our knowledge, it contained the highest proportion of intergenic regions compared to the other known Basidiomycota mitogenomes. Introns and intergenic regions were two major contributing factors to the total size of the F. filiformis mitogenome. The conserved PCG cox3 is located in an intron of another conserved PCG, nad5. This is a unique phenomenon in all known fungal mitogenomes. Gain/loss of introns was observed in cox1, nad5, and rnl. Length polymorphism was widely observed in intergenic regions. Accordingly, primers were designed as useful markers for rapid identification of F. filiformis isolates with differentiated mitogenomes. Our findings provide a basis for further studies related to variety identification and population genetics of this economically important mushroom.

The first two mitochondrial genomes for the genus Ramaria reveal mitochondrial genome evolution of Ramaria and phylogeny of Basidiomycota

In the present study, we assembled and analyzed the mitogenomes of two Ramaria species. The assembled mitogenomes of Ramaria cfr. rubripermanens and R. rubella were circularized, with sizes of 126,497 bp and 143,271 bp, respectively. Comparative mitogenome analysis showed that intron region contributed the most (contribution rate, 43.74%) to the size variations of Ramaria mitogenomes. The genetic contents, gene length, tRNAs, and codon usages of the two Ramaria mitogenomes varied greatly. In addition, the evolutionary rates of different core protein coding genes (PCGs) in Phallomycetidae mitogenomes varied. We detected large-scale gene rearrangements between Phallomycetidae mitogenomes, including gene displacement and tRNA doubling. A total of 4499 bp and 7746 bp aligned fragments were detected between the mitochondrial and nuclear genomes of R. cfr. rubripermanens and R. rubella, respectively, indicating possible gene transferring events. We further found frequent intron loss/gain and potential intron transfer events in Phallomycetidae mitogenomes during the evolution, and the mitogenomes of R. rubella contained a novel intron P44. Phylogenetic analyses using both Bayesian inference (BI) and Maximum Likelihood (ML) methods based on a combined mitochondrial gene dataset obtained an identical and well-supported phylogenetic tree for Basidiomycota, wherein R. cfr. rubripermanens and Turbinellus floccosus are sister species. This study served as the first report on mitogenomes from the genus Ramaria, which provides a basis for understanding the evolution, genetics, and taxonomy of this important fungal group.

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.

Comparative mitochondrial genome analyses reveal conserved gene arrangement but massive expansion/contraction in two closely related Exserohilum pathogens

Exserohilum turcicum and E. rostratum, two closely related fungal species, are both economically important pathogens but have quite different target hosts (specific to plants and cross-kingdom infection, respectively). In the present study, complete circular mitochondrial genomes of the two Exserohilum species were sequenced and de novo assembled, which mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). Comparative analyses indicated that these two fungi had significant mitogenomic collinearity and consistent mitochondrial gene arrangement, yet with vastly different mitogenome sizes, 264,948 bp and 64,620 bp, respectively. By contrast with the 17 introns containing 17 intronic ORFs (one-to-one) in the E. rostratum mitogenome, E. turcicum involved far more introns (70) and intronic ORFs (126), which was considered as the main contributing factors of their mitogenome expansion/contraction. Within the generally intron-rich gene cox1, a total of 18 and 10 intron position classes (Pcls) were identified separately in the two mitogenomes. Moreover, 16.16% and 10.85% ratios of intra-mitogenomic repetitive regions were detected in E. turcicum and E. rostratum, respectively. Based on the combined mitochondrial gene dataset, we established a well-supported topology of phylogeny tree of 98 ascomycetes, implying that mitogenomes may act as an effective molecular marker for fungal phylogenetic reconstruction. Our results served as the first report on mitogenomes in the genus Exserohilum, and would have significant implications in understanding the origin, evolution and pathogenic mechanisms of this fungal lineage.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the pathogenic fungus Ilyonectria destructans

Ilyonectria destructans is a pathogenic fungus causing root rot and other symptoms on trees and many crops. This paper analyses the mitochondrial genome of I. destructans and compares it with other published Nectriaceae mitogenomes. The I. destructans mitogenome appears as a circular DNA molecule of 42,895 bp and an overall GC content of 28.23%. It contains 28 protein-coding genes (15 core protein genes and 13 free-standing ORFs), two rRNAs and 27 tRNAs. The gene content and order were found to be conserved in the mitogenome of I. destructans and other Nectriaceae, although the genome size varies because of the variation in the number and length of intergenic regions and introns. For most core protein-coding genes in Nectriaceae species, Ka/Ks < 1 indicates purifying selection. Among some Nectriaceae representatives, only the rps3 gene was found under positive selection. Phylogenetic analyses based on nucleotide sequences of 15 protein-coding genes divided 45 Hypocreales species into six major clades matching the families Bionectriaceae, Cordycipitaceae, Clavicipitaceae, Ophiocordycipitaceae, Hypocreaceae and Nectriaceae. I. destructans appeared as a sister species to unidentified Ilyonectia sp., closely related to C. ilicicola, N. cinnabarina and a clad of ten Fusarium species and G. moniliformis. The complete mitogenome of I. destructans reported in the current paper will facilitate the study of epidemiology, biology, genetic diversity of the species and the evolution of family Nectriace and the Hypocreales order.

Characterization of the complete mitochondrial genome of Aspergillus terricola (Aspergillaceae, Eurotiales), isolated from soy sauce fermentation system

In the present study, the complete mitochondrial genome of Aspergillus terricola É.J. Marchal 1893 was sequenced and assembled. The mitochondrial genome of A. terricola was composed of circular DNA molecules, with a total size of 28,689 bp. The GC content of the A. terricola mitochondrial genome was 26.34%. A total of 18 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, and 26 transfer RNA (tRNA) genes were detected in the A. terricola mitochondrial genome. Phylogenetic analysis based on the combined mitochondrial gene dataset indicated that the A. terricola exhibited a close relationship with A. parasiticus.

The Gene Rearrangement, Loss, Transfer, and Deep Intronic Variation in Mitochondrial Genomes of Conidiobolus

The genus Conidiobolus s.s. was newly delimited from Conidiobolus s.l. In order to gain insight into its mitochondrial genetic background, this study sequenced six mitochondrial genomes of the genus Conidiobolus s.s. These mitogenomes were all composed of circular DNA molecules, ranging from 29,253 to 48,417 bp in size and from 26.61 to 27.90% in GC content. The order and direction for 14 core protein-coding genes (PCGs) were identical, except for the atp8 gene lost in Conidiobolus chlamydosporus, Conidiobolus polyspermus, and Conidiobolus polytocus, and rearranged in the other Conidiobolus s.s. species. Besides, the atp8 gene split the cox1 gene in Conidiobolus taihushanensis. Phylogenomic analysis based on the 14 core PCGs confirmed that all Conidiobolus s.s. species formed a monophyly in the Entomophthoromycotina lineage. The number and length of introns were the main factors contributing to mitogenomic size, and deep variations and potential transfer were detected in introns. In addition, gene transfer occurred between the mitochondrial and nuclear genomes. This study promoted the understanding of the evolution and phylogeny of the Conidiobolus s.s. genus.

Characterization of the first complete chloroplast genome of Amaranthus hybridus (Caryophyllales: Amaranthaceae) with phylogenetic implications

In the present study, the complete chloroplast genome of Amaranthus hybridus was sequenced and assembled. The complete chloroplast genome of Amaranthus hybridus is 150,709 in size, with the GC content of 36.56%. The chloroplast genome of Amaranthus hybridus contained 86 protein-coding genes (PCGs), eight ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. Phylogenetic analysis based on combined chloroplast gene dataset indicated that the Amaranthus hybridus exhibited a close relationship with A. hypochondriacus and A. caudatus.

The complete mitochondrial genome of Ophiocordyceps gracilis and its comparison with related species

In this study, the complete mitochondrial genome of O. gracilis was sequenced and assembled before being compared with related species. As the second largest mitogenome reported in the family Ophiocordycipitaceae, the mitogenome of O. gracilis (voucher OG201301) is a circular DNA molecule of 134,288 bp that contains numerous introns and longer intergenomic regions. UCA was detected as anticodon in tRNA-Sec of O. gracilis, while comparative mitogenome analysis of nine Ophiocordycipitaceae fungi indicated that the order and contents of PCGs and rRNA genes were considerably conserved and could descend from a common ancestor in Ophiocordycipitaceae. In addition, the expansion of mitochondrial organization, introns, gene length, and order of O. gracilis were determined to be similar to those of O. sinensis, which indicated common mechanisms underlying adaptive evolution in O. gracilis and O. sinensis. Based on the mitochondrial gene dataset (15 PCGs and 2 RNA genes), a close genetic relationship between O. gracilis and O. sinensis was revealed through phylogenetic analysis. This study is the first to investigate the molecular evolution, phylogenetic pattern, and genetic structure characteristics of mitogenome in O. gracilis. Based on the obtained results, the mitogenome of O. gracilis can increase understanding of the genetic diversity and evolution of cordycipitoid fungi.

The first complete mitochondrial genome of mushroom Leucoagaricus naucinus (Agaricaceae, Agaricales) and insights into its phylogeny

Leucoagaricus naucinus (Fr.) Singer is a mycorrhizal fungus widely distributed in the northern Hemisphere. In the present study, the complete mitochondrial genome of Leucoagaricus naucinus was sequenced, assembled, and annotated. The L. naucinus mitochondrial genome was composed of circular DNA molecules, with the total size of 61,434 bp. The GC content of the L. naucinus mitochondrial genome was 26.07%. A total of 30 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, and 26 transfer RNA (tRNA) genes were detected in the L. naucinus mitochondrial genome. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the L. naucinus exhibited a close relationship with Agaricus bisporus.

Characterization of the complete mitochondrial genome of Peniophora lycii (Russulales: Peniophoraceae) with its phylogenetic analysis

Peniophora lycii is a resupinate lichen-like species distributed all over the world. In the present study, we sequenced and assembled the complete mitochondrial genome of Peniophora lycii. The size of the mitochondrial genome of P. lycii was 38,296 bp, with a GC content of 25.89%. Twenty protein-coding genes, 2 ribosomal RNA genes, and 24 transfer RNA genes were identified in the mitochondrial genome of P. lycii. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the mitochondrial genome of P. lycii exhibited a close relationship with that of Heterobasidion irregulare.

The first complete mitochondrial genome of edible and medicinal fungus Chroogomphus rutilus (Gomphidiaceae, Boletales) and insights into its phylogeny

In the present study, we assembled and annotated the complete mitochondrial genome of Chroogomphus rutilus. The complete mitochondrial genome of C. rutilus was composed of circular DNA molecules, with a size of 37,508 bp. The GC content of the C. rutilus mitogenome was 22.82%. A total of 18 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, and 24 transfer RNA (tRNA) genes were detected in the C. rutilus mitogenome. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the C. rutilus exhibited a close relationship with species from the genus Rhizopogon. This study served as the first report on the complete mitochondrial genome from the family Gomphidiaceae, which will promote the understanding of phylogeny, evolution, and taxonomy of this important fungal species.

The first complete mitochondrial genome of Macalpinomyces bursus (Ustilaginales: Ustilaginaceae) and insights into its phylogeny

In the present study, the complete mitochondrial genome of Macalpinomyces bursus (Berk.) Vanky 2002 was sequenced and assembled. The complete mitochondrial genome of M. bursus was 49,024 bp in length, with the GC content of 30.4%. The M. bursus mitochondrial genome contained 27 protein-coding genes, 2 ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the M. bursus exhibited a close relationship with species from the genera Ustilago, Sporisorium, and Anthracocystis.

The first complete chloroplast genome of Fagopyrum leptopodum (Diels) Hedberg (Caryophyllales: Polygonaceae) with phylogenetic implications

In the present study, we sequenced and assembled the complete chloroplast genome of Fagopyrum leptopodum (Diels) Hedberg. The chloroplast genome of F. leptopodum was composed of 85 protein-coding genes, 8 ribosomal RNA genes, and 37 transfer RNA genes. The F. leptopodum chloroplast genome is 159,375 bp in length, with a GC content of 37.81%. Phylogenetic analysis based on the combined chloroplast gene dataset indicated that the F. leptopodum exhibited a close relationship with Fagopyrum luojishanense.

Characterization of the complete mitochondrial genome of Coprinellus micaceus, a wild saprobic mushroom in China

The Glistening Inkcap (Coprinellus micaceus) is a wild saprobic mushroom in China. In this study, we assembled and annotated its complete mitochondrial genome using raw data sequenced through Illumina NovaSeq 6000 platform (Illumina, San Diego, CA). The length of the C. micaceus mitochondrial genome is 65,450 bp with 33.05% GC content. Totally, 41 genes, including 14 protein-coding genes, 25 tRNAs and 2 rRNAs were identified in the mitochondrial genome. Phylogenetic analysis showed that the mitochondrial genome relationship between C. micaceus and Coprinopsis cinerea was the closest.

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.

Characterization and phylogenetic analysis of the complete mitochondrial genome of pathogen Trichosporon inkin (Trichosporonales: Trichosporonaceae)

In the present study, the complete mitochondrial genome of Trichosporon inkin was sequenced and assembled. The complete mitochondrial genome of T. inkin contained 22 protein-coding genes (PCG), 2 ribosomal RNA (rRNA) genes, and 24 transfer RNA (tRNA) genes. The total size of the T. inkin mitochondrial genome is 39,466 bp, with the GC content of 27.56%. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the T. inkin exhibited a close relationship with Trichosporon asahii.

Comparative Mitogenomic Analysis Reveals Dynamics of Intron Within and Between Tricholoma Species and Phylogeny of Basidiomycota

The genus of Tricholoma is a group of important ectomycorrhizal fungi. The overlapping of morphological characteristics often leads to the confusion of Tricholoma species classification. In this study, the mitogenomes of five Tricholoma species were sequenced based on the next-generation sequencing technology, including T. matsutake SCYJ1, T. bakamatsutake, T. terreum, T. flavovirens, and T. saponaceum. These five mitogenomes were all composed of circular DNA molecules, with sizes ranging from 49,480 to 103,090 bp. Intergenic sequences were considered to be the main factor contributing to size variations of Tricholoma mitogenomes. Comparative mitogenomic analysis showed that the introns of the Agaricales mitogenome experienced frequent loss/gain events. In addition, potential gene transfer was detected between the mitochondrial and nuclear genomes of the five species of Tricholoma. Evolutionary analysis showed that the rps3 gene of the Tricholoma species was under positive selection or relaxed selection in the evolutionary process. In addition, large-scale gene rearrangements were detected between some Tricholoma species. Phylogenetic analysis using the Bayesian inference and maximum likelihood methods based on a combined mitochondrial gene set yielded identical and well-supported tree topologies. This study promoted the understanding of the genetics, evolution, and phylogeny of the Tricholoma genus and related species.

Panorama of intron dynamics and gene rearrangements in the phylum Basidiomycota as revealed by the complete mitochondrial genome of Turbinellus floccosus

In the present study, the complete mitogenome of Turbinellus floccosus was sequenced, assembled, and compared with other basidiomycete mitogenomes. The mitogenome of T. floccosus consists of a circular DNA molecule, with a size of 62,846 bp. Gene arrangement analysis indicated that large-scale gene rearrangements occurred in the levels of family and genus of basidiomycete species, and the mitogenome of T. floccosus contained a unique gene order. A significant correlation between the number of introns and the mitochondrial genome size of Basidiomycota were detected (P < 0.01). A total of 896 introns were detected in the core protein-coding genes (PCGs) of 74 basidiomycete species, and the cox1 gene was the largest host gene of basidiomycete introns. Intron position class (Pcls) P383 in the cox1 gene was the most common intron in Basidiomycota, which distributed in 40 of 74 basidiomycete species. In addition, frequent intron loss/gain events were detected in basidiomycete species. More than 50% of bases around insertion sites (- 15 bp to 15 bp) of Pcls from different species were conservative, indicating site preferences of intron insertions in Basidiomycota. Further analysis showed that 76.09% of introns tended to insert downstream to a T base in Basidiomycota. Phylogenetic analysis for 74 basidiomycetes indicated mitochondrial genes are effective molecular markers for phylogeny of basidiomycetes. The study served as the first report on the mitogenome from the family Gomphaceae, which will help to understand the intron origin and evolution in Basidiomycota. KEY POINTS: • The mitogenome of Turbinellus floccosus had a unique gene arrangement. • Intron loss/gain events were detected in the 74 basidiomycete species. • Introns tend to insert downstream of a T base in basidiomycete mitogenomes.

Comparative mitochondrial genome analysis reveals intron dynamics and gene rearrangements in two Trametes species

Trametes species are efficient wood decomposers that are widespread throughout the world. Mitogenomes have been widely used to understand the phylogeny and evolution of fungi. Up to now, two mitogenomes from the Trametes genus have been revealed. In the present study, the complete mitogenomes of two novel Trametes species, Trametes versicolor and T. coccinea, were assembled and compared with other Polyporales mitogenomes. Both species contained circular DNA molecules, with sizes of 67,318 bp and 99,976 bp, respectively. Comparative mitogenomic analysis indicated that the gene number, length and base composition varied between the four Trametes mitogenomes we tested. In addition, all of the core protein coding genes in Trametes species were identified and subjected to purifying selection. The mitogenome of T. coccinea contained the largest number of introns among the four Trametes species tested, and introns were considered the main factors contributing to size variations of Polyporales. Several novel introns were detected in the Trametes species we assembled, and introns identified in Polyporales were found to undergo frequent loss/gain events. Large-scale gene rearrangements were detected between closely related Trametes species, including gene inversions, insertions, and migrations. A well-supported phylogenetic tree for 77 Basidiomycetes was obtained based on the combined mitochondrial gene set using 2 phylogenetic inference methods. The results showed that mitochondrial genes are effective molecular markers for understanding the phylogeny of Basidiomycetes. This study is the first to report the mitogenome rearrangement and intron dynamics of Trametes species, which shed light on the evolution of Trametes and other related species.

The first eleven mitochondrial genomes from the ectomycorrhizal fungal genus (Boletus) reveal intron loss and gene rearrangement

In the present study, eleven novel complete mitogenomes of Boletus were assembled and compared. The eleven complete mitogenomes were all composed of circular DNA molecules, with sizes ranging from 32,883 bp to 48,298 bp. The mitochondrial gene arrangement of Boletus varied greatly from other Boletales mitogenomes, and gene position reversal were observed frequently in the evolution of Boletus. Across the 15 core protein-coding genes (PCGs) tested, atp9 had the least and rps3 had the largest genetic distances among the eleven Boletus species, indicating varied evolution rates of core PCGs. In addition, the Ka/Ks value for nad3 gene was >1, suggesting that this gene was subject to possible positive selection pressure. Comparative mitogenomic analysis indicated that the intronic region was significantly correlated with the size of mitogenomes in Boletales. Two large-scale intron loss events were detected in the evolution of Boletus. Phylogenetic analyses based on a combined mitochondrial gene dataset yielded a well-supported (BPP ≥ 0.99; BS =100) phylogenetic tree for 72 Agaricomycetes, and the Boletus species had a close relationship with Paxillus. This study served as the first report on complete mitogenomes in Boletus, which will further promote investigations of the genetics, evolution and phylogeny of the Boletus genus.

The mitogenomes of two saprophytic Boletales species (Coniophora) reveals intron dynamics and accumulation of plasmid-derived and non-conserved genes

The order Boletales is a group of fungi with complex life styles, which include saprophytic and ectomycorrhizal mushroom-forming fungi. In the present study, the complete mitogenomes of two saprophytic Boletales species, Coniophora olivacea, and C. puteana, were assembled and compared with mitogenomes of ectomycorrhizal Boletales. Both mitogenomes comprised circular DNA molecules with sizes of 78,350 bp and 79,655 bp, respectively. Comparative mitogenomic analysis indicated that the two saprophytic Boletales species contained more plasmid-derived (7 on average) and unknown functional genes (12 on average) than the four ectomycorrhizal Boletales species previously reported. In addition, the core protein coding genes, nad2 and rps3, were found to be subjected to positive selection pressure between some Boletales species. Frequent intron gain/loss events were detected in Boletales and Basidiomycetes, and several novel intron classes were found in two Coniophora species. A total of 33 introns were detected in C. olivacea, and most were found to have undergone contraction in the C. olivacea mitogenome. Mitochondrial genes of Coniophora species were found to have undergone large-scale gene rearrangements, and the accumulation of intra-genomic repeats in the mitogenome was considered as one of the main contributing factors. Based on combined mitochondrial gene sets, we obtained a well-supported phylogenetic tree for 76 Basidiomycetes, demonstrating the utility of mitochondrial gene analysis for inferring Basidiomycetes phylogeny. The study served as the first report on the mitogenomes of the family Coniophorineae, which will help to understand the origin and evolution patterns of Boletales species with complex lifestyles.