Global gene expression and focused knockout analysis reveals genes associated with fungal fruiting body development in Neurospora crassa

Regulatory functions of homeobox domain transcription factors in fungi

Homeobox domain (HD) proteins present a crucial involvement in morphological differentiation and other functions in eukaryotes. Most HD genes encode transcription factors (TFs) that orchestrate a regulatory role in cellular and developmental decisions. In fungi, multiple studies have increased our understanding of these important HD regulators in recent years. These reports have revealed their role in fungal development, both sexual and asexual, as well as their importance in governing other biological processes in these organisms, including secondary metabolism, pathogenicity, and sensitivity to environmental stresses. Here, we provide a comprehensive review of the current knowledge on the regulatory roles of HD-TFs in fungi, with a special focus on Aspergillus species.

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.

A DEAD-box RNA helicase mediates meiotic silencing by unpaired DNA

During the sexual phase of Neurospora crassa, unpaired genes are subject to a silencing mechanism known as meiotic silencing by unpaired DNA (MSUD). MSUD targets the transcripts of an unpaired gene and utilizes typical RNA interference factors for its process. Using a reverse genetic screen, we have identified a meiotic silencing gene called sad-9, which encodes a DEAD-box RNA helicase. While not essential for vegetative growth, SAD-9 plays a crucial role in both sexual development and MSUD. Our results suggest that SAD-9, with the help of the SAD-2 scaffold protein, recruits the SMS-2 Argonaute to the perinuclear region, the center of MSUD activity.

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.

Comparative Transcriptomics of Fusarium graminearum and Magnaporthe oryzae Spore Germination Leading up To Infection

For fungal plant pathogens, the germinating spore provides the first interaction with the host. Spore germlings move across the plant surface and use diverse penetration strategies for ingress into plant surfaces. Penetration strategies include pressurized melanized appressoria, which facilitate physically punching through the plant cuticle, and nonmelanized appressoria, which penetrate with the help of enzymes or cuticular damage to breach the plant surface. Two well-studied plant pathogens, Fusarium graminearum and Magnaporthe oryzae, are typical of these two modes of penetration. We applied comparative transcriptomics to Fusarium graminearum and Magnaporthe oryzae to characterize the genetic programming of the early host-pathogen interface. Four sequential stages of development following spore localization on the plant surface, from spore swelling to appressorium formation, were sampled for each species on culture medium and on barley sheaths, and transcriptomic analyses were performed. Gene expression in the prepenetration stages in both species and under both conditions was similar. In contrast, gene expression in the final stage was strongly influenced by the environment. Appressorium formation involved the greatest number of differentially expressed genes. Laser-dissection microscopy was used to perform detailed transcriptomics of initial infection points by F. graminearum. These analyses revealed new and important aspects of early fungal ingress in this species. Expression of the trichothecene genes involved in biosynthesis of deoxynivalenol by F. graminearum implies that toxisomes are not fully functional until after penetration and indicates that deoxynivalenol is not essential for penetration under our conditions. The use of comparative gene expression of divergent fungi promises to advance highly effective targets for antifungal strategies. IMPORTANCE Fusarium graminearum and Magnaporthe oryzae are two of the most important pathogens of cereal grains worldwide. Despite years of research, strong host resistance has not been identified for F. graminearum, so other methods of control are essential. The pathogen takes advantage of multiple entry points to infect the host, including breaches in the florets due to senescence of flower parts and penetration of the weakened trichome bases to breach the epidermis. In contrast, M. oryzae directly punctures leaves that it infects, and resistant cultivars have been characterized. The threat of either pathogen causing a major disease outbreak is ever present. Comparative transcriptomics demonstrated its potential to reveal novel and effective disease prevention strategies that affect the initial stages of disease. Shedding light on the basis of this diversity of infection strategies will result in development of increasingly specific control strategies.

A moonlighting function of a chitin polysaccharide monooxygenase, CWR-1, in Neurospora crassa allorecognition

Organisms require the ability to differentiate themselves from organisms of different or even the same species. Allorecognition processes in filamentous fungi are essential to ensure identity of an interconnected syncytial colony to protect it from exploitation and disease. Neurospora crassa has three cell fusion checkpoints controlling formation of an interconnected mycelial network. The locus that controls the second checkpoint, which allows for cell wall dissolution and subsequent fusion between cells/hyphae, cwr (cell wall remodeling), encodes two linked genes, cwr-1 and cwr-2. Previously, it was shown that cwr-1 and cwr-2 show severe linkage disequilibrium with six different haplogroups present in N. crassa populations. Isolates from an identical cwr haplogroup show robust fusion, while somatic cell fusion between isolates of different haplogroups is significantly blocked in cell wall dissolution. The cwr-1 gene encodes a putative polysaccharide monooxygenase (PMO). Herein we confirm that CWR-1 is a C1-oxidizing chitin PMO. We show that the catalytic (PMO) domain of CWR-1 was sufficient for checkpoint function and cell fusion blockage; however, through analysis of active-site, histidine-brace mutants, the catalytic activity of CWR-1 was ruled out as a major factor for allorecognition. Swapping a portion of the PMO domain (V86 to T130) did not switch cwr haplogroup specificity, but rather cells containing this chimera exhibited a novel haplogroup specificity. Allorecognition to mediate cell fusion blockage is likely occurring through a protein-protein interaction between CWR-1 with CWR-2. These data highlight a moonlighting role in allorecognition of the CWR-1 PMO domain.

Differential Expression of Cell Wall Remodeling Genes Is Part of the Dynamic Phase-Specific Transcriptional Program of Conidial Germination of Trichoderma asperelloides

The nature of saprophytic and mycoparasitic hyphal growth of Trichoderma spp. has been studied extensively, yet its initiation via conidial germination in this genus is less well understood. Using near-synchronous germinating cultures of Trichoderma asperelloides, we followed the morphological progression from dormant conidia to initial polar growth to germling formation and to evidence for first branching. We found that the stage-specific transcriptional profile of T. asperelloides is one of the most dynamic described to date: transcript abundance of over 5000 genes-comprising approximately half of the annotated genome-was unremittingly reduced in the transition from dormancy to polar growth. Conversely, after the onset of germination, the transcript abundance of approximately a quarter of the genome was unremittingly elevated during the transition from elongation to initial branching. These changes are a testimony to the substantial developmental events that accompany germination. Bayesian network analysis identified several chitinase- and glucanase-encoding genes as active transcriptional hubs during germination. Furthermore, the expression of specific members of the chitin synthase and glucan elongase families was significantly increased during germination in the presence of Rhizoctonia solani-a known host of the mycoparasite-indicating that host recognition can occur during the early stages of mycoparasite development.

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.

Transcriptional Divergence Underpinning Sexual Development in the Fungal Class Sordariomycetes

Gene expression divergence through evolutionary processes is thought to be important for achieving programmed development in multicellular organisms. To test this premise in filamentous fungi, we investigated transcriptional profiles of 3,942 single-copy orthologous genes (SCOGs) in five related sordariomycete species that have morphologically diverged in the formation of their flask-shaped perithecia. We compared expression of the SCOGs to inferred gene expression levels of the most recent common ancestor of the five species, ranking genes from their largest increases to smallest increases in expression during perithecial development in each of the five species. We found that a large proportion of the genes that exhibited evolved increases in gene expression were important for normal perithecial development in Fusarium graminearum. Many of these genes were previously uncharacterized, encoding hypothetical proteins without any known functional protein domains. Interestingly, the developmental stages during which aberrant knockout phenotypes appeared largely coincided with the elevated expression of the deleted genes. In addition, we identified novel genes that affected normal perithecial development in Magnaporthe oryzae and Neurospora crassa, which were functionally and transcriptionally diverged from the orthologous counterparts in F. graminearum. Furthermore, comparative analysis of developmental transcriptomes and phylostratigraphic analysis suggested that genes encoding hypothetical proteins are generally young and transcriptionally divergent between related species. This study provides tangible evidence of shifts in gene expression that led to acquisition of novel function of orthologous genes in each lineage and demonstrates that several genes with hypothetical function are crucial for shaping multicellular fruiting bodies. IMPORTANCE The fungal class Sordariomycetes includes numerous important plant and animal pathogens. It also provides model systems for studying fungal fruiting body development, as its members develop fruiting bodies with a few well-characterized tissue types on common growth media and have rich genomic resources that enable comparative and functional analyses. To understand transcriptional divergence of key developmental genes between five related sordariomycete fungi, we performed targeted knockouts of genes inferred to have evolved significant upward shifts in expression. We found that many previously uncharacterized genes play indispensable roles at different stages of fruiting body development, which have undergone transcriptional activation in specific lineages. These novel genes are predicted to be phylogenetically young and tend to be involved in lineage- or species-specific function. Transcriptional activation of genes with unknown function seems to be more frequent than ever thought, which may be crucial for rapid adaption to changing environments for successful sexual reproduction.

Genetic Networks That Govern Sexual Reproduction in the Pezizomycotina

Sexual development in filamentous fungi is a complex process that relies on the precise control of and interaction between a variety of genetic networks and pathways. The mating-type (MAT) genes are the master regulators of this process and typically act as transcription factors, which control the expression of genes involved at all stages of the sexual cycle. In many fungi, the sexual cycle typically begins when the mating pheromones of one mating type are recognized by a compatible partner, followed by physical interaction and fertilization. Subsequently, highly specialized sexual structures are formed, within which the sexual spores develop after rounds of meiosis and mitosis. These spores are then released and germinate, forming new individuals that initiate new cycles of growth. This review provides an overview of the known genetic networks and pathways that are involved in each major stage of the sexual cycle in filamentous ascomycete fungi.

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.

Kinetic Characterization of a Putatively Chitin-Active LPMO Reveals a Preference for Soluble Substrates and Absence of Monooxygenase Activity

Enzymes known as lytic polysaccharide monooxygenases (LPMOs) are recognized as important contributors to aerobic enzymatic degradation of recalcitrant polysaccharides such as chitin and cellulose. LPMOs are remarkably abundant in nature, with some fungal species possessing more than 50 LPMO genes, and the biological implications of this diversity remain enigmatic. For example, chitin-active LPMOs have been encountered in biological niches where chitin conversion does not seem to take place. We have carried out an in-depth kinetic characterization of a putatively chitin-active LPMO from Aspergillus fumigatus (AfAA11B), which, as we show here, has multiple unusual properties, such as a low redox potential and high oxidase activity. Furthermore, AfAA11B is hardly active on chitin, while being very active on soluble oligomers of N-acetylglucosamine. In the presence of chitotetraose, the enzyme can withstand considerable amounts of H₂O₂, which it uses to efficiently and stoichiometrically convert this substrate. The unique properties of AfAA11B allowed experiments showing that it is a strict peroxygenase and does not catalyze a monooxygenase reaction. This study shows that nature uses LPMOs for breaking glycosidic bonds in non-polymeric substrates in reactions that depend on H₂O₂. The quest for the true substrates of these enzymes, possibly carbohydrates in the cell wall of the fungus or its competitors, will be of major interest.

Recombination-independent recognition of DNA homology for meiotic silencing in Neurospora crassa

The pairing of homologous chromosomes represents a critical step of meiosis in nearly all sexually reproducing species. In many organisms, pairing involves chromosomes that remain apparently intact. The mechanistic nature of homology recognition at the basis of such pairing is unknown. Using "meiotic silencing by unpaired DNA" (MSUD) as a model process, we demonstrate the existence of a cardinally different approach to DNA homology recognition in meiosis. The main advantage of MSUD over other experimental systems lies in its ability to identify any relatively short DNA fragment lacking a homologous allelic partner. Here, we show that MSUD does not rely on the canonical mechanism of meiotic recombination, yet it is promoted by REC8, a conserved component of the meiotic cohesion complex. We also show that certain patterns of interspersed homology are recognized as pairable during MSUD. Such patterns need to be colinear and must contain short tracts of sequence identity spaced apart at 21 or 22 base pairs. By using these periodicity values as a guiding parameter in all-atom molecular modeling, we discover that homologous DNA molecules can pair by forming quadruplex-based contacts with an interval of 2.5 helical turns. This process requires right-handed plectonemic coiling and additional conformational changes in the intervening double-helical segments. Our results 1) reconcile genetic and biophysical evidence for the existence of direct homologous double-stranded DNA (dsDNA)-dsDNA pairing, 2) identify a role for this process in initiating RNA interference, and 3) suggest that chromosomes can be cross-matched by a precise mechanism that operates on intact dsDNA molecules.

A global search for novel transcription factors impacting the Neurospora crassa circadian clock

Eukaryotic circadian oscillators share a common circuit architecture, a negative feedback loop in which a positive element activates the transcription of a negative one that then represses the action of the former, inhibiting its own expression. While studies in mammals and insects have revealed additional transcriptional inputs modulating the expression of core clock components, this has been less characterized in the model Neurospora crassa, where the participation of other transcriptional components impacting circadian clock dynamics remains rather unexplored. Thus, we sought to identify additional transcriptional regulators modulating the N. crassa clock, following a reverse genetic screen based on luminescent circadian reporters and a collection of transcription factors (TFs) knockouts, successfully covering close to 60% of them. Besides the canonical core clock components WC-1 and -2, none of the tested transcriptional regulators proved to be essential for rhythmicity. Nevertheless, we identified a set of 23 TFs that when absent lead to discrete, but significant, changes in circadian period. While the current level of analysis does not provide mechanistic information about how these new players modulate circadian parameters, the results of this screen reveal that an important number of light and clock-regulated TFs, involved in a plethora of processes, are capable of modulating the clockworks. This partial reverse genetic clock screen also exemplifies how the N. crassa knockout collection continues to serve as an expedite platform to address broad biological questions.

An introgressed gene causes meiotic drive in Neurospora sitophila

In order to survive, most organisms must deal with parasites. Such parasites can be other organisms or, sometimes, selfish genes found within the host genome itself. While much is known about parasitic organisms, the interaction with their hosts, and their ability to spread within and between species, much less is known about selfish genes. We here identify a selfish “spore killer” gene in the fungus Neurospora sitophila. The gene appears to have evolved within the genus but has entered the species through hybridization and introgression. We also show that the host can counteract the gene through RNA interference. These results shed light on the diversity of selfish genes in terms of origin, evolution, and host interactions.

Meiotic drive elements cause their own preferential transmission following meiosis. In fungi, this phenomenon takes the shape of spore killing, and in the filamentous ascomycete Neurospora sitophila, the Sk-1 spore killer element is found in many natural populations. In this study, we identify the gene responsible for spore killing in Sk-1 by generating both long- and short-read genomic data and by using these data to perform a genome-wide association test. We name this gene Spk-1. Through molecular dissection, we show that a single 405-nt-long open reading frame generates a product that both acts as a poison capable of killing sibling spores and as an antidote that rescues spores that produce it. By phylogenetic analysis, we demonstrate that the gene has likely been introgressed from the closely related species Neurospora hispaniola, and we identify three subclades of N. sitophila, one where Sk-1 is fixed, another where Sk-1 is absent, and a third where both killer and sensitive strain are found. Finally, we show that spore killing can be suppressed through an RNA interference-based genome defense pathway known as meiotic silencing by unpaired DNA. Spk-1 is not related to other known meiotic drive genes, and similar sequences are only found within Neurospora. These results shed light on the diversity of genes capable of causing meiotic drive, their origin and evolution, and their interaction with the host genome.

Shannon entropy as a metric for conditional gene expression in Neurospora crassa

Neurospora crassa has been an important model organism for molecular biology and genetics for over 60 years. Neurospora crassa has a complex life cycle, with over 28 distinct cell types and is capable of transcriptional responses to many environmental conditions including nutrient availability, temperature, and light. To quantify variation in N. crassa gene expression, we analyzed public expression data from 97 conditions and calculated the Shannon Entropy value for Neurospora’s approximately 11,000 genes. Entropy values can be used to estimate the variability in expression for a single gene over a range of conditions and be used to classify individual genes as constitutive or condition-specific. Shannon entropy has previously been used measure the degree of tissue specificity of multicellular plant or animal genes. We use this metric here to measure variable gene expression in a microbe and provide this information as a resource for the N. crassa research community. Finally, we demonstrate the utility of this approach by using entropy values to identify genes with constitutive expression across a wide range of conditions and to identify genes that are activated exclusively during sexual development.

Unmatched Level of Molecular Convergence among Deeply Divergent Complex Multicellular Fungi

Convergent evolution is pervasive in nature, but it is poorly understood how various constraints and natural selection limit the diversity of evolvable phenotypes. Here, we analyze the transcriptome across fruiting body development to understand the independent evolution of complex multicellularity in the two largest clades of fungi—the Agarico- and Pezizomycotina. Despite >650 My of divergence between these clades, we find that very similar sets of genes have convergently been co-opted for complex multicellularity, followed by expansions of their gene families by duplications. Over 82% of shared multicellularity-related gene families were expanding in both clades, indicating a high prevalence of convergence also at the gene family level. This convergence is coupled with a rich inferred repertoire of multicellularity-related genes in the most recent common ancestor of the Agarico- and Pezizomycotina, consistent with the hypothesis that the coding capacity of ancestral fungal genomes might have promoted the repeated evolution of complex multicellularity. We interpret this repertoire as an indication of evolutionary predisposition of fungal ancestors for evolving complex multicellular fruiting bodies. Our work suggests that evolutionary convergence may happen not only when organisms are closely related or are under similar selection pressures, but also when ancestral genomic repertoires render certain evolutionary trajectories more likely than others, even across large phylogenetic distances.

Clustering analysis of large-scale phenotypic data in the model filamentous fungus Neurospora crassa

Background

With 9730 protein-coding genes and a nearly complete gene knockout strain collection, Neurospora crassa is a major model organism for filamentous fungi. Despite this abundance of information, the phenotypes of these gene knockout mutants have not been categorized to determine whether there are broad correlations between phenotype and any genetic features.

Results

Here, we analyze data for 10 different growth or developmental phenotypes that have been obtained for 1168 N. crassa knockout mutants. Of these mutants, 265 (23%) are in the normal range, while 903 (77%) possess at least one mutant phenotype. With the exception of unclassified functions, the distribution of functional categories for genes in the mutant dataset mirrors that of the N. crassa genome. In contrast, most genes do not possess a yeast ortholog, suggesting that our analysis will reveal functions that are not conserved in Saccharomyces cerevisiae. To leverage the phenotypic data to identify pathways, we used weighted Partitioning Around Medoids (PAM) approach with 40 clusters. We found that genes encoding metabolic, transmembrane and protein phosphorylation-related genes are concentrated in subsets of clusters. Results from K-Means clustering of transcriptomic datasets showed that most phenotypic clusters contain multiple expression profiles, suggesting that co-expression is not generally observed for genes with shared phenotypes. Analysis of yeast orthologs of genes that co-clustered in MAPK signaling cascades revealed potential networks of interacting proteins in N. crassa.

Conclusions

Our results demonstrate that clustering analysis of phenotypes is a promising tool for generating new hypotheses regarding involvement of genes in cellular pathways in N. crassa. Furthermore, information about gene clusters identified in N. crassa should be applicable to other filamentous fungi, including saprobes and pathogens.

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.

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.

Abundant Perithecial Protein (APP) from Neurospora is a primitive functional analog of ocular crystallins

The eye arose during the Cambrian explosion from pre-existing proteins that would have been recruited for the formation of the specialized components of this organ, such as the transparent lens. Proteins suitable for the role of lens crystallins would need to possess unusual physical properties and the study of such earliest analogs of ocular crystallins would add to our understanding of the nature of recruitment of proteins as lens/corneal crystallins. We show that the Abundant Perithecial Protein (APP) of the fungi Neurospora and Sordaria fulfils the criteria for an early crystallin analog. The perithecia in these fungal species are phototropic, and APP accumulates at a high concentration in the neck of the pitcher-shaped perithecium. Spores are formed at the base of the perithecium, and light contributes to their maturation. The hydrodynamic properties of APP appear to exclude dimer formation or aggregation at high protein concentrations. APP is also deficient in Ca²⁺ binding, a property seen in its close homolog, the calcium-binding cell adhesion molecule (DdCAD-1) from Dictyostelium discoidum. Comparable to crystallins, APP occurs in high concentrations and seems to have dispensed with Ca²⁺ binding in exchange for greater stability. These crystallin-like attributes of APP lead us to demonstrate that it is a primitive form of ocular crystallins.

The Blue-Light Photoreceptor Sfwc-1 Gene Regulates the Phototropic Response and Fruiting-Body Development in the Homothallic Ascomycete Sordaria fimicola

Sordaria fimicola, a coprophilous ascomycete, is a homothallic fungus that can undergo sexual differentiation with cellular and morphological changes followed by multicellular tissue development to complete its sexual cycle. In this study, we identified and characterized the blue-light photoreceptor gene in S. fimicola The S. fimicola white collar-1 photoreceptor (SfWC-1) contains light-oxygen-voltage-sensing (LOV), Per-Arnt-Sim (PAS), and other conserved domains and is homologous to the WC-1 blue-light photoreceptor of Neurospora crassa The LOV domain of Sfwc-1 was deleted by homologous recombination using Agrobacterium-mediated protoplast transformation. The Sfwc-1(Δlov) mutant showed normal vegetative growth but produced less carotenoid pigment under illumination. The mutant showed delayed and less-pronounced fruiting-body formation, was defective in phototropism of the perithecial beaks, and lacked the fruiting-body zonation pattern compared with the wild type under the illumination condition. Gene expression analyses supported the light-induced functions of the Sfwc-1 gene in the physiology and developmental process of perithecial formation in S. fimicola Moreover, green fluorescent protein (GFP)-tagged SfWC-1 fluorescence signals were transiently strong upon light induction and prominently located inside the nuclei of living hyphae. Our studies focused on the putative blue-light photoreceptor in a model ascomycete and contribute to a better understanding of the photoregulatory functions and networks mediated by the evolutionarily conserved blue-light photoreceptors across diverse fungal phyla.IMPORTANCESordaria sp. has been a model for study of fruiting-body differentiation in fungi. Several environmental factors, including light, affect cellular and morphological changes during multicellular tissue development. Here, we created a light-oxygen-voltage-sensing (LOV) domain-deleted Sfwc-1 mutant to study blue-light photoresponses in Sordaria fimicola Phototropism and rhythmic zonation of perithecia were defective in the Sfwc-1(Δlov) mutant. Moreover, fruiting-body development in the mutant was reduced and also significantly delayed. Gene expression analysis and subcellular localization study further revealed the light-induced differential gene expression and cellular responses upon light stimulation in S. fimicola.

It's All in the Genes: The Regulatory Pathways of Sexual Reproduction in Filamentous Ascomycetes

Sexual reproduction in filamentous ascomycete fungi results in the production of highly specialized sexual tissues, which arise from relatively simple, vegetative mycelia. This conversion takes place after the recognition of and response to a variety of exogenous and endogenous cues, and relies on very strictly regulated gene, protein, and metabolite pathways. This makes studying sexual development in fungi an interesting tool in which to study gene-gene, gene-protein, and protein-metabolite interactions. This review provides an overview of some of the most important genes involved in this process; from those involved in the conversion of mycelia into sexually-competent tissue, to those involved in the development of the ascomata, the asci, and ultimately, the ascospores.

Metabolism and Development during Conidial Germination in Response to a Carbon-Nitrogen-Rich Synthetic or a Natural Source of Nutrition in Neurospora crassa

Fungal spores germinate and undergo vegetative growth, leading to either asexual or sexual reproductive dispersal. Previous research has indicated that among developmental regulatory genes, expression is conserved across nutritional environments, whereas pathways for carbon and nitrogen metabolism appear highly responsive-perhaps to accommodate differential nutritive processing. To comprehensively investigate conidial germination and the adaptive life history decision-making underlying these two modes of reproduction, we profiled transcription of Neurospora crassa germinating on two media: synthetic Bird medium, designed to promote asexual reproduction; and a natural maple sap medium, on which both asexual reproduction and sexual reproduction manifest. A later start to germination but faster development was observed on synthetic medium. Metabolic genes exhibited altered expression in response to nutrients-at least 34% of the genes in the genome were significantly downregulated during the first two stages of conidial germination on synthetic medium. Knockouts of genes exhibiting differential expression across development altered germination and growth rates, as well as in one case causing abnormal germination. A consensus Bayesian network of these genes indicated especially tight integration of environmental sensing, asexual and sexual development, and nitrogen metabolism on a natural medium, suggesting that in natural environments, a more dynamic and tentative balance of asexual and sexual development may be typical of N. crassa colonies.IMPORTANCE One of the most remarkable successes of life is its ability to flourish in response to temporally and spatially varying environments. Fungi occupy diverse ecosystems, and their sensitivity to these environmental changes often drives major fungal life history decisions, including the major switch from vegetative growth to asexual or sexual reproduction. Spore germination comprises the first and simplest stage of vegetative growth. We examined the dependence of this early life history on the nutritional environment using genome-wide transcriptomics. We demonstrated that for developmental regulatory genes, expression was generally conserved across nutritional environments, whereas metabolic gene expression was highly labile. The level of activation of developmental genes did depend on current nutrient conditions, as did the modularity of metabolic and developmental response network interactions. This knowledge is critical to the development of future technologies that could manipulate fungal growth for medical, agricultural, or industrial purposes.

Diverse members of the Xylariales lack canonical mating-type regions

A survey of genomes reported here for 10 isolates of Monosporascus species and an additional 25 genomes from other members of the Xylariales (representing 15 genera) available in public databases indicated that genes typically associated with MAT1-1 (mat A) or MAT1-2 (mat a) mating types are absent or have diverged greatly relative to counterparts in other Pezizomycotina. This was particularly surprising for isolates known to be homothallic, given that homothallic members of the Pezizomycotina typically possess a MAT1-1-1 (mat A-1) gene and one or both of two other closely-linked mating-type genes, MAT1-1-2 (mat A-2) and MAT1-1-3 (mat A-3), in addition to MAT1-2-1 (mat a-1). We failed to detect candidate genes for either MAT1-1-1 or MAT1-1-2 in any member of the Xylariales. Genes related to MAT1-2-1 and MAT1-1-3 are present in the genomes examined, but most appear to be orthologs of MATA_HMG (high-mobility group) genes with non-mating-type functions rather than orthologs of mating-type genes. Several MATA_HMG genes were found in genome positions that suggest they are derived from mating-type genes, but these genes are highly divergent relative to known MAT1-2-1 and MAT1-1-3 genes. The genomes examined represent substantial diversity within the order and include M. cannonballus, M. ibericus, Xylaria hypoxylon, X. striata, Daldinia eschscholzii, Eutypa lata, Rosellinia necatrix, Microdochium bolleyi and several others. We employed a number of avenues to search for homologs, including multiple BLAST approaches and examination of annotated genes adjacent to genes known to flank mating regions in other members of the Ascomycota. The results suggest that the mating regions have been lost from, or altered dramatically in, the Xylariales genomes examined and that mating and sexual development in these fungi are controlled differently than has been reported for members of the Pezizomycotina studied to date.

Mating-type factor-specific regulation of the fumagillin/pseurotin secondary metabolite supercluster in Aspergillus fumigatus

In the human pathogenic mold Aspergillus fumigatus, sexual identity is determined by the mating-type idiomorphs MAT1-1 and MAT1-2 residing at the MAT locus. Upon crossing of compatible partners, a heterothallic mating is executed to eventually form cleistothecia that contain recombinant ascospores. Given that the MAT1 gene products are DNA binding master regulators that govern this complex developmental process, we monitored the MAT1-driven transcriptomes of A. fumigatus by conditional overexpression of either MAT1 gene followed by RNA-seq analyses. Numerous genes related to the process of mating were found to be under transcriptional control, such as pheromone production and recognition. Substantial differences between the MAT1-1- and MAT1-2-driven transcriptomes could be detected by functional categorization of differentially expressed genes. Moreover, a significant and distinct impact on expression of genetic clusters of secondary metabolism became apparent, which could be verified on the product level. Unexpectedly, specific cross-regulation of the fumagillin/pseurotin supercluster was evident, thereby uncoupling its co-regulatory characteristic. These insights imply a tight interconnection of sexual development accompanied by ascosporogenesis with secondary metabolite production of a pathogenic fungus and impose evolutionary constraints that link these two fundamental aspects of the fungal lifestyle.

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.

Functional Profiling of Transcription Factor Genes in Neurospora crassa

Regulation of gene expression by DNA-binding transcription factors is essential for proper control of growth and development in all organisms. In this study, we annotate and characterize growth and developmental phenotypes for transcription factor genes in the model filamentous fungus Neurospora crassa. We identified 312 transcription factor genes, corresponding to 3.2% of the protein coding genes in the genome. The largest class was the fungal-specific Zn₂Cys₆ (C6) binuclear cluster, with 135 members, followed by the highly conserved C2H2 zinc finger group, with 61 genes. Viable knockout mutants were produced for 273 genes, and complete growth and developmental phenotypic data are available for 242 strains, with 64% possessing at least one defect. The most prominent defect observed was in growth of basal hyphae (43% of mutants analyzed), followed by asexual sporulation (38%), and the various stages of sexual development (19%). Two growth or developmental defects were observed for 21% of the mutants, while 8% were defective in all three major phenotypes tested. Analysis of available mRNA expression data for a time course of sexual development revealed mutants with sexual phenotypes that correlate with transcription factor transcript abundance in wild type. Inspection of this data also implicated cryptic roles in sexual development for several cotranscribed transcription factor genes that do not produce a phenotype when mutated.

A-to-I RNA editing is developmentally regulated and generally adaptive for sexual reproduction in Neurospora crassa

Although fungi lack adenosine deaminase acting on RNA (ADAR) enzymes, adenosine to inosine (A-to-I) RNA editing was reported recently in Fusarium graminearum during sexual reproduction. In this study, we profiled the A-to-I editing landscape and characterized its functional and adaptive properties in the model filamentous fungus Neurospora crassa A total of 40,677 A-to-I editing sites were identified, and approximately half of them displayed stage-specific editing or editing levels at different sexual stages. RNA-sequencing analysis with the Δstc-1 and Δsad-1 mutants confirmed A-to-I editing occurred before ascus development but became more prevalent during ascosporogenesis. Besides fungal-specific sequence and secondary structure preference, 63.5% of A-to-I editing sites were in the coding regions and 81.3% of them resulted in nonsynonymous recoding, resulting in a significant increase in the proteome complexity. Many genes involved in RNA silencing, DNA methylation, and histone modifications had extensive recoding, including sad-1, sms-3, qde-1, and dim-2. Fifty pseudogenes harbor premature stop codons that require A-to-I editing to encode full-length proteins. Unlike in humans, nonsynonymous editing events in N. crassa are generally beneficial and favored by positive selection. Almost half of the nonsynonymous editing sites in N. crassa are conserved and edited in Neurospora tetrasperma Furthermore, hundreds of them are conserved in F. graminearum and had higher editing levels. Two unknown genes with editing sites conserved between Neurospora and Fusarium were experimentally shown to be important for ascosporogenesis. This study comprehensively analyzed A-to-I editing in N. crassa and showed that RNA editing is stage-specific and generally adaptive, and may be functionally related to repeat induced point mutation and meiotic silencing by unpaired DNA.

An RNA Recognition Motif-Containing Protein Functions in Meiotic Silencing by Unpaired DNA

Meiotic silencing by unpaired DNA (MSUD) is a biological process that searches pairs of homologous chromosomes (homologs) for segments of DNA that are unpaired. Genes found within unpaired segments are silenced for the duration of meiosis. In this report, we describe the identification and characterization of Neurospora crassa sad-7, a gene that encodes a protein with an RNA recognition motif (RRM). Orthologs of sad-7 are found in a wide range of ascomycete fungi. In N. crassa, sad-7 is required for a fully efficient MSUD response to unpaired genes. Additionally, at least one parent must have a functional sad-7 allele for a cross to produce ascospores. Although sad-7-null crosses are barren, sad-7^Δ strains grow at a wild-type (wt) rate and appear normal under vegetative growth conditions. With respect to expression, sad-7 is transcribed at baseline levels in early vegetative cultures, at slightly higher levels in mating-competent cultures, and is at its highest level during mating. These findings suggest that SAD-7 is specific to mating-competent and sexual cultures. Although the role of SAD-7 in MSUD remains elusive, green fluorescent protein (GFP)-based tagging studies place SAD-7 within nuclei, perinuclear regions, and cytoplasmic foci of meiotic cells. This localization pattern is unique among known MSUD proteins and raises the possibility that SAD-7 coordinates nuclear, perinuclear, and cytoplasmic aspects of MSUD.

The FgSRP1 SR-protein gene is important for plant infection and pre-mRNA processing in Fusarium graminearum

The versatile functions of SR (serine/arginine-rich) proteins in pre-mRNA splicing and processing are modulated by reversible phosphorylation. Previous studies showed that FgPrp4, the only protein kinase among spliceosome components, is important for intron splicing and the FgSrp1 SR protein is phosphorylated at five conserved sites in Fusarium graminearum. In this study, we showed that the Fgsrp1 deletion mutant rarely produced conidia and caused only limited symptoms on wheat heads and corn silks. Deletion of FgSRP1 also reduced ascospore ejection and deoxynivalenol (DON) production. Interestingly, FgSRP1 had two transcript isoforms due to alternative splicing and both of them were required for its normal functions in growth and DON biosynthesis. FgSrp1 localized to the nucleus and interacted with FgPrp4 in vivo. Deletion of all four conserved phosphorylation sites but not individual ones affected the FgSRP1 function, suggesting their overlapping functions. RNA-seq analysis showed that the expression of over thousands of genes and splicing efficiency in over 140 introns were affected. Taken together, FgSRP1 is important for conidiation, and pathogenesis and alternative splicing is important for its normal functions. The FgSrp1 SR protein is likely important for pre-mRNA processing or splicing of various genes in different developmental and infection processes.

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.

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.

Genome-wide A-to-I RNA editing in fungi independent of ADAR enzymes

Yeasts and filamentous fungi do not have adenosine deaminase acting on RNA (ADAR) orthologs and are believed to lack A-to-I RNA editing, which is the most prevalent editing of mRNA in animals. However, during this study with the PUK1(FGRRES_01058) pseudokinase gene important for sexual reproduction in Fusarium graminearum, we found that two tandem stop codons, UA(1831)GUA(1834)G, in its kinase domain were changed to UG(1831)GUG(1834)G by RNA editing in perithecia. To confirm A-to-I editing of PUK1 transcripts, strand-specific RNA-seq data were generated with RNA isolated from conidia, hyphae, and perithecia. PUK1 was almost specifically expressed in perithecia, and 90% of transcripts were edited to UG(1831)GUG(1834)G. Genome-wide analysis identified 26,056 perithecium-specific A-to-I editing sites. Unlike those in animals, 70.5% of A-to-I editing sites inF. graminearum occur in coding regions, and more than two-thirds of them result in amino acid changes, including editing of 69PUK1-like pseudogenes with stop codons in ORFs.PUK1orthologs and other pseudogenes also displayed stage-specific expression and editing in Neurospora crassa and F. verticillioides Furthermore,F. graminearum differs from animals in the sequence preference and structure selectivity of A-to-I editing sites. Whereas A's embedded in RNA stems are targeted by ADARs, RNA editing inF. graminearum preferentially targets A's in hairpin loops, which is similar to the anticodon loop of tRNA targeted by adenosine deaminases acting on tRNA (ADATs). Overall, our results showed that A-to-I RNA editing occurs specifically during sexual reproduction and mainly in the coding regions in filamentous ascomycetes, involving adenosine deamination mechanisms distinct from metazoan ADARs.

Global Analysis of Predicted G Protein-Coupled Receptor Genes in the Filamentous Fungus, Neurospora crassa

G protein−coupled receptors (GPCRs) regulate facets of growth, development, and environmental sensing in eukaryotes, including filamentous fungi. The largest predicted GPCR class in these organisms is the Pth11-related, with members similar to a protein required for disease in the plant pathogen Magnaporthe oryzae. However, the Pth11-related class has not been functionally studied in any filamentous fungal species. Here, we analyze phenotypes in available mutants for 36 GPCR genes, including 20 Pth11-related, in the model filamentous fungus Neurospora crassa. We also investigate patterns of gene expression for all 43 predicted GPCR genes in available datasets. A total of 17 mutants (47%) possessed at least one growth or developmental phenotype. We identified 18 mutants (56%) with chemical sensitivity or nutritional phenotypes (11 uniquely), bringing the total number of mutants with at least one defect to 28 (78%), including 15 mutants (75%) in the Pth11-related class. Gene expression trends for GPCR genes correlated with the phenotypes observed for many mutants and also suggested overlapping functions for several groups of co-transcribed genes. Several members of the Pth11-related class have phenotypes and/or are differentially expressed on cellulose, suggesting a possible role for this gene family in plant cell wall sensing or utilization.

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.

A Comprehensive Analysis of the Transcriptomes of Marssonina brunnea and Infected Poplar Leaves to Capture Vital Events in Host-Pathogen Interactions

Background

Understanding host-pathogen interaction mechanisms helps to elucidate the entire infection process and focus on important events, and it is a promising approach for improvement of disease control and selection of treatment strategy. Time-course host-pathogen transcriptome analyses and network inference have been applied to unravel the direct or indirect relationships of gene expression alterations. However, time series analyses can suffer from absent time points due to technical problems such as RNA degradation, which limits the application of algorithms that require strict sequential sampling. Here, we introduce an efficient method using independence test to infer an independent network that is exclusively concerned with the frequency of gene expression changes.

Results

Highly resistant NL895 poplar leaves and weakly resistant NL214 leaves were infected with highly active and weakly active Marssonina brunnea, respectively, and were harvested at different time points. The independent network inference illustrated the top 1,000 vital fungus-poplar relationships, which contained 768 fungal genes and 54 poplar genes. These genes could be classified into three categories: a fungal gene surrounded by many poplar genes; a poplar gene connected to many fungal genes; and other genes (possessing low degrees of connectivity). Notably, the fungal gene M6_08342 (a metalloprotease) was connected to 10 poplar genes, particularly including two disease-resistance genes. These core genes, which are surrounded by other genes, may be of particular importance in complicated infection processes and worthy of further investigation.

Conclusions

We provide a clear framework of the interaction network and identify a number of candidate key effectors in this process, which might assist in functional tests, resistant clone selection, and disease control in the future.

Endogenous Small RNA Mediates Meiotic Silencing of a Novel DNA Transposon

Genome defense likely evolved to curtail the spread of transposable elements and invading viruses. A combination of effective defense mechanisms has been shown to limit colonization of the Neurospora crassa genome by transposable elements. A novel DNA transposon named Sly1-1 was discovered in the genome of the most widely used laboratory "wild-type" strain FGSC 2489 (OR74A). Meiotic silencing by unpaired DNA, also simply called meiotic silencing, prevents the expression of regions of the genome that are unpaired during karyogamy. This mechanism is posttranscriptional and is proposed to involve the production of small RNA, so-called masiRNAs, by proteins homologous to those involved in RNA interference-silencing pathways in animals, fungi, and plants. Here, we demonstrate production of small RNAs when Sly1-1 was unpaired in a cross between two wild-type strains. These small RNAs are dependent on SAD-1, an RNA-dependent RNA polymerase necessary for meiotic silencing. We present the first case of endogenously produced masiRNA from a novel N. crassa DNA transposable element.

Neurospora crassa female development requires the PACC and other signal transduction pathways, transcription factors, chromatin remodeling, cell-to-cell fusion, and autophagy

Using a screening protocol we have identified 68 genes that are required for female development in the filamentous fungus Neurospora crassa. We find that we can divide these genes into five general groups: 1) Genes encoding components of the PACC signal transduction pathway, 2) Other signal transduction pathway genes, including genes from the three N. crassa MAP kinase pathways, 3) Transcriptional factor genes, 4) Autophagy genes, and 5) Other miscellaneous genes. Complementation and RIP studies verified that these genes are needed for the formation of the female mating structure, the protoperithecium, and for the maturation of a fertilized protoperithecium into a perithecium. Perithecia grafting experiments demonstrate that the autophagy genes and the cell-to-cell fusion genes (the MAK-1 and MAK-2 pathway genes) are needed for the mobilization and movement of nutrients from an established vegetative hyphal network into the developing protoperithecium. Deletion mutants for the PACC pathway genes palA, palB, palC, palF, palH, and pacC were found to be defective in two aspects of female development. First, they were unable to initiate female development on synthetic crossing medium. However, they could form protoperithecia when grown on cellophane, on corn meal agar, or in response to the presence of nearby perithecia. Second, fertilized perithecia from PACC pathway mutants were unable to produce asci and complete female development. Protein localization experiments with a GFP-tagged PALA construct showed that PALA was localized in a peripheral punctate pattern, consistent with a signaling center associated with the ESCRT complex. The N. crassa PACC signal transduction pathway appears to be similar to the PacC/Rim101 pathway previously characterized in Aspergillus nidulans and Saccharomyces cerevisiae. In N. crassa the pathway plays a key role in regulating female development.