Enrichment of G4DNA and a Large Inverted Repeat Coincide in the Mitochondrial Genomes of Termitomyces

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

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

The mitochondrial genome of the bioluminescent fish Malacosteus niger Ayres, 1848 (Stomiidae, Actinopterygii) is large and complex, and contains an inverted-repeat structure

We determined the complete mitogenome sequence of the bioluminescent fish Malacosteus niger using long-read sequencing technologies. The 21,263 bp mitogenome features a complex structure with two copies of a 1198-bp inverted-repeat and a region of 2616-bp containing alternating copies of 16 and 26 bp repeat elements. Whole mitogenome phylogenies inferred from both nucleotide and amino-acid datasets place M. niger among Melanostomiinae. The need for additional complete mitogenome sequences from the subfamily Malacosteinae is discussed.

A genetic linkage map and improved genome assembly of the termite symbiont Termitomyces cryptogamus

BACKGROUND

The termite-fungus symbiosis is an ancient stable mutualism of two partners that reproduce and disperse independently. With the founding of each termite colony the symbiotic association must be re-established with a new fungus partner. Complementarity in the ability to break down plant substrate may help to stabilize this symbiosis despite horizontal symbiont transmission. An alternative, non-exclusive, hypothesis is that a reduced rate of evolution may contribute to stabilize the symbiosis, the so-called Red King Effect.

METHODS

To explore this concept, we produced the first linkage map of a species of Termitomyces, using genotyping by sequencing (GBS) of 88 homokaryotic offspring. We constructed a highly contiguous genome assembly using PacBio data and a de-novo evidence-based annotation. This improved genome assembly and linkage map allowed for examination of the recombination landscape and its potential effect on the mutualistic lifestyle.

RESULTS

Our linkage map resulted in a genome-wide recombination rate of 22 cM/Mb, lower than that of other related fungi. However, the total map length of 1370 cM was similar to that of other related fungi.

CONCLUSIONS

The apparently decreased rate of recombination is primarily due to genome expansion of islands of gene-poor repetitive sequences. This study highlights the importance of inclusion of genomic context in cross-species comparisons of recombination rate.

Horizontal transfer of tRNA genes to mitochondrial plasmids facilitates gene loss from fungal mitochondrial DNA

Fungal and plant mitochondria are known to exchange DNA with retroviral plasmids. Transfer of plasmid DNA to the organellar genome is best known and occurs through wholesale insertion of the plasmid. Less well known is the transfer of organellar DNA to plasmids, in particular tRNA genes. Presently, it is unknown whether fungal plasmids can adopt mitochondrial functions such as tRNA production through horizontal gene transfer. In this paper, we studied the exchange of DNA between fungal linear plasmids and fungal mtDNA, mainly focusing on the basidiomycete family Lyophyllaceae. We report at least six independent transfers of complete tRNA genes to fungal plasmids. Furthermore, we discovered two independent cases of loss of a tRNA gene from a fungal mitochondrial genome following transfer of such a gene to a linear mitochondrial plasmid. We propose that loss of a tRNA gene from mtDNA following its transfer to a plasmid creates a mutualistic dependency of the host mtDNA on the plasmid. We also find that tRNA genes transferred to plasmids encode codons that occur at the lowest frequency in the host mitochondrial genomes, possibly due to a higher number of unused transcripts. We discuss the potential consequences of mtDNA transfer to plasmids for both the host mtDNA and the plasmid.

Genetic population structure of the agaric Blastosporella zonata (Lyophyllaceae) reveals cryptic species and different roles for sexual and asexual spores in dispersal

Blastosporella zonata is one of the few basidiomycete fungi that produce asexual spores (conidia) on the mushroom. The role of these conidia in the fungal lifecycle is not known. We tested whether conidia are being utilized in local dispersal by looking for signatures of clonality in 21 samples from three localities separated by about three kilometres in Murillo, Colombia. To identify clonally related individuals, we sequenced three polymorphic markers at two unlinked loci (nuclear rRNA: ITS and LSU, and TEFIα) for all collections plus three herbarium samples. We identified two sets of clonally related individuals growing closely together in one of the three localities, and only one pair shared between localities. In all three localities we observed multiple non-clonally related dikaryons showing that sexual reproduction is also important. Our results indicate that the conidia on the mushroom are primarily important for local dispersal. Unexpectedly, our results also indicate two reproductively isolated populations, possibly representing cryptic biological species. Citation: Van de Peppel LJJ, Baroni TJ, Franco-Molano AE, et al. 2022. Genetic population structure of the agaric Blastosporella zonata (Lyophyllacea) reveals cryptic species and different roles for sexual and asexual spores in dispersal. Persoonia 49: 195-200. https://doi.org/10.3767/persoonia.2022.49.06.

Four new genera and six new species of lyophylloid agarics (Agaricales, Basidiomycota) from three different continents

Four new genera encompassing six new species are described in the Lyophyllaceae (Agaricales): the monotypic Australocybe from Australia; the monotypic Phaeotephrocybe from Belize; the monotypic Nigrocarnea from Laos containing the novel conidia-producing species N. radicata and Praearthromyces containing two Asian taxa, the conidia-producing P. corneri known from Malaysia and Singapore and P. griseus from Thailand, which lacks conidial production. In addition, we describe Arthromyces glabriceps, which is the first species in the genus Arthromyces that lacks conidial production on the basidiome. Alternative options for a classification of these lyophylloid taxa are discussed and based on morphological, ecological and biogeographical considerations rejected.

A chromosomal-level reference genome of the widely utilized Coccidioides posadasii laboratory strain "Silveira"

Coccidioidomycosis is a common fungal disease that is endemic to arid and semi-arid regions of both American continents. Coccidioides immitis and Coccidioides posadasii are the etiological agents of the disease, also known as Valley Fever. For several decades, the C. posadasii strain Silveira has been used widely in vaccine studies, is the source strain for production of diagnostic antigens, and is a widely used experimental strain for functional studies. In 2009, the genome was sequenced using Sanger sequencing technology, and a draft assembly and annotation were made available. In this study, the genome of the Silveira strain was sequenced using single molecule real-time sequencing PacBio technology, assembled into chromosomal-level contigs, genotyped, and the genome was reannotated using sophisticated and curated in silico tools. This high-quality genome sequencing effort has improved our understanding of chromosomal structure, gene set annotation, and lays the groundwork for identification of structural variants (e.g. transversions, translocations, and copy number variants), assessment of gene gain and loss, and comparison of transposable elements in future phylogenetic and population genomics studies.

Comparative analysis of Malassezia furfur mitogenomes and the development of a mitochondria-based typing approach

Malassezia furfur is a yeast species belonging to Malasseziomycetes, Ustilaginomycotina and Basidiomycota that is found on healthy warm-blooded animal skin, but also involved in various skin disorders like seborrheic dermatitis/dandruff and pityriasis versicolor. Moreover, Malassezia are associated with bloodstream infections, Crohn's disease and pancreatic carcinoma. Recent advances in Malassezia genomics and genetics have focused on the nuclear genome. In this work, we present the M. furfur mitochondrial (mt) genetic heterogenicity with full analysis of 14 novel and six available M. furfur mt genomes. The mitogenome analysis reveals a mt gene content typical for fungi, including identification of variable mt regions suitable for intra-species discrimination. Three of them, namely the trnK–atp6 and cox3–nad3 intergenic regions and intron 2 of the cob gene, were selected for primer design to identify strain differences. Malassezia furfur strains belonging to known genetic variable clusters, based on AFLP and nuclear loci, were assessed for their mt variation using PCR amplification and sequencing. The results suggest that these mt regions are excellent molecular markers for the typing of M. furfur strains and may provide added value to nuclear regions when assessing evolutionary relationships at the intraspecies level.

Ancestral predisposition toward a domesticated lifestyle in the termite-cultivated fungus Termitomyces

The ancestor of termites relied on gut symbionts for degradation of plant material, an association that persists in all termite families.^1,2 However, the single-lineage Macrotermitinae has additionally acquired a fungal symbiont that complements digestion of food outside the termite gut.³ Phylogenetic analysis has shown that fungi grown by these termites form a clade-the genus Termitomyces-but the events leading toward domestication remain unclear.⁴ To address this, we reconstructed the lifestyle of the common ancestor of Termitomyces using a combination of ecological data with a phylogenomic analysis of 21 related non-domesticated species and 25 species of Termitomyces. We show that the closely related genera Blastosporella and Arthromyces also contain insect-associated species. Furthermore, the genus Arthromyces produces asexual spores on the mycelium, which may facilitate insect dispersal when growing on aggregated subterranean fecal pellets of a plant-feeding insect. The sister-group relationship between Arthromyces and Termitomyces implies that insect association and asexual sporulation, present in both genera, preceded the domestication of Termitomyces and did not follow domestication as has been proposed previously. Specialization of the common ancestor of these two genera on an insect-fecal substrate is further supported by similar carbohydrate-degrading profiles between Arthromyces and Termitomyces. We describe a set of traits that may have predisposed the ancestor of Termitomyces toward domestication, with each trait found scattered in related taxa outside of the termite-domesticated clade. This pattern indicates that the origin of the termite-fungus symbiosis may not have required large-scale changes of the fungal partner.

Characterization of G4 DNA formation in mitochondrial DNA and their potential role in mitochondrial genome instability

Mitochondria possess their own genome which can be replicated independently of nuclear DNA. Mitochondria being the powerhouse of the cell produce reactive oxygen species, due to which the mitochondrial genome is frequently exposed to oxidative damage. Previous studies have demonstrated an association of mitochondrial deletions to aging and human disorders. Many of these deletions were present adjacent to non-B DNA structures. Thus, we investigate noncanonical structures associated with instability in mitochondrial genome. In silico studies revealed the presence of > 100 G-quadruplex motifs (of which 5 have the potential to form 3-plate G4 DNA), 23 inverted repeats, and 3 mirror repeats in the mitochondrial DNA (mtDNA). Further analysis revealed that among the deletion breakpoints from patients with mitochondrial disorders, majority are located at G4 DNA motifs. Interestingly, ~ 50% of the deletions were at base-pair positions 8271-8281, ~ 35% were due to deletion at 12362-12384, and ~ 12% due to deletion at 15516-15545. Formation of 3-plate G-quadruplex DNA structures at mitochondrial fragile regions was characterized using electromobility shift assay, circular dichroism (CD), and Taq polymerase stop assay. All 5 regions could fold into both intramolecular and intermolecular G-quadruplex structures in a KCl-dependent manner. G4 DNA formation was in parallel orientation, which was abolished in the presence of LiCl. The formation of G4 DNA affected both replication and transcription. Finally, immunolocalization of BG4 with MitoTracker confirmed the formation of G-quadruplex in mitochondrial genome. Thus, we characterize the formation of 5 different G-quadruplex structures in human mitochondrial region, which may contribute toward formation of mitochondrial deletions.

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

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

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

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

Large inverted repeats identified by intra-specific comparison of mitochondrial genomes provide insights into the evolution of Agrocybe aegerita

Genomic structure and content of Agrocybe aegerita mitochondrial DNA contain essential information regarding the evolution of this gourmet mushroom. In this study, eight isolates of A. aegerita were sequenced and assembled into complete mitochondrial genomes. The mtDNA of the isolate Ag0067 contained two genotypes, both of which were quadripartite architecture consisting of two identical inverted repeats, separated by a small single-copy region and a large single-copy region. The only difference was opposite directions of the small single-copy region. The mtDNAs ranged from 116,329 bp to 134,035 bp, harboring two large identical inverted repeats. Genes of plasmid-origin were present in regions flanked by inverted repeat ID2. Most of the core genes evolved at a relatively low rate, whereas five tRNA genes located in corresponding regions of Ag0002:1-14000 and Ag0002:50001-61000 showed higher diversity. A long fragment inversion (10 Kb) was suggested to have occurred during the differentiation of two main clades, leading to two different gene orders. The number and distribution of the introns varied greatly among the A. aegerita mtDNAs. Fast invasion of short insertions likely resulted in the diversity of introns as well as other non-coding regions, increasing the variation of the mtDNAs. We raised a model about the evolution of the large repeats to explain the unusual features of A. aegerita mtDNAs. This study constructed quadripartite architecture of A. aegerita mtDNAs analogous to chloroplast DNA, proposed an interconversion model of the divergent mitochondrial genotypes with large inverted repeats. The findings could increase our knowledge of fungal evolution.