Lessons from the genome sequence of Neurospora crassa: tracing the path from genomic blueprint to multicellular organism

Correlating sugar transporter expression and activities to identify transporters for an orphan sugar substrate

Filamentous fungi like Neurospora crassa are able to take up and metabolize important sugars present, for example, in agricultural and human food wastes. However, only a fraction of all putative sugar transporters in filamentous fungi has been characterized to date, and for many sugar substrates, the corresponding transporters are unknown. In N. crassa, only 14 out of the 42 putative major facilitator superfamily (MFS)-type sugar transporters have been characterized so far. To uncover this hidden potential for biotechnology, it is therefore necessary to find new strategies. By correlation of the uptake profile of sugars of interest after different induction conditions with the expression profiles of all 44 genes encoding predicted sugar transporters in N. crassa, together with an exhaustive phylogenetic analysis using sequences of characterized fungal sugar transporters, we aimed to identify transporter candidates for the tested sugars. Following this approach, we found a high correlation of uptake rates and expression strengths for many sugars with dedicated transporters, like galacturonic acid and arabinose, while the correlation is loose for sugars that are transported by several transporters due to functional redundancy. Nevertheless, this combinatorial approach allowed us to elucidate the uptake system for the disaccharide lactose, a by-product of the dairy industry, which consists of the two main cellodextrin transporters CDT-1 and CDT-2 with a minor contribution of the related transporter NCU00809. Moreover, a non-MFS transporter involved in glycerol transport was also identified. Deorphanization of sugar transporters or identification of transporters for orphan sugar substrates by correlation of uptake kinetics with transporter expression and phylogenetic information can thus provide a way to optimize the reuse of food industry by-products and agricultural wastes by filamentous fungi in order to create economic value and reduce their environmental impact. KEY POINTS: • The Neurospora crassa genome contains 30 uncharacterized putative sugar transporter genes. • Correlation of transporter expression and sugar uptake profiles can help to identify transporters for orphan sugar substrates. • CDT-1, CDT-2, and NCU00809 are key players in the transport of the dairy by-product lactose in N. crassa.

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

The sensor histidine kinase (SLN1) and acetyl-CoA carboxylase (ACC1) coordinately regulate the response of Neurospora crassa to the springtail Sinella curviseta (Collembola: Entomobryidae) attack

IMPORTANCE

Understanding the regulatory pathways by which fungi respond to environmental signals through interlinked genes provides insights into the interactions between fungi and insects. The coordinated optimization of the regulatory networks is necessary for fungi to adapt to their habitats. We demonstrated that the synergistic regulation of sensor histidine kinase (SLN1) and acetyl-CoA carboxylase (ACC1) plays a critical role in regulating the fungal response to Sinella curviseta stress. Furthermore, we found that the enhanced production of trehalose, carotenoids, and 5-MTHF plays crucial role in the resistance to the fungivore. Our results provide insights into the understanding of the adaptation of N. crassa to environmental stimuli.

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.

Histone deacetylation and cytosine methylation compartmentalize heterochromatic regions in the genome organization of Neurospora crassa

Chromosomes must correctly fold in eukaryotic nuclei for proper genome function. Eukaryotic organisms hierarchically organize their genomes, including in the fungus Neurospora crassa, where chromatin fiber loops compact into Topologically Associated Domain (TAD)-like structures formed by heterochromatic region aggregation. However, insufficient data exists on how histone post-translational modifications, including acetylation, affect genome organization. In Neurospora, the HCHC complex (comprised of the proteins HDA-1, CDP-2, HP1, and CHAP) deacetylates heterochromatic nucleosomes, as loss of individual HCHC members increases centromeric acetylation and alters the methylation of cytosines in DNA. Here, we assess if the HCHC complex affects genome organization by performing Hi-C in strains deleted of the cdp-2 or chap genes. CDP-2 loss increases intra- and inter-chromosomal heterochromatic region interactions, while loss of CHAP decreases heterochromatic region compaction. Individual HCHC mutants exhibit different patterns of histone post-translational modifications genome-wide: without CDP-2, heterochromatic H4K16 acetylation is increased, yet smaller heterochromatic regions lose H3K9 trimethylation and gain inter-heterochromatic region interactions; CHAP loss produces minimal acetylation changes but increases heterochromatic H3K9me3 enrichment. Loss of both CDP-2 and the DIM-2 DNA methyltransferase causes extensive genome disorder, as heterochromatic-euchromatic contacts increase despite additional H3K9me3 enrichment. Our results highlight how the increased cytosine methylation in HCHC mutants ensures genome compartmentalization when heterochromatic regions become hyperacetylated without HDAC activity.

The cell functions of phospholipase C-1, Ca2+/H+ exchanger-1, and secretory phospholipase A2 in tolerance to stress conditions and cellulose degradation in Neurospora crassa

We investigated the cell functions of the Ca2+ signaling genes phospholipase C-1 (plc-1), Ca2+/H+ exchanger (cpe-1), and secretory phospholipase A2 (splA2) for stress responses and cellulose utilization in Neurospora crassa. The Δplc-1, Δcpe-1, and ΔsplA2 mutants displayed increased sensitivity to the alkaline pH and reduced survival during induced thermotolerance. The ΔsplA2 mutant also exhibited hypersensitivity to the DTT-induced endoplasmic reticulum (ER) stress, increased microcrystalline cellulose utilization, increased protein secretion, and glucose accumulation in the culture supernatants. Moreover, the ΔsplA2 mutant could not grow on microcrystalline cellulose during ER stress. Furthermore, plc-1, cpe-1, and splA2 synthetically regulate the acquisition of thermotolerance induced by heat shock, responses to alkaline pH and ER stress, and utilization of cellulose and other alternate carbon sources in N. crassa. In addition, expression of the alkaline pH regulator, pac-3, and heat shock proteins, hsp60, and hsp80 was reduced in the Δplc-1, Δcpe-1, and ΔsplA2 single and double mutants. The expression of the unfolded protein response (UPR) markers grp-78 and pdi-1 was also significantly reduced in the mutants showing growth defect during ER stress. The increased cellulolytic activities of the ΔsplA2 and Δcpe-1; ΔsplA2 mutants were due to increased cbh-1, cbh-2, and endo-2 expression in N. crassa. Therefore, plc-1, cpe-1, and splA2 are involved in stress responses and cellulose utilization in N. crassa.

Manipulation of fungal cell wall integrity to improve production of fungal natural products

Fungi, as an important industrial microorganism, play an essential role in the production of natural products (NPs) due to their advantages of utilizing cheap raw materials as substrates and strong protein secretion ability. Although many metabolic engineering strategies have been adopted to enhance the biosynthetic pathway of NPs in fungi, the fungal cell wall as a natural barrier tissue is the final and key step that affects the efficiency of NPs synthesis. To date, many important progresses have been achieved in improving the synthesis of NPs by regulating the cell wall structure of fungi. In this review, we systematically summarize and discuss various strategies for modifying the cell wall structure of fungi to improve the synthesis of NPs. At first, the cell wall structure of different types of fungi is systematically described. Then, strategies to disrupt cell wall integrity (CWI) by regulating the synthesis of cell wall polysaccharides and binding proteins are summarized, which have been applied to improve the synthesis of NPs. In addition, we also summarize the studies on the regulation of CWI-related signaling pathway and the addition of exogenous components for regulating CWI to improve the synthesis of NPs. Finally, we propose the current challenges and essential strategies to usher in an era of more extensive manipulation of fungal CWI to improve the production of fungal NPs.

Calmodulin, Calcium/Calmodulin-Dependent Kinases-1 and 2 Regulate Expression of the Heat Shock Proteins for Heat Shock Tolerance and Pheromone Signaling Genes for Sexual Development in Neurospora crassa

Calmodulin (CaM) is a primary Ca2+ sensor that binds and activates numerous target proteins and regulates several cellular processes in eukaryotes. CaM is essential in Neurospora crassa; therefore, we generated a CaM mutant using repeat-induced point (RIP) mutation and investigated the cmdRIP mutant phenotypes. We also studied knockout mutants of four Ca2+/CaM kinases (camk-1, 2, 3, and 4) for their role during stress conditions and sexual development. The cmdRIP, ∆camk-1, and ∆camk-2 mutants showed reduced survival and growth rates under heat stress, oxidative stress, pH, and ER stress conditions. In addition, under the heat stress conditions, expression of the heat shock protein genes hsp70 and hsp80 was reduced in the cmdRIP, ∆camk-1, and ∆camk-2 mutants. The cmdRIP mutant was also defective in cell fusion, its vegetative hyphae could not support the fertilized wild type perithecia graft, and female sterile. Furthermore, the expression of pheromone signaling genes pre-1, pre-2, ccg-4, mfa-1, and fmf-1 was reduced in the cmdRIP, ∆camk-1, and ∆camk-2 mutants. Therefore, CaM, Ca2+/CaMK-1 and 2 are involved in the tolerance to heat stress conditions and sexual development by regulating the heat shock and pheromone response pathways, respectively, in N. crassa.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-023-01091-8.

Identification of a common secondary mutation in the Neurospora crassa knockout collection conferring a cell fusion-defective phenotype

Gene-deletion mutants represent a powerful tool to study gene function. The filamentous fungus Neurospora crassa is a well-established model organism, and features a comprehensive gene knockout strain collection. While these mutant strains have been used in numerous studies, resulting in the functional annotation of many Neurospora genes, direct confirmation of gene-phenotype relationships is often lacking, which is particularly relevant given the possibility of background mutations, sample contamination, and/or strain mislabeling. Indeed, spontaneous mutations resulting in phenotypes resembling many cell fusion mutants have long been known to occur at relatively high frequency in N. crassa, and these secondary mutations are common in the Neurospora deletion collection. The identity of these mutations, however, is largely unknown. Here, we report that the Δada-3 strain from the N. crassa knockout collection, which exhibits a cell fusion defect, harbors a secondary mutation responsible for this phenotype. Through whole-genome sequencing and genetic analyses, we found a ~30-Kb deletion in this strain affecting a known cell fusion-related gene, so/ham-1, and show that it is the absence of this gene-and not of ada-3-that underlies its cell fusion defect. We additionally found three other knockout strains harboring the same deletion, suggesting that this mutation may be common in the collection and could have impacted previous studies. Our findings provide a cautionary note and highlight the importance of proper functional validation of strains from mutant collections. We discuss our results in the context of the spread of cell fusion-defective cheater variants in N. crassa cultures. IMPORTANCE This study emphasizes the need for careful and detailed characterization of strains from mutant collections. Specifically, we found a common deletion in various strains from the Neurospora crassa gene knockout collection that results in a cell fusion-defective phenotype. This is noteworthy because this collection is known to contain background mutations-of a largely unclear nature-that produce cell fusion-defective phenotypes. Our results describe an example of such mutations, and highlight how this common genetic defect could have impacted previous studies that have used the affected strains. Furthermore, they provide a cautionary note about the use of Neurospora strains with similar phenotypes. Lastly, these findings offer additional details relevant to our understanding of the origin and spread of cell fusion-defective cheater variants in N. crassa cultures.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

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.

Insights into the Antifungal Properties of Myxobacteria Outer Membrane β-1,6-Glucanase

Fungal cell wall decomposition enzymes exhibit great potential for the development of efficient antifungal agents. However, their practical application is restricted due to incomplete understanding of the action mode. In our previous study, we identified that a novel outer membrane (OM) β-1,6-glucanase GluM is deployed by predatory myxobacteria to feed on fungi. In this work, we provide deep insights into the antifungal mechanism of β-1,6-glucanase and its potential in improving plant disease resistance. The fungal cell wall decomposition ability of GluM resulted in irregular hyphae morphology, changed chitin distribution, increased membrane permeability, and leakage of cell constituents in Magnaporthe oryzae Guy11. Under the attack pattern, the cell wall integrity pathway was activated by strain Guy11 for self-protection. GluM exhibited a distinct endo-model toward fungal cell wall; the favorite substrate of GluM toward fungal β-1,6-glucan may give reason for its efficient antifungal activity compared with Trichoderma β-1,6-glucanase. Moreover, released glucans from GluM hydrolysis of fungal cell wall functioned as an elicitor and induced rice immunity by means of jasmonic acid pathway. Based on the dual roles of antifungal properties, gluM transgenic plants conferred enhanced resistance against fungal infection.

The Neurospora crassa PCL-1 cyclin is a PHO85-1 (PGOV) kinase partner that directs the complex to glycogen metabolism and is involved in calcium metabolism regulation

Cyclins are a family of proteins characterized by possessing a cyclin box domain that mediates binding to cyclin dependent kinases (CDKs) partners. In this study, the search for a partner cyclin of the PHO85-1 CDK retrieved PCL-1 an ortholog of yeast Pcls (for Pho85 cyclins) that performs functions common to Pcls belonging to different cyclin families. We show here that PCL-1, as a typical cyclin, is involved in cell cycle control and cell progression. In addition, PCL-1 regulates glycogen metabolism; Δpcl-1 cells accumulate higher glycogen levels than wild-type cells and the glycogen synthase (GSN) enzyme is less phosphorylated and, therefore, more active in the mutant cells. Together with PHO85-1, PCL-1 phosphorylates in vitro GSN at the Ser636 amino acid residue. Modeling studies identified PHO85-1 and PCL-1 as a CDK/cyclin complex, with a conserved intermolecular region stabilized by hydrophobic and polar interactions. PCL-1 is also involved in calcium and NaCl stress response. Δpcl-1 cells are sensitive to high NaCl concentration; on the contrary, they grow better and overexpress calcium responsive genes under high calcium chloride concentration compared to the wild-type strain. The expression of the calcium-responsive CRZ-1 transcription factor is modulated by PCL-1, and this transcription factor seems to be less phosphorylated in Δpcl-1 cells since exhibits nuclear location in these cells in the absence of calcium. Our results show that PCL-1 locates at different cell regions suggesting that it may determine its activity by controlling its intracellular location and reveal an interesting functional divergence between yeast and filamentous fungus cyclins.

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

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

Mechanisms of Antifungal Properties of Metal Nanoparticles

The appearance of resistant species of fungi to the existent antimycotics is challenging for the scientific community. One emergent technology is the application of nanotechnology to develop novel antifungal agents. Metal nanoparticles (NPs) have shown promising results as an alternative to classical antimycotics. This review summarizes and discusses the antifungal mechanisms of metal NPs, including combinations with other antimycotics, covering the period from 2005 to 2022. These mechanisms include but are not limited to the generation of toxic oxygen species and their cellular target, the effect of the cell wall damage and the hyphae and spores, and the mechanisms of defense implied by the fungal cell. Lastly, a description of the impact of NPs on the transcriptomic and proteomic profiles is discussed.

Regulator of G Protein Signaling Proteins Control Growth, Development and Cellulase Production in Neurospora crassa

Heterotrimeric (αβγ) G protein signaling pathways are critical environmental sensing systems found in eukaryotic cells. Exchange of GDP for GTP on the Gα subunit leads to its activation. In contrast, GTP hydrolysis on the Gα is accelerated by Regulator of G protein Signaling (RGS) proteins, resulting in a return to the GDP-bound, inactive state. Here, we analyzed growth, development and extracellular cellulase production in strains with knockout mutations in the seven identified RGS genes (rgs-1 to rgs-7) in the filamentous fungus, Neurospora crassa. We compared phenotypes to those of strains with either knockout mutations or expressing predicted constitutively activated, GTPase-deficient alleles for each of the three Gα subunit genes (gna-1^Q204L, gna-2^Q205L or gna-3^Q208L). Our data revealed that six RGS mutants have taller aerial hyphae than wild type and all seven mutants exhibit reduced asexual sporulation, phenotypes shared with strains expressing the gna-1^Q204L or gna-3^Q208L allele. In contrast, Δrgs-1 and Δrgs-3 were the only RGS mutants with a slower growth rate phenotype, a defect in common with gna-1^Q204L strains. With respect to female sexual development, Δrgs-1 possessed defects most similar to gna-3^Q208L strains, while those of Δrgs-2 mutants resembled strains expressing the gna-1^Q204L allele. Finally, we observed that four of the seven RGS mutants had significantly different extracellular cellulase levels relative to wild type. Of interest, the Δrgs-2 mutant had no detectable activity, similar to the gna-3^Q208L strain. In contrast, the Δrgs-1 and Δrgs-4 mutants and gna-1^Q204L and gna-2^Q205L strains exhibited significantly higher cellulase activity than wild type. With the exception of sexual development, our results demonstrate the greatest number of genetic interactions between rgs-1 and gna-1 and rgs-2 and gna-3 in N. crassa.

Trichoderma - genomes and genomics as treasure troves for research towards biology, biotechnology and agriculture

The genus Trichoderma is among the best studied groups of filamentous fungi, largely because of its high relevance in applications from agriculture to enzyme biosynthesis to biofuel production. However, the physiological competences of these fungi, that led to these beneficial applications are intriguing also from a scientific and ecological point of view. This review therefore summarizes recent developments in studies of fungal genomes, updates on previously started genome annotation efforts and novel discoveries as well as efforts towards bioprospecting for enzymes and bioactive compounds such as cellulases, enzymes degrading xenobiotics and metabolites with potential pharmaceutical value. Thereby insights are provided into genomes, mitochondrial genomes and genomes of mycoviruses of Trichoderma strains relevant for enzyme production, biocontrol and mycoremediation. In several cases, production of bioactive compounds could be associated with responsible genes or clusters and bioremediation capabilities could be supported or predicted using genome information. Insights into evolution of the genus Trichoderma revealed large scale horizontal gene transfer, predominantly of CAZyme genes, but also secondary metabolite clusters. Investigation of sexual development showed that Trichoderma species are competent of repeat induced point mutation (RIP) and in some cases, segmental aneuploidy was observed. Some random mutants finally gave away their crucial mutations like T. reesei QM9978 and QM9136 and the fertility defect of QM6a was traced back to its gene defect. The Trichoderma core genome was narrowed down to 7000 genes and gene clustering was investigated in the genomes of multiple species. Finally, recent developments in application of CRISPR/Cas9 in Trichoderma, cloning and expression strategies for the workhorse T. reesei as well as the use genome mining tools for bioprospecting Trichoderma are highlighted. The intriguing new findings on evolution, genomics and physiology highlight emerging trends and illustrate worthwhile perspectives in diverse fields of research with Trichoderma.

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.

HDA-2-Containing Complex Is Required for Activation of Catalase-3 Expression in Neurospora crassa

It is essential for aerobic organisms to maintain the homeostasis of intracellular reactive oxygen species (ROS) for survival and adaptation to the environment. In line with other eukaryotes, the catalase of Neurospora crassa is an important enzyme for clearing ROS, and its expression is tightly regulated by the growth phase and various oxidative stresses. Our study reveals that, in N. crassa, histone deacetylase 2 (HDA-2) and its catalytic activity positively regulate the expression of the catalase-3 (cat-3) gene. HDA-2, SIF-2, and SNT-1 may form a subcomplex with such a regulation role. As expected, deletion of HDA-2 or SIF-2 subunit increased acetylation levels of histone H4, indicating that loss of HDA-2 complex fails to deacetylate H4 at the cat-3 locus. Furthermore, loss of HDA-2 or its catalytic activity led to dramatic decreases of TFIIB and RNA polymerase II (RNAP II) recruitment at the cat-3 locus and also resulted in high deposition of H2A.Z at the promoter and transcription start site (TSS) regions of the cat-3 gene. Collectively, this study strongly demonstrates that the HDA-2-containing complex activates the transcription of the cat-3 gene by facilitating preinitiation complex (PIC) assembly and antagonizing the inhibition of H2A.Z at the cat-3 locus through H4 acetylation. IMPORTANCE Clearance of reactive oxygen species (ROS) is critical to the survival of aerobic organisms. In the model filamentous fungus Neurospora crassa, catalase-3 (cat-3) expression is activated in response to H₂O₂-induced ROS stress. We found that histone deacetylase 2 (HDA-2) positively regulates cat-3 transcription in N. crassa; this is widely divergent from the classical repressive role of most histone deacetylases. Like HDA-2, the SIF-2 or SNT-1 subunit of HDA-2-containing complex plays a positive role in cat-3 transcription. Furthermore, we also found that HDA-2-containing complex provides an appropriate chromatin environment to facilitate PIC assembly and to antagonize the inhibition role of H2A.Z at the cat-3 locus through H4 acetylation. Taken together, our results establish a mechanism for how the HDA-2-containing complex regulates transcription of the cat-3 gene in N. crassa.

Transcription Factor Mavib-1 Negatively Regulates Conidiation by Affecting Utilization of Carbon and Nitrogen Source in Metarhizium acridum

Conidium is the main infection unit and reproductive unit of pathogenic fungi. Exploring the mechanism of conidiation and its regulation contributes to understanding the pathogenicity of pathogenic fungi. Vib-1, a transcription factor, was reported to participate in the conidiation process. However, the regulation mechanism of Vib-1 in conidiation is still unclear. In this study, we analyzed the function of Vib-1 and its regulation mechanism in conidiation through knocking out and overexpression of Vib-1 in entomopathogenic fungus Metarhizium acridum. Results showed that the colonial growth of Mavib-1 disruption mutant (ΔMavib-1) was significantly decreased, and conidiation was earlier compared to wild type (WT), while overexpression of Mavib-1 led to a delayed conidiation especially when carbon or nitrogen sources were insufficient. Overexpression of Mavib-1 resulted in a conidiation pattern shift from microcycle conidiation to normal conidiation on nutrient-limited medium. These results indicated that Mavib-1 acted as a positive regulator in vegetative growth and a negative regulator in conidiation by affecting utilization of carbon and nitrogen sources in M. acridum. Transcription profile analysis demonstrated that many genes related to carbon and nitrogen source metabolisms were differentially expressed in ΔMavib-1 and OE strains compared to WT. Moreover, Mavib-1 affects the conidial germination, tolerance to UV-B and heat stresses, cell wall integrity, conidial surface morphology and conidial hydrophobicity in M. acridum. These findings unravel the regulatory mechanism of Mavib-1 in fungal growth and conidiation, and enrich the knowledge to conidiation pattern shift of filamentous fungi.

Development of an Efficient C-to-T Base-Editing System and Its Application to Cellulase Transcription Factor Precise Engineering in Thermophilic Fungus Myceliophthora thermophila

Myceliophthora thermophila is a thermophilic fungus with great potential in biorefineries and biotechnology. The base editor is an upgraded version of the clustered regularly interspaced short palindromic repeats (CRISPR)-dependent genome-editing tool that introduces precise point mutations without causing DNA double-strand breaks (DSBs) and has been used in various organisms but rarely in filamentous fungi, especially thermophilic filamentous fungi. Here, for the first time, we constructed three cytosine base editors (CBEs) in M. thermophila, namely, evolved apolipoprotein B mRNA-editing enzyme catalytic subunit 1 (APOBEC1) cytosine base editor 4 max (Mtevo-BE4max), bacteriophage Mu Gam protein cytosine base editor 4 max (MtGAM-BE4max), and evolved CDA1 deaminase cytosine base editor (Mtevo-CDA1), and efficiently inactivated genes by precisely converting three codons (CAA, CAG, and CGA) into stop codons without DSB formation. The Mtevo-CDA1 editor with up to 92.6% editing efficiency is a more suitable tool for cytosine base editing in thermophilic fungi. To investigate the function of each motif of the cellulase transcription factor M. thermophila CLR-2 (MtCLR-2), we used the Mtevo-CDA1 editor. The fungal-specific motif of MtCLR-2 was found to be strongly involved in cellulase secretion, conidium formation, hyphal branching, and colony formation. Mutation of the fungus-specific motif caused significant defects in these characteristics. Thus, we developed an efficient thermophilic fungus-compatible base-editing system that could also be used for genetic engineering in other relevant filamentous fungi.

IMPORTANCE A CRISPR/Cas-based base-editing approach has been developed to introduce point mutations without inducing double-strand breaks (DSBs) and attracted substantial academic and industrial interest. Our study developed the deaminase-cytosine base-editing system to efficiently edit three target genes, amdS, cre-1, and the essential cellulase regulator gene Mtclr-2, in Myceliophthora thermophila. A variety of point mutations in the target loci of the DNA-binding domain and fungus-specific motif of M. thermophila CLR-2 (MtCLR-2) were successfully generated via our base editor Mtevo-CDA1 to elucidate its function. Here, we show that the DNA-binding domain of MtCLR-2 is important for the fungal response to cellulose conditions, while its fungus-specific motif is involved in fungal growth. These findings indicate that our base editor can be an effective tool for elucidating the functions of motifs of target genes in filamentous fungi and for metabolic engineering in the field of synthetic biology.

Diversity of biologically active secondary metabolites in the ascomycete order Sordariales

Ascomycetes belonging to the order Sordariales are a well-known reservoir of secondary metabolites with potential beneficial applications. Species of the Sordariales are ubiquitous, and they are commonly found in soils and in lignicolous, herbicolous, and coprophilous habitats. Some of their species have been used as model organisms in modern fungal biology or were found to be prolific producers of potentially useful secondary metabolites. However, the majority of sordarialean species are poorly studied. Traditionally, the classification of the Sordariales has been mainly based on morphology of the ascomata, ascospores, and asexual states, characters that have been demonstrated to be homoplastic by modern taxonomic studies based on multi-locus phylogeny. Herein, we summarize for the first time relevant information about the available knowledge on the secondary metabolites and the biological activities exerted by representatives of this fungal order, as well as a current outlook of the potential opportunities that the recent advances in omic tools could bring for the discovery of secondary metabolites in this order.

Heat shock proteins and the calcineurin-crz1 signaling regulate stress responses in fungi

The heat shock proteins (Hsps) act as a molecular chaperone to stabilize client proteins involved in various cell functions in fungi. Hsps are classified into different families such as HSP90, HSP70, HSP60, HSP40, and small HSPs (sHsps). Hsp90, a well-studied member of the Hsp family proteins, plays a role in growth, cell survival, and pathogenicity in fungi. Hsp70 and sHsps are involved in the development, tolerance to stress conditions, and drug resistance in fungi. Hsp60 is a mitochondrial chaperone, and Hsp40 regulates fungal ATPase machinery. In this review, we describe the cell functions, regulation, and the molecular link of the Hsps with the calcineurin-crz1 calcium signaling pathway for their role in cell survival, growth, virulence, and drug resistance in fungi and related organisms.

The ACF chromatin-remodeling complex is essential for Polycomb repression

Establishing and maintaining appropriate gene repression is critical for the health and development of multicellular organisms. Histone H3 lysine 27 (H3K27) methylation is a chromatin modification associated with repressed facultative heterochromatin, but the mechanism of this repression remains unclear. We used a forward genetic approach to identify genes involved in transcriptional silencing of H3K27-methylated chromatin in the filamentous fungus Neurospora crassa. We found that the N. crassa homologs of ISWI (NCU03875) and ACF1 (NCU00164) are required for repression of a subset of H3K27-methylated genes and that they form an ACF chromatin-remodeling complex. This ACF complex interacts with chromatin throughout the genome, yet association with facultative heterochromatin is specifically promoted by the H3K27 methyltransferase, SET-7. H3K27-methylated genes that are upregulated when iswi or acf1 are deleted show a downstream shift of the +1 nucleosome, suggesting that proper nucleosome positioning is critical for repression of facultative heterochromatin. Our findings support a direct role of the ACF complex in Polycomb repression.

Not all Is SET for Methylation: Evolution of Eukaryotic Protein Methyltransferases

Dynamic posttranslational modifications to canonical histones that constitute the nucleosome (H2A, H2B, H3, and H4) control all aspects of enzymatic transactions with DNA. Histone methylation has been studied heavily for the past 20 years, and our mechanistic understanding of the control and function of individual methylation events on specific histone arginine and lysine residues has been greatly improved over the past decade, driven by excellent new tools and methods. Here, we will summarize what is known about the distribution and some of the functions of protein methyltransferases from all major eukaryotic supergroups. The main conclusion is that protein, and specifically histone, methylation is an ancient process. Many taxa in all supergroups have lost some subfamilies of both protein arginine methyltransferases (PRMT) and the heavily studied SET domain lysine methyltransferases (KMT). Over time, novel subfamilies, especially of SET domain proteins, arose. We use the interactions between H3K27 and H3K36 methylation as one example for the complex circuitry of histone modifications that make up the "histone code," and we discuss one recent example (Paramecium Ezl1) for how extant enzymes that may resemble more ancient SET domain KMTs are able to modify two lysine residues that have divergent functions in plants, fungi, and animals. Complexity of SET domain KMT function in the well-studied plant and animal lineages arose not only by gene duplication but also acquisition of novel DNA- and histone-binding domains in certain subfamilies.

Complete Genome Sequences and Genome-Wide Characterization of Trichoderma Biocontrol Agents Provide New Insights into their Evolution and Variation in Genome Organization, Sexual Development, and Fungal-Plant Interactions

Trichoderma spp. represent one of the most important fungal genera to mankind and in natural environments. The genus harbors prolific producers of wood-decaying enzymes, biocontrol agents against plant pathogens, plant-growth-promoting biofertilizers, as well as model organisms for studying fungal-plant-plant pathogen interactions. Pursuing highly accurate, contiguous, and chromosome-level reference genomes has become a primary goal of fungal research communities. Here, we report the chromosome-level genomic sequences and whole-genome annotation data sets of four strains used as biocontrol agents or biofertilizers (Trichoderma virens Gv29-8, Trichoderma virens FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1). Our results provide comprehensive categorization, correct positioning, and evolutionary detail of both nuclear and mitochondrial genomes, including telomeres, AT-rich blocks, centromeres, transposons, mating-type loci, nuclear-encoded mitochondrial sequences, as well as many new secondary metabolic and carbohydrate-active enzyme gene clusters. We have also identified evolutionarily conserved core genes contributing to plant-fungal interactions, as well as variations potentially linked to key behavioral traits such as sex, genome defense, secondary metabolism, and mycoparasitism. The genomic resources we provide herein significantly extend our knowledge not only of this economically important fungal genus, but also fungal evolution and basic biology in general. IMPORTANCE Telomere-to-telomere and gapless reference genome assemblies are necessary to ensure that all genomic variants are studied and discovered, including centromeres, telomeres, AT-rich blocks, mating type loci, biosynthetic, and metabolic gene clusters. Here, we applied long-range sequencing technologies to determine the near-completed genome sequences of four widely used biocontrol agents or biofertilizers: Trichoderma virens Gv29-8 and FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1. Like those of three Trichoderma reesei wild isolates [QM6a, CBS999.97(MAT1-1) and CBS999.97(MAT1-2)] we reported previously, these four biocontrol agent genomes each contain seven nuclear chromosomes and a circular mitochondrial genome. Substantial intraspecies and intragenus diversities are also discovered, including single nucleotide polymorphisms, chromosome shuffling, as well as genomic relics derived from historical transposition events and repeat-induced point (RIP) mutations.

Fungal Lysine Deacetylases in Virulence, Resistance, and Production of Small Bioactive Compounds

The growing number of immunocompromised patients begs for efficient therapy strategies against invasive fungal infections. As conventional antifungal treatment is increasingly hampered by resistance to commonly used antifungals, development of novel therapy regimens is required. On the other hand, numerous fungal species are industrially exploited as cell factories of enzymes and chemicals or as producers of medically relevant pharmaceuticals. Consequently, there is immense interest in tapping the almost inexhaustible fungal portfolio of natural products for potential medical and industrial applications. Both the pathogenicity and production of those small metabolites are significantly dependent on the acetylation status of distinct regulatory proteins. Thus, classical lysine deacetylases (KDACs) are crucial virulence determinants and important regulators of natural products of fungi. In this review, we present an overview of the members of classical KDACs and their complexes in filamentous fungi. Further, we discuss the impact of the genetic manipulation of KDACs on the pathogenicity and production of bioactive molecules. Special consideration is given to inhibitors of these enzymes and their role as potential new antifungals and emerging tools for the discovery of novel pharmaceutical drugs and antibiotics in fungal producer strains.

Revisiting the Neurospora crassa mitochondrial genome

The mitochondrial genome of Neurospora crassa has been less studied than its nuclear counterpart, yet it holds great potential for understanding the diversity and evolution of this important fungus. Here we describe a new mitochondrial DNA (mtDNA) complete sequence of a N. crassa wild type strain. The genome with 64 839 bp revealed 21 protein-coding genes and several hypothetical open reading frames with no significant homology to any described gene. Five large repetitive regions were identified across the genome, including partial or complete genes. The largest repeated region holds a partial nd2 section that was also detected in Neurospora intermedia, suggesting a rearrangement that occurred before the N. crassa speciation. Interestingly, N. crassa has a palindrome adjacent to the partial nd2 repeated region possibly related to the genomic rearrangement, which is absent in N. intermedia. Finally, we compared the sequences of the three available N. crassa complete mtDNAs and found low levels of intraspecific variability. Most differences among strains were due to small indels in noncoding regions. The revisiting of the N. crassa mtDNA forms the basis for future studies on mitochondrial genome organization and variability.

Chromatin accessibility profiling in Neurospora crassa reveals molecular features associated with accessible and inaccessible chromatin

BACKGROUND

Regulation of chromatin accessibility and transcription are tightly coordinated processes. Studies in yeast and higher eukaryotes have described accessible chromatin regions, but little work has been done in filamentous fungi.

RESULTS

Here we present a genome-scale characterization of accessible chromatin regions in Neurospora crassa, which revealed characteristic molecular features of accessible and inaccessible chromatin. We present experimental evidence of inaccessibility within heterochromatin regions in Neurospora, and we examine features of both accessible and inaccessible chromatin, including the presence of histone modifications, types of transcription, transcription factor binding, and relative nucleosome turnover rates. Chromatin accessibility is not strictly correlated with expression level. Accessible chromatin regions in the model filamentous fungus Neurospora are characterized the presence of H3K27 acetylation and commonly associated with pervasive non-coding transcription. Conversely, methylation of H3 lysine-36 catalyzed by ASH1 is correlated with inaccessible chromatin within promoter regions.

CONCLUSIONS

In N. crassa, H3K27 acetylation is the most predictive histone modification for open chromatin. Conversely, our data show that H3K36 methylation is a key marker of inaccessible chromatin in gene-rich regions of the genome. Our data are consistent with an expanded role for H3K36 methylation in intergenic regions of filamentous fungi compared to the model yeasts, S. cerevisiae and S. pombe, which lack homologs of the ASH1 methyltransferase.

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.

IMITATION SWITCH is required for normal chromatin structure and gene repression in PRC2 target domains

Polycomb Group (PcG) proteins are part of an epigenetic cell memory system that plays essential roles in multicellular development, stem cell biology, X chromosome inactivation, and cancer. In animals, plants, and many fungi, Polycomb Repressive Complex 2 (PRC2) catalyzes trimethylation of histone H3 lysine 27 (H3K27me3) to assemble transcriptionally repressed facultative heterochromatin. PRC2 is structurally and functionally conserved in the model fungus Neurospora crassa, and recent work in this organism has generated insights into PRC2 control and function. To identify components of the facultative heterochromatin pathway, we performed a targeted screen of Neurospora deletion strains lacking individual ATP-dependent chromatin remodeling enzymes. We found the Neurospora homolog of IMITATION SWITCH (ISW) is critical for normal transcriptional repression, nucleosome organization, and establishment of typical histone methylation patterns in facultative heterochromatin domains. We also found that stable interaction between PRC2 and chromatin depends on ISW. A functional ISW ATPase domain is required for gene repression and normal H3K27 methylation. ISW homologs interact with accessory proteins to form multiple complexes with distinct functions. Using proteomics and molecular approaches, we identified three distinct Neurospora ISW-containing complexes. A triple mutant lacking three ISW accessory factors and disrupting multiple ISW complexes led to widespread up-regulation of PRC2 target genes and altered H3K27 methylation patterns, similar to an ISW-deficient strain. Taken together, our data show that ISW is a key component of the facultative heterochromatin pathway in Neurospora, and that distinct ISW complexes perform an apparently overlapping role to regulate chromatin structure and gene repression at PRC2 target domains.

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.

A novel antisense long non-coding RNA participates in asexual and sexual reproduction by regulating the expression of GzmetE in Fusarium graminearum

Fusarium graminearum is an important worldwide pathogen that causes Fusarium head blight in wheat, barley, maize and other grains. LncRNAs play important roles in many biological processes, but little is known about their functions and mechanisms in filamentous fungi. Here, we report that a natural antisense RNA, GzmetE-AS, is transcribed from the opposite strand of GzmetE. GzmetE encodes a homoserine O-acetyltransferase, which is important for sexual development and plant infection. The expression of GzmetE-AS was increased significantly during the conidiation stage, while GzmetE was upregulated in the late stage of sexual reproduction. Overexpression of GzmetE-AS inhibited the transcription of GzmetE. In contrast, the expression of GzmetE was significantly increased in GzmetE-AS transcription termination strain GzmetE-AS-T. Furthermore, GzmetE-AS-T produced more perithecia and facilitated the ascospore discharge, resembling the phenotype of GzmetE overexpressing strains. However, overexpression of GzmetE-AS in ∆dcl1/2 strain cannot inhibit the expression of GzmetE, and the GzmetE nat-siRNA is also significantly reduced in ∆dcl1/2 mutant. Taken together, we have identified a novel antisense lncRNA GzmetE-AS, which is involved in asexual and sexual reproduction by regulating its antisense gene GzmetE through RNAi pathway. Our findings reveal that the lncRNA plays critical roles in the development of F. graminearum.

Gene diversity explains variation in biological features of insect killing fungus, Beauveria bassiana

Beauveria bassiana is a species complex whose isolates show considerable natural genetic variability. However, little is known about how this genetic diversity affects the fungus performance. Herein, we characterized the diversity of genes involved in various mechanisms of the infective cycle of 42 isolates that have different growth rates, thermotolerance and virulence. The analysed genes showed general genetic diversity measured as non-synonymous changes (NSC) and copy number variation (CNV), with most of them being subjected to positive episodic diversifying selection. Correlation analyses between NSC or CNV and the isolate virulence, thermotolerance and growth rate revealed that various genes shaped the biological features of the fungus. Lectin-like, mucin signalling, Biotrophy associated and chitinase genes NSCs correlated with the three biological features of B. bassiana. In addition, other genes (i.e. DNA photolyase and cyclophilin B) that had relatively conserved sequences, had variable CNs across the isolates which were correlated with the variability of either virulence or thermotolerance of B. bassiana isolates. The data obtained is important for a better understanding of population structure, ecological and potential impact when isolates are used as mycoinsecticides and can justify industrialization of new isolates.

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.

Distribution, evolution and expression of GATA-TFs provide new insights into their functions in light response and fruiting body development of Tolypocladium guangdongense

Background

Fungal GATA-type transcription factors (GATA-TFs) are a class of transcriptional regulators involved in various biological processes. However, their functions are rarely analyzed systematically, especially in edible or medicinal fungi, such as Tolypocladium guangdongense, which has various medicinal and food safety properties with a broad range of potential applications in healthcare products and the pharmaceutical industry.

Methods

GATA-TFs in T. guangdongense (TgGATAs) were identified using InterProScan. The type, distribution, and gene structure of TgGATAs were analyzed by genome-wide analyses. A phylogenetic tree was constructed to analyze their evolutionary relationships using the neighbor-joining (NJ) method. To explore the functions of GATA-TFs, conserved domains were analyzed using MEME, and cis-elements were predicted using the PlantCARE database. In addition, the expression patterns of TgGATAs under different light conditions and developmental stages were studied using qPCR.

Results

Seven TgGATAs were identified. They were randomly distributed on four chromosomes and contained one to four exons. Phylogenetic analysis indicated that GATA-TFs in each subgroup are highly conserved, especially for GATA1 to GATA5. Intron distribution analyses suggested that GATA1 and GATA3 possessed the most conserved gene structures. Light treatments induced the expression levels of TgGATA1 and TgGATA5-7, but the expression levels varied depending on the duration of illumination. The predicted protein structures indicate that TgGATA1 and TgGATA2 possess typical light-responsive domains and may function as photoreceptors to regulate downstream biological processes. TgGATA3 and TgGATA5 may be involved in nitrogen metabolism and siderophore biosynthesis, respectively. TgGATA6 and TgGATA7 possess unique Zn finger loop sequences, suggesting that they may have special functions. Furthermore, gene expression analysis indicated that TgGATA1 (WC1) was notably involved in mycelial color transformation, while other genes were involved in fruiting body development to some extent. These results provide valuable information to further explore the mechanisms through which TgGATAs are regulated during fruiting body development.

Normal Patterns of Histone H3K27 Methylation Require the Histone Variant H2A.Z in Neurospora crassa

Neurospora crassa contains a minimal Polycomb repression system, which provides rich opportunities to explore Polycomb-mediated repression across eukaryotes and enables genetic studies that can be difficult in plant and animal systems. Polycomb Repressive Complex 2 is a multi-subunit complex that deposits mono-, di-, and trimethyl groups on lysine 27 of histone H3, and trimethyl H3K27 is a molecular marker of transcriptionally repressed facultative heterochromatin. In mouse embryonic stem cells and multiple plant species, H2A.Z has been found to be colocalized with H3K27 methylation. H2A.Z is required for normal H3K27 methylation in these experimental systems, though the regulatory mechanisms are not well understood. We report here that Neurospora crassa mutants lacking H2A.Z or SWR-1, the ATP-dependent histone variant exchanger, exhibit a striking reduction in levels of H3K27 methylation. RNA-sequencing revealed downregulation of eed, encoding a subunit of PRC2, in an hH2Az mutant compared to wild type, and overexpression of EED in a ΔhH2Az;Δeed background restored most H3K27 methylation. Reduced eed expression leads to region-specific losses of H3K27 methylation, suggesting that differential dependence on EED concentration is critical for normal H3K27 methylation at certain regions in the genome.

Calcium signaling is involved in diverse cellular processes in fungi

Calcium (Ca²⁺) is a universal signalling molecule of life. The Ca²⁺ signalling is an evolutionarily conserved process from prokaryotes to eukaryotes. Ca²⁺ at high concentration is deleterious to the cell; therefore, cell maintains a low resting level of intracellular free Ca²⁺ concentration ([Ca²⁺]_c). The resting [Ca²⁺]_c is tightly regulated, and a transient increase of the [Ca²⁺]_c initiates a signalling cascade in the cell. Ca²⁺ signalling plays an essential role in various processes, including growth, development, reproduction, tolerance to stress conditions, and virulence in fungi. In this review, we describe the evolutionary aspects of Ca²⁺ signalling and cell functions of major Ca²⁺ signalling proteins in different fungi.

Polycomb Repression without Bristles: Facultative Heterochromatin and Genome Stability in Fungi

Genome integrity is essential to maintain cellular function and viability. Consequently, genome instability is frequently associated with dysfunction in cells and associated with plant, animal, and human diseases. One consequence of relaxed genome maintenance that may be less appreciated is an increased potential for rapid adaptation to changing environments in all organisms. Here, we discuss evidence for the control and function of facultative heterochromatin, which is delineated by methylation of histone H3 lysine 27 (H3K27me) in many fungi. Aside from its relatively well understood role in transcriptional repression, accumulating evidence suggests that H3K27 methylation has an important role in controlling the balance between maintenance and generation of novelty in fungal genomes. We present a working model for a minimal repressive network mediated by H3K27 methylation in fungi and outline challenges for future research.

Quantification of Neurospora crassa mitochondrial DNA using quantitative real-time PCR

The filamentous fungus Neurospora crassa is a popular model organism used in a wide range of biochemical and genetic studies and vastly used in mitochondrial research. Despite the relevance of mitochondria in N. crassa biology, no method for quantification of mitochondrial DNA (mtDNA) is currently available. Quantitative real-time PCR (qPCR) is a powerful tool, with a wide range of applications, and has been used for the quantification of nucleic acids in humans and a few other species. Here we present a new qPCR assay for relative quantification of N. crassa mtDNA. Three sets of qPCR primers targeting different regions of the mitochondrial genome were tested for mtDNA quantification. The qPCR was successfully validated in N. crassa strains from different geographical locations, representing the vast genetic diversity of this species, and knockout mutant strains. Moreover the assay proved to be efficient in templates with varied amounts of mitochondria, obtained through different DNA extraction methods. The qPCR performed well in all tested samples revealing a higher amount of mtDNA than nuclear DNA in all cases. This technique will facilitate the characterization of mtDNA of N. crassa in future studies and can be used as a tool to validate methods of mitochondria isolation. SIGNIFICANCE AND IMPACT OF THE STUDY: The standardization of quantitative real-time PCR (qPCR) techniques is essential to enable and facilitate future comparisons. Neurospora crassa is a model organism with a lot of potential in different fields of study. Here we use N. crassa to develop and establish an assay to quantify mitochondrial DNA using qPCR. We tested strains with different geographical background and our data demonstrated the usefulness of this assay to quantify mitochondrial DNA in N. crassa. This technique can be useful in a wide variety of applications and in different types of studies.

Two dominant selectable markers for genetic manipulation in Neurospora crassa

Neurospora crassa is an excellent model fungus for studies on molecular genetics, biochemistry, physiology, and molecular cell biology. Along with the rapid progress of Neurospora research, new tools facilitating more efficient and accurate genetic analysis are in high demand. Here, we tested whether the dominant selective makers widely used in yeasts are applicable in N. crassa. Among them, we found that the strains of N. crassa are sensitive to the aminoglycoside antibiotics, G418 and nourseothricin. 1000 μg/mL of G418 or 50 μg/mL of nourseothricin is sufficient to inhibit Neurospora growth completely. When the neomycin phosphotransferase gene (neo) used in mammalian cells is expressed, N. crassa shows potent resistance to G418. This establishes G418-resistant marker as a dominant selectable marker to use in N. crassa. Similarly, when the nourseothricin acetyltransferase gene (nat) from Streptomyces noursei is induced by qa-2 promoter in the presence of quinic acid (QA), N. crassa shows potent resistance to nourseothricin. When nat is constitutively expressed by full-length or truncated versions of the promoter from the N. crassa cfp gene (NCU02193), or by the trpC promoter of Aspergillus nidulans, the growth of N. crassa in the presence of nourseothricin is proportional to the expression levels of Nat. Finally, these two markers are used to knock-out wc-2 or al-1 gene from the N. crassa genome. The successful development of these two markers in this study expands the toolbox for N. crassa and very likely for other filamentous fungi as well.

WHI-2 Regulates Intercellular Communication via a MAP Kinase Signaling Complex

The formation of the fungal mycelial network is facilitated by somatic cell fusion of germinating asexual spores (or germlings). Neurospora crassa germlings in close proximity display chemotropic growth that is dependent upon an intracellular network of mitogen-activated protein kinase (MAPK) signaling cascades. Approximately 80 genes involved in intercellular communication and fusion have been identified, including three mutants with similar morphological phenotypes: Δwhi-2, Δcsp-6, and Δamph-1. Here we show that WHI-2 localizes to the cell periphery and regulates endocytosis, mitochondrial organization, sporulation, and cell fusion. WHI-2 was required to transduce signals through a conserved MAPK pathway (NRC-1/MEK-2/MAK-2) and target transcription factors (PP-1/ADV-1). The amph-1 locus encodes a Bin/Amphiphysin/Rvs domain-containing protein and mis-expression of whi-2 compensated for the cell fusion and endocytosis deficiencies of a Δamph-1 mutant. The csp-6 locus encodes a haloacid dehalogenase phosphatase whose activity was essential for cell fusion. Although fusion-deficient with themselves, cells that lacked whi-2, csp-6, or amph-1 showed a low frequency of chemotropic interactions with wild type cells. We hypothesize that WHI-2 could be important for signal perception during chemotropic interactions via a role in endocytosis.

Comparative transcriptome analysis revealed genes involved in the fruiting body development of Ophiocordyceps sinensis

Ophiocordyceps sinensis is a highly valued fungus that has been used as traditional Asian medicine. This fungus is one of the most important sources of income for the nomadic populations of the Tibetan Plateau. With global warming and excessive collection, the wild O. sinensis resources declined dramatically. The cultivation of O. sinensis hasn’t been fully operational due to the unclear genetic basis of the fruiting body development. Here, our study conducted pairwise comparisons between transcriptomes acquired from different growth stages of O. sinensis including asexual mycelium (CM), developing fruiting body (DF) and mature fruiting body (FB). All RNA-Seq reads were aligned to the genome of O. sinensis CO18 prior to comparative analyses. Cluster analysis showed that the expression profiles of FB and DF were highly similar compared to CM. Alternative splicing analysis (AS) revealed that the stage-specific splicing genes may have important functions in the development of fruiting body. Functional enrichment analyses showed that differentially expressed genes (DEGs) were enriched in protein synthesis and baseline metabolism during fruiting body development, indicating that more protein and energy might be required for fruiting body development. In addition, some fruiting body development-associated genes impacted by ecological factors were up-regulated in FB samples, such as the nucleoside diphosphate kinase gene (ndk), β subunit of the fatty acid synthase gene (cel-2) and the superoxide dismutase gene (sod). Moreover, the expression levels of several cytoskeletons genes were significantly altered during all these growth stages, suggesting that these genes play crucial roles in both vegetative growth and the fruiting body development. Quantitative PCR (qPCR) was used to validate the gene expression profile and the results supported the accuracy of the RNA-Seq and DEGs analysis. Our study offers a novel perspective to understand the underlying growth stage-specific molecular differences and the biology of O. sinensis fruiting body development.

G-protein-coupled Receptors in Fungi

G-protein-coupled receptors (GPCRs) are the largest family of transmembrane receptors in fungi. These receptors have an important role in the transduction of extracellular signals into intracellular sites in response to diverse stimuli. They enable fungi to coordinate cell function and metabolism, thereby promoting their survival and propagation, and sense certain fundamentally conserved elements, such as nutrients, pheromones, and stress, for adaptation to their niches, environmental stresses, and host environment, causing disease and pathogen virulence. This chapter highlights the role of GPCRs in fungi in coordinating cell function and metabolism. Fungal cells sense the molecular interactions between extracellular signals. Their respective sensory systems are described here in detail.

A Comparison of the Neurospora and Drosophila Clocks

In Neurospora and other fungi, the protein frequency (FRQ) is an integral part and a negative element in the fungal circadian oscillator. In Drosophila and many other higher organisms, the protein period (PER) is an integral part and a negative element of their circadian oscillator. Employing bioinformatic techniques, such as BLAST, CLUSTAL, and MEME (Multiple Em for Motif Elicitation), 11 regions (sequences) of potential similarity were found between the fungal FRQ and the Drosophila PER. Many of these FRQ regions are conserved in many fungal FRQ(s). Many of these PER regions are conserved in many insects. In addition, these regions are also of biological significance since mutations in these regions lead to changes in the circadian clock of Neurospora and Drosophila. Many of these regions of similarity between FRQ and PER are also conserved between the Drosophila PER and the mouse PER (mPER2). This suggests conserved and important regions for all 3 proteins and a common ancestor, possibly in those amoeba, such as Capsaspora, that sits at the base of the phylogenetic tree where fungi and animals diverged. Two additional examples of a possible common ancestor between Neurospora and Drosophila were found. One, the white collar (WC-1) protein of Neurospora and the Drosophila PER, shows significant similarity in its Per/Arnt/Sim (PAS) motifs to the PAS motif of an ARNT-like protein found in the amoeba, Capsaspora. Two, both of the positive elements in each system (i.e., WC-1 in Neurospora and cycle [CYC] in Drosophila), show significant similarity to this Capsaspora ARNT protein. A discussion of these findings centers on the long-time debate about the origins of the many different clock systems (i.e., independent evolution or common ancestor as well as to the question of how new genes are formed).

The bZIP transcription factor FpAda1 is essential for fungal growth and conidiation in Fusarium pseudograminearum

Fusarium pseudograminearum is an important pathogen of Fusarium crown rot and Fusarium head blight, which is able to infect wheat and barley worldwide, causing great economic losses. Transcription factors (TFs) of the basic leucine zipper (bZIP) protein family control important processes in all eukaryotes. In this study, we identified a gene, designated FpAda1, encoding a bZIP TF in F. pseudograminearum. The homolog of FpAda1 is also known to affect hyphal growth in Neurospora crassa. Deletion of FpAda1 in F. pseudograminearum resulted in defects in hyphal growth, mycelial branching and conidia formation. Pathogenicity assays showed that virulence of the Δfpada1 mutant was dramatically decreased on wheat coleoptiles and barley leaves. However, wheat coleoptile inoculation assay showed that Δfpada1 could penetrate and proliferate in wheat cells. Moreover, the FpAda1 was required for abnormal nuclear morphology in conidia and transcription of FpCdc2 and FpCdc42. Taken together, these results indicate that FpAda1 is an important transcription factor involved in growth and development in F. pseudograminearum.