Search

DNA N6-adenine methylation (6mA) has recently gained importance as an epigenetic modification in eukaryotes. Its function in lineages with high levels, such as early-diverging fungi (EDF), is of particular interest. Here, we investigated the biological significance and evolutionary implications of 6mA in EDF, which exhibit divergent evolutionary patterns in 6mA usage. The analysis of two Mucorales species displaying extreme 6mA usage reveals that species with high 6mA levels show symmetric methylation enriched in highly expressed genes. In contrast, species with low 6mA levels show mostly asymmetric 6mA. Interestingly, transcriptomic regulation throughout development and in response to environmental cues is associated with changes in the 6mA landscape. Furthermore, we identify an EDF-specific methyltransferase, likely originated from endosymbiotic bacteria, as responsible for asymmetric methylation, while an MTA-70 methylation complex performs symmetric methylation. The distinct phenotypes observed in the corresponding mutants reinforced the critical role of both types of 6mA in EDF.

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

Papiliotrema laurentii 5307AH was isolated from an aircraft polymer-coated surface. The genome size is 19,510,785 bp with a G + C content of 56%. The genome harbors genes encoding oxygenases, cutinases, lipases, and enzymes for styrene degradation, all of which could play a critical role in survival on xenobiotic surfaces.

Modification of lignin in feedstocks via genetic engineering aims to reduce biomass recalcitrance to facilitate efficient conversion processes. These improvements can be achieved by expressing exogenous enzymes that interfere with native biosynthetic pathways responsible for the production of the lignin precursors. In-planta expression of a 3-dehydroshikimate dehydratase (QsuB) in poplar trees reduced lignin content and altered their monomer composition, which enabled higher yields of sugars after cell wall polysaccharide hydrolysis. Understanding how plants respond to such genetic modifications at the transcriptional and metabolic levels is needed to facilitate further improvement and field deployment. In this work, we amassed fundamental knowledge on lignin-modified QsuB poplar using RNA-seq and metabolomics. The data clearly demonstrate that changes in gene expression and metabolite abundance can occur in a strict spatiotemporal fashion, revealing tissue-specific responses in the xylem, phloem, or periderm. In the poplar line that exhibits the strongest reduction in lignin, we found that 3% of the transcripts had altered expression levels and ~19% of the detected metabolites had differential abundance in the xylem from older stems. Changes affect predominantly the shikimate and phenylpropanoid pathways as wells as secondary cell wall metabolism, and result in significant accumulation of hydroxybenzoates derived from protocatechuate and salicylate.

Given the multitude of extracellular enzymes at their disposal, many of which are designed to degrade natures polymers (lignin, cutin, cellulose, etc.), fungi are adept at targeting synthetic polyesters with similar chemical composition. Microbial-influenced deterioration of xenobiotic polymeric surfaces is an area of interest for material scientists as these are important for the conservation of the underlying structural materials. Here, we describe the isolation and characterization of the Papiliotrema laurentii 5307AH (P. laurentii) cutinase, Plcut1. P. laurentii is basidiomycete yeast with the ability to disperse Impranil-DLN (Impranil), a colloidal polyester polyurethane, in agar plates. To test whether the fungal factor involved in this clearing was a secreted enzyme, we screened the ability of P. laurentii culture supernatants to disperse Impranil. Using size exclusion chromatography (SEC), we isolated fractions that contained Impranil-clearing activity. These fractions harbored a single ~22 kD band, which was excised and subjected to peptide sequencing. Homology searches using the peptide sequences identified, revealed that the protein Papla1 543643 (Plcut1) displays similarities to serine esterase and cutinase family of proteins. Biochemical assays using recombinant Plcut1 confirmed that this enzyme has the capability to hydrolyze Impranil, soluble esterase substrates, and apple cutin. Finally, we confirmed the presence of the Plcut1 in culture supernatants using a custom antibody that specifically recognizes this protein. The work shown here supports a major role for the Plcut1 in the fungal degradation of natural polyesters and xenobiotic polymer surfaces.IMPORTANCEFungi play a vital role in the execution of a broad range of biological processes that drive ecosystem function through production of a diverse arsenal of enzymes. However, the universal reactivity of these enzymes is a current problem for the built environment and the undesired degradation of polymeric materials in protective coatings. Here, we report the identification and characterization of a hydrolase from Papiliotrema laurentii 5307AH, an aircraft-derived fungal isolate found colonizing a biodeteriorated polymer-coated surface. We show that P. laurentii secretes a cutinase capable of hydrolyzing soluble esters as well as ester-based compounds forming solid surface coatings. These findings indicate that this fungus plays a significant role in biodeterioration through the production of a cutinase adept at degrading ester-based polymers, some of which form the backbone of protective surface coatings. The work shown here provides insights into the mechanisms employed by fungi to degrade xenobiotic polymers.

Mutant populations are crucial for functional genomics and discovering novel traits for crop breeding. Sorghum, a drought and heat-tolerant C4 species, requires a vast, large-scale, annotated, and sequenced mutant resource to enhance crop improvement through functional genomics research. Here, we report a sorghum large-scale sequenced mutant population with 9.5 million ethyl methane sulfonate (EMS)-induced mutations that covered 98% of sorghum's annotated genes using inbred line BTx623. Remarkably, a total of 610 320 mutations within the promoter and enhancer regions of 18 000 and 11 790 genes, respectively, can be leveraged for novel research of cis-regulatory elements. A comparison of the distribution of mutations in the large-scale mutant library and sorghum association panel (SAP) provides insights into the influence of selection. EMS-induced mutations appeared to be random across different regions of the genome without significant enrichment in different sections of a gene, including the 5' UTR, gene body, and 3'-UTR. In contrast, there were low variation density in the coding and UTR regions in the SAP. Based on the Ka /Ks value, the mutant library (~1) experienced little selection, unlike the SAP (0.40), which has been strongly selected through breeding. All mutation data are publicly searchable through SorbMutDB (https://www.depts.ttu.edu/igcast/sorbmutdb.php) and SorghumBase (https://sorghumbase.org/). This current large-scale sequence-indexed sorghum mutant population is a crucial resource that enriched the sorghum gene pool with novel diversity and a highly valuable tool for the Poaceae family, that will advance plant biology research and crop breeding.

Favolaschia claudopus, a wood-inhabiting basidiomycete of the Mycenaceae family, is considered an invasive species that has recently spread from Oceania to Europe. The CIRM-BRFM 2984 strain of this fungus was originally isolated from a basidiome collected from the fallen limb of a decayed oak tree in Southwest France. The genome sequence of this strain shared characteristics with other Mycenaceae species, including a large genome size and enriched content of protein-coding genes. The genome sequence provided here will facilitate further investigation on the factors that contribute to the successful global dissemination of F. claudopus.

The rice receptor kinase XA21 confers broad-spectrum resistance to Xanthomonas oryzae pv. oryzae (Xoo), the causal agent of rice bacterial blight disease. To investigate the relationship between the expression level of XA21 and resulting resistance, we generated independent HA-XA21 transgenic rice lines accumulating the XA21 immune receptor fused with an HA epitope tag. Whole-genome sequence analysis identified the T-DNA insertion sites in sixteen independent T0 events. Through quantification of the HA-XA21 protein and assessment of the resistance to Xoo strain PXO99 in six independent transgenic lines, we observed that XA21-mediated resistance is dose dependent. In contrast, based on the four agronomic traits quantified in these experiments, yield is unlikely to be affected by the expression level of HA-XA21. These findings extend our knowledge of XA21-mediated defense and contribute to the growing number of well-defined genomic landing pads in the rice genome that can be targeted for gene insertion without compromising yield.

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

BackgroundCost-effective production of biofuels from lignocellulose requires the fermentation of D-xylose. Many yeast species within and closely related to the genera Spathaspora and Scheffersomyces (both of the order Serinales) natively assimilate and ferment xylose. Other species consume xylose inefficiently, leading to extracellular accumulation of xylitol. Xylitol excretion is thought to be due to the different cofactor requirements of the first two steps of xylose metabolism. Xylose reductase (XR) generally uses NADPH to reduce xylose to xylitol, while xylitol dehydrogenase (XDH) generally uses NAD+ to oxidize xylitol to xylulose, creating an imbalanced redox pathway. This imbalance is thought to be particularly consequential in hypoxic or anoxic environments.ResultsWe screened the growth of xylose-fermenting yeast species in high and moderate aeration and identified both ethanol producers and xylitol producers. Selected species were further characterized for their XR and XDH cofactor preferences by enzyme assays and gene expression patterns by RNA-Seq. Our data revealed that xylose metabolism is more redox balanced in some species, but it is strongly affected by oxygen levels. Under high aeration, most species switched from ethanol production to xylitol accumulation, despite the availability of ample oxygen to accept electrons from NADH. This switch was followed by decreases in enzyme activity and the expression of genes related to xylose metabolism, suggesting that bottlenecks in xylose fermentation are not always due to cofactor preferences. Finally, we expressed XYL genes from multiple Scheffersomyces species in a strain of Saccharomyces cerevisiae. Recombinant S. cerevisiae expressing XYL1 from Scheffersomyces xylosifermentans, which encodes an XR without a cofactor preference, showed improved anaerobic growth on xylose as the primary carbon source compared to S. cerevisiae strain expressing XYL genes from Scheffersomyces stipitis.ConclusionCollectively, our data do not support the hypothesis that xylitol accumulation occurs primarily due to differences in cofactor preferences between xylose reductase and xylitol dehydrogenase; instead, gene expression plays a major role in response to oxygen levels. We have also identified the yeast Sc. xylosifermentans as a potential source for genes that can be engineered into S. cerevisiae to improve xylose fermentation and biofuel production.

Plant lignocellulosic biomass, i.e. secondary cell walls of plants, is a vital alternative source for bioenergy. However, the acetylation of xylan in secondary cell walls impedes the conversion of biomass to biofuels. Previous studies have shown that REDUCED WALL ACETYLATION (RWA) proteins are directly involved in the acetylation of xylan but the regulatory mechanism of RWAs is not fully understood. In this study, we demonstrate that overexpression of a Populus trichocarpa PtRWA-C gene increases the level of xylan acetylation and increases the lignin content and S/G ratio, ultimately yielding poplar woody biomass with reduced saccharification efficiency. Furthermore, through gene coexpression network and expression quantitative trait loci (eQTL) analysis, we found that PtRWA-C was regulated not only by the secondary cell wall hierarchical regulatory network but also by an AP2 family transcription factor HARDY (HRD). Specifically, HRD activates PtRWA-C expression by directly binding to the PtRWA-C promoter, which is also the cis-eQTL for PtRWA-C. Taken together, our findings provide insights into the functional roles of PtRWA-C in xylan acetylation and consequently saccharification and shed light on synthetic biology approaches to manipulate this gene and alter cell wall properties. These findings have substantial implications for genetic engineering of woody species, which could be used as a sustainable source of biofuels, valuable biochemicals, and biomaterials.

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.

Lipomyces tetrasporous is an oleaginous yeast that can utilize a variety of plant-based sugars. It accumulates lipids during growth on lignocellulosic biomass hydrolysates. We present the annotated genome sequence of L. tetrasporous NRRL Y-64009 to aid in its development as a platform organism for producing lipids and lipid-based bioproducts.

Brown-rot fungi lack many enzymes associated with complete wood degradation, such as lignin-attacking peroxidases, and have developed alternative mechanisms for rapid wood breakdown. To identify the effects of culture conditions and wood substrates on gene expression, we grew Fibroporia radiculosa in submerged cultures containing Wiley milled wood (5 days) and solid wood wafers (30 days), using aspen, pine, and spruce as a substrate. The comparative analysis revealed that wood species had a limited effect on the transcriptome:

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

The fungal genus Armillaria contains necrotrophic pathogens and some of the largest terrestrial organisms that cause tremendous losses in diverse ecosystems, yet how they evolved pathogenicity in a clade of dominantly non-pathogenic wood degraders remains elusive. Here we show that Armillaria species, in addition to gene duplications and de novo gene origins, acquired at least 1,025 genes via 124 horizontal gene transfer events, primarily from Ascomycota. Horizontal gene transfer might have affected plant biomass degrading and virulence abilities of Armillaria, and provides an explanation for their unusual, soft rot-like wood decay strategy. Combined multi-species expression data revealed extensive regulation of horizontally acquired and wood-decay related genes, putative virulence factors and two novel conserved pathogenicity-induced small secreted proteins, which induced necrosis in planta. Overall, this study details how evolution knitted together horizontally and vertically inherited genes in complex adaptive traits of plant biomass degradation and pathogenicity in important fungal pathogens.

Black yeasts are polyextremotolerant fungi that contain high amounts of melanin in their cell wall and maintain a primar yeast form. These fungi grow in xeric, nutrient depletes environments which implies that they require highly flexible metabolisms and have been suggested to contain the ability to form lichen-like mutualisms with nearby algae and bacteria. However, the exact ecological niche and interactions between these fungi and their surrounding community are not well understood. We have isolated 2 novel black yeasts from the genus Exophiala that were recovered from dryland biological soil crusts. Despite notable differences in colony and cellular morphology, both fungi appear to be members of the same species, which has been named Exophiala viscosa (i.e. E. viscosa JF 03-3 Goopy and E. viscosa JF 03-4F Slimy). A combination of whole genome sequencing, phenotypic experiments, and melanin regulation experiments have been performed on these isolates to fully characterize these fungi and help decipher their fundamental niche within the biological soil crust consortium. Our results reveal that E. viscosa is capable of utilizing a wide variety of carbon and nitrogen sources potentially derived from symbiotic microbes, can withstand many forms of abiotic stresses, and excretes melanin which can potentially provide ultraviolet resistance to the biological soil crust community. Besides the identification of a novel species within the genus Exophiala, our study also provides new insight into the regulation of melanin production in polyextremotolerant fungi.

Fungi have evolved over millions of years and their species diversity is predicted to be the second largest on the earth. Fungi have cross-kingdom interactions with many organisms that have mutually shaped their evolutionary trajectories. Zygomycete fungi hold a pivotal position in the fungal tree of life and provide important perspectives on the early evolution of fungi from aquatic to terrestrial environments. Phylogenomic analyses have found that zygomycete fungi diversified into two separate clades, the Mucoromycota which are frequently associated with plants and Zoopagomycota that are commonly animal-associated fungi. Genetic elements that contributed to the fitness and divergence of these lineages may have been shaped by the varied interactions these fungi have had with plants, animals, bacteria, and other microbes. To investigate this, we performed comparative genomic analyses of the two clades of zygomycetes in the context of Kingdom Fungi, benefiting from our generation of a new collection of zygomycete genomes, including nine produced for this study. We identified lineage-specific genomic content that may contribute to the disparate biology observed in these zygomycetes. Our findings include the discovery of undescribed diversity in CotH, a Mucormycosis pathogenicity factor, which was found in a broad set of zygomycetes. Reconciliation analysis identified multiple duplication events and an expansion of CotH copies throughout the Mucoromycotina, Mortierellomycotina, Neocallimastigomycota, and Basidiobolus lineages. A kingdom-level phylogenomic analysis also identified new evolutionary relationships within the subphyla of Mucoromycota and Zoopagomycota, including supporting the sister-clade relationship between Glomeromycotina and Mortierellomycotina and the placement of Basidiobolus as sister to other Zoopagomycota lineages.

Lignocellulose forms plant cell walls, and its three constituent polymers, cellulose, hemicellulose and lignin, represent the largest renewable organic carbon pool in the terrestrial biosphere. Insights into biological lignocellulose deconstruction inform understandings of global carbon sequestration dynamics and provide inspiration for biotechnologies seeking to address the current climate crisis by producing renewable chemicals from plant biomass. Organisms in diverse environments disassemble lignocellulose, and carbohydrate degradation processes are well defined, but biological lignin deconstruction is described only in aerobic systems. It is currently unclear whether anaerobic lignin deconstruction is impossible because of biochemical constraints or, alternatively, has not yet been measured. We applied whole cell-wall nuclear magnetic resonance, gel-permeation chromatography and transcriptome sequencing to interrogate the apparent paradox that anaerobic fungi (Neocallimastigomycetes), well-documented lignocellulose degradation specialists, are unable to modify lignin. We find that Neocallimastigomycetes anaerobically break chemical bonds in grass and hardwood lignins, and we further associate upregulated gene products with the observed lignocellulose deconstruction. These findings alter perceptions of lignin deconstruction by anaerobes and provide opportunities to advance decarbonization biotechnologies that depend on depolymerizing lignocellulose.

Deciduous woody plants like poplar (Populus spp.) have seasonal bud dormancy. It has been challenging to simultaneously delay the onset of bud dormancy in the fall and advance bud break in the spring, as bud dormancy, and bud break were thought to be controlled by different genetic factors. Here, we demonstrate that heterologous expression of the REVEILLE1 gene (named AaRVE1) from Agave (Agave americana) not only delays the onset of bud dormancy but also accelerates bud break in poplar in field trials. AaRVE1 heterologous expression increases poplar biomass yield by 166% in the greenhouse. Furthermore, we reveal that heterologous expression of AaRVE1 increases cytokinin contents, represses multiple dormancy-related genes, and up-regulates bud break-related genes, and that AaRVE1 functions as a transcriptional repressor and regulates the activity of the DORMANCY-ASSOCIATED PROTEIN 1 (DRM1) promoter. Our findings demonstrate that AaRVE1 appears to function as a regulator of bud dormancy and bud break, which has important implications for extending the growing season of deciduous trees in frost-free temperate and subtropical regions to increase crop yield.

Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes. We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia-Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase (lecsl) and γ-glutamyl transpeptidase (leggt), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula. Paralogs of both genes that are unique to Lentinula (lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes. The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.

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.

Low-cost plant substrates, such as soybean hulls, are used for various industrial applications. Filamentous fungi are important producers of Carbohydrate Active enZymes (CAZymes) required for the degradation of these plant biomass substrates. CAZyme production is tightly regulated by several transcriptional activators and repressors. One such transcriptional activator is CLR-2/ClrB/ManR, which has been identified as a regulator of cellulase and mannanase production in several fungi. However, the regulatory network governing the expression of cellulase and mannanase encoding genes has been reported to differ between fungal species. Previous studies showed that Aspergillus niger ClrB is involved in the regulation of (hemi-)cellulose degradation, although its regulon has not yet been identified. To reveal its regulon, we cultivated an A. niger ΔclrB mutant and control strain on guar gum (a galactomannan-rich substrate) and soybean hulls (containing galactomannan, xylan, xyloglucan, pectin and cellulose) to identify the genes that are regulated by ClrB. Gene expression data and growth profiling showed that ClrB is indispensable for growth on cellulose and galactomannan and highly contributes to growth on xyloglucan in this fungus. Therefore, we show that A. niger ClrB is crucial for the utilization of guar gum and the agricultural substrate, soybean hulls. Moreover, we show that mannobiose is most likely the physiological inducer of ClrB in A. niger and not cellobiose, which is considered to be the inducer of N. crassa CLR-2 and A. nidulans ClrB.

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

In the North-Central United States, lowland ecotype switchgrass can increase yield by up to 50% compared with locally adapted but early flowering cultivars. However, lowland ecotypes are not winter tolerant. The mechanism for winter damage is unknown but previously has been associated with late flowering time. This study investigated heading date (measured for two years) and winter survivorship (measured for three years) in a multi-generation population generated from two winter-hardy lowland individuals and diverse southern lowland populations. Sequencing data (311,776 markers) from 1,306 individuals were used to evaluate genome-wide trait prediction through cross-validation and progeny prediction (n = 52). Genetic variance for heading date and winter survivorship was additive with high narrow-sense heritability (0.64 and 0.71, respectively) and reliability (0.68 and 0.76, respectively). The initial negative correlation between winter survivorship and heading date degraded across generations (F1r = -0.43, pseudo-F2r = -0.28, pseudo-F2 progeny r = -0.15). Within-family predictive ability was moderately high for heading date and winter survivorship (0.53 and 0.52, respectively). A multi-trait model did not improve predictive ability for either trait. Progeny predictive ability was 0.71 for winter survivorship and 0.53 for heading date. These results suggest that lowland ecotype populations can obtain sufficient survival rates in the northern United States with two or three cycles of effective selection. Despite accurate genomic prediction, naturally occurring winter mortality successfully isolated winter tolerant genotypes and appears to be an efficient method to develop high-yielding, cold-tolerant switchgrass cultivars.

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

Pleurotus ostreatus is a white-rot fungus that can degrade lignin in a preferential manner using a variety of extracellular enzymes, including manganese and versatile peroxidases (encoded by the vp1-3 and mnp1-6 genes, respectively). This fungus also secretes a family of structurally related small secreted proteins (SSPs) encoded by the ssp1-6 genes. Using RNA sequencing (RNA-seq), we determined that ssp4 and ssp6 are the predominant members of this gene family that were expressed by P. ostreatus during the first three weeks of growth on wheat straw. Downregulation of ssp4 in a strain harboring an ssp RNAi construct (KDssp1) was then confirmed, which, along with an increase in ssp6 transcript levels, coincided with reduced lignin degradation and the downregulation of vp2 and mnp1. In contrast, we observed an increase in the expression of genes related to pectin and side-chain hemicellulose degradation, which was accompanied by an increase in extracellular pectin-degrading capacity. Genome-wide comparisons between the KDssp1 and the wild-type strains demonstrated that ssp silencing conferred accumulated changes in gene expression at the advanced cultivation stages in an adaptive rather than an inductive mode of transcriptional response. Based on co-expression networking, crucial gene modules were identified and linked to the ssp knockdown genotype at different cultivation times. Based on these data, as well as previous studies, we propose that P. ostreatus SSPs have potential roles in modulating the lignocellulolytic and pectinolytic systems, as well as a variety of fundamental biological processes related to fungal growth and development.

Macrocystis pyrifera (giant kelp), is a brown macroalga of great ecological importance as a primary producer and structure-forming foundational species that provides habitat for hundreds of species. It has many commercial uses (e.g. source of alginate, fertilizer, cosmetics, feedstock). One of the limitations to exploiting giant kelp's economic potential and assisting in giant kelp conservation efforts is a lack of genomic tools like a high quality, contiguous reference genome with accurate gene annotations. Reference genomes attempt to capture the complete genomic sequence of an individual or species, and importantly provide a universal structure for comparison across a multitude of genetic experiments, both within and between species. We assembled the giant kelp genome of a haploid female gametophyte de novo using PacBio reads, then ordered contigs into chromosome level scaffolds using Hi-C. We found the giant kelp genome to be 537 MB, with a total of 35 scaffolds and 188 contigs. The assembly N50 is 13,669,674 with GC content of 50.37%. We assessed the genome completeness using BUSCO, and found giant kelp contained 94% of the BUSCO genes from the stramenopile clade. Annotation of the giant kelp genome revealed 25,919 genes. Additionally, we present genetic variation data based on 48 diploid giant kelp sporophytes from three different Southern California populations that confirms the population structure found in other studies of these populations. This work resulted in a high-quality giant kelp genome that greatly increases the genetic knowledge of this ecologically and economically vital species.

Plant biomass is one of the most abundant renewable carbon sources, which holds great potential for replacing current fossil-based production of fuels and chemicals. In nature, fungi can efficiently degrade plant polysaccharides by secreting a broad range of carbohydrate-active enzymes (CAZymes), such as cellulases, hemicellulases, and pectinases. Due to the crucial role of plant biomass-degrading (PBD) CAZymes in fungal growth and related biotechnology applications, investigation of their genomic diversity and transcriptional dynamics has attracted increasing attention. In this project, we systematically compared the genome content of PBD CAZymes in six taxonomically distant species, Aspergillus niger, Aspergillus nidulans, Penicillium subrubescens, Trichoderma reesei, Phanerochaete chrysosporium, and Dichomitus squalens, as well as their transcriptome profiles during growth on nine monosaccharides. Considerable genomic variation and remarkable transcriptomic diversity of CAZymes were identified, implying the preferred carbon source of these fungi and their different methods of transcription regulation. In addition, the specific carbon utilization ability inferred from genomics and transcriptomics was compared with fungal growth profiles on corresponding sugars, to improve our understanding of the conversion process. This study enhances our understanding of genomic and transcriptomic diversity of fungal plant polysaccharide-degrading enzymes and provides new insights into designing enzyme mixtures and metabolic engineering of fungi for related industrial applications.

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Filamentous fungi degrade complex plant material to its monomeric building blocks, which have many biotechnological applications. Transcription factors play a key role in plant biomass degradation, but little is known about their interactions in the regulation of polysaccharide degradation. Here, we deepened the knowledge about the storage polysaccharide regulators AmyR and InuR in Aspergillus niger. AmyR controls starch degradation, while InuR is involved in sucrose and inulin utilization. In our study, the phenotypes of A. niger parental, ΔamyR, ΔinuR and ΔamyRΔinuR strains were assessed in both solid and liquid media containing sucrose or inulin as carbon source to evaluate the roles of AmyR and InuR and the effect of culture conditions on their functions. In correlation with previous studies, our data showed that AmyR has a minor contribution to sucrose and inulin utilization when InuR is active. In contrast, growth profiles and transcriptomic data showed that the deletion of amyR in the ΔinuR background strain resulted in more pronounced growth reduction on both substrates, mainly evidenced by data originating from solid cultures. Overall, our results show that submerged cultures do not always reflect the role of transcription factors in the natural growth condition, which is better represented on solid substrates. Importance: The type of growth has critical implications in enzyme production by filamentous fungi, a process that is controlled by transcription factors. Submerged cultures are the preferred setups in laboratory and industry and are often used for studying the physiology of fungi. In this study, we showed that the genetic response of A. niger to starch and inulin was highly affected by the culture condition, since the transcriptomic response obtained in a liquid environment did not fully match the behavior of the fungus in a solid environment. These results have direct implications in enzyme production and would help industry choose the best approaches to produce specific CAZymes for industrial purposes.

ObjectiveAdaptive laboratory evolution (ALE) is an effective approach to study the evolution behavior of bacterial cultures and to select for strains with desired metabolic features. In this study, we explored the possibility of evolving Thermotoga sp. strain RQ7 for cellulose-degrading abilities.ResultsWild type RQ7 strain was subject to a series of transfers over six and half years with cellulose filter paper as the main and eventually the sole carbon source. Each transfer was accompanied with the addition of 50 μg of Caldicellulosiruptor saccharolyticus DSM 8903 genomic DNA. A total of 331 transfers were completed. No cellulose degradation was observed with the RQ7 cultures. Thirty three (33) isolates from six time points were sampled and sequenced. Nineteen (19) of the 33 isolates were unique, and the rest were duplicated clones. None of the isolates acquired C. saccharolyticus DNA, but all accumulated small-scale mutations throughout their genomes. Sequence analyses revealed 35 mutations that were preserved throughout the generations and another 15 mutations emerged near the end of the study. Many of the affected genes participate in phosphate metabolism, substrate transport, stress response, sensory transduction, and gene regulation.

Wildfires drastically impact the soil environment, altering the soil organic matter, forming pyrolyzed compounds, and markedly reducing the diversity of microorganisms. Pyrophilous fungi, especially the species from the orders Pezizales and Agaricales, are fire-responsive fungal colonizers of post-fire soil that have historically been found fruiting on burned soil and thus may encode mechanisms of processing these compounds in their genomes. Pyrophilous fungi are diverse. In this work, we explored this diversity and sequenced six new genomes of pyrophilous Pezizales fungi isolated after the 2013 Rim Fire near Yosemite Park in California, USA: Pyronema domesticum, Pyronema omphalodes, Tricharina praecox, Geopyxis carbonaria, Morchella snyderi, and Peziza echinospora. A comparative genomics analysis revealed the enrichment of gene families involved in responses to stress and the degradation of pyrolyzed organic matter. In addition, we found that both protein sequence lengths and G + C content in the third base of codons (GC3) in pyrophilous fungi fall between those in mesophilic/nonpyrophilous and thermophilic fungi. A comparative transcriptome analysis of P. domesticum under two conditions - growing on charcoal, and during sexual development - identified modules of genes that are co-expressed in the charcoal and light-induced sexual development conditions. In addition, environmental sensors such as transcription factors STE12, LreA, LreB, VosA, and EsdC were upregulated in the charcoal condition. Taken together, these results highlight genomic adaptations of pyrophilous fungi and indicate a potential connection between charcoal tolerance and fruiting body formation in P. domesticum.

Crassulacean acid metabolism (CAM) plants exhibit elevated drought and heat tolerance compared to C3 and C4 plants through an inverted pattern of day/night stomatal closure and opening for CO2 assimilation. However, the molecular responses to water-deficit conditions remain unclear in obligate CAM species. In this study, we presented genome-wide transcription sequencing analysis using leaf samples of an obligate CAM species Kalanchoë fedtschenkoi under moderate and severe drought treatments at two-time points of dawn (2-h before the start of light period) and dusk (2-h before the dark period). Differentially expressed genes were identified in response to environmental drought stress and a whole genome wide co-expression network was created as well. We found that the expression of CAM-related genes was not regulated by drought stimuli in K. fedtschenkoi. Our comparative analysis revealed that CAM species (K. fedtschenkoi) and C3 species (Arabidopsis thaliana, Populus deltoides 'WV94') share some common transcriptional changes in genes involved in multiple biological processes in response to drought stress, including ABA signaling and biosynthesis of secondary metabolites.

The brown rot fungus Fomitopsis pinicola efficiently depolymerizes wood cellulose via the combined activities of oxidative and hydrolytic enzymes. Mass spectrometric analyses of culture filtrates identified specific proteins, many of which were differentially regulated in response to substrate composition.

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.

Fungi survive in diverse ecological niches by secreting proteins and other molecules into the environment to acquire food and interact with various biotic and abiotic stressors. Fungal secretome content is, therefore, believed to be tightly linked to fungal ecologies. We sampled 132 genomes from the early-diverging terrestrial fungal lineage zygomycetes (Mucoromycota and Zoopagomycota) and characterized their secretome composition. Our analyses revealed that phylogeny played an important role in shaping the secretome composition of zygomycete fungi with trophic mode contributing a smaller amount. Reconstruction of the evolution of secreted digestive enzymes revealed lineage-specific expansions, indicating that Mucoromycota and Zoopagomycota followed different trajectories early in their evolutionary history. We identified the presence of multiple pathogenicity-related proteins in the lineages known as saprotrophs, suggesting that either the ecologies of these fungi are incompletely known, and/or that these pathogenicity-related proteins have important functions associated with saprotrophic ecologies, both of which invite further investigation.

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

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.

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.

The halotolerant and osmotolerant yeast Zygosaccharomyces rouxii can produce multiple volatile compounds and has the ability to grow on lignocellulosic hydrolysates. We report the annotated genome sequence of Z. rouxii NRRL Y-64007 to support its development as a platform organism for biofuel and bioproduct production.

Members of the fungal genus Phyllosticta can colonize a variety of plant hosts, including several Citrus species such as Citrus sinensis (orange), Citrus limon (lemon), and Citrus maxima (pomelo). Some Phyllosticta species have the capacity to cause disease, such as Citrus Black Spot, while others have only been observed as endophytes. Thus far, genomic differences underlying lifestyle adaptations of Phyllosticta species have not yet been studied. Furthermore, the lifestyle of Phyllosticta citrichinaensis is ambiguous, as it has been described as a weak pathogen but Koch's postulates may not have been established and the presence of this species was never reported to cause any crop or economic losses. Here, we examined the genomic differences between pathogenic and endophytic Phyllosticta spp. colonizing Citrus and specifically aimed to elucidate the lifestyle of Phyllosticta citrichinaensis. We found several genomic differences between species of different lifestyles, including groups of genes that were only present in pathogens or endophytes. We also observed that species, based on their carbohydrate active enzymes, group independent of their phylogenetic association, and this clustering correlated with trophy prediction. Phyllosticta citrichinaensis shows an intermediate lifestyle, sharing genomic and phenotypic attributes of both pathogens and endophytes. We thus present the first genomic comparison of multiple citrus-colonizing pathogens and endophytes of the genus Phyllosticta, and therefore provide the basis for further comparative studies into the lifestyle adaptations within this genus.

IntroductionMushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi.ObjectivesThis study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution.MethodsWe analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution.ResultsWe detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation.ConclusionAdaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

Efficient utilization of agro-industrial waste, such as sugar beet pulp, is crucial for the bio-based economy. The fungus Aspergillus niger possesses a wide array of enzymes that degrade complex plant biomass substrates, and several regulators have been reported to play a role in their production. The role of the regulators GaaR, AraR, and RhaR in sugar beet pectin degradation has previously been reported. However, genetic regulation of the degradation of sugar beet pulp has not been assessed in detail. In this study, we generated a set of single and combinatorial deletion mutants targeting the pectinolytic regulators GaaR, AraR, RhaR, and GalX as well as the (hemi-)cellulolytic regulators XlnR and ClrB to address their relative contribution to the utilization of sugar beet pulp. We show that A. niger has a flexible regulatory network, adapting to the utilization of (hemi-)cellulose at early timepoints when pectin degradation is impaired.

The current impetus towards a sustainable bio-based economy has accelerated research to better understand the mechanisms through which filamentous fungi convert plant biomass, a valuable feedstock for biotechnological applications. Several transcription factors have been reported to control the polysaccharide degradation and metabolism of the resulting sugars in fungi. However, little is known about their individual contributions, interactions and crosstalk. D-galactose is a hexose sugar present mainly in hemicellulose and pectin in plant biomass. Here, we study D-galactose conversion by Aspergillus niger and describe the involvement of the arabinanolytic and xylanolytic activators AraR and XlnR, in addition to the D-galactose-responsive regulator GalX. Our results deepen the understanding of the complexity of the filamentous fungal regulatory network for plant biomass degradation and sugar catabolism, and facilitate the generation of more efficient plant biomass-degrading strains for biotechnological applications.