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

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

Introduction

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

Methods

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

Results and Discussion

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

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.

The effect of short-term fallowing on the microbial communities in forest soil cultivated with ginseng: Preliminary research

Background

Continuous cultivation of ginseng crops in fixed plots can lead to disease outbreaks, yield losses and replanting failures. Fallow periods can help restore soil health and increase the sustainability of agricultural systems; however, taking land out of production for extended periods is often not feasible. Short-term fallow periods could restore soil health, but few studies have examined the effects of short-term fallow treatment on the health of soil in ginseng fields.

Methods

In this preliminary study, we used metagenomic analysis to assess changes in the abundance of major ginseng pathogens and soil health overall following a short-term fallow period in a region in the Changbai Mountains. A sample from a forest plot (Hx0ks), was compared to a sample from a field where ginseng was previously cultivated and then had been left fallow for two years (Hx2), and a sample from a field that had been fallow for two years and was subsequently replanted with ginseng (Clsd).

Results

Soil that was fallow for two years, and then replanted with ginseng, showed reduced nutrient content and lower diversity of soil bacterial and fungal communities than soil that remained fallow. Candidatus Solibacter (5%) and Rhizomicrobium (3%) were the most abudant bacterial genera in Hx2. Rhizomicrobium (4%) and Gemmatimonas (3%) were the most abundant bacterial genera in Clsd. Mortierella (22%) and Peziza (12%) dominated the fungal community in Hx2. Lecanicillium (38%) and Mortierella (13%) dominated the fungal community in Clsd. Fallow periods also increased the functional diversity of soil as predicted by PICRUSt and decreased the relative abundance of the pathogenic fungi.

Conclusions

Preliminary findings were consistent with the hypothesis that fallow management in ginseng cultivation can improve soil microbial community structure and function and reduces the number of plant pathogens; however, testing this hypothesis will require replicated plots.

Comparative Mitogenomics of Two Sympatric Catfishes of Exostoma (Siluriformes: Sisoridae) from the Lower Yarlung Tsangpo River and Its Application for Phylogenetic Consideration

The genus Exostoma is a group of stenotopic and rheophilic glyptosternine catfishes distributed in South and Southeast Asia. So far, comprehensive studies on mitogenomics referring to this genus are very scarce. In this study, we first sequenced and annotated the complete mitochondrial genomes of Exostoma tibetanum and Exostoma tenuicaudatum—two sympatric congeners from the lower Yarlung Tsangpo River, Tibet, China. The mitogenomes of both species contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one light-strand origin of replication, and one control region, with lengths of 16,528 bp and 16,533 bp, respectively. The mitogenome architecture, nucleotide composition, and codon usage of protein-coding genes were almost identical between the two Exostoma species, although some estimated parameters varied. Phylogenetic analysis strongly supported the monophyly of Exostoma in the subfamily Glyptosternae, and Exostoma tibetanum had the closest relationship to Exostoma tenuicaudatum. The divergence time estimation demonstrated that these two species diverged approximately 1.51 Ma during the early Pleistocene, which was speculated to be triggered by the river system changes caused by the uplift of the southeastern Tibetan Plateau. Selection pressure analyses indicated that all protein-coding genes of Exostoma species underwent a strong purifying selection, while minority positive sites from NADH dehydrogenase complex genes were detected. These findings are expected to promote our understanding of the molecular phylogeny of the genus Exostoma and provide valuable mitogenomic resources for the subfamily Glyptosternae

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.

The mitochondrial genome of a wild edible mushroom, Russula rosea

Russula rosea is a common wild edible ectomycorrhizal fungus, which is widely distributed all over the world. We assembled the complete mitochondrial genome of R. rosea with the total length was 54177 bp and the GC content of 22.34%. It contains a total of 57 genes, including 14 standard protein-coding genes, one conserved ribosomal protein S3 gene (rps3), two rRNA genes, 24 tRNA genes, 15 free-standing open reading frames (ORFs) and one DNA polymerase gene (dpo). Mitochondrial genome found a close evolutionary relationship between Russula rosea and Russula lepida, which was helpful to study the genetic evolutionary relationship of edible fungi.

Genetic Differentiation and Widespread Mitochondrial Heteroplasmy among Geographic Populations of the Gourmet Mushroom Thelephora ganbajun from Yunnan, China

The mitochondrial genomes are generally considered non-recombining and homoplasmic in nature. However, our previous study provided the first evidence of extensive and stable mitochondrial heteroplasmy in natural populations of the basidiomycete fungus Thelephora ganbajun from Yunnan province, China. The heteroplasmy was characterized by the presence of two types of introns residing at adjacent but different sites in the cytochrome oxidase subunits I (cox1) gene within an individual strain. However, the frequencies of these two introns among isolates from different geographical populations and the implications for the genetic structure in natural populations have not been investigated. In this study, we analyzed DNA sequence variation at the internal transcribed spacer (ITS) regions of the nuclear ribosomal RNA gene cluster among 489 specimens from 30 geographic locations from Yunnan and compared that variation with distribution patterns of the two signature introns in the cox1 gene that are indicative of heteroplasmy in this species. In our samples, evidence for gene flow, abundant genetic diversity, and genotypic uniqueness among geographic samples in Yunnan were revealed by ITS sequence variation. While there was insignificant positive correlation between geographic distance and genetic differentiation among the geographic samples based on ITS sequences, a moderate significant correlation was found between ITS sequence variation, geographical distance of sampling sites, and distribution patterns of the two heteroplasmic introns in the cox1 gene. Interestingly, there was a significantly negative correlation between the copy numbers of the two co-existing introns. We discussed the implications of our results for a better understanding of the spread of stable mitochondrial heteroplasmy, mito-nuclear interactions, and conservation of this important gourmet mushroom.

First two mitochondrial genomes for the order Filobasidiales reveal novel gene rearrangements and intron dynamics of Tremellomycetes

In the present study, two mitogenomes from the Filobasidium genus were assembled and compared with other Tremellomycetes mitogenomes. The mitogenomes of F. wieringae and F. globisporum both comprised circular DNA molecules, with sizes of 27,861 bp and 71,783 bp, respectively. Comparative mitogenomic analysis revealed that the genetic contents, tRNAs, and codon usages of the two Filobasidium species differed greatly. The sizes of the two Filobasidium mitogenomes varied greatly with the introns being the main factor contributing to mitogenome expansion in F. globisporum. Positive selection was observed in several protein-coding genes (PCGs) in the Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina species, including cob, cox2, nad2, and rps3 genes. Frequent intron loss/gain events were detected to have occurred during the evolution of the Tremellomycetes mitogenomes, and the mitogenomes of 17 species from Agaricomycotina, Pucciniomycotina, and Ustilaginomycotina have undergone large-scale gene rearrangements. Phylogenetic analyses based on Bayesian inference and the maximum likelihood methods using a combined mitochondrial gene set generated identical and well-supported phylogenetic trees, wherein Filobasidium species had close relationships with Trichosporonales species. This study, which is the first report on mitogenomes from the order Filobasidiales, provides a basis for understanding the genomics, evolution, and taxonomy of this important fungal group.

Characterization of the complete mitochondrial genome of Dioszegia changbaiensis (Tremellales: Bulleribasidiaceae) with phylogenetic implications

In this study, the complete mitochondrial genome of Dioszegia changbaiensis we sequenced and assembled by the next-generation sequencing. The complete mitochondrial genome of Dioszegia changbaiensis contained 22 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. The total length of the Dioszegia changbaiensis mitochondrial genome is 34,853 bp, and the GC content of the mitochondrial genome is 41.88%. Phylogenetic analysis based on a combined mitochondrial gene dataset indicated that the mitochondrial genome of Dioszegia changbaiensis exhibited a close relationship with that of Hannaella oryzae.

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.

Genomic and Transcriptomic Insight of Giant Sclerotium Formation of Wood-Decay Fungi

Many fungi form persistent and dormant sclerotia with compact hardened mycelia during unfavorable circumstances. While most of these sclerotia are small in size, Wolfiporia cocos, a wood-decay fungus, grows into giant sclerotia, which are mainly composed of polysaccharides of linear (1→3)-β-D-glucans. To explore the underlying mechanism of converting sophisticated wood polysaccharides for biosynthesis of highly homogenized glucans in W. cocos, we sequenced and assembled the genome of a cultivated W. cocos strain (WCLT) in China. The 62-Mb haploid genome contains 44.2% repeat sequences, of which, 48.0% are transposable elements (TEs). Contrary to the genome of W. cocos from North America, WCLT has independently undergone a partial genome duplication (PGD) event. The large-scale TE insertion and PGD occurrence overlapped with an archeological Pleistocene stage of low oxygen and high temperature, and these stresses might have induced the differences in sclerotium due to geographical distribution. The wood decomposition enzymes, as well as sclerotium-regulator kinases, aquaporins, and highly expanded gene families such as NAD-related families, together with actively expressed 1,3-β-glucan synthase for sclerotium polysaccharides, all have contributed to the sclerotium formation and expansion. This study shall inspire further exploration on how fungi convert wood into simple glucans in the sclerotium of W. cocos.

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.

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 mitochondrial genome of the grape powdery mildew pathogen Erysiphe necator is intron rich and exhibits a distinct gene organization

Powdery mildews are notorious fungal plant pathogens but only limited information exists on their genomes. Here we present the mitochondrial genome of the grape powdery mildew fungus Erysiphe necator and a high-quality mitochondrial gene annotation generated through cloning and Sanger sequencing of full-length cDNA clones. The E. necator mitochondrial genome consists of a circular DNA sequence of 188,577 bp that harbors a core set of 14 protein-coding genes that are typically present in fungal mitochondrial genomes, along with genes encoding the small and large ribosomal subunits, a ribosomal protein S3, and 25 mitochondrial-encoded transfer RNAs (mt-tRNAs). Interestingly, it also exhibits a distinct gene organization with atypical bicistronic-like expression of the nad4L/nad5 and atp6/nad3 gene pairs, and contains a large number of 70 introns, making it one of the richest in introns mitochondrial genomes among fungi. Sixty-four intronic ORFs were also found, most of which encoded homing endonucleases of the LAGLIDADG or GIY-YIG families. Further comparative analysis of five E. necator isolates revealed 203 polymorphic sites, but only five were located within exons of the core mitochondrial genes. These results provide insights into the organization of mitochondrial genomes of powdery mildews and represent valuable resources for population genetic and evolutionary studies.

Evolutionary Insights Into Two Widespread Ectomycorrhizal Fungi (Pisolithus) From Comparative Analysis of Mitochondrial Genomes

The genus Pisolithus is a group of global ectomycorrhizal fungi. The characterizations of Pisolithus mitochondrial genomes have still been unknown. In the present study, the complete mitogenomes of two Pisolithus species, Pisolithus microcarpus, and Pisolithus tinctorius, were assembled and compared with other Boletales mitogenomes. Both Pisolithus mitogenomes comprised circular DNA molecules with sizes of 43,990 and 44,054 bp, respectively. Comparative mitogenomic analysis showed that the rps3 gene differentiated greatly between Boletales species, and this gene may be subjected to strong pressure of positive selection between some Boletales species. Several plasmid-derived genes and genes with unknown functions were detected in the two Pisolithus mitogenomes, which needs further analysis. The two Pisolithus species show a high degree of collinearity, which may represent the gene arrangement of the ancestors of ectomycorrhizal Boletales species. Frequent intron loss/gain events were detected in Boletales and basidiomycetes, and intron P717 was only detected in P. tinctorius out of the eight Boletales mitogenomes tested. We reconstructed phylogeny of 79 basidiomycetes based on combined mitochondrial gene dataset, and obtained well-supported phylogenetic topologies. This study served as the first report on the mitogenomes of the family Pisolithaceae, which will promote the understanding of the evolution of Pisolithus species.

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.

Comparative Analysis of Mitochondrial Genome Features among Four Clonostachys Species and Insight into Their Systematic Positions in the Order Hypocreales

The mycoparasite fungi of Clonostachys have contributed to the biological control of plant fungal disease and nematodes. The Clonostachys fungi strains were isolated from Ophiocordyceps highlandensis, Ophiocordycepsnigrolla and soil, which identified as Clonostachyscompactiuscula, Clonostachysrogersoniana, Clonostachyssolani and Clonostachys sp. To explore the evolutionary relationship between the mentioned species, the mitochondrial genomes of four Clonostachys species were sequenced and assembled. The four mitogenomes consisted of complete circular DNA molecules, with the total sizes ranging from 27,410 bp to 42,075 bp. The GC contents, GC skews and AT skews of the mitogenomes varied considerably. Mitogenomic synteny analysis indicated that these mitogenomes underwent gene rearrangements. Among the 15 protein-coding genes within the mitogenomes, the nad4L gene exhibited the least genetic distance, demonstrating a high degree of conservation. The selection pressure analysis of these 15 PCGs were all below 1, indicating that PCGs were subject to purifying selection. Based on protein-coding gene calculation of the significantly supported topologies, the four Clonostachys species were divided into a group in the phylogenetic tree. The results supplemented the database of mitogenomes in Hypocreales order, which might be a useful research tool to conduct a phylogenetic analysis of Clonostachys. Additionally, the suitable molecular marker was significant to study phylogenetic relationships in the Bionectriaceae family.

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.

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.

The 287,403 bp Mitochondrial Genome of Ectomycorrhizal Fungus Tuber calosporum Reveals Intron Expansion, tRNA Loss, and Gene Rearrangement

In the present study, the mitogenome of Tuber calosporum was assembled and analyzed. The mitogenome of T. calosporum comprises 15 conserved protein-coding genes, two rRNA genes, and 14 tRNAs, with a total size of 287,403 bp. Fifty-eight introns with 170 intronic open reading frames were detected in the T. calosporum mitogenome. The intronic region occupied 69.41% of the T. calosporum mitogenome, which contributed to the T. calosporum mitogenome significantly expand relative to most fungal species. Comparative mitogenomic analysis revealed large-scale gene rearrangements occurred in the mitogenome of T. calosporum, involving gene relocations and position exchanges. The mitogenome of T. calosporum was found to have lost several tRNA genes encoding for cysteine, aspartate, histidine, etc. In addition, a pair of fragments with a total length of 32.91 kb in both the nuclear and mitochondrial genomes of T. calosporum was detected, indicating possible gene transfer events. A total of 12.83% intragenomic duplications were detected in the T. calosporum mitogenome. Phylogenetic analysis based on mitochondrial gene datasets obtained well-supported tree topologies, indicating that mitochondrial genes could be reliable molecular markers for phylogenetic analyses of Ascomycota. This study served as the first report on mitogenome in the family Tuberaceae, thereby laying the groundwork for our understanding of the evolution, phylogeny, and population genetics of these important ectomycorrhizal fungi.

Mitogenomes of Two Phallus Mushroom Species Reveal Gene Rearrangement, Intron Dynamics, and Basidiomycete Phylogeny

Phallus indusiatus and Phallus echinovolvatus are edible bamboo mushrooms with pharmacological properties. We sequenced, assembled, annotated, and compared the mitogenomes of these species. Both mitogenomes were composed of circular DNA molecules, with sizes of 89,139 and 50,098 bp, respectively. Introns were the most important factor in mitogenome size variation within the genus Phallus. Phallus indusiatus, P. echinovolvatus, and Turbinellus floccosus in the subclass Phallomycetidae have conservative gene arrangements. Large-scale gene rearrangements were observed in species representing 42 different genera of Basidiomycetes. A variety of intron position classes were found in the 44 Basidiomycete species analyzed. A novel group II intron from the P. indusiatus mitogenome was compared with other fungus species containing the same intron, and we demonstrated that the insertion sites of the intron had a base preference. Phylogenetic analyses based on combined gene datasets yielded well-supported Bayesian posterior probability (BPP = 1) topologies. This indicated that mitochondrial genes are reliable molecular markers for analyzing the phylogenetic relationships of the Basidiomycetes. This is the first study of the mitogenome of the genus Phallus, and it increases our understanding of the population genetics and evolution of bamboo mushrooms and related species.

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

In the present study, the complete mitogenome of Clavaria fumosa, was sequenced, assembled, and compared. The complete mitogenome of C. fumosa is 256,807 bp in length and is the largest mitogenomes among all Basidiomycota mitogenomes reported. Comparative mitogenomic analysis indicated that the C. fumosa mitogenome contained the most introns and intronic ORFs among all fungal mitogenomes. Large intergenic regions, intronic regions, accumulation of repeat sequences and plasmid-derived genes together promoted the size expansion of the C. fumosa mitogenome. In addition, the rps3 gene was found subjected to positive selection between some Agaricales species. We found frequent intron gain/loss events in Agaricales mitogenomes, and four novel intron classes were detected in the C. fumosa mitogenome. Large-scale gene rearrangements were found occurred in Agaricales species and the C. fumosa mitogenome had a unique gene arrangement which differed from other Agaricales species. Phylogenetic analysis for 76 Basidiomycetes based on combined mitochondrial gene sets indicated that mitochondrial genes could be used as effective molecular markers for reconstructing evolution of Basidiomycota. The study served as the first report on the mitogenomes of the family Clavariaceae, which will promote the understanding of the genetics, evolution and taxonomy of C. fumosa and related species.

The complete mitochondrial genome of medicinal fungus Taiwanofungus camphoratus reveals gene rearrangements and intron dynamics of Polyporales

Taiwanofungus camphoratus is a highly valued medicinal mushroom that is endemic to Taiwan, China. In the present study, the mitogenome of T. camphoratus was assembled and compared with other published Polyporales mitogenomes. The T. camphoratus mitogenome was composed of circular DNA molecules, with a total size of 114,922 bp. Genome collinearity analysis revealed large-scale gene rearrangements between the mitogenomes of Polyporales, and T. camphoratus contained a unique gene order. The number and classes of introns were highly variable in 12 Polyporales species we examined, which proved that numerous intron loss or gain events occurred in the evolution of Polyporales. The Ka/Ks values for most core protein coding genes in Polyporales species were less than 1, indicating that these genes were subject to purifying selection. However, the rps3 gene was found under positive or relaxed selection between some Polyporales species. Phylogenetic analysis based on the combined mitochondrial gene set obtained a well-supported topology, and T. camphoratus was identified as a sister species to Laetiporus sulphureus. This study served as the first report on the mitogenome in the Taiwanofungus genus, which will provide a basis for understanding the phylogeny and evolution of this important fungus.

Characterization of the complete mitochondrial genome of Sterigmatomyces hyphaenes (Agaricostilbales: Agaricostilbaceae) and implications for its phylogeny

In this study, the complete mitochondrial genome of Sterigmatomyces hyphaenes was sequenced by the next-generation sequencing. The complete mitochondrial genome of S. hyphaenes contained 17 protein-coding genes (PCG), 2 ribosomal RNA (rRNA) genes, and 23 transfer RNA (tRNA) genes. The total size of the S. hyphaenes mitochondrial genome is 26,198 bp, and the GC content of the mitochondrial genome is 42.08%. Phylogenetic analysis based on the combined mitochondrial gene dataset indicated that the mitochondrial genome of S. hyphaenes exhibited a close relationship with that of Rhodotorula mucilaginosa.

The First Mitochondrial Genome for Geastrales (Sphaerobolus stellatus) Reveals Intron Dynamics and Large-Scale Gene Rearrangements of Basidiomycota

In this study, the mitogenome of artillery fungus, Sphaerobolus stellatus, was assembled and compared with other Basidiomycota mitogenomes. The Sphaerobolus stellatus mitogenome was composed of circular DNA molecules, with a total size of 152,722 bp. Accumulation of intergenic and intronic sequences contributed to the Sphaerobolus stellatus mitogenome becoming the fourth largest mitogenome among Basidiomycota. We detected large-scale gene rearrangements in Basidiomycota mitogenomes, and the Sphaerobolus stellatus mitogenome contains a unique gene order. The quantity and position classes of intron varied between 75 Basidiomycota species we tested, indicating frequent intron loss/gain events occurred in the evolution of Basidiomycota. A novel intron position classes (P1281) was detected in the Sphaerobolus stellatus mitogenome, without any homologous introns from other Basidiomycota species. A pair of fragments with a total length of 9.12 kb in both the nuclear and mitochondrial genomes of Sphaerobolus stellatus was detected, indicating possible gene transferring events. Phylogenetic analysis based on the combined mitochondrial gene set obtained well-supported tree topologies (Bayesian posterior probabilities ≥ 0.99; bootstrap values ≥98). This study served as the first report on the mitogenome from the order Geastrales, which will promote the understanding of the phylogeny, population genetics, and evolution of the artillery fungus, Sphaerobolus stellatus.

Characterization of the complete mitochondrial genome of biocontrol yeast Sporobolomyces sp. (Sporidiobolales: Sporidiobolaceae) with phylogenetic analysis

In this study, we obtained the complete mitochondrial genome of Sporobolomyces sp. using next-generation sequencing. The complete mitochondrial genome of Sporobolomyces sp. contained 15 protein-coding genes (PCG), two ribosomal RNA (rRNA) genes, and 25 transfer RNA (tRNA) genes. The total length of the Sporobolomyces sp. mitochondrial genome is 26,430 bp, and the GC content of the mitochondrial genome is 39.32%. Phylogenetic analysis based on combined mitochondrial gene dataset indicated that the mitochondrial genome of Sporobolomyces sp. exhibited a close relationship with that of Rhodotorula mucilaginosa.

Characterization of the complete mitochondrial genome of Meira sp. (Exobasidiales: Brachybasidiaceae) with phylogenetic analysis

In this study, the complete mitochondrial genome of Meira sp. was sequenced and assembled. The complete mitochondrial genome of Meira sp. has 15 protein-coding (PCG) genes, 2 ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. The mitochondrial genome of Meira sp. has a total size of 23,353 bp, with the base composition as follows: A (30.62%), T (32.82%), G (17.84%) and C (18.73%). Phylogenetic analysis indicated that the mitogenome of Meira sp. exhibited a close relationship with the mitogenome of Jaminaea angkorensis.

Comparative mitogenome analysis of two ectomycorrhizal fungi (Paxillus) reveals gene rearrangement, intron dynamics, and phylogeny of basidiomycetes

In this study, the mitogenomes of two Paxillus species were assembled, annotated and compared. The two mitogenomes of Paxillus involutus and P. rubicundulus comprised circular DNA molecules, with the size of 39,109 bp and 41,061 bp, respectively. Evolutionary analysis revealed that the nad4L gene had undergone strong positive selection in the two Paxillus species. In addition, 10.64 and 36.50% of the repetitive sequences were detected in the mitogenomes of P. involutus and P. rubicundulus, respectively, which might transfer between mitochondrial and nuclear genomes. Large-scale gene rearrangements and frequent intron gain/loss events were detected in 61 basidiomycete species, which revealed large variations in mitochondrial organization and size in Basidiomycota. In addition, the insertion sites of the basidiomycete introns were found to have a base preference. Phylogenetic analysis of the combined mitochondrial gene set gave identical and well-supported tree topologies, indicating that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships of Basidiomycota. This study is the first report on the mitogenomes of Paxillus, which will promote a better understanding of their contrasted ecological strategies, molecular evolution and phylogeny of these important ectomycorrhizal fungi and related basidiomycete species.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Clavulina sp. (Cantharellales: Clavulinaceae)

In this study, the complete mitochondrial genome of Clavulina sp. was sequenced and assembled. The complete mitochondrial genome of Clavulina sp. contains 20 protein-coding (PCG) genes, 2 ribosomal RNA (rRNA) genes, and 25 transfer RNA (tRNA) genes. The total size of the Clavulina sp. complete mitochondrial genome is 31,816 bp, with the GC content of 27.72%. Phylogenetic analysis indicated that the mitochondrial genome of Clavulina sp. exhibited a close relationship with that of the genus Cantharellus.

Complete Mitogenomes of Three Carangidae (Perciformes) Fishes: Genome Description and Phylogenetic Considerations

Carangidae are ecologically and economically important marine fish. The complete mitogenomes of three Carangidae species (Alectis indicus, Decapterus tabl, and Alepes djedaba) were sequenced, characterized, and compared with 29 other species of the family Carangidae in this study. The length of the three mitogenomes ranged from 16,530 to 16,610 bp, and the structures included 2 rRNA genes (12S rRNA and 16S rRNA), 1 control region (a non-coding region), 13 protein-coding genes, and 22 tRNA genes. Among the 22 tRNA genes, only tRNA-Ser (GCT) was not folded into a typical cloverleaf secondary structure and had no recognizable DHU stem. The full-length sequences and protein-coding genes (PCGs) of the mitogenomes of the three species all had obvious AT biases. The majority of the AT-skew and GC-skew values of the PCGs among the three species were negative, demonstrating bases T and C were more plentiful than A and G. Analyses of Ka/Ks and overall p-genetic distance demonstrated that ATP8 showed the highest evolutionary rate and COXI/COXII were the most conserved genes in the three species. The phylogenetic tree based on PCGs sequences of mitogenomes using maximum likelihood and Bayesian inference analyses showed that three clades were divided corresponding to the subfamilies Caranginae, Naucratinae, and Trachinotinae. The monophyly of each superfamily was generally well supported. The divergence time analyses showed that Carangidae evolved during three geological periods, the Cretaceous, Paleogene, and Neogene. A. indicus began to differentiate from other species about 27.20 million years ago (Mya) in the early Miocene, while D. tabl (21.25 Mya) and A. djedaba (14.67 Mya) differentiated in the middle Oligocene.

Mitochondrial genome characterization and phylogenetic analysis of bird schistosome Trichobilharzia szidati

In the present study, the complete mitochondrial genome of T. szidati was assembled by next generation sequencing (NGS). We found that the complete mitochondrial genome of T.szidati is 14, 303 bp in length and consists of 3023 (21.1%) adenine, 1153 (8.1%) cytosine, 3432 (24.0%) guanosine and 6695 (46.8%) thymine. The genome contains 12 conserved core protein-coding genes (atp6, cox1, cox2, cox3, nad1, nad2, nad3, nad4, nad4L, nad5, nad6, CYTB), 21 tRNA genes, 2 rRNA genes and 1 D-loop region. Phylogenetic analysis showed that T. szidati has a close relationship with T. regent. Knowledge of mitochondrial genome of T. szidati could provide useful information for the further studies of evolutionary biology, epidemiology and species identification.

Comparative Mitogenome Analysis Reveals Mitochondrial Genome Differentiation in Ectomycorrhizal and Asymbiotic Amanita Species

In this present study, we assembled and analyzed the mitogenomes of two asymbiotic and six ectomycorrhizal Amanita species based on next-generation sequencing data. The size of the eight Amanita mitogenomes ranged from 37,341 to 137,428 bp, and we considered introns to be one of the main factors contributing to the size variation of Amanita. The introns of the cox1 gene experienced frequent gain/loss events in Amanita; and the intron position class cox1P386 was lost in the six ectomycorrhizal Amanita species. In addition, ectomycorrhizal Amanita species had more repetitive sequences and fewer intergenic sequences than asymbiotic Amanita species in their mitogenomes. Large-scale gene rearrangements were detected in the Amanita species we tested, including gene displacements and inversions. On the basis of the combined mitochondrial gene set, we reconstructed the phylogenetic relationships of 66 Basidiomycetes. The six ectomycorrhizal Amanita species were of single origin, and the two saprophytic Amanita species formed two distinct clades. This study is the first to elucidate the functions of the mitogenome in the evolution and ecological adaptation of Amanita species.

Characterization and phylogenetic analysis of the complete mitochondrial genome of a basidiomycetous yeast Cystobasidium sp. (Cystobasidiales: Cystobasidiaceae)

In the present study, the complete mitochondrial genome of a basidiomycetous yeast Cystobasidium sp. was assembled and obtained. The mitochondrial genome of Cystobasidium sp. contains 16 protein-coding genes, 2 ribosomal RNA genes (rRNA), and 24 transfer RNA (tRNA) genes. The complete mitogenome of Cystobasidium sp. has a total length of 24,914 bp, with the base composition as follows: A (30.82), T (32.88%), C (18.37%) and G (17.93%). The Cystobasidium sp. mitogenome exhibited a close relationship with the mitogenome of Microbotryum cf. violaceum, M. lychnidis-dioicae, and Rhodotorula mucilaginosa.

Population Genomic Analysis Reveals a Highly Conserved Mitochondrial Genome in Fusarium asiaticum

Fusarium asiaticum is one of the pivotal members of the Fusarium graminearum species complex (FGSC) causing Fusarium head blight (FHB) on wheat, barley and rice in large parts of Asia. Besides resulting in yield losses, FHB also causes the accumulation of mycotoxins such as nivalenol (NIV) and deoxynivalenol (DON). The aim of this study was to conduct population studies on F. asiaticum from Southern China through mitochondrial genome analyses. All strains were isolated from wheat or rice from several geographic areas in seven provinces in Southern China. Based on geographic location and host, 210 isolates were selected for next generation sequencing, and their mitogenomes were assembled by GRAbB and annotated to explore the mitochondrial genome variability of F. asiaticum. The F. asiaticum mitogenome proves extremely conserved and variation is mainly caused by absence/presence of introns harboring homing endonuclease genes. These variations could be utilized to develop molecular markers for track and trace of migrations within and between populations. This study illustrates how mitochondrial introns can be used as markers for population genetic analysis. SNP analysis demonstrate the occurrence of mitochondrial recombination in F. asiaticum as was previously found for F. oxysporum and implied for F. graminearum. Furthermore, varying degrees of genetic diversity and recombination showed a high association with different geographic regions as well as with cropping systems. The mitogenome of F. graminearum showed a much higher SNP diversity while the interspecies intron variation showed no evidence of gene flow between the two closely related and sexual compatible species.

Comparative Analyses of Mitochondrial Genomes Provide Evolutionary Insights Into Nematode-Trapping Fungi

Predatory fungi in Orbiliaceae (Ascomycota) have evolved a diversity of trapping devices that enable them to trap and kill nematodes, other small animals, and protozoans. These trapping devices include adhesive hyphae, adhesive knobs, adhesive networks, constricting rings, and non-constricting rings. Their diversity and practical importance have attracted significant attention from biologists, making them excellent model organisms for studying adaptative evolution and as biological control agents against parasitic nematodes. The putative origins and evolutionary relationships among these carnivorous fungi have been investigated using nuclear protein-encoding genes, but their patterns of mitogenome relationships and divergences remain unknown. Here we analyze and compare the mitogenomes of 12 fungal strains belonging to eight species, including six species representing all four types of nematode trapping devices and two from related but non-predatory fungi. All 12 analyzed mitogenomes were of circular DNA molecules, with lengths ranging from 146,101 bp to 280,699 bp. Gene synteny analysis revealed several gene rearrangements and intron transfers among the mitogenomes. In addition, the number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among these mitogenomes. The increased number and total size of introns were the main contributors to the length differences among the mitogenomes. Phylogenetic analyses of the protein-coding genes indicated that mitochondrial and nuclear genomes evolved at different rates, and signals of positive selection were found in several genes involved in energy metabolism. Our study provides novel insights into the evolution of nematode-trapping fungi and shall facilitate further investigations of this ecologically and agriculturally important group of fungi.

Mitochondrial genome characterization and phylogenetic analysis of Blastocladiella sp. (Blastocladiales: Blastocladiaceae)

In the present study, the complete mitochondrial genome of an early diverging fungus Blastocladiella sp. was assembled by the next-generation sequencing. The complete mitochondrial genome of Blastocladiella sp. is 33, 800 bp in length and consists of 11,620 (34.38%) adenine, 5,047 (14.93%) cytosine, 6,025 (17.83%) guanosine and 11,108 (32.86%) thymine. The genome contains 19 protein-coding genes, 24 tRNA genes and 2 rRNA genes. Phylogenetic analysis based on the combined mitochondrial gene set showed that Blastocladiella sp. has a close relationship with Allomyces macrogynus and Blastocladiella emersonii.

Mitochondrial genome characterization and phylogenetic analysis of arbuscular mycorrhizal fungus Rhizophagus sp

In the present study, the complete mitochondrial genome of Rhizophagus sp. was assembled by the next-generation sequencing. We found that the complete mitochondrial genome of Rhizophagus sp. is 50,449 bp in length and consists of 14,741 (29.22%) adenine, 9427 (18.69%) cytosine, 9248 (18.33%) guanosine, and 17,033 (33.76%) thymine. The genome contains 24 conserved core protein-coding genes, 25 tRNA genes, and 2 rRNA genes. Phylogenetic analysis based on the combined mitochondrial gene set showed that Rhizophagus sp. has a close relationship with Rhizophagus fasciculatus, Glomus irregular, and G. intratadices.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Conidiobolus sp. (Entomophthorales: Ancylistaceae)

In the present study, we presented the complete mitochondrial genome of an entomophthoroid fungus Conidiobolus sp. The mitogenome of Conidiobolus sp. has a total length of 26,612 bp, with the base composition as follows: A (44.22%), T (27.10%), C (10.99%) and G (17.68%). The mitogenome contains 19 protein-coding genes, 2 ribosomal RNA genes (rRNA), and 23 transfer RNA (tRNA) genes. The taxonomic status of the Conidiobolus sp. mitogenome exhibited a close relationship with the mitogenome of Conidiobolus heterosporus.

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.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the human pathogenic fungus Cryptococcus sp. (Tremellales: Cryptococcaceae)

In the present study, the complete mitogenome of Cryptococcus sp. were sequenced and assembled. The complete mitogenome of Cryptococcus sp. was composed of circular DNA molecules, with a total length of 30,029 bp. The base composition of this mitochondrial genome is as follows: A (31.94%), T (34.89%), G (15.97%), and C (17.21%). The mitogenome contains 20 protein-coding genes, 2 ribosomal RNA genes (rRNA), and 22 transfer RNA (tRNA) genes. Phylogenetic analysis showed that the mitogenome of the Cryptococcus sp. exhibited a closest relationship with Cryptococcus gattii.

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

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

Comparative Mitochondrial Genome Analysis of Two Ectomycorrhizal Fungi (Rhizopogon) Reveals Dynamic Changes of Intron and Phylogenetic Relationships of the Subphylum Agaricomycotina

In the present study, we assembled and compared two mitogenomes from the Rhizopogon genus. The two mitogenomes of R. salebrosus and R. vinicolor comprised circular DNA molecules, with the sizes of 66,704 bp and 77,109 bp, respectively. Comparative mitogenome analysis indicated that the length and base composition of protein coding genes (PCGs), rRNA genes and tRNA genes varied between the two species. Large fragments aligned between the mitochondrial and nuclear genomes of both R. salebrosus (43.41 kb) and R. vinicolor (12.83 kb) indicated that genetic transfer between mitochondrial and nuclear genomes has occurred over evolutionary time of Rhizopogon species. Intronic regions were found to be the main factors contributing to mitogenome expansion in R. vinicolor. Variations in the number and type of introns in the two mitogenomes indicated that frequent intron loss/gain events occurred during the evolution of Rhizopogon species. Phylogenetic analyses based on Bayesian inference (BI) and Maximum likelihood (ML) methods using a combined mitochondrial gene set yielded identical and well-supported tree topologies, wherein Rhizopogon species showed close relationships with Agaricales species. This is the first study of mitogenomes within the genus Rhizopogon, and it provides a basis for understanding the evolution and differentiation of mitogenomes from the ectomycorrhizal fungal genus.

The complete mitochondrial genome of a wild edible mushroom, Russula griseocarnosa

Russula griseocarnosa is a wild edible ectomycorrhizal mushroom in southern China. In this study, we assembled the complete mitochondrial genome of R. griseocarnosa. Its total length was 60995 bp with a GC content of 21% and contained a total of 52 genes, including 14 standard protein-coding genes, two rRNA genes, 21 tRNA genes and 15 free-standing open reading frames (ORFs). Phylogenetic analysis reflected that the evolutionary processes between R. griseocarnosa and some agaricomycetes.

RETRACTED ARTICLE: Characterization of the complete mitochondrial genome of an important edible fungus Auricularia Polytricha

We, the Editor, authors and Publisher of Mitochondrial DNA Part B: Resources, have retracted the following article:

Qiang Li, Xin Jin, Zuqin Chen, Chuan Xiong, Ping Li & Wenli Huang (2019) ‘Characterization of the complete mitochondrial genome of an important edible fungus Auricularia Polytricha’ 4:1, 1107–1108, DOI https://doi.org/10.1080/23802359.2019.1586494.

This article has been retracted at the request of the authors as, subsequent to publication, it has been identified that the species reported in the article is incorrect. The species should be Suillus sp. and not Auricularia polytricha as reported. The article has therefore been retracted as the data reported in the article are inaccurate.

We have been informed in our decision-making by our policy on publishing ethics and integrity and the COPE guidelines on retractions.

The retracted article will remain online to maintain the scholarly record, but it will be digitally watermarked on each page as “Retracted”.