Fungal mitochondrial genomes and genetic polymorphisms

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

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

Phyllosticta paracitricarpa is synonymous with the EU quarantine fungus P. citricarpa based on phylogenomic analyses

Phyllosticta citricarpa is an important citrus-pathogen and a quarantine organism in the European Union. Its recently described relative, P. paracitricarpa, is very closely related and not listed as a quarantine organism. P. paracitricarpa is very difficult to distinguish from P. citricarpa, since its morphological features overlap and the barcoding gene sequences that were originally used to delimit them as distinct species have a low number of species-specific polymorphisms that have subsequently been shown to overlap between the two clades. Therefore, we performed extensive genomic analyses to determine whether the genetic variation between P. citricarpa and P. paracitricarpa strains should be considered to represent infraspecific variation within P. citricarpa, or whether it is indicative of distinct species. Using a phylogenomic analysis with 3,000 single copy ortholog genes and whole-genome comparisons, we determined that the variation between P. citricarpa and P. paracitricarpa can be considered as infraspecies variation within P. citricarpa. We also determined the level of variation in mitochondrial assemblies of several Phyllosticta species and concluded there are only minimal differences between the assemblies of P. citricarpa and P. paracitricarpa. Thus, using several orthogonal approaches, we here demonstrate that variation within the nuclear and mitochondrial genomes of other Phyllosticta species is larger than variation between genomes obtained from P. citricarpa and P. paracitricarpa strains. Thus, P. citricarpa and P. paracitricarpa should be considered as conspecific.

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.

Exploring Mitochondrial Heterogeneity and Evolutionary Dynamics in Thelephora ganbajun through Population Genomics

Limited exploration in fungal mitochondrial genetics has uncovered diverse inheritance modes. The mitochondrial genomes are inherited uniparentally in the majority of sexual eukaryotes, our discovery of persistent mitochondrial heterogeneity within the natural population of the basidiomycete fungus Thelephora ganbajun represents a significant advance in understanding mitochondrial inheritance and evolution in eukaryotes. Here, we present a comprehensive analysis by sequencing and assembling the complete mitogenomes of 40 samples exhibiting diverse cox1 heterogeneity patterns from various geographical origins. Additionally, we identified heterogeneous variants in the nad5 gene, which, similar to cox1, displayed variability across multiple copies. Notably, our study reveals a distinct prevalence of introns and homing endonucleases in these heterogeneous genes. Furthermore, we detected potential instances of horizontal gene transfer involving homing endonucleases. Population genomic analyses underscore regional variations in mitochondrial genome composition among natural samples exhibiting heterogeneity. Thus, polymorphisms in heterogeneous genes, introns, and homing endonucleases significantly influence mitochondrial structure, structural variation, and evolutionary dynamics in this species. This study contributes valuable insights into mitochondrial genome architecture, population dynamics, and the evolutionary implications of mitochondrial heterogeneity in sexual eukaryotes.

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.

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.

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made towards understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

Mitochondrial genome fragmentation is correlated with increased rates of molecular evolution

While mitochondrial genome content and organization is quite diverse across all Eukaryotes, most bilaterian animal mitochondrial genomes (mitogenomes) exhibit highly conserved gene content and organisation, with genes typically encoded on a single circular chromosome. However, many species of parasitic lice (Insecta: Phthiraptera) are among the notable exceptions, having mitogenomes fragmented into multiple circular chromosomes. To better understand the process of mitogenome fragmentation, we conducted a large-scale genomic study of a major group of lice, Amblycera, with extensive taxon sampling. Analyses of the evolution of mitogenome structure across a phylogenomic tree of 90 samples from 53 genera revealed evidence for multiple independent origins of mitogenome fragmentation, some inferred to have occurred less than five million years ago. We leveraged these many independent origins of fragmentation to compare the rates of DNA substitution and gene rearrangement, specifically contrasting branches with fragmented and non-fragmented mitogenomes. We found that lineages with fragmented mitochondrial genomes had significantly higher rates of mitochondrial sequence evolution. In addition, lineages with fragmented mitochondrial genomes were more likely to have mitogenome gene rearrangements than those with single-chromosome mitochondrial genomes. By combining phylogenomics and mitochondrial genomics we provide a detailed portrait of mitogenome evolution across this group of insects with a remarkably unstable mitogenome structure, identifying processes of molecular evolution that are correlated with mitogenome fragmentation.

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.

Mitochondrial genome diversity across the subphylum Saccharomycotina

Introduction

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available.

Methods

By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades.

Results

Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements.

Discussion

As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

Mitogenome Variations in a Global Population of Aspergillus fumigatus

Aspergillus fumigatus is a ubiquitous, critical priority human fungal pathogen. Despite its clinical importance, there is limited knowledge regarding the variations of the genome within mitochondria, the powerhouse organelle within eukaryotic cells. In this study, we leveraged publicly available, raw, whole genome sequence data isolates from 1939 to investigate the variations in the mitochondrial genomes of A. fumigatus. These isolates were isolated from 22 countries on six continents, as well as from outer space and from within the International Space Station. In total, our analysis revealed 39 mitochondrial single nucleotide polymorphisms (mtSNPs) within this global sample, and, together, these 39 mtSNPs grouped the 1939 isolates into 79 mitochondrial multilocus genotypes (MLGs). Among the 79 MLGs, 39 were each distributed in at least two countries and 30 were each shared by at least two continents. The two most frequent MLGs were also broadly distributed: MLG11 represented 420 isolates from 11 countries and four continents and while MLG79 represented 418 isolates from 18 countries and five continents, consistent with long-distance dispersals of mitogenomes. Our population genetic analyses of the mtSNPs revealed limited differentiation among continental populations, but highly variable genetic differences among national populations, largely due to localized clonal expansions of different MLGs. Phylogenetic analysis and Discriminant Analysis of Principal Components of mtSNPs suggested the presence of at least three mitogenome clusters. Linkage disequilibrium, Index of Association, and phylogenetic incompatibility analyses collectively suggested evidence for mitogenome recombination in natural populations of A. fumigatus. In addition, sequence read depth analyses revealed an average ratio of ~20 mitogenomes per nuclear genome in this global population, but the ratios varied among strains within and between certain geographic populations. Together, our results suggest evidence for organelle dynamics, genetic differentiation, recombination, and both widespread and localized clonal expansion of the mitogenomes in the global A. fumigatus population.

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.

Mitochondrial Genome Diversity across the Subphylum Saccharomycotina

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available. By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades. Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements. As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

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

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

De novo mitochondrial genome sequencing of Cladonia subulata and phylogenetic analysis with other dissimilar species

In this study, the complete mitochondrial genome of Cladonia subulata (L.) FH Wigg was sequenced and assembled and then compared with those of other Cladonia species. The mitogenome of Cladonia subulata, the type species of Cladonia, consisted of a circular DNA molecule of 58,895 bp 44 genes (15 protein-coding genes, 2 rRNA genes, and 27 tRNA genes). The base composition had shown an obvious AT preference, and all 27 tRNA genes formed a typical clover structure. Comparison with other 7 Cladonia species indicated that the duplication/loss of tRNAs had occurred during evolution, and introns appeared to explain the variation in cox1 genes in Cladonia, the mitochondrial genome tends to be generally conservative and local dynamic changes. Repeat sequences were mainly located in gene intervals, which were mainly distributed among intergenic spacers and may cause rearrangement of the mitogenome. The phylogenetic results showed that Cladonia subulata and C. polycarpoides were assigned to the Cladonia Subclade. The results add to the available mitochondrial genome sequence information of Cladonia subulata, provide basic data for the systematic development, resource protection, and genetic diversity research in Cladonia subulata, and also provide theoretical support for further genomic research of lichens.

Uniparental Inheritance and Recombination as Strategies to Avoid Competition and Combat Muller's Ratchet among Mitochondria in Natural Populations of the Fungus Amanita phalloides

Uniparental inheritance of mitochondria enables organisms to avoid the costs of intracellular competition among potentially selfish organelles. By preventing recombination, uniparental inheritance may also render a mitochondrial lineage effectively asexual and expose mitochondria to the deleterious effects of Muller's ratchet. Even among animals and plants, the evolutionary dynamics of mitochondria remain obscure, and less is known about mitochondrial inheritance among fungi. To understand mitochondrial inheritance and test for mitochondrial recombination in one species of filamentous fungus, we took a population genomics approach. We assembled and analyzed 88 mitochondrial genomes from natural populations of the invasive death cap Amanita phalloides, sampling from both California (an invaded range) and Europe (its native range). The mitochondrial genomes clustered into two distinct groups made up of 57 and 31 mushrooms, but both mitochondrial types are geographically widespread. Multiple lines of evidence, including negative correlations between linkage disequilibrium and distances between sites and coalescent analysis, suggest low rates of recombination among the mitochondria (ρ = 3.54 × 10-4). Recombination requires genetically distinct mitochondria to inhabit a cell, and recombination among A. phalloides mitochondria provides evidence for heteroplasmy as a feature of the death cap life cycle. However, no mushroom houses more than one mitochondrial genome, suggesting that heteroplasmy is rare or transient. Uniparental inheritance emerges as the primary mode of mitochondrial inheritance, even as recombination appears as a strategy to alleviate Muller's ratchet.

Insight into the evolutionary and domesticated history of the most widely cultivated mushroom Agaricus bisporus via mitogenome sequences of 361 global strains

Agaricus bisporus is the most widely cultivated edible mushroom in the world with a only around three hundred years known history of cultivation. Therefore, it represents an ideal organism not only to investigate the natural evolutionary history but also the understanding on the evolution going back to the early era of domestication. In this study, we generated the mitochondrial genome sequences of 352 A. bisporus strains and 9 strains from 4 closely related species around the world. The population mitogenomic study revealed all A. bisporus strains can be divided into seven clades, and all domesticated cultivars present only in two of those clades. The molecular dating analysis showed this species origin in Europe on 4.6 Ma and we proposed the main dispersal routes. The detailed mitogenome structure studies showed that the insertion of the plasmid-derived dpo gene caused a long fragment (MIR) inversion, and the distributions of the fragments of dpo gene were strictly in correspondence with these seven clades. Our studies also showed A. bisporus population contains 30 intron distribution patterns (IDPs), while all cultivars contain only two IDPs, which clearly exhibit intron loss compared to the others. Either the loss occurred before or after domestication, that could suggest that the change facilitates their adaptation to the cultivated environment.

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.

Macroevolutionary diversity of traits and genomes in the model yeast genus Saccharomyces

Species is the fundamental unit to quantify biodiversity. In recent years, the model yeast Saccharomyces cerevisiae has seen an increased number of studies related to its geographical distribution, population structure, and phenotypic diversity. However, seven additional species from the same genus have been less thoroughly studied, which has limited our understanding of the macroevolutionary events leading to the diversification of this genus over the last 20 million years. Here, we show the geographies, hosts, substrates, and phylogenetic relationships for approximately 1,800 Saccharomyces strains, covering the complete genus with unprecedented breadth and depth. We generated and analyzed complete genome sequences of 163 strains and phenotyped 128 phylogenetically diverse strains. This dataset provides insights about genetic and phenotypic diversity within and between species and populations, quantifies reticulation and incomplete lineage sorting, and demonstrates how gene flow and selection have affected traits, such as galactose metabolism. These findings elevate the genus Saccharomyces as a model to understand biodiversity and evolution in microbial eukaryotes.

Phylogenetic divergences in brown rot fungal pathogens of Monilinia species from a worldwide collection: inferences based on the nuclear versus mitochondrial genes

BACKGROUND

Phylogenetic analyses for plant pathogenic fungi explore many questions on diversities, relationships, origins, and divergences of populations from different sources such as species, host, and geography. This information is highly valuable, especially from a large global sampling, to understand the evolutionary paths of the pathogens worldwide. Monilinia fructicola and M. laxa are two important fungal pathogens of stone fruits that cause the widespread disease commonly known as brown rot. Three nuclear genes (Calmodulin, SDHA, TEF1α) and three mitochondrial genes (Cytochrome_b, NAD2, and NAD5) of the two pathogen species from a worldwide collection including five different countries from four different continents were studied in this work.

RESULTS

Both Maximum Likelihood and Bayesian approaches were applied to the data sets, and in addition, Maximum Parsimony based approaches were used for the regions having indel polymorphisms. Calmodulin, SDHA, NAD2, and NAD5 regions were found phylogenetically informative and utilized for phylogenetics of Monilinia species for the first time. Each gene region presented a set of haplotypes except Cytochrome_b, which was monomorphic. According to this large collection of two Monilinia species around the world, M. fructicola showed more diversity than M. laxa, a result that should be carefully considered, as M. fructicola is known to be a quarantine pathogen. Moreover, the other two mitochondrial genes (NAD2 and NAD5) did not have any substitution type mutations but presented an intron indel polymorphism indicating the contribution of introns as well as mobile introns to the fungal diversity and evolution. Based on the concatenated gene sets, nuclear DNA carries higher mutations and uncovers more phylogenetic clusters in comparison to the mitochondrial DNA-based data for these fungal species.

CONCLUSIONS

This study provides the most comprehensive knowledge on the phylogenetics of both nuclear and mitochondrial genes of two prominent brown rot pathogens, M. fructicola and M. laxa. Based on the regions used in this study, the nuclear genes resolved phylogenetic branching better than the mitochondrial genes and discovered new phylogenetic lineages for these species.

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.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

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 importance of complete and high-quality genome sequences in Aspergillus niger research

The possibility to sequence the entire genome of an organism revolutionized the fields of biology and biotechnology. The first genome sequence of the important filamentous fungus Aspergillus niger was obtained in 2007, 11 years after the release of the first eukaryotic genome sequence. From that moment, genomics of A. niger has seen major progresses, facilitated by the advances in the sequencing technologies and in the methodologies for gene function prediction. However, there are still challenges to face when trying to obtain complete genomes, equipped with all the repetitive sequences that they contain and without omitting the mitochondrial sequences. The aim of this perspective article is to discuss the current status of A. niger genomics and draw attention to the open challenges that the fungal community should address to move research of this important fungus forward.

Significance of CYP3A4 ∗ 1G and OPRM1 A118G Polymorphisms in Postoperative Sufentanil Analgesia in Women of Different Ethnicities

Objective. To investigate the association between CYP3A4 $\ast$ 1G and OPRM1A118G gene polymorphisms and postoperative analgesia with sufentanil in women of Zhuang ethnicity from western Guangxi, China. Methods. Forty-eight Chinese Zhuang women who underwent elective myomectomy under general anesthesia in our hospital from January 2019 to December 2020 were selected, and another 47 Chinese Han patients in the same period were selected as the control subjects. CYP3A4 $\ast$ 1G and OPRM1 A118G gene polymorphisms as well as sedation and pain scores at different time points after surgery were compared between the two groups of patients to analyze the relationship between the degree of pain and dosage of sufentanil and to analyze the effect of gene polymorphisms on the occurrence of adverse reactions. Results. The frequencies of $\ast$ 1/ $\ast$ 1G and $\ast$ 1G/ $\ast$ 1G genotypes, allele $1\ast \text{G}$ of CYP3A4 $\ast$ 1G and genotypes AA, and allele A of OPRM1 A118G in Zhuang patients were lower than those in Han patients, while the frequencies of $\ast$ 1/ $\ast$ 1, allele $\ast$ 1 of CYP3A4 $\ast$ 1G and genotypes AG, genotypes GG, and allele G of OPRM1 A118G were higher in Zhuang women. There was no significant difference in the Ramsay and VAS scores between the two groups at different time points after surgery, but the sufentanil use in Zhuang patients was higher than in Han patients at different time points after surgery. In addition, sufentanil use was highest in Zhuang patients with the $\ast$ 1/ $\ast$ 1 genotype of the CYP3A4 $\ast$ 1G gene. No significant difference was found in the incidence of adverse reactions during analgesia between the two groups. Conclusion. CYP3A4 $\ast$ 1G could be associated with postoperative sufentanil analgesia in Zhuang patients in western Guangxi and should be considered when developing personalized analgesia regimens.

Mitochondrial Genome Contributes to the Thermal Adaptation of the Oomycete Phytophthora infestans

As a vital element of climate change, elevated temperatures resulting from global warming present new challenges to natural and agricultural sustainability, such as ecological disease management. Mitochondria regulate the energy production of cells in responding to environmental fluctuation, but studying their contribution to the thermal adaptation of species is limited. This knowledge is needed to predict future disease epidemiology for ecology conservation and food security. Spatial distributions of the mitochondrial genome (mtDNA) in 405 Phytophthora infestans isolates originating from 15 locations were characterized. The contribution of MtDNA to thermal adaptation was evaluated by comparative analysis of mtDNA frequency and intrinsic growth rate, relative population differentiation in nuclear and mtDNA, and associations of mtDNA distribution with local geography climate conditions. Significant variation in frequency, intrinsic growth rate, and spatial distribution was detected in mtDNA. Population differentiation in mtDNA was significantly higher than that in the nuclear genome, and spatial distribution of mtDNA was strongly associated with local climatic conditions and geographic parameters, particularly air temperature, suggesting natural selection caused by a local temperature is the main driver of the adaptation. Dominant mtDNA grew faster than the less frequent mtDNA. Our results provide useful insights into the evolution of pathogens under global warming. Given its important role in biological functions and adaptation to local air temperature, mtDNA intervention has become an increasing necessity for future disease management. To secure ecological integrity and food production under global warming, a synergistic study on the interactive effect of changing temperature on various components of biological and ecological functions of mitochondria in an evolutionary frame is urgently needed.

The mitochondrial genome of a wild edible mushroom, Russula rosea

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

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.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts belong to diverse clades within the kingdom of fungi and are important to human everyday life. This work investigates the evolutionary relationships among these yeasts from a mitochondrial (mt) genomic perspective. A comparative study of 155 yeast mt genomes representing all major phylogenetic lineages of Saccharomycotina was performed, including genome size and content variability, intron and intergenic regions’ diversity, genetic code alterations, and syntenic variation. Findings from this study suggest that mt genome size diversity is the result of a ceaseless random process, mainly based on genetic recombination and intron mobility. Gene order analysis revealed conserved syntenic units and many occurring rearrangements, which can be correlated with major evolutionary events as shown by the phylogenetic analysis of the concatenated mt protein matrix. For the first time, molecular dating indicated a slower mt genome divergence rate in the early stages of yeast evolution, in contrast with a faster rate in the late evolutionary stages, compared to their nuclear time divergence. Genetic code reassignments of mt genomes are a perpetual process happening in many different parallel evolutionary steps throughout the evolution of Saccharomycotina. Overall, this work shows that phylogenetic studies based on the mt genome of yeasts highlight major evolutionary events.

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 analyses of complete Peronosporaceae (Oomycota) mitogenome sequences - insights into structural evolution and phylogeny

Members of the Peronosporaceae (Oomycota, Chromista), which currently consists of 25 genera and approximately 1000 recognised species, are responsible for disease on a wide range of plant hosts. Molecular phylogenetic analyses over the last two decades have improved our understanding of evolutionary relationships within Peronosporaceae. To date, 16 numbered and three named clades have been recognised; it is clear from these studies that the current taxonomy does not reflect evolutionary relationships. Whole organelle genome sequences are an increasingly important source of phylogenetic information, and in this study we present comparative and phylogenetic analyses of mitogenome sequences from 15 of the 19 currently recognized clades of Peronosporaceae, including 44 newly assembled sequences. Our analyses suggest strong conservation of mitogenome size and gene content across Peronosporaceae but, as previous studies have suggested, limited conservation of synteny. Specifically, we identified 28 distinct syntenies amongst the 71 examined isolates. Moreover, 19 of the isolates contained inverted or direct repeats, suggesting repeated sequences may be more common than previously thought. In terms of phylogenetic relationships, our analyses of 34 concatenated mitochondrial gene sequences resulted in a topology that was broadly consistent with previous studies. However, unlike previous studies concatenated mitochondrial sequences provided strong support for higher level relationships within the family.

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.

The First Mitochondrial Genome of Ciborinia camelliae and Its Position in the Sclerotiniaceae Family

Ciborinia camelliae is the causal agent of camellia flower blight (CFB). It is a hemibiotrophic pathogen, inoperculate Discomycete of the family Sclerotiniaceae. It shows host and organ specificity infecting only flowers of species belonging to the genus Camellia, causing serious damage to the ornamental component of the plant. In this work, the first mitochondrial genome of Ciborinia camellia is reported. The mitogenome was obtained by combining Illumina short read and Nanopore long read technology. To resolve repetitive elements, specific primers were designed and used for Sanger sequencing. The manually curated mitochondrial DNA (mtDNA) of the Italian strain DSM 112729 is a circular sequence of 114,660 bp, with 29.6% of GC content. It contains two ribosomal RNA genes, 33 transfer RNAs, one RNase P gene, and 62 protein-coding genes. The latter include one gene coding for a ribosomal protein (rps3) and the 14 typical proteins involved in the oxidative metabolism. Moreover, a partial mtDNA assembled from a contig list was obtained from the deposited genome assembly of a New Zealand strain of C. camelliae. The present study contributes to understanding the mitogenome arrangement and the evolution of this phytopathogenic fungus in comparison to other Sclerotiniaceae species and confirms the usefulness of mitochondrial analysis to define phylogenetic positioning of this newly sequenced species.

Sensitivity of the U.S. Wheat Powdery Mildew Population to Quinone Outside Inhibitor Fungicides and Determination of the Complete Blumeria graminis f. sp. tritici Cytochrome b Gene

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici, is managed primarily with cultivar resistance and foliar fungicides. Quinone outside inhibitors (QoIs), which target the mitochondrial cytochrome b (cytb) gene, are one of the two main fungicide classes used on wheat. While European populations of B. graminis f. sp. tritici are widely insensitive to QoIs, largely because of the cytb mutation G143A, the QoI sensitivity of the U.S. B. graminis f. sp. tritici population had never been evaluated despite years of QoI use on U.S. wheat. A total of 381 B. graminis f. sp. tritici isolates from 15 central and eastern U.S. states were screened for sensitivity to QoI fungicides pyraclostrobin and picoxystrobin. A modest range of sensitivities was observed, with maximum resistance factors of 11.2 for pyraclostrobin and 5.3 for picoxystrobin. The F129L, G137R, and G143A cytb mutations were not detected in the U.S. B. graminis f. sp. tritici population, nor were mutations identified in the PEWY loop, a key part of the Qo site. Thus, no genetic basis for the observed quantitative variation in QoI sensitivity of U.S. B. graminis f. sp. tritici was identified. Isolate sporulation was weakly negatively associated with reduced QoI sensitivity, suggesting a fitness cost. In the course of the study, the complete B. graminis f. sp. tritici cytb gene sequence was determined for the first time in the isolate 96224 v. 3.16 reference genome. Contrary to previous reports, the gene has an intron that appears to belong to intron group II, which is unusual in fungi. The study was the first QoI sensitivity screening of a large, geographically diverse set of U.S. B. graminis f. sp. tritici isolates, and while the population as a whole remains relatively sensitive, some quantitative loss of efficacy was observed.

Population genomic analyses reveal evidence for limited recombination in the superbug Candida auris in nature

Candida auris is a recently emerged, multidrug-resistant pathogenic yeast capable of causing a diversity of human infections worldwide. Genetic analyses based on whole-genome sequences have clustered strains in this species into five divergent clades, with each clade containing limited genetic variation and one of two mating types, MTLa or MTLα. The patterns of genetic variations suggest simultaneous emergence and clonal expansion of multiple clades of this pathogen across the world. At present, it is unclear whether recombination has played any role during the evolution of C. auris. In this study, we analyzed patterns of associations among single nucleotide polymorphisms in both the nuclear and the mitochondrial genomes of 1,285 strains to investigate potential signatures of recombination in natural C. auris populations. Overall, we found that polymorphisms in the nuclear and mitochondrial genomes clustered the strains similarly into the five clades, consistent with a lack of evidence for recombination among the clades after their divergence. However, variable percentages of SNP pairs showed evidence of phylogenetic incompatibility and linkage equilibrium among samples in both the nuclear and the mitochondrial genomes, with the percentages higher in the total population than those within individual clades. Our results are consistent with limited but greater frequency of recombination before the divergence of the clades than afterwards. SNPs at loci related to antifungal resistance showed frequencies of recombination similar to or lower than those observed for SNPs in other parts of the genome. Together, though very limited, evidence for the observed recombination for both before and after the divergence of the clades suggests the possibility for continuous genetic exchange in natural populations of this important yeast pathogen.

Mitonuclear Genetic Interactions in the Basidiomycete Heterobasidion parviporum Involve a Non-conserved Mitochondrial Open Reading Frame

The mitochondrial and nuclear genomes of Eukaryotes are inherited separately and consequently follow distinct evolutionary paths. Nevertheless, the encoding of many mitochondrial proteins by the nuclear genome shows the high level of integration they have reached, which makes mitonuclear genetic interactions all the more conceivable. For each species, natural selection has fostered the evolution of coadapted alleles in both genomes, but a population-wise divergence of such alleles could lead to important phenotypic variation, and, ultimately, to speciation. In this study in the Basidiomycete Heterobasidion parviporum, we have investigated the genetic basis of phenotypic variation among laboratory-designed heterokaryons carrying the same pair of haploid nuclei, but a different mitochondrial genome. Radial growth rate data of thirteen unrelated homokaryotic parents and of their heterokaryotic offspring were combined with SNP data extracted from parental genome sequences to identify nuclear and mitochondrial loci involved in mitonuclear interactions. Two nuclear loci encoding mitochondrial proteins appeared as best candidates to engage in a genetic interaction affecting radial growth rate with a non-conserved mitochondrial open reading frame of unknown function and not reported apart from the Russulales order of Basidiomycete fungi. We believe our approach could be useful to investigate several important traits of fungal biology where mitonuclear interactions play a role, including virulence of fungal pathogens.

Comparative Mitogenomics of Fungal Species in Stachybotryaceae Provides Evolutionary Insights into Hypocreales

Stachybotrys chartarum is one of the world’s ten most feared fungi within the family Stachybotryaceae, although to date, not a single mitogenome has been documented for Stachybotryaceae. Herein, six mitogenomes of four different species in Stachybotryaceae are newly reported. The S. chartarum mitogenome was 30.7 kb in length and contained two introns (one each in rnl and cox1). A comparison of the mitogenomes of three different individuals of S. chartarum showed few nucleotide variations and conservation of gene content/order and intron insertion. A comparison of the mitogenomes of four different Stachybotryaceae species (Memnoniella echinata, Myrothecium inundatum, S. chartarum, and S. chlorohalonata), however, revealed variations in intron insertion, gene order/content, and nad2/nad3 joining pattern. Further investigations on all Hypocreales species with available mitogenomes showed greater variabilities in gene order (six patterns) and nad2/nad3 joining pattern (five patterns) although a dominant pattern always existed in each case. Ancestral state estimation showed that in each case the dominant pattern was always more ancestral than those rare patterns. Phylogenetic analyses based on mitochondrion-encoded genes supported the placement of Stachybotryaceae in Hypocreales. The crown age of Stachybotryaceae was estimated to be approximately the Early Cretaceous (141–142 Mya). This study greatly promotes our understanding of the evolution of fungal species in Hypocreales.

Mitochondrial Genomics of Six Cacao Pathogens From the Basidiomycete Family Marasmiaceae

Thread blight disease has recently been described as an emerging disease on cacao (Theobroma cacao) in Ghana. In Ghana, thread blight disease is caused by multiple species of the Marasmiaceae family: Marasmius tenuissimus, M. crinis-equi, M. palmivorus, and Marasmiellus scandens. Interestingly, two additional members of the Marasmiaceae; Moniliophthora roreri (frosty pod rot) and Moniliophthora perniciosa (witches’ broom disease), are major pathogens of cacao in the Western hemisphere. It is important to accurately characterize the genetic relationships among these economically important species in support of their disease management. We used data from Illumina NGS-based genome sequencing efforts to study the mitochondrial genomes (mitogenomes) of the four cacao thread blight associated pathogens from Ghana and compared them with published mitogenomes of Mon. roreri and Mon. perniciosa. There is a remarkable interspecies variation in mitogenome size within the six cacao-associated Marasmiaceae species, ranging from 43,121 to 109,103 bp. The differences in genome lengths are primarily due to the number and lengths of introns, differences in intergenic space, and differences in the size and numbers of unidentified ORFs (uORF). Among seven M. tenuissimus mitogenomes sequenced, there is variation in size and sequence pointing to divergent evolution patterns within the species. The intronic regions show a high degree of sequence variation compared to the conserved sequences of the 14 core genes. The intronic ORFs identified, regardless of species, encode GIY-YIG or LAGLIDADG domain-containing homing endonuclease genes. Phylogenetic relationships using the 14 core proteins largely mimic the phylogenetic relationships observed in gene order patterns, grouping M. tenuissimus with M. crinis-equi, and M. palmivorus with Mon. roreri and Mon. perniciosa, leaving Mar. scandens as an outlier. The results from this study provide evidence of independent expansion/contraction events and sequence diversification in each species and establish a foundation for further exploration of the evolutionary trajectory of the fungi in Marasmiaceae family.

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.

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

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

Dynamic evolution of eukaryotic mitochondrial and nuclear genomes: a case study in the gourmet pine mushroom Tricholoma matsutake

Fungi, as eukaryotic organisms, contain two genomes, the mitochondrial genome and the nuclear genome, in their cells. How the two genomes evolve and correlate to each other is debated. Herein, taking the gourmet pine mushroom Tricholoma matsutake as an example, we performed comparative mitogenomic analysis using samples collected from diverse locations and compared the evolution of the two genomes. The T. matsutake mitogenome encodes 49 genes and is rich of repetitive and non-coding DNAs. Six genes were invaded by up to 11 group I introns, with one cox1 intron cox1P372 showing presence/absence dynamics among different samples. Bioinformatic analyses suggested limited or no evidence of mitochondrial heteroplasmy. Interestingly, hundreds of mitochondrial DNA fragments were found in the nuclear genome, with several larger than 500 nt confirmed by PCR assays and read count comparisons, indicating clear evidence of transfer of mitochondrial DNA into the nuclear genome. Nuclear DNA of T. matsutake showed a higher mutation rate than mitochondrial DNA. Furthermore, we found evidence of incongruence between phylogenetic trees derived from mitogenome and nuclear DNA sequences. Together, our results reveal the dynamic genome evolution of the gourmet pine mushroom.

Characterization of mitogenomes from four Mucorales species and insights into pathogenicity

BACKGROUND

Mucorales, as one major order of Zygomycetes fungi, can infect human beings and cause serious consequence. We have noticed the pathogenicity of Mucorales is closely related to energy metabolism, while mitochondria play the role of energy factories in almost all biological activities.

METHODS

Virulence of M irregularis, M hiemalis, L corymbifera and R arrhizus were verified in Galleria mellonella larvae, as well as mitochondrial gene copies analysed with RT-qPCR. Mitogenomes of the four Mucorales species were sequenced based on illumina NovaSeq technology to study their characteristic features and functional regions.

RESULTS

Variant virulence of M irregularis, M hiemalis, L corymbifera and R arrhizu were verified by clinical retrospective data and our G mellonella infection models, also copies of mitochondrial genes indicated the significant associations with pathogenicity. A total of 274.18 clean reads were generated to be assembled; the complete mitogenomes of the four Mucorales species were obtained with totally different length. After the genomes annotated and compared, M irregularis was found more similar with M hiemalis than those of L corymbifera and R arrhizus, especially the small (rrns) and large (rrnl) subunits of mitochondrial ribosomal RNA (rRNA) genes. The GC content, ncRNAs and the distribution of the SNPs and InDels were also compared, and the GC content rate of fungi seems to be related to the fungal thermal adaptability. In addition, linear mitogenomes of these four Mucorales showed diverse arrangements of orf genes and directionality of some conserved gene elements.

CONCLUSION

This study uncovered the pathogenicity variances among the four Mucorales species and the relationship between their mitogenomic features and clinical pathogenicity. Further studies like spatial structure of mitochondrial genomes and the comprehensive analysis of transcription regulation are needed.

Revisiting the Neurospora crassa mitochondrial genome

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

Mitochondrial genomes of the human pathogens Coccidioides immitis and Coccidioides posadasii

Fungal mitochondrial genomes encode genes involved in crucial cellular processes, such as oxidative phosphorylation and mitochondrial translation, and the molecule has been used as a molecular marker for population genetics studies. Coccidioides immitis and C. posadasii are endemic fungal pathogens that cause coccidioidomycosis in arid regions across both American continents. To date, approximately 150 Coccidioides isolates have been sequenced to infer patterns of variation in nuclear genomes. However, less attention has been given to the mitochondrial genomes of Coccidioides. In this report, we describe the assembly and annotation of mitochondrial reference genomes for two representative strains of C. posadasii and C. immitis, as well as assess population variation among 77 selected genomes. The sizes of the circular-mapping molecules are 68.2 Kb in C. immitis and 75.1 Kb in C. posadasii. We identify 14 mitochondrial protein-coding genes common to most fungal mitochondria, which are largely syntenic across different populations and species of Coccidioides. Both Coccidioides species are characterized by a large number of group I and II introns, harboring twice the number of elements as compared to closely related Onygenales. The introns contain complete or truncated ORFs with high similarity to homing endonucleases of the LAGLIDADG and GIY-YIG families. Phylogenetic comparisons of mitochondrial and nuclear genomes show extensive phylogenetic discordance suggesting that the evolution of the two types of genetic material is not identical. This work represents the first assessment of mitochondrial genomes among isolates of both species of Coccidioides, and provides a foundation for future functional work.

Evidence for Persistent Heteroplasmy and Ancient Recombination in the Mitochondrial Genomes of the Edible Yellow Chanterelles From Southwestern China and Europe

Mitochondrial genes and genomes have patterns of inheritance that are distinctly different from those of nuclear genes and genomes. In nature, the mitochondrial genomes in eukaryotes are generally considered non-recombining and homoplasmic. If heteroplasmy and recombination exist, they are typically very limited in both space and time. Here we show that mitochondrial heteroplasmy and recombination may not be limited to a specific population nor exit only transiently in the basidiomycete Cantharellus cibarius and related species. These edible yellow chanterelles are an ecologically very important group of fungi and among the most prominent wild edible mushrooms in the Northern Hemisphere. At present, very little is known about the genetics and population biology of these fungia cross large geographical distances. Our study here analyzed a total of 363 specimens of edible yellow chanterelles from 24 geographic locations in Yunnan in southwestern China and six geographic locations in five countries in Europe. For each mushroom sample, we obtained the DNA sequences at two genes, one in the nuclear genome and one in the mitochondrial genome. Our analyses of the nuclear gene, translation elongation factor 1-alpha (tef-1) and the DNA barcode of C. cibarius and related species, suggested these samples belong to four known species and five potential new species. Interestingly, analyses of the mitochondrial ATP synthase subunit 6 (atp6) gene fragment revealed evidence of heteroplasmy in two geographic samples in Yunnan and recombination within the two new putative species in Yunnan. Specifically, all four possible haplotypes at two polymorphic nucleotide sites within the mitochondrial atp6 gene were found distributed across several geographic locations in Yunnan. Furthermore, these four haplotypes were broadly distributed across multiple phylogenetic clades constructed based on nuclear tef-1 sequences. Our results suggest that heteroplasmy and mitochondrial recombination might have happened repeatedly during the evolution of the yellow chanterelles. Together, our results suggest that the edible yellow chanterelles represent an excellent system from which to study the evolution of mitochondrial-nuclear genome relationships.

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.

Genome reduction and relaxed selection is associated with the transition to symbiosis in the basidiomycete genus Podaxis

Insights into the genomic consequences of symbiosis for basidiomycete fungi associated with social insects remain sparse. Capitalizing on viability of spores from centuries-old herbarium specimens of free-living, facultative, and specialist termite-associated Podaxis fungi, we obtained genomes of 10 specimens, including two type species described by Linnaeus >240 years ago. We document that the transition to termite association was accompanied by significant reductions in genome size and gene content, accelerated evolution in protein-coding genes, and reduced functional capacities for oxidative stress responses and lignin degradation. Functional testing confirmed that termite specialists perform worse under oxidative stress, while all lineages retained some capacity to cleave lignin. Mitochondrial genomes of termite associates were significantly larger; possibly driven by smaller population sizes or reduced competition, supported by apparent loss of certain biosynthetic gene clusters. Our findings point to relaxed selection that mirrors genome traits observed among obligate endosymbiotic bacteria of many insects.

Recent and Ongoing Horizontal Transfer of Mitochondrial Introns Between Two Fungal Tree Pathogens

Two recently introduced fungal plant pathogens (Ceratocystis lukuohia and Ceratocystis huliohia) are responsible for Rapid ‘ōhi‘a Death (ROD) in Hawai‘i. Despite being sexually incompatible, the two pathogens often co-occur in diseased ‘ōhi‘a sapwood, where genetic interaction is possible. We sequenced and annotated 33 mitochondrial genomes of the two pathogens and related species, and investigated 35 total Ceratocystis mitogenomes. Ten mtDNA regions [one group I intron, seven group II introns, and two autonomous homing endonuclease (HE) genes] were heterogeneously present in C. lukuohia mitogenomes, which were otherwise identical. Molecular surveys with specific primers showed that the 10 regions had uneven geographic distribution amongst populations of C. lukuohia. Conversely, identical orthologs of each region were present in every studied isolate of C. huliohia regardless of geographical origin. Close relatives of C. lukuohia lacked or, rarely, had few and dissimilar orthologs of the 10 regions, whereas most relatives of C. huliohia had identical or nearly identical orthologs. Each region included or worked in tandem with HE genes or reverse transcriptase/maturases that could facilitate interspecific horizontal transfers from intron-minus to intron-plus alleles. These results suggest that the 10 regions originated in C. huliohia and are actively moving to populations of C. lukuohia, perhaps through transient cytoplasmic contact of hyphal tips (anastomosis) in the wound surface of ‘ōhi‘a trees. Such contact would allow for the transfer of mitochondria followed by mitochondrial fusion or cytoplasmic exchange of intron intermediaries, which suggests that further genomic interaction may also exist between the two pathogens.

Pan-Mitogenomics Approach Discovers Diversity and Dynamism in the Prominent Brown Rot Fungal Pathogens

Monilinia fructicola and Monilinia laxa species are the most destructive and economically devastating fungal plant pathogens causing brown rot disease on stone and pome fruits worldwide. Mitochondrial genomes (mitogenomes) play critical roles influencing the mechanisms and directions of the evolution of fungal pathogens. The pan-mitogenomics approach predicts core and accessory regions of the mitochondrial genomes and explains the gain or loss of variation within and between species. The present study is a fungal pan-mitogenome of M. fructicola (N = 8) and M. laxa (N = 8) species. The completely sequenced and annotated mitogenomes showed high variability in size within and between the species. The mitogenomes of M. laxa were larger, ranging from 178,351 to 179,780bp, than the mitogenomes of M. fructicola, ranging from 158,607 to 167,838bp. However, size variation within the species showed that M. fructicola isolates were more variable in the size range than M. laxa isolates. All the mitogenomes included conserved mitochondrial genes, as well as variable regions including different mobile introns encoding homing endonucleases or maturase, non-coding introns, and repetitive elements. The linear model analysis supported the hypothesis that the mitogenome size expansion is due to presence of variable (accessory) regions. Gene synteny was mostly conserved among all samples, with the exception for order of the rps3 in the mitogenome of one isolate. The mitogenomes presented AT richness; however, A/T and G/C skew varied among the mitochondrial genes. The purifying selection was detected in almost all the protein-coding genes (PCGs) between the species. However, cytochrome b was the only gene showing a positive selection signal among the total samples. Combined datasets of amino acid sequences of 14 core mitochondrial PCGs and rps3 obtained from this study together with published mitochondrial genome sequences from some other species from Heliotales were used to infer a maximum likelihood (ML) phylogenetic tree. ML tree indicated that both Monilinia species highly diverged from each other as well as some other fungal species from the same order. Mitogenomes harbor much information about the evolution of fungal plant pathogens, which could be useful to predict pathogenic life strategies.