In Neurospora crassa, the transcription factor COL-26 functions as a regulator of glucose signaling and metabolism. Its loss leads to resistance to carbon catabolite repression. Here, we report that COL-26 is necessary for the expression of amylolytic genes in N. crassa and is required for the utilization of maltose and starch. Additionally, the Δcol-26 mutant shows growth defects on preferred carbon sources, such as glucose, an effect that was alleviated if glutamine replaced ammonium as the primary nitrogen source. This rescue did not occur when maltose was used as a sole carbon source. Transcriptome and metabolic analyses of the Δcol-26 mutant relative to its wild type parental strain revealed that amino acid and nitrogen metabolism, the TCA cycle and GABA shunt were adversely affected. Phylogenetic analysis showed a single col-26 homolog in Sordariales, Ophilostomatales, and the Magnaporthales, but an expanded number of col-26 homologs in other filamentous fungal species. Deletion of the closest homolog of col-26 in Trichoderma reesei, bglR, resulted in a mutant with similar preferred carbon source growth deficiency, and which was alleviated if glutamine was the sole nitrogen source, suggesting conservation of COL-26 and BglR function. Our finding provides novel insight into the role of COL-26 for utilization of starch and in integrating carbon and nitrogen metabolism for balanced metabolic activities for optimal carbon and nitrogen distribution.

The herbivore digestive tract is home to a complex community of anaerobic microbes that work together to break down lignocellulose. These microbiota are an untapped resource of strains, pathways and enzymes that could be applied to convert plant waste into sugar substrates for green biotechnology. We carried out more than 400 parallel enrichment experiments from goat faeces to determine how substrate and antibiotic selection influence membership, activity, stability and chemical productivity of herbivore gut communities. We assembled 719 high-quality metagenome-assembled genomes (MAGs) that are unique at the species level. More than 90% of these MAGs are from previously unidentified herbivore gut microorganisms. Microbial consortia dominated by anaerobic fungi outperformed bacterially dominated consortia in terms of both methane production and extent of cellulose degradation, which indicates that fungi have an important role in methane release. Metabolic pathway reconstructions from MAGs of 737 bacteria, archaea and fungi suggest that cross-domain partnerships between fungi and methanogens enabled production of acetate, formate and methane, whereas bacterially dominated consortia mainly produced short-chain fatty acids, including propionate and butyrate. Analyses of carbohydrate-active enzyme domains present in each anaerobic consortium suggest that anaerobic bacteria and fungi employ mostly complementary hydrolytic strategies. The division of labour among herbivore anaerobes to degrade plant biomass could be harnessed for industrial bioprocessing.

Centromeres are rapidly evolving across eukaryotes, despite performing a conserved function to ensure high-fidelity chromosome segregation. CENP-A chromatin is a hallmark of a functional centromere in most organisms. Due to its critical role in kinetochore architecture, the loss of CENP-A is tolerated in only a few organisms, many of which possess holocentric chromosomes. Here, we characterize the consequence of the loss of CENP-A in the fungal kingdom. Mucor circinelloides, an opportunistic human pathogen, lacks CENP-A along with the evolutionarily conserved CENP-C but assembles a monocentric chromosome with a localized kinetochore complex throughout the cell cycle. Mis12 and Dsn1, two conserved kinetochore proteins, were found to co-localize to a short region, one in each of nine large scaffolds, composed of an ∼200-bp AT-rich sequence followed by a centromere-specific conserved motif that echoes the structure of budding yeast point centromeres. Resembling fungal regional centromeres, these core centromere regions are embedded in large genomic expanses devoid of genes yet marked by Grem-LINE1s, a novel retrotransposable element silenced by the Dicer-dependent RNAi pathway. Our results suggest that these hybrid features of point and regional centromeres arose from the absence of CENP-A, thus defining novel mosaic centromeres in this early-diverging fungus.

Sphaerosporella brunnea is a pioneer ectomycorrhizal fungus with facultative saprophytic capacities. Here, we sequenced the genome of S. brunnea strain Sb_GMNB300, which is estimated at 51.6 Mb in size with 872 assembled contigs accounting for 12,597 predicted coding genes. This genome will be useful for comparative studies of Pezizales ectomycorrhizal symbioses.

Anaerobic fungi (Neocallimastigomycota) are emerging non-model hosts for biotechnology due to their wealth of biomass-degrading enzymes, yet tools to engineer these fungi have not yet been established. Here, we show that the anaerobic gut fungi have the most GC depleted genomes among 443 sequenced organisms in the fungal kingdom, which has ramifications for heterologous expression of genes as well as for emerging CRISPR-based genome engineering approaches. Comparative genomic analyses suggest that anaerobic fungi may contain cellular machinery to aid in sexual reproduction, yet a complete mating pathway was not identified. Predicted proteomes of the anaerobic fungi also contain an unusually large fraction of proteins with homopolymeric amino acid runs consisting of five or more identical consecutive amino acids. In particular, threonine runs are especially enriched in anaerobic fungal carbohydrate active enzymes (CAZymes) and this, together with a high abundance of predicted N-glycosylation motifs, suggests that gut fungal CAZymes are heavily glycosylated, which may impact heterologous production of these biotechnologically useful enzymes. Finally, we present a codon optimization strategy to aid in the development of genetic engineering tools tailored to these early-branching anaerobic fungi.

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or "specialized" metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

Filamentous fungi, such as Neurospora crassa, are very efficient in deconstructing plant biomass by the secretion of an arsenal of plant cell wall-degrading enzymes, by remodeling metabolism to accommodate production of secreted enzymes, and by enabling transport and intracellular utilization of plant biomass components. Although a number of enzymes and transcriptional regulators involved in plant biomass utilization have been identified, how filamentous fungi sense and integrate nutritional information encoded in the plant cell wall into a regulatory hierarchy for optimal utilization of complex carbon sources is not understood. Here, we performed transcriptional profiling of N. crassa on 40 different carbon sources, including plant biomass, to provide data on how fungi sense simple to complex carbohydrates. From these data, we identified regulatory factors in N. crassa and characterized one (PDR-2) associated with pectin utilization and one with pectin/hemicellulose utilization (ARA-1). Using in vitro DNA affinity purification sequencing (DAP-seq), we identified direct targets of transcription factors involved in regulating genes encoding plant cell wall-degrading enzymes. In particular, our data clarified the role of the transcription factor VIB-1 in the regulation of genes encoding plant cell wall-degrading enzymes and nutrient scavenging and revealed a major role of the carbon catabolite repressor CRE-1 in regulating the expression of major facilitator transporter genes. These data contribute to a more complete understanding of cross talk between transcription factors and their target genes, which are involved in regulating nutrient sensing and plant biomass utilization on a global level.

Environmental DNA surveys reveal that most fungal diversity represents uncultured species. We sequenced the genomes of eight uncultured species across the fungal tree of life using a new single-cell genomics pipeline. We show that, despite a large variation in genome and gene space recovery from each single amplified genome (SAG), ≥90% can be recovered by combining multiple SAGs. SAGs provide robust placement for early-diverging lineages and infer a diploid ancestor of fungi. Early-diverging fungi share metabolic deficiencies and show unique gene expansions correlated with parasitism and unculturability. Single-cell genomics holds great promise in exploring fungal diversity, life cycles and metabolic potential.

Some of the most unique and compelling survival strategies in the natural world are fixed in isolated species1. To date, molecular insight into these ancient adaptations has been limited, as classic experimental genetics has focused on interfertile individuals in populations2. Here we use a new mapping approach, which screens mutants in a sterile interspecific hybrid, to identify eight housekeeping genes that underlie the growth advantage of Saccharomyces cerevisiae over its distant relative Saccharomyces paradoxus at high temperature. Pro-thermotolerance alleles at these mapped loci were required for the adaptive trait in S. cerevisiae and sufficient for its partial reconstruction in S. paradoxus. The emerging picture is one in which S. cerevisiae improved the heat resistance of multiple components of the fundamental growth machinery in response to selective pressure. Our study lays the groundwork for the mapping of genotype to phenotype in clades of sister species across Eukarya.

Cytochrome P450s form an important group of enzymes involved in xenobiotics degradation and metabolism, both primary and secondary. These enzymes are also useful in industry as biotechnological tools for bioconversion and a few are reported to be involved in pathogenicity. Trichoderma spp. are widely used in industry and agriculture and are known for their biosynthetic potential of a large number of secondary metabolites. For realising the full biosynthetic potential of an organism, it is important to do a genome-wide annotation and cataloguing of these enzymes.

Lignocellulose biomass is known as a recalcitrant material towards enzymatic hydrolysis, increasing the process cost in biorefinery. In nature, filamentous fungi naturally degrade lignocellulose, using an arsenal of hydrolytic and oxidative enzymes. Assessment of enzyme hydrolysis efficiency generally relies on the yield of glucose for a given biomass. To better understand the markers governing recalcitrance to enzymatic degradation, there is a need to enlarge the set of parameters followed during deconstruction.

Next to d-glucose, the pentoses l-arabinose and d-xylose are the main monosaccharide components of plant cell wall polysaccharides and are therefore of major importance in biotechnological applications that use plant biomass as a substrate. Pentose catabolism is one of the best-studied pathways of primary metabolism of Aspergillus niger, and an initial outline of this pathway with individual enzymes covering each step of the pathway has been previously established. However, although growth on l-arabinose and/or d-xylose of most pentose catabolic pathway (PCP) single deletion mutants of A. niger has been shown to be negatively affected, it was not abolished, suggesting the involvement of additional enzymes. Detailed analysis of the single deletion mutants of the known A. niger PCP genes led to the identification of additional genes involved in the pathway. These results reveal a high level of complexity and redundancy in this pathway, emphasizing the need for a comprehensive understanding of metabolic pathways before entering metabolic engineering of such pathways for the generation of more efficient fungal cell factories.

Microalgae efficiently convert sunlight into lipids and carbohydrates, offering bio-based alternatives for energy and chemical production. Improving algal productivity and robustness against abiotic stress requires a systems level characterization enabled by functional genomics. Here, we characterize a halotolerant microalga Scenedesmus sp. NREL 46B-D3 demonstrating peak growth near 25 °C that reaches 30 g/m2/day and the highest biomass accumulation capacity post cell division reported to date for a halotolerant strain. Functional genomics analysis revealed that genes involved in lipid production, ion channels and antiporters are expanded and expressed. Exposure to temperature stress shifts fatty acid metabolism and increases amino acids synthesis. Co-expression analysis shows that many fatty acid biosynthesis genes are overexpressed with specific transcription factors under cold stress. These and other genes involved in the metabolic and regulatory response to temperature stress can be further explored for strain improvement.

The Tremellomycetes are a species-rich group within the basidiomycete fungi; however, most analyses of this group to date have focused on pathogenic Cryptococcus species within the order Tremellales. Recent genome-assisted studies of other Tremellomycetes have identified interesting features with respect to biotechnological applications as well as the evolution of genes involved in mating and sexual development. Here, we report genome sequences of two strains of Filobasidium floriforme, a species from the order Filobasidiales, which branches basally to the Tremellales, Trichosporonales, and Holtermanniales. The assembled genomes of strains CBS6241 and CBS6242 are 27.4 Mb and 26.4 Mb in size, respectively, with 8314 and 7695 predicted protein-coding genes. Overall sequence identity at nucleic acid level between the strains is 97%. Among the predicted genes are pheromone precursor and pheromone receptor genes as well as two genes encoding homedomain (HD) transcription factors, which are predicted to be part of the mating type (MAT) locus. Sequence analysis indicates that CBS6241 and CBS6242 carry different alleles for both the pheromone/receptor genes as well as the HD transcription factors. Orthology inference identified 1482 orthogroups exclusively found in F. floriforme, some of which were involved in carbohydrate transport and metabolism. Subsequent CAZyme repertoire characterization identified 267 and 247 enzymes for CBS6241 and CBS6242, respectively, the second highest number of CAZymes among the analyzed Tremellomycete species. In addition, F. floriforme contains five CAZymes absent in other species and several plant-cell-wall degrading CAZymes with the highest copy number in Tremellomycota, indicating the biotechnological potential of this species.

The mutualistic ectomycorrhizal (ECM) fungal genus Pisolithus comprises 19 species defined to date which colonize the roots of >50 hosts worldwide suggesting that substantial genomic and functional evolution occurred during speciation. To better understand this intra-genus variation, we undertook a comparative multi-omic study of nine Pisolithus species sampled from North America, South America, Asia, and Australasia. We found that there was a small core set of genes common to all species (13%), and that these genes were more likely to be significantly regulated during symbiosis with a host than accessory or species-specific genes. Thus, the genetic "toolbox" foundational to the symbiotic lifestyle in this genus is small. Transposable elements were located significantly closer to gene classes including effector-like small secreted proteins (SSPs). Poorly conserved SSPs were more likely to be induced by symbiosis, suggesting that they may be a class of protein that tune host specificity. The Pisolithus gene repertoire is characterized by divergent CAZyme profiles when compared with other fungi, both symbiotic and saprotrophic. This was driven by differences in enzymes associated with symbiotic sugar processing, although metabolomic analysis suggest that neither copy number nor expression of these genes is sufficient to predict sugar capture from a host plant or its metabolism in fungal hyphae. Our results demonstrate that intra-genus genomic and functional diversity within ECM fungi is greater than previously thought, underlining the importance of continued comparative studies within the fungal tree of life to refine our focus on pathways and evolutionary processes foundational to this symbiotic lifestyle.

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host-pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype-phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

Filamentous fungi produce a wide variety of enzymes in order to efficiently degrade plant cell wall polysaccharides. The production of these enzymes is controlled by transcriptional regulators, which also control the catabolic pathways that convert the released monosaccharides. Two transcriptional regulators, GalX and GalR, control d-galactose utilization in the model filamentous fungus Aspergillus nidulans, while the arabinanolytic regulator AraR regulates l-arabinose catabolism. d-Galactose and l-arabinose are commonly found together in polysaccharides, such as arabinogalactan, xylan and rhamnogalacturonan I. Therefore, the catabolic pathways that convert d-galactose and l-arabinose are often also likely to be active simultaneously. In this study, we investigated the interaction between GalX, GalR and AraR in d-galactose and l-arabinose catabolism. For this, we generated single, double and triple mutants of the three regulators, and analysed their growth and enzyme and gene expression profiles. Our results clearly demonstrated that GalX, GalR and AraR co-regulate d-galactose catabolism in A. nidulans. GalX has a prominent role on the regulation of genes of d-galactose oxido-reductive pathway, while AraR can compensate for the absence of GalR and/or GalX.

Wood-decaying fungi of the class Agaricomycetes (phylum Basidiomycota) are saprotrophs that break down lignocellulose and play an important role in nutrient recycling. They secrete a wide range of extracellular plant cell wall degrading enzymes that break down cellulose, hemicellulose, and lignin, the main building blocks of plant biomass. Although the production of these enzymes is regulated mainly at the transcriptional level, no activating regulators have been identified in any wood-decaying fungus in the class Agaricomycetes. We studied the regulation of cellulase expression in the wood-decaying fungus Schizophyllum commune. Comparative genomics and transcriptomics on two wild isolates revealed a Zn2Cys6-type transcription factor gene (roc1) that was highly upregulated during growth on cellulose, compared to glucose. It is only conserved in the class Agaricomycetes. A roc1 knockout strain showed an inability to grow on medium with cellulose as sole carbon source, and growth on cellobiose and xylan (other components of wood) was inhibited. Growth on non-wood-related carbon sources was not inhibited. Cellulase gene expression and enzyme activity were reduced in the Δroc1 strain. ChIP-Seq identified 1474 binding sites of the Roc1 transcription factor. Promoters of genes involved in lignocellulose degradation were enriched with these binding sites, especially those of LPMO (lytic polysaccharide monooxygenase) CAZymes, indicating that Roc1 directly regulates these genes. A conserved motif was identified as the binding site of Roc1, which was confirmed by a functional promoter analysis. Together, Roc1 is a key regulator of cellulose degradation and the first identified in wood-decaying fungi in the phylum Basidiomycota. IMPORTANCE Wood-degrading fungi in the phylum Basidiomycota play a crucial role in nutrient recycling by breaking down all components of wood. Fungi have evolved transcriptional networks that regulate expression of wood-degrading enzymes, allowing them to prioritize one nutrient source over another. However, to date all these transcription factors have been identified in the phylum Ascomycota, which is only distantly related to the phylum Basidiomycota. Here, we identified the transcription factor Roc1 as a key regulator of cellulose degradation in the mushroom-forming and wood-degrading fungus Schizophyllum commune. Roc1 is highly conserved in the phylum Basidiomycota. Using comparative genomics, transcriptomics, ChIP-Seq and promoter analysis we have identified direct targets of Roc1, as well as other aspects of the transcriptional response to cellulose.

The genome of entomopathogenic fungus Tolypocladium inflatum Gams encodes 43 putative biosynthetic gene clusters for specialized metabolites, although genotype-phenotype linkages have been reported only for the cyclosporins and fumonisins. T. inflatum was cultured in defined minimal media, supplemented with or without one of nine different amino acids. Acquisition of LC-MS/MS data for molecular networking and manual analysis facilitated annotation of putative known and unknown metabolites. These data led us to target a family of peptaibols and guided the isolation and purification of tolypocladamide H (1), which showed modest antibacterial activity and toxicity to mammalian cells at micromolar concentrations. HRMS/MS, NMR, and advanced Marfey's analysis were used to assign the structure of 1 as a peptaibol containing 4-[(E)-2-butenyl]-4-methyl-l-threonine (Bmt), a hallmark structural motif of the cyclosporins. LC-MS detection of homologous tolypocladamide metabolites and phylogenomic analyses of peptaibol biosynthetic genes in other cultured Tolypocladium species allowed assignment of a putative tolypocladamide nonribosomal peptide synthetase gene.

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.