Gene expression differences among three Neurospora species reveal genes required for sexual reproduction in Neurospora crassa

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

A mitochondrial NAD/NADH kinase governs fungal virulence through an oxidative stress response and arginine biosynthesis in Fusarium graminearum

NADP/NADPH plays an indispensable role in cellular metabolism, serving as a pivotal cofactor in numerous enzymatic processes involved in anabolic pathways, antioxidant defense, and the biosynthesis of essential cellular components. NAD/NADH kinases (NADKs) phosphorylate NAD/NADH, constituting the sole de novo synthetic pathway for NADP/NADPH generation. Despite the pivotal role of NADP/NADPH in cellular functions, the physiological role of NADK remains largely unexplored in filamentous fungi. In this study, we identified three putative NADKs in Fusarium graminearum-FgNadk1, FgNadk2, and FgNadk3-responsible for NAD/NADH phosphorylation. NADK-mediated formation of intracellular NADPH proved crucial for vegetative growth, sexual reproduction, and virulence. Specifically, FgNadk2, the mitochondrial NADK, played a role in oxidative stress resistance and the maintenance of mitochondrial reactive oxygen species levels. Moreover, the deletion of FgNADK2 resulted in arginine auxotrophy, contributing to the reduced fungal virulence. These findings underscore the necessity of mitochondrial NADK in fungal virulence in F. graminearum, revealing its involvement in mitochondrial redox homeostasis and the arginine biosynthetic pathway. This study provides critical insights into the interconnectedness of metabolic pathways essential for fungal growth, stress response, and pathogenicity.

Lineage-specific genes are clustered with HET-domain genes and respond to environmental and genetic manipulations regulating reproduction in Neurospora

Lineage-specific genes (LSGs) have long been postulated to play roles in the establishment of genetic barriers to intercrossing and speciation. In the genome of Neurospora crassa, most of the 670 Neurospora LSGs that are aggregated adjacent to the telomeres are clustered with 61% of the HET-domain genes, some of which regulate self-recognition and define vegetative incompatibility groups. In contrast, the LSG-encoding proteins possess few to no domains that would help to identify potential functional roles. Possible functional roles of LSGs were further assessed by performing transcriptomic profiling in genetic mutants and in response to environmental alterations, as well as examining gene knockouts for phenotypes. Among the 342 LSGs that are dynamically expressed during both asexual and sexual phases, 64% were detectable on unusual carbon sources such as furfural, a wildfire-produced chemical that is a strong inducer of sexual development, and the structurally-related furan 5-hydroxymethyl furfural (HMF). Expression of a significant portion of the LSGs was sensitive to light and temperature, factors that also regulate the switch from asexual to sexual reproduction. Furthermore, expression of the LSGs was significantly affected in the knockouts of adv-1 and pp-1 that regulate hyphal communication, and expression of more than one quarter of the LSGs was affected by perturbation of the mating locus. These observations encouraged further investigation of the roles of clustered lineage-specific and HET-domain genes in ecology and reproduction regulation in Neurospora, especially the regulation of the switch from the asexual growth to sexual reproduction, in response to dramatic environmental conditions changes.

Origins of lineage-specific elements via gene duplication, relocation, and regional rearrangement in Neurospora crassa

The origin of new genes has long been a central interest of evolutionary biologists. However, their novelty means that they evade reconstruction by the classical tools of evolutionary modelling. This evasion of deep ancestral investigation necessitates intensive study of model species within well-sampled, recently diversified, clades. One such clade is the model genus Neurospora, members of which lack recent gene duplications. Several Neurospora species are comprehensively characterized organisms apt for studying the evolution of lineage-specific genes (LSGs). Using gene synteny, we documented that 78% of Neurospora LSG clusters are located adjacent to the telomeres featuring extensive tracts of non-coding DNA and duplicated genes. Here, we report several instances of LSGs that are likely from regional rearrangements and potentially from gene rebirth. To broadly investigate the functions of LSGs, we assembled transcriptomics data from 68 experimental data points and identified co-regulatory modules using Weighted Gene Correlation Network Analysis, revealing that LSGs are widely but peripherally involved in known regulatory machinery for diverse functions. The ancestral status of the LSG mas-1, a gene with roles in cell-wall integrity and cellular sensitivity to antifungal toxins, was investigated in detail alongside its genomic neighbours, indicating that it arose from an ancient lysophospholipase precursor that is ubiquitous in lineages of the Sordariomycetes. Our discoveries illuminate a "rummage region" in the N. crassa genome that enables the formation of new genes and functions to arise via gene duplication and relocation, followed by fast mutation and recombination facilitated by sequence repeats and unconstrained non-coding sequences.

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.

Secondary Metabolism Gene Clusters Exhibit Increasingly Dynamic and Differential Expression during Asexual Growth, Conidiation, and Sexual Development in Neurospora crassa

Secondary metabolite clusters (SMCs) encode the machinery for fungal toxin production. However, understanding their function and analyzing their products requires investigation of the developmental and environmental conditions in which they are expressed. Gene expression is often restricted to specific and unexamined stages of the life cycle. Therefore, we applied comparative genomics analyses to identify SMCs in Neurospora crassa and analyzed extensive transcriptomic data spanning nine independent experiments from diverse developmental and environmental conditions to reveal their life cycle-specific gene expression patterns. We reported 20 SMCs comprising 177 genes-a manageable set for investigation of the roles of SMCs across the life cycle of the fungal model N. crassa-as well as gene sets coordinately expressed in 18 predicted SMCs during asexual and sexual growth under three nutritional and two temperature conditions. Divergent activity of SMCs between asexual and sexual development was reported. Of 126 SMC genes that we examined for knockout phenotypes, al-2 and al-3 exhibited phenotypes in asexual growth and conidiation, whereas os-5, poi-2, and pmd-1 exhibited phenotypes in sexual development. SMCs with annotated function in mating and crossing were actively regulated during the switch between asexual and sexual growth. Our discoveries call for attention to roles that SMCs may play in the regulatory switches controlling mode of development, as well as the ecological associations of those developmental stages that may influence expression of SMCs. IMPORTANCE Secondary metabolites (SMs) are low-molecular-weight compounds that often mediate interactions between fungi and their environments. Fungi enriched with SMs are of significant research interest to agriculture and medicine, especially from the aspects of pathogen ecology and environmental epidemiology. However, SM clusters (SMCs) that have been predicted by comparative genomics alone have typically been poorly defined and insufficiently functionally annotated. Therefore, we have investigated coordinate expression in SMCs in the model system N. crassa, and our results suggest that SMCs respond to environmental signals and to stress that are associated with development. This study examined SMC regulation at the level of RNA to integrate observations and knowledge of these genes in various growth and development conditions, supporting combining comparative genomics and inclusive transcriptomics to improve computational annotation of SMCs. Our findings call for detailed study of the function of SMCs during the asexual-sexual switch, a key, often-overlooked developmental stage.

Cross-species transcriptomics identifies core regulatory changes differentiating the asymptomatic asexual and virulent sexual life cycles of grass-symbiotic Epichloë fungi

Fungi from the genus Epichloë form systemic endobiotic infections of cool season grasses, producing a range of host-protective natural products in return for access to nutrients. These infections are asymptomatic during vegetative host growth, with associations between asexual Epichloë spp. and their hosts considered mutualistic. However, the sexual cycle of Epichloë spp. involves virulent growth, characterized by the envelopment and sterilization of a developing host inflorescence by a dense sheath of mycelia known as a stroma. Microscopic analysis of stromata revealed a dramatic increase in hyphal propagation and host degradation compared with asymptomatic tissues. RNAseq was used to identify differentially expressed genes in asymptomatic vs stromatized tissues from 3 diverse Epichloë-host associations. Comparative analysis identified a core set of 135 differentially expressed genes that exhibited conserved transcriptional changes across all 3 associations. The core differentially expressed genes more strongly expressed during virulent growth encode proteins associated with host suppression, digestion, adaptation to the external environment, a biosynthetic gene cluster, and 5 transcription factors that may regulate Epichloë stroma formation. An additional 5 transcription factor encoding differentially expressed genes were suppressed during virulent growth, suggesting they regulate mutualistic processes. Expression of biosynthetic gene clusters for natural products that suppress herbivory was universally suppressed during virulent growth, and additional biosynthetic gene clusters that may encode production of novel host-protective natural products were identified. A comparative analysis of 26 Epichloë genomes found a general decrease in core differentially expressed gene conservation among asexual species, and a specific decrease in conservation for the biosynthetic gene cluster expressed during virulent growth and an unusual uncharacterized gene.

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.

A putative transcription factor LFC1 negatively regulates development and yield of winter mushroom

Basidioma is the fruiting body of mushroom species. The deep understanding on the mechanism of basidioma development is valuable for mushroom breeding and cultivation. From winter mushroom (Flammulina velutipes), one of the top five industrially cultivated mushrooms, a novel putative Zn(II)2Cys6 transcription factor LFC1 with negative regulatory function in basidioma development was identified. The transcript level of lfc1 was dramatically decreased during basidioma development. Neither overexpression nor knockdown of lfc1 affected hyphal vegetative growth. However, knockdown of lfc1 could promote basidioma development and shorten cultivation time by 2 days, while overexpression of lfc1 delayed the optimal harvest time by 3 days. In the lfc1 knockdown strain, in which the lfc1 expression was reduced by 72%, mushroom yield and biological efficiency could be increased at least by 24%. Knockdown of lfc1 did not affect the shape of caps but significantly increased basidioma length and number, while its overexpression did not affect basidioma length but dramatically reduced basidioma number. In addition, rather than producing basidiomata with round caps as in wild type, the caps of basidiomata in the lfc1 overexpression mutants were significantly larger and the cap edge was wrinkled. RNA-seq analysis revealed that 455 genes had opposite transcriptional responses to lfc1 overexpression and knockdown. Some of them were previously reported as genes involved in basidioma development, including 3 hydrophobin encoding genes, 2 lectin encoding genes, FVFD16, an Eln2 ortholog encoding gene, and 3 genes encoding membrane components. As LFC1 homologs are widely present in mushroom species, lfc1 can be useful in mushroom breeding.Key Points• A novel transcription factor LFC1 negatively regulates fruiting in winter mushroom• LFC1 regulated transcription of more than 400 genes.• Reduction of LFC1 expression could shorten cultivation time and increase yield.• lfc1 could be a potentially useful reference gene for mushroom breeding.

Pathological, Morphological, Cytogenomic, Biochemical and Molecular Data Support the Distinction between Colletotrichum cigarro comb. et stat. nov. and Colletotrichum kahawae

The genus Colletotrichum has witnessed tremendous variations over the years in the number of species recognized, ranging from 11 to several hundreds. Host-specific fungal species, once the rule, are now the exception, with polyphagous behavior regarded as normal in this genus. The species Colletotrichum kahawae was created to accommodate the pathogens that have the unique ability to infect green developing coffee berries causing the devastating Coffee Berry Disease in Africa, but its close phylogenetic relationship to a polyphagous group of fungi in the C. gloeosporioides species complex led some researchers to regard these pathogens as members of a wider species. In this work we combine pathological, morphological, cytogenomic, biochemical, and molecular data of a comprehensive set of phylogenetically-related isolates to show that the Coffee Berry Disease pathogen forms a separate species, C. kahawae, and also to assign the closely related fungi, previously in C. kahawae subsp. cigarro, to a new species, C. cigarro comb. et stat. nov. This taxonomic clarification provides an opportunity to link phylogeny and functional biology, and additionally enables a much-needed tool for plant pathology and agronomy, associating exclusively C. kahawae to the Coffee Berry Disease pathogen.

A WD40 Protein Encoding Gene Fvcpc2 Positively Regulates Mushroom Development and Yield in Flammulina velutipes

Ascomycota and Basidiomycota are two closely related phyla and fungi in two phyla share some common morphological developmental process during fruiting body formation. In Neurospora crassa, the Gβ-like protein CPC-2 with a seven-WD40 repeat domain was previously reported. By transforming CPC-2 ortholog encoding genes, from 7 different fungal species across Ascomycota and Basidiomycota, into the cpc-2 deletion mutant of N. crassa, we demonstrate that all tested CPC-2 ortholog genes were able to complement the defects of the cpc-2 deletion mutant in sexual development, indicating that CPC-2 proteins from Ascomycota and Basidiomycota have the similar cellular function. Using Flammulina velutipes as a model system for mushroom species, the CPC-2 ortholog FvCPC2 was characterized. Fvcpc2 increased transcription during fruiting body development. Knockdown of Fvcpc2 by RNAi completely impaired fruiting body formation. In three Fvcpc2 knockdown mutants, transcriptional levels of genes encoding adenylate cyclase and protein kinase A catalytic subunit were significantly lower and colony growth became slower than wild type. The addition of cAMP or the PKA-activator 8-Bromo-cAMP into the medium restored the Fvcpc2 knockdown mutants to the wild-type colony growth phenotype, suggesting that the involvement of cAMP production in the regulatory mechanisms of FvCPC2. Knockdown of Fvcpc2 also weakened transcriptional responses to sexual development induction by some genes related to fruiting body development, including 4 jacalin-related lectin encoding genes, 4 hydrophobin encoding genes, and 3 functionally-unknown genes, suggesting the participation of these genes in the mechanisms by which FvCPC2 regulates fruiting body development. All three Fvcpc2 overexpression strains displayed increased mushroom yield and shortened cultivation time compared to wild type, suggesting that Fvcpc2 can be a promising reference gene for Winter Mushroom breeding. Since the orthologs of FvCPC2 were highly conserved and specifically expressed during fruiting body development in different edible mushrooms, genes encoding FvCPC2 orthologs in other mushroom species also have potential application in breeding.

Integrative Activity of Mating Loci, Environmentally Responsive Genes, and Secondary Metabolism Pathways during Sexual Development of Chaetomium globosum

Fungal diversity has amazed evolutionary biologists for decades. One societally important aspect of this diversity manifests in traits that enable pathogenicity. The opportunistic pathogen Chaetomium globosum is well adapted to a high-humidity environment and produces numerous secondary metabolites that defend it from predation. Many of these chemicals can threaten human health. Understanding the phases of the C. globosum life cycle in which these products are made enables better control and even utilization of this fungus. Among its intriguing traits is that it both is self-fertile and lacks any means of propagule-based asexual reproduction. By profiling genome-wide gene expression across the process of sexual reproduction in C. globosum and comparing it to genome-wide gene expression in the model filamentous fungus N. crassa and other closely related fungi, we revealed associations among mating-type genes, sexual developmental genes, sexual incompatibility regulators, environmentally responsive genes, and secondary metabolic pathways.

The origins and maintenance of the rich fungal diversity have been longstanding issues in evolutionary biology. To investigate how differences in expression regulation contribute to divergences in development and ecology among closely related species, transcriptomes were compared between Chaetomium globosum, a homothallic pathogenic fungus thriving in highly humid ecologies, and Neurospora crassa, a heterothallic postfire saprotroph. Gene expression was quantified in perithecia at nine distinct morphological stages during nearly synchronous sexual development. Unlike N. crassa, expression of all mating loci in C. globosum was highly correlated. Key regulators of the initiation of sexual development in response to light stimuli—including orthologs of N. crassasub-1, sub-1-dependent gene NCU00309, and asl-1—showed regulatory dynamics matching between C. globosum and N. crassa. Among 24 secondary metabolism gene clusters in C. globosum, 11—including the cochliodones biosynthesis cluster—exhibited highly coordinated expression across perithecial development. C. globosum exhibited coordinately upregulated expression of histidine kinases in hyperosmotic response pathways—consistent with gene expression responses to high humidity we identified in fellow pathogen Fusarium graminearum. Bayesian networks indicated that gene interactions during sexual development have diverged in concert with the capacities both to reproduce asexually and to live a self-compatible versus self-incompatible life cycle, shifting the hierarchical roles of genes associated with conidiation and heterokaryon incompatibility in N. crassa and C. globosum. This divergence supports an evolutionary history of loss of conidiation due to unfavorable combinations of heterokaryon incompatibility in homothallic species.

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

Many filamentous ascomycetes develop three-dimensional fruiting bodies for production and dispersal of sexual spores. Fruiting bodies are among the most complex structures differentiated by ascomycetes; however, the molecular mechanisms underlying this process are insufficiently understood. Previous comparative transcriptomics analyses of fruiting body development in different ascomycetes suggested that there might be a core set of genes that are transcriptionally regulated in a similar manner across species. Conserved patterns of gene expression can be indicative of functional relevance, and therefore such a set of genes might constitute promising candidates for functional analyses. In this study, we have sequenced the genome of the Pezizomycete Ascodesmis nigricans, and performed comparative transcriptomics of developing fruiting bodies of this fungus, the Pezizomycete Pyronema confluens, and the Sordariomycete Sordaria macrospora With only 27 Mb, the A. nigricans genome is the smallest Pezizomycete genome sequenced to date. Comparative transcriptomics indicated that gene expression patterns in developing fruiting bodies of the three species are more similar to each other than to nonsexual hyphae of the same species. An analysis of 83 genes that are upregulated only during fruiting body development in all three species revealed 23 genes encoding proteins with predicted roles in vesicle transport, the endomembrane system, or transport across membranes, and 13 genes encoding proteins with predicted roles in chromatin organization or the regulation of gene expression. Among four genes chosen for functional analysis by deletion in S. macrospora, three were shown to be involved in fruiting body formation, including two predicted chromatin modifier genes.

The Genetics and Biochemistry of Cell Wall Structure and Synthesis in Neurospora crassa, a Model Filamentous Fungus

This review discusses the wealth of information available for the N. crassa cell wall. The basic organization and structure of the cell wall is presented and how the wall changes during the N. crassa life cycle is discussed. Over forty cell wall glycoproteins have been identified by proteomic analyses. Genetic and biochemical studies have identified many of the key enzymes needed for cell wall biogenesis, and the roles these enzymes play in cell wall biogenesis are discussed. The review includes a discussion of how the major cell wall components (chitin, β-1,3-glucan, mixed β-1,3-/ β-1,4- glucans, glycoproteins, and melanin) are synthesized and incorporated into the cell wall. We present a four-step model for how cell wall glycoproteins are covalently incorporated into the cell wall. In N. crassa, the covalent incorporation of cell wall glycoproteins into the wall occurs through a glycosidic linkage between lichenin (a mixed β-1,3-/β-1,4- glucan) and a "processed" galactomannan that has been attached to the glycoprotein N-linked oligosaccharides. The first step is the addition of the galactomannan to the N-linked oligosaccharide. Mutants affected in galactomannan formation are unable to incorporate glycoproteins into their cell walls. The second step is carried out by the enzymes from the GH76 family of α-1,6-mannanases, which cleave the galactomannan to generate a processed galactomannan. The model suggests that the third and fourth steps are carried out by members of the GH72 family of glucanosyltransferases. In the third step the glucanosyltransferases cleave lichenin and generate enzyme/substrate intermediates in which the lichenin is covalently attached to the active site of the glucanosyltransferases. In the final step, the glucanosyltransferases attach the lichenin onto the processed galactomannans, which creates new glycosidic bonds and effectively incorporates the glycoproteins into the cross-linked cell wall glucan/chitin matrix.

Hiding in plain sight: New virus genomes discovered via a systematic analysis of fungal public transcriptomes

The distribution and diversity of RNA viruses in fungi is incompletely understood due to the often cryptic nature of mycoviral infections and the focused study of primarily pathogenic and/or economically important fungi. As most viruses that are known to infect fungi possess either single-stranded or double-stranded RNA genomes, transcriptomic data provides the opportunity to query for viruses in diverse fungal samples without any a priori knowledge of virus infection. Here we describe a systematic survey of all transcriptomic datasets from fungi belonging to the subphylum Pezizomycotina. Using a simple but effective computational pipeline that uses reads discarded during normal RNA-seq analyses, followed by identification of a viral RNA-dependent RNA polymerase (RdRP) motif in de novo assembled contigs, 59 viruses from 44 different fungi were identified. Among the viruses identified, 88% were determined to be new species and 68% are, to our knowledge, the first virus described from the fungal species. Comprehensive analyses of both nucleotide and inferred protein sequences characterize the phylogenetic relationships between these viruses and the known set of mycoviral sequences and support the classification of up to four new families and two new genera. Thus the results provide a deeper understanding of the scope of mycoviral diversity while also increasing the distribution of fungal hosts. Further, this study demonstrates the suitability of analyzing RNA-seq data to facilitate rapid discovery of new viruses.

The Zn(II)2Cys6-Type Transcription Factor ADA-6 Regulates Conidiation, Sexual Development, and Oxidative Stress Response in Neurospora crassa

Conidiation and sexual development are critical for reproduction, dispersal and better-adapted survival in many filamentous fungi. The Neurospora crassa gene ada-6 encodes a Zn(II)2Cys6-type transcription factor, whose deletion resulted in reduced conidial production and female sterility. In this study, we confirmed the positive contribution of ada-6 to conidiation and sexual development by detailed phenotypic characterization of its deletion mutant and the complemented mutant. To understand the regulatory mechanisms of ADA-6 in conidiation and sexual development, transcriptomic profiles generated by RNA-seq from the Δada-6 mutant and wild type during conidiation and sexual development were compared. During conidial development, differential expressed genes (DEGs) between the Δada-6 mutant and wild type are mainly involved in oxidation-reduction process and single-organism metabolic process. Several conidiation related genes are positively regulated by ADA-6, including genes that positively regulate conidiation (fluffy and acon-3), and genes preferentially expressed during conidial development (eas, con-6, con-8, con-10, con-13, pcp-1, and NCU9357), as the expression of these genes were lower in the Δada-6 mutant compared to wild type during conidial development. Phenotypic observation of deletion mutants for other genes with unknown function down-regulated by ada-6 deletion revealed that deletion mutants for four genes (NCU00929, NCU05260, NCU00116, and NCU04813) produced less conidia than wild type. Deletion of ada-6 resulted in female sterility, which might be due to that ADA-6 affects oxidation-reduction process and transmembrane transport process, and positively regulates the transcription of pre-2, poi-2, and NCU05832, three key genes participating in sexual development. In both conidiation and the sexual development process, ADA-6 regulates the transcription of cat-3 and other genes participating in reactive oxygen species production according to RNA-seq data, indicating a role of ADA-6 in oxidative stress response. This was further confirmed by the results that deletion of ada-6 led to hypersensitivity to oxidants H₂O₂ and menadione. Together, these results proved that ADA-6, as a global regulator, plays a crucial role in conidiation, sexual development, and oxidative stress response of N. crassa.

Altering Neurospora crassa MOB2A exposes its functions in development and affects its interaction with the NDR kinase COT1

The Neurospora crassa Mps One Binder (MOB) proteins MOB2A and MOB2B physically interact with the Nuclear Dbf2 Related (NDR) kinase COT1 and have been shown to have overlapping functions in various aspects of asexual development. Here, we identified two N. crassa MOB2A residues, Tyr117 and Tyr119, which are potentially phosphorylated. Using phosphomimetic mob-2a mutants we have been able to establish that apart from their previously described roles, MOB2A/B are involved in additional developmental processes. Enhanced conidial germination, accompanied by conidial agglutination, in the phosphomimetic mutants indicated that MOB2A is a negative regulator of germination. Thick-section imaging of perithecia revealed slow maturation and a lack of asci alignment in the mutant strains demonstrating a role for MOB2A in sexual development. We demonstrate that even though MOB2A and MOB2B have some overlapping functions, MOB2B cannot compensate for the roles MOB2A has in conidiation and germination. Altering Tyr residues 117 and 119 impaired the physical interactions between MOB2A and COT1, most likely contributing to some of the observed effects. As cot-1 and the phosphomimetic mutants share an extragenic suppressor (gul-1), we concluded that at least some of the effects imposed by altering Tyr117 and Tyr119 are mediated by the NDR kinase.

Light sensing by opsins and fungal ecology: NOP-1 modulates entry into sexual reproduction in response to environmental cues

Understanding the genetic basis of the switch from asexual to sexual lifestyles in response to sometimes rapid environmental changes is one of the major challenges in fungal ecology. Light appears to play a critical role in the asexual-sexual switch-but fungal genomes harbour diverse light sensors. Fungal opsins are homologous to bacterial green-light-sensory rhodopsins, and their organismal functions in fungi have not been well understood. Three of these opsin-like proteins were widely distributed across fungal genomes, but homologs of the Fusarium opsin-like protein CarO were present only in plant-associated fungi. Key amino acids, including potential retinal binding sites, functionally diverged on the phylogeny of opsins. This diversification of opsin-like proteins could be correlated with life history-associated differences among fungi in their expression and function during morphological development. In Neurospora crassa and related species, knockout of the opsin NOP-1 led to a phenotype in the regulation of the asexual-sexual switch, modulating response to both light and oxygen conditions. Sexual development commenced early in ∆nop-1 strains cultured in unsealed plates under constant blue and white light. Furthermore, comparative transcriptomics showed that the expression of nop-1 is light-dependent and that the ∆nop-1 strain abundantly expresses genes involved in oxidative stress response, genes enriched in NAD/NADP binding sites, genes with functions in proton transmembrane movement and catalase activity, and genes involved in the homeostasis of protons. Based on these observations, we contend that light and oxidative stress regulate the switch via light-responsive and ROS pathways in model fungus N. crassa and other fungi.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

Comparative Transcriptome Analyses in Zymoseptoria tritici Reveal Significant Differences in Gene Expression Among Strains During Plant Infection

Zymoseptoria tritici is an ascomycete fungus that causes Septoria tritici blotch, a globally distributed foliar disease on wheat. Z. tritici populations are highly polymorphic and exhibit significant quantitative variation for virulence. Despite its importance, the genes responsible for quantitative virulence in this pathogen remain largely unknown. We investigated the expression profiles of four Z. tritici strains differing in virulence in an experiment conducted under uniform environmental conditions. Transcriptomes were compared at four different infection stages to characterize the regulation of gene families thought to be involved in virulence and to identify new virulence factors. The major components of the fungal infection transcriptome showed consistent expression profiles across strains. However, strain-specific regulation was observed for many genes, including some encoding putative virulence factors. We postulate that strain-specific regulation of virulence factors can determine the outcome of Z. tritici infections. We show that differences in gene expression may be major determinants of virulence variation among Z. tritici strains, adding to the already known contributions to virulence variation based on differences in gene sequence and gene presence/absence polymorphisms.

The transcriptional landscape of basidiosporogenesis in mature Pisolithus microcarpus basidiocarp

Background

Pisolithus microcarpus (Cooke & Massee) G. Cunn is a gasteromycete that produces closed basidiocarps in symbiosis with eucalypts and acacias. The fungus produces a complex basidiocarp composed of peridioles at different developmental stages and an upper layer of basidiospores free of the hyphae and ready for wind dispersal upon the rupture of the basidiocarp pellis. During basidiosporogenesis, a process that takes place inside the basidiocarp peridioles, a conspicuous reserve of fatty acids is present throughout development. While several previous studies have described basidiosporogenesis inside peridioles, very little is known about gene expression changes that may occur during this part of the fungal life cycle. The objective of this work was to analyze gene transcription during peridiole and basidiospore development, while focusing specifically on cell cycle progression and lipid metabolism.

Results

Throughout different developmental stages of the peridioles we analyzed, 737 genes were regulated between adjacent compartments (>5 fold, FDR-corrected p-value < 0.05) corresponding to 3.49% of the genes present in the P. microcarpus genome. We identified three clusters among the regulated genes which showed differential expression between the peridiole developmental stages and the basidiospores. During peridiole development, transcripts for proteins involved in cellular processes, signaling, and information storage were detected, notably those for coding transcription factors, DNA polymerase subunits, DNA repair proteins, and genes involved in chromatin structure. For both internal embedded basidiospores (hereto referred to as “Internal spores”, IS) and external free basidiospores (hereto referred to as “Free spores”, FS), upregulated transcripts were found to involve primary metabolism, particularly fatty acid metabolism (FA). High expression of transcripts related to β-oxidation and the glyoxylate shunt indicated that fatty acids served as a major carbon source for basidiosporogenesis.

Conclusion

Our results show that basidiocarp formation in P. microcarpus involves a complex array of genes that are regulated throughout peridiole development. We identified waves of transcripts with coordinated regulation and identified transcription factors which may play a role in this regulation. This is the first work to describe gene expression patterns during basidiocarp formation in an ectomycorrhizal gasteromycete fungus and sheds light on genes that may play important roles in the developmental process.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3545-5) contains supplementary material, which is available to authorized users.

A-to-I RNA editing independent of ADARs in filamentous fungi

ADAR mediated A-to-I RNA editing is thought to be unique to animals and occurs mainly in the non-coding regions. Recently filamentous fungi such as Fusarium graminearum were found to lack orthologs of animal ADARs but have stage-specific A-to-I editing during sexual reproduction. Unlike animals, majority of editing sites are in the coding regions and often result in missense and stop loss changes in fungi. Furthermore, whereas As in RNA stems are targeted by animal ADARs, RNA editing in fungi preferentially targets As in hairpin loops, implying that fungal RNA editing involves mechanisms related to editing of the anticodon loop by ADATs. Identification and characterization of fungal adenosine deaminases and their stage-specific co-factors may be helpful to understand the evolution of human ADARs. Fungi also can be used to study biological functions of missense and stop loss RNA editing events in eukaryotic organisms.

The Fast-Evolving phy-2 Gene Modulates Sexual Development in Response to Light in the Model Fungus Neurospora crassa

Rapid responses to changes in incident light are critical to the guidance of behavior and development in most species. Phytochrome light receptors in particular play key roles in bacterial physiology and plant development, but their functions and regulation are less well understood in fungi. Nevertheless, genome-wide expression measurements provide key information that can guide experiments that reveal how genes respond to environmental signals and clarify their role in development. We performed functional genomic and phenotypic analyses of the two phytochromes in Neurospora crassa, a fungal model adapted to a postfire environment that experiences dramatically variable light conditions. Expression of phy-1 and phy-2 was low in early sexual development and in the case of phy-2 increased in late sexual development. Under light stimulation, strains with the phytochromes deleted exhibited increased expression of sexual development-related genes. Moreover, under red light, the phy-2 knockout strain commenced sexual development early. In the evolution of phytochromes within ascomycetes, at least two duplications have occurred, and the faster-evolving phy-2 gene has frequently been lost. Additionally, the three key cysteine sites that are critical for bacterial and plant phytochrome function are not conserved within fungal phy-2 homologs. Through the action of phytochromes, transitions between asexual and sexual reproduction are modulated by light level and light quality, presumably as an adaptation for fast asexual growth and initiation of sexual reproduction of N. crassa in exposed postfire ecosystems.

Environmental signals, including light, play critical roles in regulating fungal growth and pathogenicity, and balance of asexual and sexual reproduction is critical in fungal pathogens’ incidence, virulence, and distribution. Red light sensing by phytochromes is well known to play critical roles in bacterial physiology and plant development. Homologs of phytochromes were first discovered in the fungal model Neurospora crassa and then subsequently in diverse other fungi, including many plant pathogens. Our study investigated the evolution of red light sensors in ascomycetes and confirmed—using the model fungus Neurospora crassa—their roles in modulating the asexual-sexual reproduction balance in fungi. Our findings also provide a key insight into one of the most poorly understood aspects of fungal biology, suggesting that further study of the function of phytochromes in fungi is critical to reveal the genetic basis of the asexual-sexual switch responsible for fungal growth and distribution, including diverse and destructive plant pathogens.

A Tool Set for the Genome-Wide Analysis of Neurospora crassa by RT-PCR

Neurospora crassa is an important model organism for filamentous fungi as well as for circadian biology and photobiology. Although the community-accumulated tool set for the molecular analysis of Neurospora is extensive, two components are missing: (1) dependable reference genes whose level of expression are relatively constant across light/dark cycles and as a function of time of day and (2) a catalog of primers specifically designed for real-time PCR (RT-PCR). To address the first of these we have identified genes that are optimal for use as reference genes in RT-PCR across a wide range of expression levels; the mRNA/transcripts from these genes have potential for use as reference noncycling transcripts outside of Neurospora. In addition, we have generated a genome-wide set of RT-PCR primers, thereby streamlining the analysis of gene expression. In validation studies these primers successfully identified target mRNAs arising from 70% (34 of 49) of all tested genes and from all (28) of the moderately to highly expressed tested genes.

Functional Analysis of Developmentally Regulated Genes chs7 and sec22 in the Ascomycete Sordaria macrospora

During sexual development, filamentous ascomycetes form complex, three-dimensional fruiting bodies for the generation and dispersal of spores. In previous studies, we identified genes with evolutionary conserved expression patterns during fruiting body formation in several fungal species. Here, we present the functional analysis of two developmentally up-regulated genes, chs7 and sec22, in the ascomycete Sordaria macrospora. The genes encode a class VII (division III) chitin synthase and a soluble N-ethylmaleimide-sensitive-factor attachment protein receptor (SNARE) protein, respectively. Deletion mutants of chs7 had normal vegetative growth and were fully fertile but showed sensitivity toward cell wall stress. Deletion of sec22 resulted in a reduced number of ascospores and in defects in ascospore pigmentation and germination, whereas vegetative growth was normal in the mutant. A SEC22-EGFP fusion construct under control of the native sec22 promoter and terminator regions was expressed during different stages of sexual development. Expression of several development-related genes was deregulated in the sec22 mutant, including three genes involved in melanin biosynthesis. Our data indicate that chs7 is dispensable for fruiting body formation in S. macrospora, whereas sec22 is required for ascospore maturation and germination and thus involved in late stages of sexual development.