Frequent heteroplasmy and recombination in the mitochondrial genomes of the basidiomycete mushroom Thelephora ganbajun

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.

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.

Molecular phylogeny and morphology reveal four new species of Thelephora (Thelephorales, Basidiomycota) from subtropical China, closely related to T. ganbajun

The genus of Thelephora is a group of cosmopolitan ectomycorrhizal fungi with basidiocarps of morphological diversity that has an extremely scarce species reported from the forest ecosystem in China. In this study, phylogenetic analyses of Thelephora species from subtropical China were carried out based on multiple loci including the internal transcribed spacer (ITS) regions, the large subunit of nuclear ribosomal RNA gene (nLSU), and the small subunit of mitochondrial rRNA gene (mtSSU). Maximum likelihood and Bayesian analyses were used to construct the phylogenetic tree. The phylogenetic positions of four new species, Th. aquila, Th. glaucoflora, Th. nebula, and Th. pseudoganbajun, were revealed based on morphological and molecular evidence. Molecular analyses demonstrated that the four new species were closely related to Th. ganbajun and formed a clade with robust support in the phylogenetic tree. Regarding morphology, they share some common morphological characteristics, including flabelliform to imbricate pilei, generative hyphae more or less covered by crystals, and subglobose to irregularly lobed basidiospores (5–8 × 4–7 μm) with tuberculate ornamentation. These new species are described and illustrated and are compared to similar morphological or phylogenetically related species. A key to the new and allied species from China is provided.

Genome Assembly and Genetic Traits of the Pleuromutilin-Producer Clitopilus passeckerianus DSM1602

The gilled mushroom Clitopilus passeckerianus (Entolomataceae, Agaricales, Basidiomycota) is well known to produce the terpenoid pleuromutilin, which is the biotechnological basis for medically important antibiotics such as lefamulin and retapamulin. Their unique mode of action and good tolerance entails an increasing demand of pleuromutilin-derived antibiotics in veterinary and human health care. Surprisingly, despite their pharmaceutical importance, no genome sequence is available of any pleuromutilin-producing fungus. Here, we present the high-quality draft genome sequence of the pleuromutilin-producer C. passeckerianus DSM1602 including functional genome annotation. More precisely, we employed a hybrid assembly strategy combining Illumina sequencing and Nanopore sequencing to assemble the mitochondrial genome as well as the nuclear genome. In accordance with the dikaryotic state of the fungus, the nuclear genome has a diploid character. Interestingly, the mitochondrial genome appears duplicated. Bioinformatic analysis revealed a versatile secondary metabolism with an emphasis on terpenoid biosynthetic enzymes in C. passeckerianus and also in related strains. Two alleles of biosynthetic gene clusters for pleuromutilin were found in the genome of C. passeckerianus. The pleuromutilin genes were reassembled with yeast-specific elements for heterologous expression in Saccharomyces cerevisiae. Our work lays the foundation for metabolic strain engineering towards higher yields of the valuable compound pleuromutilin.

A new mushroom hyperaccumulator: Cadmium and arsenic in the ectomycorrhizal basidiomycete Thelephora penicillata

Fruit-bodies of six Thelephora species (Fungi, Basidiomycota, Thelephoraceae) were analyzed for their trace element concentrations. In Thelephora penicillata, extremely high concentrations of Cd and As were found, followed by highly elevated concentrations of Cu and Zn. The highest accumulation ability was found for Cd with a mean concentration of 1.17 ± 0.37 g kg^-1 (dry mass) in fruit-bodies collected from 20 unpolluted sites; the mean As concentration was 0.878 ± 0.242 g kg^-1. Furthermore, striking accumulation of Se (923 ± 28 mg kg^-1) was found in one sample of T. vialis and elevated concentrations of S were detected in T. palmata (19.6 ± 5.9 g kg^-1). The analyzed Thelephora species were sequenced and, based on the Maximum Likelihood phylogenetic analysis (ITS rDNA) of the genus, possible other Thelephora (hyper)accumulators were predicted on the basis of their phylogenetic relationship with the discovered (hyper)accumulators. The striking ability of T. penicillata to accumulate simultaneously Cd, As, Cu, and Zn has no parallel in the Fungal Kingdom and raises the question of a biological importance of metal(loid) hyperaccumulation in mushrooms.

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

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

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.

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.

Revealing the high variability on nonconserved core and mobile elements of Austropuccinia psidii and other rust mitochondrial genomes

Mitochondrial genomes are highly conserved in many fungal groups, and they can help characterize the phylogenetic relationships and evolutionary biology of plant pathogenic fungi. Rust fungi are among the most devastating diseases for economically important crops around the world. Here, we report the complete sequence and annotation of the mitochondrial genome of Austropuccinia psidii (syn. Puccinia psidii), the causal agent of myrtle rust. We performed a phylogenomic analysis including the complete mitochondrial sequences from other rust fungi. The genome composed of 93.299 bp has 73 predicted genes, 33 of which encoded nonconserved proteins (ncORFs), representing almost 45% of all predicted genes. A. psidii mtDNA is one of the largest rust mtDNA sequenced to date, most likely due to the abundance of ncORFs. Among them, 33% were within intronic regions of diverse intron groups. Mobile genetic elements invading intron sequences may have played significant roles in size but not shaping of the rust mitochondrial genome structure. The mtDNAs from rust fungi are highly syntenic. Phylogenetic inferences with 14 concatenated mitochondrial proteins encoded by the core genes placed A. psidii according to phylogenetic analysis based on 18S rDNA. Interestingly, cox1, the gene with the greatest number of introns, provided phylogenies not congruent with the core set. For the first time, we identified the proteins encoded by three A. psidii ncORFs using proteomics analyses. Also, the orf208 encoded a transmembrane protein repressed during in vitro morphogenesis. To the best of our knowledge, we presented the first report of a complete mtDNA sequence of a member of the family Sphaerophragmiacea.

Distribution of Zinc in Mycelial Cells and Antioxidant and Anti-Inflammatory Activities of Mycelia Zinc Polysaccharides from Thelephora ganbajun TG-01

This study demonstrates that Thelephora ganbajun had a strong ability to absorb zinc, and zinc can be compartmentally stored in the small vesicles and mainly accumulated in the form of zinc-enriched polysaccharides (zinc content was 25.0 ± 1.27 mg/g). Mycelia zinc polysaccharides (MZPS) and its fractions were isolated. The main fraction (MZPS-2) with the highest antioxidant activity in vitro was composed of mannose : galacturonic acid : glucose : galactose in a molar ratio of 61.19 : 1 : 39.67 : 48.67, with a weight-averaged molecular weight of 5.118 × 105 Da. MZPS-2 had both α-pyranose and β-pyranose configuration and had a triple helical conformation. By establishing zebrafish models, we found that MZPS-2 can significantly scavenge free radicals, reduce the generation of reactive oxygen species caused by inflammation, and inhibit the recruitment of neutrophils toward the injury site. Therefore, MZPS-2 exhibited antioxidant and anti-inflammatory effects and can be used as a zinc supplement with specific biological activities to alleviate zinc deficiency complications, such as chronic oxidative stress or inflammation.

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

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

Evidence of Extensive Intraspecific Noncoding Reshuffling in a 169-kb Mitochondrial Genome of a Basidiomycetous Fungus

Comparative genomics of fungal mitochondrial genomes (mitogenomes) have revealed a remarkable pattern of rearrangement between and within major phyla owing to horizontal gene transfer and recombination. The role of recombination was exemplified at a finer evolutionary time scale in basidiomycetes group of fungi as they display a diversity of mitochondrial DNA inheritance patterns. Here, we assembled mitogenomes of six species from the Hymenochaetales order of basidiomycetes and examined 59 mitogenomes from 2 genetic lineages of Phellinus noxius. Gene order is largely collinear, while intergene regions are major determinants of mitogenome size variation. Substantial sequence divergence was found in shared introns consistent with high horizontal gene transfer frequency observed in yeasts, but we also identified a rare case where an intron was retained in five species since speciation. In contrast to the hyperdiversity observed in nuclear genomes of Phellinus noxius, mitogenomes’ intraspecific polymorphisms at protein-coding sequences are extremely low. Phylogeny network based on introns revealed turnover as well as exchange of introns between two lineages. Strikingly, some strains harbor a mosaic origin of introns from both lineages. Analysis of intergenic sequence indicated substantial differences between and within lineages, and an expansion may be ongoing as a result of exchange between distal intergenes. These findings suggest that the evolution in mitochondrial DNAs is usually lineage specific but chimeric mitotypes are frequently observed, thus capturing the possible evolutionary processes shaping mitogenomes in a basidiomycete. The large mitogenome sizes reported in various basidiomycetes appear to be a result of interspecific reshuffling of intergenes.

Notes, outline and divergence times of Basidiomycota

The Basidiomycota constitutes a major phylum of the kingdom Fungi and is second in species numbers to the Ascomycota. The present work provides an overview of all validly published, currently used basidiomycete genera to date in a single document. An outline of all genera of Basidiomycota is provided, which includes 1928 currently used genera names, with 1263 synonyms, which are distributed in 241 families, 68 orders, 18 classes and four subphyla. We provide brief notes for each accepted genus including information on classification, number of accepted species, type species, life mode, habitat, distribution, and sequence information. Furthermore, three phylogenetic analyses with combined LSU, SSU, 5.8s, rpb1, rpb2, and ef1 datasets for the subphyla Agaricomycotina, Pucciniomycotina and Ustilaginomycotina are conducted, respectively. Divergence time estimates are provided to the family level with 632 species from 62 orders, 168 families and 605 genera. Our study indicates that the divergence times of the subphyla in Basidiomycota are 406–430 Mya, classes are 211–383 Mya, and orders are 99–323 Mya, which are largely consistent with previous studies. In this study, all phylogenetically supported families were dated, with the families of Agaricomycotina diverging from 27–178 Mya, Pucciniomycotina from 85–222 Mya, and Ustilaginomycotina from 79–177 Mya. Divergence times as additional criterion in ranking provide additional evidence to resolve taxonomic problems in the Basidiomycota taxonomic system, and also provide a better understanding of their phylogeny and evolution.

Flip-flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae, Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein-coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip-flop organization. The flip-flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70-gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data-rich approach (i.e., organelle genome data) from broader taxon sampling.

Clonorchis sinensis and Clonorchiasis: The Relevance of Exploring Genetic Variation

Parasitic trematodes (flukes) cause substantial mortality and morbidity in humans. The Chinese liver fluke, Clonorchis sinensis, is one of the most destructive parasitic worms in humans in China, Vietnam, Korea and the Russian Far East. Although C. sinensis infection can be controlled relatively well using anthelmintics, the worm is carcinogenic, inducing cholangiocarcinoma and causing major suffering in ~15 million people in Asia. This chapter provides an account of C. sinensis and clonorchiasis research-covering aspects of biology, epidemiology, pathogenesis and immunity, diagnosis, treatment and control, genetics and genomics. It also describes progress in the area of molecular biology (genetics, genomics, transcriptomics and proteomics) and highlights challenges associated with comparative genomics and population genetics. It then reviews recent advances in the sequencing and characterisation of the mitochondrial and nuclear genomes for a Korean isolate of C. sinensis and summarises salient comparative genomic work and the implications thereof. The chapter concludes by considering how advances in genomic and informatics will enable research on the genetics of C. sinensis and related parasites, as well as the discovery of new fluke-specific intervention targets.