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

The Mitogenomic Characterization and Phylogenetic Analysis of the Plant Pathogen Phyllosticta yuccae

Phyllosticta yuccae is an important plant pathogen causing leaf spot disease in Yucca gigantea Lem. It is imperative to note that the amount of information available about the mitogenome of this subject is severely limited. This must be addressed immediately, as it is crucial to our understanding and progress in this field. To better understand the mitogenomic characteristics of P. yuccae, we conducted its sequencing by MGISEQ. Afterwards, the mitogenome was assembled and annotated. The mitogenomic characteristics and phylogenetic placement of the P. yuccae strain KUMCC 6213 were analyzed. The study revealed that the mitogenome of P. yuccae is a circular DNA molecule, consisting of 178,540 base pairs. It contains a total of 64 genes, including 14 protein-coding genes (PCGs), 26 transfer RNA genes (tRNA), 2 ribosomal RNA genes (rRNA), and 22 open reading frame genes (ORF), accounting for 80.98% of the total size. Repetitive sequences accounted for 15.42% of the mitogenome. The analysis of codon usage indicated that the codon UUA was the most commonly utilized, whereas the amino acid Leu was the most frequently employed. A comparative analysis of mitogenomes between P. yuccae and Macrophomina phaseolina (Tassi) Goid. showed notable variations in the position and size of gene clusters, with cox1, nad4, and nad4L genes exhibiting relatively low conservation. Phylogenetic analysis based on the 14 PCGs revealed that P. yuccae has the closest genetic relationship with M. phaseolina (Botryosphaeriaceae, Botryosphaeriales). This study first reports the mitogenome of P. yuccae and validates its phylogenetic placement. The findings enhance the knowledge of mitogenomes in Botryosphaeriales, offering novel perspectives on the genetics and evolution of the plant pathogen P. yuccae. This is crucial for the accurate prevention and management of leaf spot disease in Y. gigantea.

The First Two Complete Mitochondrial Genomes for the Subfamily Meligethinae (Coleoptera: Nitidulidae) and Implications for the Higher Phylogeny of Nitidulidae

The phylogenetic status of the family Nitidulidae and its sister group relationship remain controversial. Also, the status of the subfamily Meligethinae is not fully understood, and previous studies have been mainly based on morphology, molecular fragments, and biological habits, rather than the analysis of the complete mitochondrial genome. Up to now, there has been no complete mitochondrial genome report of Meligethinae. In this study, the complete mitochondrial genomes of Meligethinus tschungseni and Brassicogethes affinis (both from China) were provided, and they were compared with the existing complete mitochondrial genomes of Nitidulidae. The phylogenetic analysis among 20 species of Coleoptera was reconstructed via PhyloBayes analysis and Maximum likelihood (ML) analysis, respectively. The results showed that the full lengths of Meligethinus tschungseni and Brassicogethes affinis were 15,783 bp and 16,622 bp, and the AT contents were 77% and 76.7%, respectively. Each complete mitochondrial genome contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and a control region (A + T-rich region). All the PCGs begin with the standard start codon ATN (ATA, ATT, ATG, ATC). All the PCGs terminate with a complete terminal codon, TAA or TAG, except cox1, cox2, nad4, and nad5, which terminate with a single T. Furthermore, all the tRNAs have a typical clover-leaf secondary structure except trnS1, whose DHU arm is missing in both species. The two newly sequenced species have different numbers and lengths of tandem repeat regions in their control regions. Based on the genetic distance and Ka/Ks analysis, nad6 showed a higher variability and faster evolutionary rate. Based on the available complete mitochondrial genomes, the results showed that the four subfamilies (Nitidulinae, Meligethinae, Carpophilinae, Epuraeinae) of Nitidulidae formed a monophyletic group and further supported the sister group relationship of Nitidulidae + Kateretidae. In addition, the taxonomic status of Meligethinae and the sister group relationship between Meligethinae and Nitidulinae (the latter as currently circumscribed) were also preliminarily explored.

First complete mitogenome of Massarineae and its contribution to phylogenetic implications in Pleosporales

Endophytic fungi play an important role in the growth and development of traditional Chinese medicine plants. We isolated a strain of Acrocalymma vagum from the endophytic fungi of the traditional Chinese plants Paris. To accurately identify this endophytic fungal species of interest, we sequenced the mitochondrial genome of A. vagum, which is the first discovered mitochondrial genome in Massarineae. The A. vagum mitochondrial genome consists of a 35,079-bp closed circular DNA molecule containing 36 genes. Then, we compared the general sequence characteristics of A. vagum with those of Pleosporales, and the second structure of the 22 tRNAs was predicted. The phylogenetic relationship of A. vagum was constructed using two different data sets (protein-coding genes and amino acids). The phylogenetic tree shows that A. vagum is located at the root of Pleosporales. The analysis of introns shows that the number of introns increases with the increase in branch length. The results showed that monophyly was confirmed for all families in Pleosporales except for Pleosporaceae. A. vagum is an ancient species in the Pleosporales, and Pleosporaceae may require further revision. In Pleosporales, the number of introns is positively correlated with branch length, providing data for further study on the origin of introns.

The First Five Mitochondrial Genomes for the Family Nidulariaceae Reveal Novel Gene Rearrangements, Intron Dynamics, and Phylogeny of Agaricales

The family Nidulariaceae, consisting of five genera including Cyathus, is a unique group of mushrooms commonly referred to as bird's nest fungi due to their striking resemblance to bird's nests. These mushrooms are considered medicinal mushrooms in Chinese medicine and have received attention in recent years for their anti-neurodegenerative properties. However, despite the interest in these mushrooms, very little is known about their mitochondrial genomes (mitogenomes). This study is the first comprehensive investigation of the mitogenomes of five Nidulariaceae species with circular genome structures ranging in size from 114,236 bp to 129,263 bp. Comparative analyses based on gene content, gene length, tRNA, and codon usage indicate convergence within the family Nidulariaceae and heterogeneity within the order Agaricales. Phylogenetic analysis based on a combined mitochondrial conserved protein dataset provides a well-supported phylogenetic tree for the Basidiomycetes, which clearly demonstrates the evolutionary relationships between Nidulariaceae and other members of Agaricales. Furthermore, phylogenetic inferences based on four different gene sets reveal the stability and proximity of evolutionary relationships within Agaricales. These results reveal the uniqueness of the family Nidulariaceae and its similarity to other members of Agaricales; provide valuable insights into the origin, evolution, and genetics of Nidulariaceae species; and enrich the fungal mitogenome resource. This study will help to expand the knowledge and understanding of the mitogenomes in mushrooms.

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.

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 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.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes

Introduction

Codon usage bias is a prevalent phenomenon observed across various species and genes. However, the specific attributes of codon usage in the mitochondrial genome of Ganoderma species remain unknown.

Methods

In this study, we investigated the codon bias of 12 mitochondrial core protein-coding genes (PCGs) in 9 Ganoderma species, including 13 Ganoderma strains.

Results

The codons of all Ganoderma strains showed a preference for ending in A/T. Additionally, correlations between codon base composition and the codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) were identified, demonstrating the impact of base composition on codon bias. Various base bias indicators were found to vary between or within Ganoderma strains, including GC3s, the CAI, the CBI, and the FOP. The results also revealed that the mitochondrial core PCGs of Ganoderma have an average effective number of codons (ENC) lower than 35, indicating strong bias toward certain codons. Evidence from neutrality plot and PR2-bias plot analysis indicates that natural selection is a major factor affecting codon bias in Ganoderma. Additionally, 11 to 22 optimal codons (ΔRSCU>0.08 and RSCU>1) were identified in 13 Ganoderma strains, with GCA, AUC, and UUC being the most widely used optimal codons in Ganoderma. By analyzing the combined mitochondrial sequences and relative synonymous codon usage (RSCU) values, the genetic relationships between or within Ganoderma strains were determined, indicating variations between them. Nevertheless, RSCU-based analysis illustrated the intra- and interspecies relationships of certain Ganoderma species.

Discussion

This study deepens our insight into the synonymous codon usage characteristics, genetics, and evolution of this important fungal group.

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.

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.

The complete mitochondrial genome of the edible mushroom Grifola frondosa

The culinary-medicinal mushroom Grifola frondosa is widely cultivated in East Asia. In this study, the complete mitochondrial genome of G. frondosa was determined using Illumina sequencing. The circular molecule was 197,486 bp in length with a content of 25.01% GC, which was one of the largest mitochondrial genomes in the order Polyporales. A total of 39 known genes encoding 13 common mitochondrial genes, 24 tRNA genes, 1 ribosomal protein s3 gene (rps3), and 1 DNA polymerase gene (dpo) were predicted in this genome. The phylogenetic analysis showed that G. frondosa clustered together with Sparassis crispa, Laetiporus sulphureus, Wolfiporia cocos, and Taiwanofungus camphoratus. The complete mitochondrial genome reported here may provide new insight into genetic information and evolution for further studies.

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 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.

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.

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 chloroplast genome of Hordeum vulgare L. var. trifurcatum with phylogenetic analysis

In the present study, the complete chloroplast genome of Hordeum vulgare L. var. trifurcatum was sequenced, assembled and compared with closely related species. The chloroplast genome of Hordeum vulgare L. var. trifurcatum was composed of 84 protein-coding genes (PCG), 8 ribosomal RNA (rRNA) genes, and 38 transfer RNA (tRNA) genes. The Hordeum vulgare L. var. trifurcatum chloroplast genome is 136,485 bp in size, with the GC content of 38.32%. Phylogenetic analysis based on the combined chloroplast gene dataset indicated that the Hordeum vulgare L. var. trifurcatum exhibited a close relationship with Hordeum vulgare subsp. spontaneum and Hordeum vulgare subsp. vulgare.

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.

Novel gene rearrangement in the mitochondrial genome of Siliqua minima (Bivalvia, Adapedonta) and phylogenetic implications for Imparidentia

Siliqua minima (Gmelin, 1791) is an important economic shellfish species belonging to the family Pharidae. To date, the complete mitochondrial genome of only one species in this family (Sinonovacula constricta) has been sequenced. Research on the Pharidae family is very limited; to improve the evolution of this bivalve family, we sequenced the complete mitochondrial genome of S. minima by next-generation sequencing. The genome is 17,064 bp in length, consisting of 12 protein-coding genes (PCGs), 22 transfer RNA genes (tRNA), and two ribosomal RNA genes (rRNA). From the rearrangement analysis of bivalves, we found that the gene sequences of bivalves greatly variable among species, and with closer genetic relationship, the more consistent of the gene arrangement is higher among the species. Moreover, according to the gene arrangement of seven species from Adapedonta, we found that gene rearrangement among families is particularly obvious, while the gene order within families is relatively conservative. The phylogenetic analysis between species of the superorder Imparidentia using 12 conserved PCGs. The S. minima mitogenome was provided and will improve the phylogenetic resolution of Pharidae 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.

Comparative Mitogenomic Analysis of Two Cuckoo Bees (Apoidea: Anthophila: Megachilidae) with Phylogenetic Implications

Bees (Hymenoptera, Apoidea and Anthophila) are distributed worldwide and considered the primary pollinators of angiosperm. Megachilidae is one of the largest families of Anthophila. In this study, two complete mitogenomes of cuckoo bees in Megachilidae, namely Coelioxys fenestrata and Euaspis polynesia, were amplified and sequenced, with a length of 17,004 bp (C. fenestrata) and 17,682 bp (E. polynesia). The obtained results show that 37 mitogenomic genes and one putative control region were conserved within Hymenoptera. Truncated stop codon T was found in the cox3 gene of E. polynesia. The secondary structure of small (rrnS) and large (rrnL) rRNA subunits contained three domains (28 helices) and five domains (44 helices) conserved within Hymenoptera, respectively. Compared with ancestral gene order, gene rearrangement events included local inversion and gene shuffling. In order to reveal the phylogenetic position of cuckoo bees, we performed phylogenetic analysis. The results supported that all families of Anthophila were monophyletic, the tribe-level relationship of Megachilidae was Osmiini + (Anthidiini + Megachilini) and Coelioxys fenestrata was clustered to the Megachile genus, which was more closely related to Megachile sculpturalis and Megachile strupigera than Euaspis polynesia.

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.

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.

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.

Characterization of Nuclear and Mitochondrial Genomes of Two Tobacco Endophytic Fungi Leptosphaerulina chartarum and Curvularia trifolii and Their Contributions to Phylogenetic Implications in the Pleosporales

The symbiont endophytic fungi in tobacco are highly diverse and difficult to classify. Here, we sequenced the genomes of Curvularia trifolii and Leptosphaerulina chartarum isolated from tobacco plants. Finally, 41.68 Mb and 37.95 Mb nuclear genomes were sequenced for C. trifolii and L. chartarum with the scaffold N50, accounting for 638.94 Kb and 284.12 Kb, respectively. Meanwhile, we obtained 68,926 bp and 59,100 bp for their mitochondrial genomes. To more accurately classify C. trifolii and L. chartarum, we extracted seven nuclear genes and 12 mitochondrial genes from these two genomes and their closely related species. The genes were then used for calculation of evolutionary rates and for phylogenetic analysis. Results showed that it was difficult to achieve consistent results using a single gene due to their different evolutionary rates, while the phylogenetic trees obtained by combining datasets showed stable topologies. It is, therefore, more accurate to construct phylogenetic relationships for endophytic fungi based on multi-gene datasets. This study provides new insights into the distribution and characteristics of endophytic fungi in tobacco.