High variability of mitochondrial gene order among fungi

Uncovering the role of mitochondrial genome in pathogenicity and drug resistance in pathogenic fungi

Fungal infections are becoming more prevalent globally, particularly affecting immunocompromised populations, such as people living with HIV, organ transplant recipients and those on immunomodulatory therapy. Globally, approximately 6.55 million people are affected by invasive fungal infections annually, leading to serious health consequences and death. Mitochondria are membrane-bound organelles found in almost all eukaryotic cells and play an important role in cellular metabolism and energy production, including pathogenic fungi. These organelles possess their own genome, the mitochondrial genome, which is usually circular and encodes proteins essential for energy production. Variation and evolutionary adaptation within and between species' mitochondrial genomes can affect mitochondrial function, and consequently cellular energy production and metabolic activity, which may contribute to pathogenicity and drug resistance in certain fungal species. This review explores the link between the mitochondrial genome and mechanisms of fungal pathogenicity and drug resistance, with a particular focus on Cryptococcus neoformans and Candida albicans. These insights deepen our understanding of fungal biology and may provide new avenues for developing innovative therapeutic strategies.

Complete mitogenome and phylogenetic analysis of Brassica oleracea L. var. italica Plenck

In this study, we assembled and annotated the mitochondrial genome (mitogenome) of B. oleracea L. var. italica Plenck. This mitogenome was found to span 219 964 bp, with a GC content of 45.25%. It comprised 61 genes, including 35 protein-coding, 23 transfer RNA, and three ribosomal RNA genes. Among these, only 11 genes contained introns. Codon preference analysis revealed a bias toward codons ending in A/U bases. A phylogenetic analysis demonstrated a close relationship between B. oleracea L. var. italica Plenck, B. oleracea L. botrytis, and B. oleracea var. capitata. This reference mitogenome provides a basis for research on genetic conservation, phylogenetic relationships, and molecular breeding strategies among members of the Brassica genus.

Description and characterization of the small mitochondrial genome of Trichoderma cerinum (Hypocreales, Hypocreaceae) and its evolutionary perspectives

Trichoderma Persoon (1794) is a genus of fungus found in soils and decaying wood all over the world. We present the assembly and annotation of the mitochondrial genome of Trichoderma cerinum Bissett, Kubicek & Szakacs (2003) and an assessment of the phylogeny of the group, discussing the loss and gain of shared genes in the evolutionary history of fungi. We downloaded the raw data of T. cerinum from the NCBI database and exported it to the Galaxy Europe platform, where we performed the mitogenome assembly using the NOVOplasty tool. We used three tools for annotation. The phylogeny was conducted with 12 Trichoderma species and the T. cerinum. Fusarium oxysporum was used as an outgroup. We got a circularized mitochondrial genome of 26,696 bp, with 15 protein-coding genes, 25 tRNAs, two rRNAs, two endonuclease sequences, and an orf40. This species contains only one intron in the cob gene. In our phylogenetic reconstruction, T. cerinum was recovered as a sister group of a clade containing Trichoderma lixii, Trichoderma afroharzianum, Trichoderma simmonsii, and Trichoderma harzianum.

Evolutionary Relationships of the Phytophthora 1a Subclade Species Based on Complete Mitogenomes, and Novel Markers for Their Differentiation

Phytophthora is a genus of oomycetes that includes many aggressive pathogens capable of devastating farmlands and forests worldwide. Among the oldest and most well-known species, P. cactorum exhibits morphological and genetic similarities to other homothallic species within subclade 1a, which complicates the understanding of their evolutionary relationships. This study primarily compared seven P. cactorum strains from diverse origins with three closely related species in subclade 1a, utilizing mitogenome sequences for analysis. The circular mitogenomes of the four species were nearly identical in size and comprised 38 protein-coding genes (PCGs), 25 transfer ribonucleic acid genes, and 2 ribosomal RNA genes. The mitochondrial genomes exhibited a higher percentage of A/T compared to G/C content. The majority of AT-skew and GC-skew values among the 38 PCGs were positive, with the AT-skew demonstrating a more pronounced bias than the GC-skew. The Ka/Ks ratios revealed that 35 PCGs underwent significant purifying selection. Although the AliGROOVE analysis indicated notable similarities among the subclade 1a species, four PCGs exhibited significantly higher pairing frequency compared to the complete mitogenome. The results from the phylogenetic analysis aligned with the pairwise genetic distances, indicating that P. cactorum is more closely related to P. pseudotsugae than to P. hedraiandra. Furthermore, we found that the nad9 gene is informative to differentiate closely related Phytophthora species within subclade 1a, akin to the cox1 gene.

Characterisation and comparative analysis of mitochondrial genomes of false, yellow, black and blushing morels provide insights on their structure and evolution

Morchella species have considerable significance in terrestrial ecosystems, exhibiting a range of ecological lifestyles along the saprotrophism-to-symbiosis continuum. However, the mitochondrial genomes of these ascomycetous fungi have not been thoroughly studied, thereby impeding a comprehensive understanding of their genetic makeup and ecological role. In this study, we analysed the mitogenomes of 30 Morchellaceae species, including yellow, black, blushing and false morels. These mitogenomes are either circular or linear DNA molecules with lengths ranging from 217 to 565 kbp and GC content ranging from 38% to 48%. Fifteen core protein-coding genes, 28-37 tRNA genes and 3-8 rRNA genes were identified in these Morchellaceae mitogenomes. The gene order demonstrated a high level of conservation, with the cox1 gene consistently positioned adjacent to the rnS gene and cob gene flanked by apt genes. Some exceptions were observed, such as the rearrangement of atp6 and rps3 in Morchellaimportuna and the reversed order of atp6 and atp8 in certain morel mitogenomes. However, the arrangement of the tRNA genes remains conserved. We additionally investigated the distribution and phylogeny of homing endonuclease genes (HEGs) of the LAGLIDADG (LAGs) and GIY-YIG (GIYs) families. A total of 925 LAG and GIY sequences were detected, with individual species containing 19-48HEGs. These HEGs were primarily located in the cox1, cob, cox2 and nad5 introns and their presence and distribution displayed significant diversity amongst morel species. These elements significantly contribute to shaping their mitogenome diversity. Overall, this study provides novel insights into the phylogeny and evolution of the Morchellaceae.

Assembly and comparative analysis of the complete mitochondrial genome of Daedaleopsissinensis (Polyporaceae, Basidiomycota), contributing to understanding fungal evolution and ecological functions

Daedaleopsissinensis is a crucial wood-decaying fungus with significant lignocellulose-degrading ability, which plays a vital role in the material cycle and energy flow of forest ecosystems. However, the mitochondrial genome of D.sinensis has not yet been revealed. In the present study, the complete mitochondrial genome of D.sinensis was assembled and compared with related species. The mitochondrial genome spans 69,155 bp and has a GC content of 25.0%. It comprises 15 protein-coding genes (PCGs), 26 transfer RNA genes, two ribosomal RNA genes and one DNA polymerase gene (dpo). Herein, we characterised and analysed the codon preferences, variation and evolution of PCGs, repeats, intron dynamics, as well as RNA editing events in the D.sinensis mitochondrial genome. Further, a phylogenetic analysis of D.sinensis and the other 86 Basidiomycota species was performed using mitochondrial genome data. The results revealed that four species, D.confragosa, D.sinensis, D.nitida and Fomesfomentarius, were grouped in a closely-related cluster with high support values, indicating that a close phylogenetic relationship existed between Daedaleopsis and Fomes. This study reported on the initial assembly and annotation of the mitochondrial genome of D.sinensis, which greatly improved the knowledge of the fungus. These results contribute to the limited understanding of the mitochondrial repository of wood-decaying fungi, thereby laying the foundation for subsequent research on fungal evolution and ecological functions.

Global transcriptome changes during growth of a novel Penicillium coffeae isolate on the wheat stripe rust fungus, Puccinia striiformis f. sp. tritici

Wheat stripe rust caused by the fungus Puccinia striiformis f. sp. tritici (Pst) is currently the most destructive disease of wheat. The major control methods which include the deployment of resistant wheat cultivars and application of chemical fungicides are losing efficiency as the fungus evolves. Natural antagonists of Pst may be an avenue for alternative and environmentally sustainable control of the disease in the field. Here we describe a novel fungus found growing on Pst pustules. We identified the fungus as a novel isolate of the plant endophyte Penicillium coffeae. We present a high-quality reference genome and a comparative transcriptomic analysis used to investigate how the fungus deploys its genes during growth amongst Pst spores. The gene content of the P. coffeae ANU01 genome is suggestive of a generalist that makes use of diverse substrates. An abundance of genes related to lipid, amino acid and carbohydrate metabolism indicate that P. coffeae ANU01 has evolved the ability to exploit nutrient stores in Pst urediniospores. P. coffeae ANU01 deploys a number of biosynthetic gene clusters during growth on Pst spores, potentially to inhibit urediniospores germination and halt defence responses. A number of genes encoding carbohydrate active enzymes are also highly upregulated, suggesting targeting and degradation of Pst urediniospores structures. Alongside carbohydrates, P. coffeae ANU01 appears to target spore lipids as a nutrient source, secreting several highly upregulated lipases. Our findings broaden the understanding of growth associated with rust spores as an evolutionary strategy and provide insight into the genes potentially required for this process.

Frequent genetic exchanges revealed by a pan-mitogenome graph of a fungal plant pathogen

Mitochondria are present in almost all eukaryotic lineages. The mitochondrial genomes (mitogenomes) evolve separately from nuclear genomes, and they can therefore provide relevant insights into the evolution of their host species. Fusarium oxysporum is a major fungal plant pathogen that is assumed to reproduce clonally. However, horizontal chromosome transfer between strains can occur through heterokaryon formation, and recently, signs of sexual recombination have been observed. Similarly, signs of recombination in F. oxysporum mitogenomes challenged the prevailing assumption of clonal reproduction in this species. Here, we construct, to our knowledge, the first fungal pan-mitogenome graph of nearly 500 F. oxysporum mitogenome assemblies to uncover the variation and evolution. In general, the gene order of fungal mitogenomes is not well conserved, yet the mitogenome of F. oxysporum and related species are highly colinear. We observed two strikingly contrasting regions in the F. oxysporum pan-mitogenome, comprising a highly conserved core mitogenome and a long variable region (6-16 kb in size), of which we identified three distinct types. The pan-mitogenome graph reveals that only five intron insertions occurred in the core mitogenome and that the long variable regions drive the difference between mitogenomes. Moreover, we observed that their evolution is neither concurrent with the core mitogenome nor with the nuclear genome. Our large-scale analysis of long variable regions uncovers frequent recombination between mitogenomes, even between strains that belong to different taxonomic clades. This challenges the common assumption of incompatibility between genetically diverse F. oxysporum strains and provides new insights into the evolution of this fungal species.IMPORTANCEInsights into plant pathogen evolution is essential for the understanding and management of disease. Fusarium oxysporum is a major fungal pathogen that can infect many economically important crops. Pathogenicity can be transferred between strains by the horizontal transfer of pathogenicity chromosomes. The fungus has been thought to evolve clonally, yet recent evidence suggests active sexual recombination between related isolates, which could at least partially explain the horizontal transfer of pathogenicity chromosomes. By constructing a pan-genome graph of nearly 500 mitochondrial genomes, we describe the genetic variation of mitochondria in unprecedented detail and demonstrate frequent mitochondrial recombination. Importantly, recombination can occur between genetically diverse isolates from distinct taxonomic clades and thus can shed light on genetic exchange between fungal strains.

Assembly, Annotation, and Comparative Analysis of Mitochondrial Genomes in Trichoderma

Trichoderma is a widely studied ascomycete fungal genus, including more than 400 species. However, genetic information on Trichoderma is limited, with most species reporting only DNA barcodes. Mitochondria possess their own distinct DNA that plays a pivotal role in molecular function and evolution. Here, we report 42 novel mitochondrial genomes (mitogenomes) combined with 18 published mitogenomes of Trichoderma. These circular mitogenomes exhibit sizes of 26,276-94,608 bp, typically comprising 15 core protein-coding genes (PCGs), 2 rRNAs, and 16-30 tRNAs; however, the number of endonucleases and hypothetical proteins encoded in the introns of PCGs increases with genome size enlargement. According to the result of phylogenetic analysis of the whole mitogenome, these strains diverged into six distinct evolutionary branches, supported by the phylogeny based on 2830 single-copy nuclear genes. Comparative analysis revealed that dynamic Trichoderma mitogenomes exhibited variations in genome size, gene number, GC content, tRNA copy, and intron across different branches. We identified three mutation hotspots near the regions encoding nad3, cox2, and nad5 that caused major changes in the mitogenomes. Evolutionary analysis revealed that atp9, cob, nad4L, nad5, and rps3 have been influenced by positive selection during evolution. This study provides a valuable resource for exploring the important roles of the genetic and evolutionary dynamics of Trichoderma mitogenome in the adaptive evolution of biocontrol fungi.

Morphology, phylogeny, and mitogenomics reveal a new entomopathogenic fungus, Blackwellomyces changningensis (Hypocreales, Clavicipitaceae), from southwestern China

Two species of Blackwellomyces (Clavicipitaceae, Hypocreales, Ascomycota) were discovered during an investigation of the diversity of entomopathogenic fungi. A new fungus and one known fungal species that were gathered from Yunnan Province were described in this study. Blackwellomyces changningensis, sp. nov. was described using morphology and phylogenetic evidence from 14 mitochondrial protein-coding gene (PCG) data sets (atp6, atp8, atp9, cob, cox1, cox2, cox3, nad1, nad2, nad3, nad4, nad5, nad6, and nad4L) and six nuclear genes (ITS [ITS4 and ITS5], 18S nuc rDNA [18S], 28S nuc rDNA [28S], tef1-α, rpb1, and rpb2). B. changningensis were found parasitic to Lepidoptera larvae and to produce filiform ascospores with septations. The asexual conidia were ovoid to ellipsoid in shape. Phylogenetic analysis and morphological observations concurred that the fungus belonged to a different species within the genus of Blackwellomyces. The foundation for further taxonomic, genetic, and evolutionary biological studies of the genus Blackwellomyces was laid by this work.

Comparative mitochondrial genomics of Thelebolaceae in Antarctica: insights into their extremophilic adaptations and evolutionary dynamics

Species of Antarctomyces and Thelebolus (Thelebolaceae), primarily found in Antarctic environments, exhibit psychrophilic adaptations, yet their mitochondrial genomes have not been extensively studied. Furthermore, few studies have compared the mitochondrial genomes of psychrophilic, psychrotrophic, and mesophilic fungi. After successful sequencing and assembly, this study annotated the mitochondrial genomes of Antarctomyces psychrotrophicus CPCC 401038 and Thelebolus microsporus CPCC 401041. We also performed a comparative analysis with the previously characterized mitochondrial genomes of psychrotrophic and mesophilic fungi. The analysis revealed that nad4L was the most conserved gene across the mitochondrial genomes, characterized by its synonymous and non-synonymous substitution rates (Ks and Ka), genetic distance, and GC content and skew within the protein-coding genes (PCGs). Additionally, the mitochondrial genomes of psychrophilic and psychrotrophic fungi showed a higher proportion of protein-coding regions and a lower GC content compared to those of mesophilic fungi, underscoring the genetic basis of cold adaptation. Phylogenetic analyses based on these mitochondrial genes also confirmed the phylogenetic relationships of Thelebolaceae in the class Leotiomycetes. These findings advance our understanding of the phylogenetic relationships and evolutionary dynamics within the family Thelebolaceae, highlighting how different environmental temperatures influence fungal mitochondrial genomic structure and adaptation.

Comparative Mitogenomics Analysis Revealed Evolutionary Divergence among Purpureocillium Species and Gene Arrangement and Intron Dynamics of Ophiocordycipitaceae

The species of Purpureocillium are cosmopolitan and multitrophic fungi that can infect a wide range of invertebrate hosts. This study reports the mitogenome of P. atypicola, a specialized spider pathogenic fungus. The 112,465 bp mitogenome encoded genes typically found in fungal mitogenomes, and a total of 52 introns inserted into seven genes. A comparison with three other Purpureocillium species revealed significant differences in length and intron number, primarily due to intron variation; however, there was no dynamic variation in the introns of the cox1 gene within the same species of the Purpureocillium genus. Different mitochondrial protein-coding genes showed variable degrees of genetic differentiation among these species, but they were all under purifying selection. Additionally, frequent intron loss or gain events were detected to have occurred during the evolution of the Ophiocordycipitaceae mitogenomes, yet the gene arrangement remains conserved. A phylogenetic analysis of the combined mitochondrial gene set gave identical and well-supported tree topologies. The estimated age of the crown of Ophiocordycipitaceae and Purpureocillium were around the Early Cretaceous period (127 Mya) and Late Cretaceous period (83 Mya), respectively. The results of this study advance our understanding of the genomics, evolution, and taxonomy of this important fungal group.

The complete mitochondrial genome of the banana pathogen Fusarium oxysporum f. sp. cubense M5

We report the complete mitochondrial genome of a causal agent of banana fusarium wilt isolated in Mexico. The whole set of genes encoding proteins related to respiration and ATP synthesis, rRNAs, tRNAs are enlisted. Two open reading frames of unknown function conserved in Fusarium oxysporum were also identified.

Mitogenomics, phylogeny and morphology reveal two new entomopathogenic species of Ophiocordyceps (Ophiocordycipitaceae, Hypocreales) from south-western China

Ophiocordyceps encompasses over 300 species, demonstrating a wide range of morphological features, hosts and habitats within its species diversity. In this study, two novel species in Ophiocordyceps were revealed parasitising Hepialidae larva buried in soil. Ophiocordycepsalbastroma was morphologically characterised by white stromata, solitary and cylindrical conidiogenous cells and smooth ovoid or ellipsoidal conidia. Ophiocordycepsnigristroma was characterised by woody and dark brown stromata, monophialidic, swollen base and lageniform conidiogenous cells and smooth fusiform or oval conidia. The two new species formed a separate clade, respectively, based on the phylogenetic analyses of a combined dataset including nrSSU, nrLSU, rpb1, rpb2, and tef-1α, as well as a dataset of mitochondrial 14 protein coding genes (PCGs). They were all closely grouped with O.sinensis. The mitochondrial genomes of them were first reported. Their mitogenomes were all typical of circular molecules, with positive AT and GC skew, similar GC content, similar genetic composition, similar codon usage and conservative gene positions. However, the length of the mitogenomes varied. Changes in the length of the genes were the leading cause of changes in the length of mitochondrial genome of Ophiocordyceps. The discovery and identification of new Ophiocordyceps species and analysis their mitochondrial genomes may serve as foundations for phylogeny and diversity research within the genus Ophiocordyceps.

Comparative mitogenomic analysis of Sporisorium reilianum f. sp. zeae suggests recombination events during its evolutionary history

Introduction

Modern understanding of the concept of genetic diversity must include the study of both nuclear and organellar DNA, which differ greatly in terms of their structure, organization, gene content and distribution. This study comprises an analysis of the genetic diversity of the smut fungus Sporisorium reilianum f. sp. zeae from a mitochondrial perspective.

Methods

Whole-genome sequencing data was generated from biological samples of S. reilianum collected from different geographical regions. Multiple sequence alignment and gene synteny analysis were performed to further characterize genetic diversity in the context of mitogenomic polymorphisms.

Results

Mitochondria of strains collected in China contained unique sequences. The largest unique sequence stretch encompassed a portion of cox1, a mitochondrial gene encoding one of the subunits that make up complex IV of the mitochondrial electron transport chain. This unique sequence had high percent identity to the mitogenome of the related species Sporisorium scitamineum and Ustilago bromivora.

Discussion

The results of this study hint at potential horizontal gene transfer or mitochondrial genome recombination events during the evolutionary history of basidiomycetes. Additionally, the distinct polymorphic region detected in the Chinese mitogenome provides the ideal foundation to develop a diagnostic method to discern between mitotypes and enhance knowledge on the genetic diversity of this organism.

Determining Gene Order Patterns in the Suillus and Boletales through Comparative Analysis of Their Mitogenomes

Suillus is one of the most important genera of ectomycorrhizal fungi. As a model for studying host specificity, its molecular fragments and nuclear genome have been analyzed. However, its mitochondrial genome has not yet been reported. In this study, we assembled five mitogenomes of Suillus and analyzed and compared their basic characteristics. Owing to the large number of introns as well as intergenic regions, the mitogenomic lengths of species of Suillus were greater than those of other species of Boletales. We identified two main patterns of gene order arrangement in the members of the order Boletales. The Ka/Ks values of 15 protein-coding genes were <1 for the mitochondrial genes of 39 Boletales species, indicating their conserved evolution. Phylogenetic trees, reconstructed using the mitogenomes, indicated that the genus Suillus was monophyletic. Phylogenetic results based on the internal transcribed spacer region and mitogenome were used to confirm the distribution of Suillus placidus in China. The results showed that the mitogenome was superior in distinguishing species compared with a single molecular fragment. This is the first study to investigate the mitogenome of Suillus, enriching the mitogenome information and providing basic data for the phylogeny, resource conservation, and genetic diversity of this genus.

De novo assembly and comprehensive analysis of the mitochondrial genome of Taxus wallichiana reveals different repeats mediate recombination to generate multiple conformations

Taxus plants are the exclusive source of paclitaxel, an anticancer drug with significant medicinal and economic value. Interspecies hybridization and gene introgression during evolution have obscured distinctions among Taxus species, complicating their phylogenetic classification. While the chloroplast genome of Taxus wallichiana, a widely distributed species in China, has been sequenced, its mitochondrial genome (mitogenome) remains uncharacterized.We sequenced and assembled the T. wallichiana mitogenome using BGI short reads and Nanopore long reads, facilitating comparisons with other gymnosperm mitogenomes. The T. wallichiana mitogenome spanning 469,949 bp, predominantly forms a circular configuration with a GC content of 50.51%, supplemented by 3 minor configurations mediated by one pair of LRs and two pairs of IntRs. It includes 32 protein-coding genes, 7 tRNA genes, and 3 rRNA genes, several of which exist in multiple copies.We detailed the mitogenome's structure, codon usage, RNA editing, and sequence migration between organelles, constructing a phylogenetic tree to elucidate evolutionary relationships. Unlike typical gymnosperm mitochondria, T. wallichiana shows no evidence of mitochondrial-plastid DNA transfer (MTPT), highlighting its unique genomic architecture. Synteny analysis indicated extensive genomic rearrangements in T. wallichiana, likely driven by recombination among abundant repetitive sequences. This study offers a high-quality T. wallichiana mitogenome, enhancing our understanding of gymnosperm mitochondrial evolution and supporting further cultivation and utilization of Taxus species.

Comparative analysis of mitochondrial genomes in Ceratocystis fimbriata complex across diverse hosts

The decline ofAcacia mangiumWilld. in Malaysia, especially in Sabah since 2010, is primarily due to Ceratocystiswilt and canker disease (CWCD) caused by theCeratocystis fimbriataEllis & Halst. complex. This study was aimed to investigate the mitochondrial genome architecture of two differentC. fimbriatacomplex isolates from Malaysia: one fromA. mangiumin Pahang (FRIM1162) and another fromEucalyptus pellitain Sarawak (FRIM1441). This research employed Next-Generation Sequencing (NGS) to contrast genomes from diverse hosts with nine additional mitochondrial sequences, identifying significant genetic diversity and mutational hotspots in the mitochondrial genome alignment. The mitochondrial genome-based phylogenetic analysis revealed a significant genetic relationship between the studied isolates and theC. fimbriatacomplex in the South American Subclade, indicating that theC. fimbriatacomplex discovered in Malaysia isC. manginecans. The comparative mitochondrial genome demonstrates the adaptability of the complex due to mobile genetic components and genomic rearrangements in the studiedfungal isolates. This research enhances our knowledge of the genetic diversity and evolutionary patterns within theC. fimbriatacomplex, aiding in a deeper understanding of fungal disease development and host adaption processes. The acquired insights are crucial for creating specific management strategies for CWCD, improving the overall understanding of fungal disease evolution and control.

Comparative Analysis of the Mitochondrial Genome Sequences of Diaporthe longicolla (syn. Phomopsis longicolla) Isolates Causing Phomopsis Seed Decay in Soybean

Diaporthe longicolla (syn. Phomopsis longicolla) is an important seed-borne fungal pathogen and the primary cause of Phomopsis seed decay (PSD) in soybean. PSD is one of the most devastating seed diseases, reducing soybean seed quality and yield worldwide. As part of a genome sequencing project on the fungal Diaporthe-Phomopsis complex, draft genomes of eight D. longicolla isolates were sequenced and assembled. Sequences of mitochondrial genomes were extracted and analyzed. The circular mitochondrial genomes ranged from 52,534 bp to 58,280 bp long, with a mean GC content of 34%. A total of 14 core protein-coding genes, 23 tRNA, and 2 rRNA genes were identified. Introns were detected in the genes of atp6, cob, cox1, cox2, cox3, nad1, nad2, nad5, and rnl. Three isolates (PL7, PL10, and PL185E) had more introns than other isolates. Approximately 6.4% of the mitochondrial genomes consist of repetitive elements. Moreover, 48 single-nucleotide polymorphisms (SNPs) and were identified. The mitochondrial genome sequences of D. longicolla will be useful to further study the molecular basis of seed-borne pathogens causing seed diseases, investigate genetic variation among isolates, and develop improved control strategies for Phomopsis seed decay of soybean.

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

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

The complete mitochondrial genome of Penicillium expansum: Insights into the fungal evolution and phylogeny

Penicillium expansum is an important phytopathogenic fungus that causes blue mold disease. In this study, the novel mitochondrial genome of P. expansum was sequenced, assembled, annotated, and compared with the previously published Penicillium mitogenomes. P. expansum mitogenome is composed of circular DNA molecules with a genome size of 25,496 bp. It encodes 16 protein-encoding genes (PCGs), two rRNA genes, and 25 tRNA genes. Comparative analysis with six other Penicillium species revealed that gene length, GC content, AT skew, and GC skew were variable among the core protein-coding genes. The Penicillium species' gene synteny analysis identified several gene rearrangements. Among the core 15 PCGs, atp8 had the lowest K2P genetic distance, which shows that this gene is highly conserved. The Ka/Ks value of most PCGs was less than 1, which shows that these genes have undergone purifying selection. Phylogenetic analysis based on 14 concatenated core mitochondrial genes revealed that P. expansum shares a close relationship with P. solitum. This study served as a first report on the complete mitochondrial genome of P. expansum and its comparative analysis that will contribute to population genetics and rapid evolutionary studies among Penicillium species.

Highly Reactive Group I Introns Ubiquitous in Pathogenic Fungi

Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.

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

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

The first two whole mitochondrial genomes for the genus Dactylis species: assembly and comparative genomics analysis

BACKGROUND

Orchardgrass (Dactylis glomerata L.), a perennial forage, has the advantages of rich leaves, high yield, and good quality and is one of the most significant forage for grassland animal husbandry and ecological management in southwest China. Mitochondrial (mt) genome is one of the major genetic systems in plants. Studying the mt genome of the genus Dactylis could provide more genetic information in addition to the nuclear genome project of the genus.

RESULTS

In this study, we sequenced and assembled two mitochondrial genomes of Dactylis species of D. glomerata (597, 281 bp) and D. aschersoniana (613, 769 bp), based on a combination of PacBio and Illumina. The gene content in the mitochondrial genome of D. aschersoniana is almost identical to the mitochondrial genome of D. glomerata, which contains 22-23 protein-coding genes (PCGs), 8 ribosomal RNAs (rRNAs) and 30 transfer RNAs (tRNAs), while D. glomerata lacks the gene encoding the Ribosomal protein (rps1) and D. aschersoniana contains one pseudo gene (atp8). Twenty-three introns were found among eight of the 30 protein-coding genes, and introns of three genes (nad 1, nad2, and nad5) were trans-spliced in Dactylis aschersoniana. Further, our mitochondrial genome characteristics investigation of the genus Dactylis included codon usage, sequences repeats, RNA editing and selective pressure. The results showed that a large number of short repetitive sequences existed in the mitochondrial genome of D. aschersoniana, the size variation of two mitochondrial genomes is due largely to the presence of a large number of short repetitive sequences. We also identified 52-53 large fragments that were transferred from the chloroplast genome to the mitochondrial genome, and found that the similarity was more than 70%. ML and BI methods used in phylogenetic analysis revealed that the evolutionary status of the genus Dactylis.

CONCLUSIONS

Thus, this study reveals the significant rearrangements in the mt genomes of Pooideae species. The sequenced Dactylis mt genome can provide more genetic information and improve our evolutionary understanding of the mt genomes of gramineous plants.

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 mitogenomes characterization of two Peltigera species (Peltigera elisabethae and Peltigera polydactylon) and comparative mitogenomic analyses of six Peltigera

In the present study, the mitochondrial genomes of Peltigera elisabethae and P. polydactylon were sequenced and assembled. The two mitogenomes were composed of circular DNA molecules, with sizes of 64,034 bp and 59,208 bp, respectively. Comparative analysis showed that the genome size, GC content, GC skew, and AT skew varied between the two mitochondrial genomes. In codon analysis, phenylalanine (Phe), isoleucine (Ile), and leucine (Leu) were most frequently used in six Peltigera genomes. Evolutionary analysis showed that all 14 protein-coding genes (PCGs) were subject to purifying selection in the six Peltigera species. Regarding gene rearrangement, the PCGs of Peltigera had the same gene sequence and gene content, and a few intron sequences and spacer sequences were rearranged in Peltigera. In the phylogenetic analysis, we used Bayesian and ML methods to construct a phylogenetic tree. Two phylogenetic trees with consistent topology with high support indicate that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships. The present study enriches the mitochondrial genome data of Peltigera and promotes further understanding of the genetics and evolution of the Peltigera genus.

The mitogenomes of Leptographium aureum, Leptographium sp., and Grosmannia fruticeta: expansion by introns

Introduction

Many members of the Ophiostomatales are of economic importance as they are bark-beetle associates and causative agents for blue stain on timber and in some instances contribute towards tree mortality. The taxonomy of these fungi has been challenging due to the convergent evolution of many traits associated with insect dispersal and a limited number of morphological characters that happen to be highly pleomorphic. This study examines the mitochondrial genomes for three members of Leptographium sensu lato [Leptographium aureum (also known as Grosmannia aurea), Grosmannia fruticeta (also known as Leptographium fruticetum), and Leptographium sp. WIN(M)1376)].

Methods

Illumina sequencing combined with gene and intron annotations and phylogenetic analysis were performed.

Results

Sequence analysis showed that gene content and gene synteny are conserved but mitochondrial genome sizes were variable: G. fruticeta at 63,821 bp, Leptographium sp. WIN(M)1376 at 81,823 bp and L. aureum at 104,547 bp. The variation in size is due to the number of introns and intron-associated open reading frames. Phylogenetic analysis of currently available mitochondrial genomes for members of the Ophiostomatales supports currently accepted generic arrangements within this order and specifically supports the separation of members with Leptographium-like conidiophores into two genera, with L. aureum grouping with Leptographium and G. fruticeta aligning with Grosmannia.

Discussion

Mitochondrial genomes are promising sequences for resolving evolutionary relationships within the Ophiostomatales.

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.

Co-evolution of large inverted repeats and G-quadruplex DNA in fungal mitochondria may facilitate mitogenome stability: the case of Malassezia

Mitogenomes are essential due to their contribution to cell respiration. Recently they have also been implicated in fungal pathogenicity mechanisms. Members of the basidiomycetous yeast genus Malassezia are an important fungal component of the human skin microbiome, linked to various skin diseases, bloodstream infections, and they are increasingly implicated in gut diseases and certain cancers. In this study, the comparative analysis of Malassezia mitogenomes contributed to phylogenetic tree construction for all species. The mitogenomes presented significant size and gene order diversity which correlates to their phylogeny. Most importantly, they showed the inclusion of large inverted repeats (LIRs) and G-quadruplex (G4) DNA elements, rendering Malassezia mitogenomes a valuable test case for elucidating the evolutionary mechanisms responsible for this genome diversity. Both LIRs and G4s coexist and convergently evolved to provide genome stability through recombination. This mechanism is common in chloroplasts but, hitherto, rarely found in mitogenomes.

Functional annotation and comparative analysis of four Botrytis cinerea mitogenomes reported from Punjab, Pakistan

Botrytis cinerea is one of the top phytopathogenic fungus which ubiquitously cause grey mold on a variety of horticultural plants. The mechanism of respiration in the fungus occurs within the mitochondria. Mitogenomes serve as a key molecular marker for the investigation of fungal evolutionary patterns. This study aimed at the complete assembly, characterization, and comparative relationship of four mitogenomes of Botrytis cinerea strains including Kst5C, Kst14A, Kst32B, Kst33A, respectively. High throughput sequencing of four mitogenomes allowed the full assembly and annotation of these sequences. The total genome length of these 4 isolates Kst5C Kst14A, Kst32B, Kst33A was 69,986 bp, 77,303 bp, 76,204 bp and 55, 226 bp respectively. The distribution of features represented 2 ribosomal RNA genes,14 respiration encoding proteins, 1 mitochondrial ribosomal protein-encoding gene, along with varying numbers of transfer RNA genes, protein-coding genes, mobile intronic regions and homing endonuclease genes including LAGLIDADG and GIY-YIG domains were found in all four mitogenomes. The comparative analyses performed also decipher significant results for four mitogenomes among fungal isolates included in the study. This is the first report on the detailed annotation of mitogenomes as a proof for investigation of variation patterns present with in the B. cinerea causing grey mold on strawberries in Pakistan. This study will also contribute to the rapid evolutionary analysis and population patterns present among Botrytis cinerea.

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.

Morphology, phylogeny, mitogenomics and metagenomics reveal a new entomopathogenic fungus Ophiocordycepsnujiangensis (Hypocreales, Ophiocordycipitaceae) from Southwestern China

Ophiocordyceps contains the largest number of Cordyceps sensu lato, various species of which are of great medicinal value. In this study, a new entomopathogenic fungus, Ophiocordycepsnujiangensis, from Yunnan in southwestern China, was described using morphological, phylogenetic, and mitogenomic evidence, and its fungal community composition was identified. It was morphologically characterized by a solitary, woody, and dark brown stromata, smooth-walled and septate hyphae, solitary and gradually tapering conidiogenous cells with plenty of warty protrusions, and oval or fusiform conidia (6.4-11.2 × 3.7-6.4 µm) with mucinous sheath. The phylogenetic location of O.nujiangensis was determined based on the Bayesian inference (BI) and the maximum likelihood (ML) analyses by concatenating nrSSU, nrLSU, tef-1a, rpb1, and rpb2 datasets, and ten mitochondrial protein-coding genes (PCGs) datasets (atp6, atp9, cob, cox2, nad1, nad2, nad3, nad4, nad4L, and nad5). Phylogenetic analyses revealed that O.nujiangensis belonged to the Hirsutellasinensis subclade within the Hirsutella clade of Ophiocordyceps. And O.nujiangensis was phylogenetically clustered with O.karstii, O.liangshanensis, and O.sinensis. Simultaneously, five fungal phyla and 151 fungal genera were recognized in the analysis of the fungal community of O.nujiangensis. The fungal community composition differed from that of O.sinensis, and differences in the microbial community composition of closely related species might be appropriate as further evidence for taxonomy.

Microheliella maris possesses the most gene-rich mitochondrial genome in Diaphoretickes

The mitochondrial genomes are very diverse, but their evolutionary history is unclear due to the lack of efforts to sequence those of protists (unicellular eukaryotes), which cover a major part of the eukaryotic tree. Cryptista comprises cryptophytes, goniomonads, kathablepharids, and Palpitomonas bilix, and their mitochondrial genomes (mt-genomes) are characterized by various gene contents, particularly the presence/absence of an ancestral (bacterial) system for the cytochrome c maturation system. To shed light on mt-genome evolution in Cryptista, we report the complete mt-genome of Microheliella maris, which was recently revealed to branch at the root of Cryptista. The M. maris mt-genome was reconstructed as a circular mapping chromosome of 61.2 kbp with a pair of inverted repeats (12.9 kbp) and appeared to be the most gene-rich among the mt-genomes of the members of Diaphoretickes (a mega-scale eukaryotic assembly including Archaeplastida, Cryptista, Haptista, and SAR) studied so far, carrying 53 protein-coding genes. With this newly sequenced mt-genome, we inferred and discussed the evolution of the mt-genome in Cryptista and Diaphoretickes.

Trichoderma - genomes and genomics as treasure troves for research towards biology, biotechnology and agriculture

The genus Trichoderma is among the best studied groups of filamentous fungi, largely because of its high relevance in applications from agriculture to enzyme biosynthesis to biofuel production. However, the physiological competences of these fungi, that led to these beneficial applications are intriguing also from a scientific and ecological point of view. This review therefore summarizes recent developments in studies of fungal genomes, updates on previously started genome annotation efforts and novel discoveries as well as efforts towards bioprospecting for enzymes and bioactive compounds such as cellulases, enzymes degrading xenobiotics and metabolites with potential pharmaceutical value. Thereby insights are provided into genomes, mitochondrial genomes and genomes of mycoviruses of Trichoderma strains relevant for enzyme production, biocontrol and mycoremediation. In several cases, production of bioactive compounds could be associated with responsible genes or clusters and bioremediation capabilities could be supported or predicted using genome information. Insights into evolution of the genus Trichoderma revealed large scale horizontal gene transfer, predominantly of CAZyme genes, but also secondary metabolite clusters. Investigation of sexual development showed that Trichoderma species are competent of repeat induced point mutation (RIP) and in some cases, segmental aneuploidy was observed. Some random mutants finally gave away their crucial mutations like T. reesei QM9978 and QM9136 and the fertility defect of QM6a was traced back to its gene defect. The Trichoderma core genome was narrowed down to 7000 genes and gene clustering was investigated in the genomes of multiple species. Finally, recent developments in application of CRISPR/Cas9 in Trichoderma, cloning and expression strategies for the workhorse T. reesei as well as the use genome mining tools for bioprospecting Trichoderma are highlighted. The intriguing new findings on evolution, genomics and physiology highlight emerging trends and illustrate worthwhile perspectives in diverse fields of research with Trichoderma.

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.

Intraspecific comparison of mitochondrial genomes reveals the evolution in medicinal fungus Ganoderma lingzhi

Several mitogenomes of the genus Ganoderma have been assembled, but intraspecific comparisons of mitogenomes in Ganoderma lingzhi have not been reported. In this study, 19 G. lingzhi mitogenomes were assembled and analyzed combined with three mitogenomes of G. lingzhi from GenBank in term of the characteristics, evolution, and phylogeny. The results showed that the mitogenomes of the G. lingzhi strains are closed circular ranging from 49.23 kb to 68.37 kb. The genetic distance, selective pressure, and base variation indicate that the 14 common protein coding genes were highly conserved. The differences in introns, open reading frames, and repetitive sequences in the mitogenome were the main factors leaded to the variations in mitogenome. The introns were horizontally transferred in mitogenomes, and the differences between introns in the same insertion, which were primarily caused by the repetitive sequence, showed that the introns may be under degeneration. Besides, the frequent insertion and deletion of introns showed an evolutionary rate faster than protein coding genes. Phylogenetic analysis showed that the G. lingzhi strains gathered with high support, and those with the same intron distribution law had closer clustering relationships.

Intraspecies variation of the mitochondrial genome: An evaluation for phylogenetic approaches based on the conventional choices of genes and segments on mitogenome

Intraspecies nucleotide sequence variation is a key to understanding the evolutionary history of a species, such as the geographic distribution and population structure. To date, numerous phylogenetic and population genetics studies have been conducted based on the sequences of a gene or an intergenic region on the mitochondrial genome (mtDNA), such as cytochrome c oxidase subunits or the D-loop. To evaluate the credibility of the usage of such 'classic' markers, we compared the phylogenetic inferences based on the analyses of the partial and entire mtDNA sequences. Importantly, the phylogenetic reconstruction based on the short marker sequences did not necessarily reproduce the tree topologies based on the analyses of the entire mtDNA. In addition, analyses on the datasets of various organisms revealed that the analyses based on the classic markers yielded phylogenetic trees with poor confidence in all tested cases compared to the results based on full-length mtDNA. These results demonstrated that phylogenetic analyses based on complete mtDNA sequences yield more insightful results compared to those based on mitochondrial genes and segments. To ameliorate the shortcomings of the classic markers, we identified a segment of mtDNA that may be used as an 'approximate marker' to closely reproduce the phylogenetic inference obtained from the entire mtDNA in the case of mammalian species, which can be utilized to design amplicon-seq-based studies. Our study demonstrates the importance of the choice of mitochondrial markers for phylogenetic analyses and proposes a novel approach to choosing appropriate markers for mammalian mtDNA that reproduces the phylogenetic inferences obtained from full-length mtDNA.

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 compact mitogenome of Ceratocystiopsis pallidobrunnea

Ceratocystiopsis is a fungal genus that has been assigned to the Ophiostomatales, fungi known for their association with various bark beetles and other arthropods. The mitochondrial genome of Ceratocystiopsis pallidobrunnea has been characterized and compared with other members of the genus Ceratocystiopsis and Ophiostomatales. At 29 022 bp, the mitogenome of C. pallidobrunnea is the smallest reported so far for this genus. Gene arrangement was observed to be conserved for this group of fungi and mitogenome variation appears to be mostly due to the absence and presence of introns. The long-term goal is to apply mitogenomes to resolve taxonomic issues within the Ophiostomatales and within the various genera that comprise the Ophiostomataceae.

The mitogenome of Urnula craterium

Urnula craterium (Schwein.) Fr. (1851) has been reported from North America, Europe, and Asia, and can be a pathogen on various hardwood species. In this study we investigated the mitochondrial genome of U. craterium. The biology and taxonomy of this fungus is poorly studied and there are no mitogenomes currently available for any member of the Sarcosomataceae (Order Pezizales). The complete mitogenome of U. craterium comprises 43 967 bps and encodes 14 protein-coding genes, a complete set of tRNAs and rRNA genes. A novel feature of the mitogenome is the presence of a single subunit DNA polymerase coding region that is typically associated with linear invertron-type plasmids. The mitogenome may offer insights into the evolution of mitogenomes among members of the Pezizales with regards to gene content and order, mobile elements, and genome sizes.

Population genetics and microevolution of clinical Candida glabrata reveals recombinant sequence types and hyper-variation within mitochondrial genomes, virulence genes, and drug targets

Candida glabrata is the second most common etiological cause of worldwide systemic candidiasis in adult patients. Genome analysis of 68 isolates from 8 hospitals across Scotland, together with 83 global isolates, revealed insights into the population genetics and evolution of C. glabrata. Clinical isolates of C. glabrata from across Scotland are highly genetically diverse, including at least 19 separate sequence types that have been recovered previously in globally diverse locations, and 1 newly discovered sequence type. Several sequence types had evidence for ancestral recombination, suggesting transmission between distinct geographical regions has coincided with genetic exchange arising in new clades. Three isolates were missing MATα1, potentially representing a second mating type. Signatures of positive selection were identified in every sequence type including enrichment for epithelial adhesins thought to facilitate fungal adhesin to human epithelial cells. In patent microevolution was identified from 7 sets of recurrent cases of candidiasis, revealing an enrichment for nonsynonymous and frameshift indels in cell surface proteins. Microevolution within patients also affected epithelial adhesins genes, and several genes involved in drug resistance including the ergosterol synthesis gene ERG4 and the echinocandin target FKS1/2, the latter coinciding with a marked drop in fluconazole minimum inhibitory concentration. In addition to nuclear genome diversity, the C. glabrata mitochondrial genome was particularly diverse, with reduced conserved sequence and conserved protein-encoding genes in all nonreference ST15 isolates. Together, this study highlights the genetic diversity within the C. glabrata population that may impact virulence and drug resistance, and 2 major mechanisms generating this diversity: microevolution and genetic exchange/recombination.

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.

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

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

Comparative analysis of the complete mitochondrial genomes of four cordyceps fungi

Cordyceps is a large group of entomogenous, medicinally important fungi. In this study, we sequenced, assembled, and annotated the entire mitochondrial genome of Ophiocordyceps xuefengensis, in addition to comparing it against other three complete cordyceps mitogenomes that were previously published. Comparative analysis indicated that the four complete mitogenomes are all composed of circular DNA molecules, although their sizes significantly differ due to high variability in intron and intergenic region sizes in the Ophiocordyceps sinensis and O. xuefengensis mitogenomes. All mitogenomes contain 14 conserved genes and two ribosomal RNA genes, but varying numbers of tRNA introns. The Ka/Ks ratios for all 14 PCGs and rps3 were all less than 1, indicating that these genes have been subject to purifying selection. Phylogenetic analysis was conducted using concatenated amino acid and nucleotide sequences of the 14 PCGs and rps3 using two different methods (Maximum Likelihood and Bayesian analysis), revealing highly supported relationships between O. xuefengensis and other Ophiocordyceps species, in addition to a close relationship with O. sinensis. Further, the analyses indicated that cox1 and rps3 play important roles in population differentiation. These mitogenomes will allow further study of the population genetics, taxonomy, and evolutionary biology of medicinally important cordyceps species.

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

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

Mitochondrial Transcription of Entomopathogenic Fungi Reveals Evolutionary Aspects of Mitogenomes

Entomopathogenic fungi and more specifically genera Beauveria and Metarhizium have been exploited for the biological control of pests. Genome analyses are important to understand better their mode of action and thus, improve their efficacy against their hosts. Until now, the sequences of their mitochondrial genomes were studied, but not at the level of transcription. Except of yeasts and Neurospora crassa, whose mt gene transcription is well described, in all other Ascomycota, i.e., Pezizomycotina, related information is extremely scarce. In this work, mt transcription and key enzymes of this function were studied. RT-PCR experiments and Northern hybridizations reveal the transcriptional map of the mt genomes of B. bassiana and M. brunneum species. The mt genes are transcribed in six main transcripts and undergo post-transcriptional modifications to create single gene transcripts. Promoters were determined in both mt genomes with a comparative in silico analysis, including all known information from other fungal mt genomes. The promoter consensus sequence is 5'-ATAGTTATTAT-3' which is in accordance with the definition of the polycistronic transcripts determined with the experiments described above. Moreover, 5'-RACE experiments in the case of premature polycistronic transcript nad1-nad4-atp8-atp6 revealed the 5' end of the RNA transcript immediately after the in silico determined promoter, as also found in other fungal species. Since several conserved elements were retrieved from these analyses compared to the already known data from yeasts and N. crassa, the phylogenetic analyses of mt RNA polymerase (Rpo41) and its transcriptional factor (Mtf1) were performed in order to define their evolution. As expected, it was found that fungal Rpo41 originate from the respective polymerase of T7/T3 phages, while the ancestor of Mtf1 is of alpha-proteobacterial origin. Therefore, this study presents insights about the fidelity of the mt single-subunit phage-like RNA polymerase during transcription, since the correct identification of mt promoters from Rpo41 requires an ortholog to bacterial sigma factor, i.e., Mtf1. Thus, a previously proposed hypothesis of a phage infected alpha-proteobacterium as the endosymbiotic progenitor of mitochondrion is confirmed in this study and further upgraded by the co-evolution of the bacterial (Mtf1) and viral (Rpo41) originated components in one functional unit.