Comparative Genomics and Transcriptomics To Analyze Fruiting Body Development in Filamentous Ascomycetes

Evidence for the Existence of Mating Subtypes Within the Schizophyllum commune: Mating Behavior and Genetic Divergence

Schizophyllum commune, a Basidiomycota fungus with a tetrapolar mating system, serves as a key model for studying sexual reproduction. In this study, two distinct mating subtypes (I and II) were identified in strain 20R-7-ZF01, isolated from subseafloor sediment, which exhibited eight different mating interaction phenotypes. Intra-subtypes exhibited colony-symmetric tetrapolar interactions (G1), whereas inter-subtype crosses yielded colony-asymmetric phenotypes (G2) and a reduced number of fruiting bodies. Nuclear migration analysis revealed that both subtypes follow the same sexual reproductive process, suggesting functional similarities despite the different reproductive outcomes. Gene silencing of mating-type loci identified the genes bbp2-9 and bbp2-7 within the B locus as key factors in determining mating subtype identity. Additionally, a similar pattern of mating subtype differentiation was observed in five other S. commune strains from both subseafloor and terrestrial environments. These findings highlight the genetic diversity within S. commune, challenge the classical understanding of fungal mating systems, and provide new insights into the genetic evolutionary mechanisms governing fungi with tetrapolar mating systems.

The Fungal and Protist Community as Affected by Tillage, Crop Residue Burning and N Fertilizer Application

The bacterial community in soil is often affected by agricultural practices, but how they affect protists and fungi is less documented. Soil from treatments that combined different N fertilizer application rates, tillage and crop residue management was sampled from a field trial started by the International Maize and Wheat Improvement Center (CIMMYT) at the 'Campo Experimental Norman E. Borlaug' (CENEB) in the Yaqui Valley in the northwest of Mexico in the early 1990s, and the fungal and protist community determined. Tillage, residue burning, and N fertilizer application had no significant effect on the fungal and protists alpha diversity expressed as Hill numbers and no significant effect on the fungal and protist community structure considering all species. The relative abundance of plant pathogens and undefined saprotrophs as determined with FUNGuildR increased significantly with tillage, while that of dung-plant and dung-soil saprotroph, and plant pathogens by burning (P < 0.05). It was found that the protists and fungal community structures were not altered by different agricultural practices, but some fungal guilds were, i.e., plant pathogens and saprotrophs, which might affect soil organic matter decomposition, nutrient cycling and crop growth.

VE-1 regulation of MAPK signaling controls sexual development in Neurospora crassa

Sexual reproduction in fungi allows genetic recombination and increases genetic diversity, allowing adaptation and survival. The velvet complex is a fungal-specific protein assembly that regulates development, pathogenesis, and secondary metabolism in response to environmental cues, such as light. In Neurospora crassa, this complex comprises VE-1, VE-2, and LAE-1. Deletion of ve-1 or ve-2, but not lae-1, leads to increased conidiation (asexual spore formation) and reduced sexual development. Mutants lacking ve-1 and/or ve-2 are female sterile and male fertile, indicating that a VE-1/VE-2 complex regulates the development of female structures. During sexual development, we observed differential regulation of 2,117 genes in dark and 4,364 genes in light between the wild type and the ∆ve-1 strain. The pheromone response and cell wall integrity pathways were downregulated in the ∆ve-1 mutant, especially in light. Additionally, we found reduced levels of both total and phosphorylated MAK-1 and MAK-2 kinases. In vitro experiments demonstrated the binding of VE-1 and VE-2 to the promoters of mak-1 and mak-2, suggesting a direct regulatory role of VE-1/VE-2 in the transcriptional control of MAPK genes to regulate sexual development. Deletion of the photosensor gene white-collar 1 prevented the light-dependent inhibition of sexual development in the ∆ve-1 mutant by increasing transcription of the pheromone response and cell wall integrity pathway genes to the levels in the dark. Our results support the proposal that the regulation of the MAP kinase pathways by the VE-1/VE-2 complex is a key element in transcriptional regulation that occurs during sexual development.

IMPORTANCE

Sexual reproduction generates new gene combinations and novel phenotypic traits and facilitates evolution. Induction of sexual development in fungi is often regulated by environmental conditions, such as the presence of light and nutrients. The velvet protein complex coordinates internal cues and environmental signals to regulate development. We have found that VE-1, a component of the velvet complex, regulates transcription during sexual development in the fungus Neurospora crassa. VE-1 regulates the transcription of many genes, including those involved in mitogen-activated protein kinase (MAPK) signaling pathways that are essential in the regulation of sexual development, and regulates the activity of the MAPK pathway. Our findings provide valuable insights into how fungi respond to environmental signals and integrate them into their reproductive processes.

Whole-genome sequencing and evolutionary analysis of the wild edible mushroom, Morchella eohespera

Morels (Morchella, Ascomycota) are an extremely desired group of edible mushrooms with worldwide distribution. Morchella eohespera is a typical black morel species, belonging to the Elata clade of Morchella species. The biological and genetic studies of this mushroom are rare, largely hindering the studies of molecular breeding and evolutionary aspects. In this study, we performed de novo sequencing and assembly of the M. eohespera strain m200 genome using the third-generation nanopore sequencing platform. The whole-genome size of M. eohespera was 53.81 Mb with a contig N50 of 1.93 Mb, and the GC content was 47.70%. A total of 9,189 protein-coding genes were annotated. Molecular dating showed that M. eohespera differentiated from its relative M. conica at ~19.03 Mya (million years ago) in Burdigalian. Evolutionary analysis showed that 657 gene families were contracted and 244 gene families expanded in M. eohespera versus the related morel species. The non-coding RNA prediction results showed that there were 336 tRNAs, 76 rRNAs, and 45 snRNAs in the M. eohespera genome. Interestingly, there was a high degree of repetition (20.93%) in the M. eohespera genome, and the sizes of long interspersed nuclear elements, short interspersed nuclear elements, and long terminal repeats were 0.83 Mb, 0.009 Mb, and 4.56 Mb, respectively. Additionally, selection pressure analysis identified that a total of 492 genes in the M. eohespera genome have undergone signatures of positive selection. The results of this study provide new insights into the genome evolution of M. eohespera and lay the foundation for in-depth research into the molecular biology of the genus Morchella in the future.

Mutations in Podospora anserina MCM1 and VelC Trigger Spontaneous Development of Barren Fruiting Bodies

The ascomycete Podospora anserina is a heterothallic filamentous fungus found mainly on herbivore dung. It is commonly used in laboratories as a model system, and its complete life cycle lasting eight days is well mastered in vitro. The main objective of our team is to understand better the global process of fruiting body development, named perithecia, induced normally in this species by fertilization. Three allelic mutants, named pfd3, pfd9, and pfd23 (for "promoting fruiting body development") obtained by UV mutagenesis, were selected in view of their abilities to promote barren perithecium development without fertilization. By complete genome sequencing of pfd3 and pfd9, and mutant complementation, we identified point mutations in the mcm1 gene as responsible for spontaneous perithecium development. MCM1 proteins are MADS box transcription factors that control diverse developmental processes in plants, metazoans, and fungi. We also identified using the same methods a mutation in the VelC gene as responsible for spontaneous perithecium development in the vacua mutant. The VelC protein belongs to the velvet family of regulators involved in the control of development and secondary metabolite production. A key role of MCM1 and VelC in coordinating the development of P. anserina perithecia with gamete formation and fertilization is highlighted.

Adaptive advantages of restorative RNA editing in fungi for resolving survival-reproduction trade-offs

RNA editing in various organisms commonly restores RNA sequences to their ancestral state, but its adaptive advantages are debated. In fungi, restorative editing corrects premature stop codons in pseudogenes specifically during sexual reproduction. We characterized 71 pseudogenes and their restorative editing in Fusarium graminearum, demonstrating that restorative editing of 16 pseudogenes is crucial for germ tissue development in fruiting bodies. Our results also revealed that the emergence of premature stop codons is facilitated by restorative editing and that premature stop codons corrected by restorative editing are selectively favored over ancestral amino acid codons. Furthermore, we found that ancestral versions of pseudogenes have antagonistic effects on reproduction and survival. Restorative editing eliminates the survival costs of reproduction caused by antagonistic pleiotropy and provides a selective advantage in fungi. Our findings highlight the importance of restorative editing in the evolution of fungal complex multicellularity and provide empirical evidence that restorative editing serves as an adaptive mechanism enabling the resolution of genetic trade-offs.

Pezizomycetes Genomes Reveal Diverse P450 Complements Characteristic of Saprotrophic and Ectomycorrhizal Lifestyles

Cytochrome P450 monooxygenases (CYPs/P450s) are heme proteins that play a role in organisms' primary and secondary metabolism. P450s play an important role in organism adaptation since lifestyle influences P450 composition in their genome. This phenomenon is well-documented in bacteria but less so in fungi. This study observed this phenomenon where diverse P450 complements were identified in saprophytic and ectomycorrhizal Pezizomycetes. Genome-wide data mining, annotation, and phylogenetic analysis of P450s in 19 Pezizomycetes revealed 668 P450s that can be grouped into 153 P450 families and 245 P450 subfamilies. Only four P450 families, namely, CYP51, CYP61, CYP5093, and CYP6001, are conserved across 19 Pezizomycetes, indicating their important role in these species. A total of 5 saprophyte Pezizomycetes have 103 P450 families, whereas 14 ectomycorrhizal Pezizomycetes have 89 P450 families. Only 39 P450 families were common, and 50 and 64 P450 families, respectively, were unique to ectomycorrhizal and saprophytic Pezizomycetes. These findings suggest that the switch from a saprophytic to an ectomycorrhizal lifestyle led to both the development of diverse P450 families as well as the loss of P450s, which led to the lowest P450 family diversity, despite the emergence of novel P450 families in ectomycorrhizal Pezizomycetes.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

The Role of Chromatin and Transcriptional Control in the Formation of Sexual Fruiting Bodies in Fungi

Fungal fruiting bodies are complex, three-dimensional structures that arise from a less complex vegetative mycelium. Their formation requires the coordinated action of many genes and their gene products, and fruiting body formation is accompanied by major changes in the transcriptome. In recent years, numerous transcription factor genes as well as chromatin modifier genes that play a role in fruiting body morphogenesis were identified, and through research on several model organisms, the underlying regulatory networks that integrate chromatin structure, gene expression, and cell differentiation are becoming clearer. This review gives a summary of the current state of research on the role of transcriptional control and chromatin structure in fruiting body development. In the first part, insights from transcriptomics analyses are described, with a focus on comparative transcriptomics. In the second part, examples of more detailed functional characterizations of the role of chromatin modifiers and/or transcription factors in several model organisms (Neurospora crassa, Aspergillus nidulans, Sordaria macrospora, Coprinopsis cinerea, and Schizophyllum commune) that have led to a better understanding of regulatory networks at the level of chromatin structure and transcription are discussed.

Genetic Diversity of Phymatotrichopsis omnivora Based on Mating Type and Microsatellite Markers Reveals Heterothallic Mating System

Phymatotrichopsis omnivora is a member of Pezizomycetes and causes root rot disease on a broad range of dicotyledonous plants. Using recently generated draft genome sequence data from four P. omnivora isolates, we developed simple sequence repeat (SSR) markers and identified both mating type genes (MAT1-1-1 and MAT1-2-1) in this fungus. To understand the genetic diversity of P. omnivora isolates (n = 43) and spore mats (n = 29) collected from four locations (Oklahoma, Texas, Arizona, and Mexico) and four host crops (cotton, alfalfa, peach, and soybean), we applied 24 SSR markers and showed that of the 72 P. omnivora isolates and spore mats tested, 41 were distinct genotypes. Furthermore, the developed SSR markers did not show cross-transferability to other close relatives of P. omnivora in the class Pezizomycetes. A multiplex PCR detecting both mating type idiomorphs and a reference gene (TUB2) was developed to screen P. omnivora isolates. Based on the dataset we tested, P. omnivora is a heterothallic fungus with both mating types present in the United States in a ratio close to 1:1. We tested P. omnivora spore mats obtained from spatially distinct disease rings that developed in a center-pivot alfalfa field and showed that both mating types can be present not only in the same field but also within a single spore mat. This study shows that P. omnivora has the genetic toolkit for generating sexually diverse progeny, providing impetus for future studies that focus on identifying sexual morphs in nature.

Functional characterization of the developmental genes asm2, asm3, and spt3 required for fruiting body formation in the filamentous ascomycete Sordaria macrospora

The formation of fruiting bodies is one of the most complex developmental processes in filamentous ascomycetes. It requires the development of sexual structures that give rise to meiosporangia (asci) and meiotic spores (ascospores) as well as surrounding structures for protection and dispersal of the spores. Previous studies have shown that these developmental processes are accompanied by significant changes of the transcriptome, and comparative transcriptomics of different fungi as well as the analysis of transcriptome changes in developmental mutants have aided in the identification of differentially regulated genes that are themselves involved in regulating fruiting body development. In previous analyses, we used transcriptomics to identify the genes asm2 and spt3, which result in developmental phenotypes when deleted in Sordaria macrospora. In this study, we identified another gene, asm3, required for fruiting body formation, and performed transcriptomics analyses of Δasm2, Δasm3, and Δspt3. Deletion of spt3, which encodes a subunit of the SAGA complex, results in a block at an early stage of development and drastic changes in the transcriptome. Deletion mutants of asm2 and asm3 are able to form fruiting bodies, but have defects in ascospore maturation. Transcriptomics analysis of fruiting bodies revealed a large overlap in differentially regulated genes in Δasm2 and Δasm3 compared to the wild type. Analysis of nuclear distribution during ascus development showed that both mutants undergo meiosis and postmeiotic divisions, suggesting that the transcriptomic and morphological changes might be related to defects in the morphogenesis of structural features of the developing asci and ascospores.

mRNA-seq and miRNA-seq profiling analyses reveal molecular mechanisms regulating induction of fruiting body in Ophiocordyceps sinensis

Ophiocordyceps sinensis has been a source of valuable materials in traditional Asian medicine for over two thousand years. With recent global warming and overharvest, however, the availability of these wild fungi has decreased dramatically. While fruiting body of O. sinensis has been artificially cultivated, the molecular mechanisms that govern the induction of fruiting body at the transcriptional and post-transcriptional levels are unclear. In this study, we carried out both mRNA and small RNA sequencing to identify crucial genes and miRNA-like RNAs (milRNAs) involved in the development of fruiting body. A total of 2875 differentially expressed genes (DEGs), and 71 differentially expressed milRNAs (DEMs) were identified among the mycoparasite complex, the sclerotium (ST) and the fruiting body stage. Functional enrichment and Gene Set Enrichment Analysis indicated that the ST had increased oxidative stress and energy metabolism and that mitogen-activated protein kinase signaling might induce the formation of fruiting body. Integrated analysis of DEGs and DEMs revealed that n_os_milR16, n_os_milR21, n_os_milR34, and n_os_milR90 could be candidate milRNAs that regulate the induction of fruiting body. This study provides transcriptome-wide insight into the molecular basis of fruiting body formation in O. Sinensis and identifies potential candidate genes for improving induction rate.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Genetic variation and function: revealing potential factors associated with microbial phenotypes

Innovations in sequencing technology have generated voluminous microbial and host genomic data, making it possible to detect these genetic variations and analyze the function influenced by them. Recently, many studies have linked such genetic variations to phenotypes through association or comparative analysis, which have further advanced our understanding of multiple microbial functions. In this review, we summarized the application of association analysis in microbes like Mycobacterium tuberculosis, focusing on screening of microbial genetic variants potentially associated with phenotypes such as drug resistance, pathogenesis and novel drug targets etc.; reviewed the application of additional comparative genomic or transcriptomic methods to identify genetic factors associated with functions in microbes; expanded the scope of our study to focus on host genetic factors associated with certain microbes or microbiome and summarized the recent host genetic variations associated with microbial phenotypes, including susceptibility and load after infection of HIV, presence/absence of different taxa, and quantitative traits of microbiome, and lastly, discussed the challenges that may be encountered and the apparent or potential viable solutions. Gene-function analysis of microbe and microbiome is still in its infancy, and in order to unleash its full potential, it is necessary to understand its history, current status, and the challenges hindering its development.

Golden Gate vectors for efficient gene fusion and gene deletion in diverse filamentous fungi

The cloning of plasmids can be time-consuming or expensive. Yet, cloning is a prerequisite for many standard experiments for the functional analysis of genes, including the generation of deletion mutants and the localization of gene products. Here, we provide Golden Gate vectors for fast and easy cloning of gene fusion as well as gene deletion vectors applicable to diverse fungi. In Golden Gate cloning, restriction and ligation occur simultaneously in a one-pot reaction. Our vector set contains recognition sites for the commonly used type IIS restriction endonuclease BsaI. We generated plasmids for C- as well as N-terminal tagging with GFP, mRFP and 3xFLAG. For gene deletion, we provide five different donor vectors for selection marker cassettes. These include standard cassettes for hygromycin B, nourseothricin and phleomycin resistance genes as well as FLP/FRT-based marker recycling cassettes for hygromycin B and nourseothricin resistance genes. To make cloning most feasible, we provide robust protocols, namely (1) an overview of cloning procedures described in this paper, (2) specific Golden Gate reaction protocols and (3) standard primers for cloning and sequencing of plasmids and generation of deletion cassettes by PCR and split-marker PCR. We show that our vector set is applicable for the biotechnologically relevant Penicillium chrysogenum and the developmental model system Sordaria macrospora. We thus expect these vectors to be beneficial for other fungi as well. Finally, the vectors can easily be adapted to organisms beyond the kingdom fungi.

Sordaria macrospora: 25 years as a model organism for studying the molecular mechanisms of fruiting body development

Abstract

Fruiting bodies are among the most complex multicellular structures formed by fungi, and the molecular mechanisms that regulate their development are far from understood. However, studies with a number of fungal model organisms have started to shed light on this developmental process. One of these model organisms is Sordaria macrospora, a filamentous ascomycete from the order Sordariales. This fungus has been a genetic model organism since the 1950s, but its career as a model organism for molecular genetics really took off in the 1990s, when the establishment of a transformation protocol, a mutant collection, and an indexed cosmid library provided the methods and resources to start revealing the molecular mechanisms of fruiting body development. In the 2000s, “omics” methods were added to the S. macrospora tool box, and by 2020, 58 developmental genes have been identified in this fungus. This review gives a brief overview of major method developments for S. macrospora, and then focuses on recent results characterizing different processes involved in regulating development including several regulatory protein complexes, autophagy, transcriptional and chromatin regulation, and RNA editing.

Key points

•Sordaria macrospora is a model system for analyzing fungal fruiting body development.

•More than 100 developmental mutants are available for S. macrospora.

•More than 50 developmental genes have been characterized in S. macrospora.