Search

10 p.-6 fig., 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., This research was supported by the NSF awards DEB-1456588 to D.H. and DEB-1442113 to A.R., the Department of Energy award 0000238431 to D.H., the GENOBIOREF (BIO2017-86559-R) project of the Spanish Ministry of Science and Innovation (cofinanced by FEDER funds) to F.J.R.-D., S.C., and A.T.M., grant PID2019-103901GB-I00 of the Spanish Ministry of Science and Innovation to P.F., and the NIH award R56AI146096 to A.R. J.L.S. and A.R. were supported by the HHMI and the Burroughs Wellcome Fund. The work (proposal: 10.46936/10.25585/60001048) conducted by the US Department of Energy JGI (https://ror.org/04xm1d337), a DOE Office of Science User Facility, was supported by the Office of Science of the US Department of Energy operated under contract no. DE-AC02-05CH11231.

Microbial biosynthetic gene clusters (BGCs) encoding secondary metabolites are thought to impact a plethora of biologically mediated environmental processes, yet their discovery and functional characterization in natural microbiomes remains challenging. Here we describe deep long-read sequencing and assembly of metagenomes from biological soil crusts, a group of soil communities that are rich in BGCs. Taking advantage of the unusually long assemblies produced by this approach, we recovered nearly 3,000 BGCs for analysis, including 712 full-length BGCs. Functional exploration through metatranscriptome analysis of a 3-day wetting experiment uncovered phylum-specific BGC expression upon activation from dormancy, elucidating distinct roles and complex phylogenetic and temporal dynamics in wetting processes. For example, a pronounced increase in BGC transcription occurs at night primarily in cyanobacteria, implicating BGCs in nutrient scavenging roles and niche competition. Taken together, our results demonstrate that long-read metagenomic sequencing combined with metatranscriptomic analysis provides a direct view into the functional dynamics of BGCs in environmental processes and suggests a central role of secondary metabolites in maintaining phylogenetically conserved niches within biocrusts., Marc Van Goethem et al. combine short- and long-read sequencing from biocrust samples to report nearly 3,000 biosynthetic gene clusters (BGCs) encoding microbial secondary metabolites. Their results demonstrate the advantage of integrated metatranscriptomic and long-read metagenomic sequencing in analyzing BGCs and provides insight into the role of secondary metabolites in maintaining phylogenetic niches in biocrusts.

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.

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

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

A number of crop wild relatives can tolerate extreme stressed to a degree outside the range observed in their domesticated relatives. However, it is unclear whether or how the molecular mechanisms employed by these species can be translated to domesticated crops. PaspalumPaspalum vaginatumis a self-incompatible and multiply stress-tolerant wild relative of maize and sorghum. Here we describe the sequencing and pseudomolecule level assembly of a vegetatively propagated accession ofP. vaginatum. Phylogenetic analysis based on 6,151 single-copy syntenic orthologous conserved in 6 related grass species placed paspalum as an outgroup of the maize-sorghum clade demonstrating paspalum as their closest sequenced wild relative. In parallel metabolic experiments, paspalum, but neither maize nor sorghum, exhibited significant increases in trehalose when grown under nutrient-deficit conditions. Inducing trehalose accumulation in maize, imitating the metabolic phenotype of paspalum, resulting in autophagy dependent increases in biomass accumulation.

Sex dimorphism and gene expression were studied in developing catkins in 159 F2 individuals from the bioenergy crop Salix purpurea, and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination in Salix purpurea is likely different from that observed in the related genus Populus. This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

Sex dimorphism and gene expression were studied in developing catkins in 159 F2individuals from the bioenergy cropSalix purpurea, and potential mechanisms and pathways for regulating sex development were explored. Differential expression, eQTL, bisulfite sequencing, and network analysis were used to characterize sex dimorphism, detect candidate master regulator genes, and identify pathways through which the sex determination region (SDR) may mediate sex dimorphism. Eleven genes are presented as candidates for master regulators of sex, supported by gene expression and network analyses. These include genes putatively involved in hormone signaling, epigenetic modification, and regulation of transcription. eQTL analysis revealed a suite of transcription factors and genes involved in secondary metabolism and floral development that were predicted to be under direct control of the sex determination region. Furthermore, data from bisulfite sequencing and small RNA sequencing revealed strong differences in expression between males and females that would implicate both of these processes in sex dimorphism pathways. These data indicate that the mechanism of sex determination inSalix purpureais likely different from that observed in the related genusPopulus. This further demonstrates the dynamic nature of SDRs in plants, which involves a multitude of mechanisms of sex determination and a high rate of turnover.

The zoosporic obligate endoparasites, Olpidium, hold a pivotal position to the reconstruction of the flagellum loss in fungi, one of the key morphological transitions associated with the colonization of land by the early fungi. We generated genome and transcriptome data from non-axenic zoospores of Olpidium bornovanus and used a metagenome approach to extract phylogenetically informative fungal markers. Our phylogenetic reconstruction strongly supported Olpidium as the closest zoosporic relative of the non-flagellated terrestrial fungi. Super-alignment analyses resolved Olpidium as sister to the non-flagellated terrestrial fungi, whereas a super-tree approach recovered different placements of Olpidium, but without strong support. Further investigations detected little conflicting signal among the sampled markers but revealed a potential polytomy in early fungal evolution associated with the branching order among Olpidium, Zoopagomycota and Mucoromycota. The branches defining the evolutionary relationships of these lineages were characterized by short branch lengths and low phylogenetic content and received equivocal support for alternative phylogenetic hypotheses from individual markers. These nodes were marked by important morphological innovations, including the transition to hyphal growth and the loss of flagellum, which enabled early fungi to explore new niches and resulted in rapid and temporally concurrent Precambrian diversifications of the ancestors of several phyla of fungi.

The zoosporic obligate endoparasites, Olpidium, hold a pivotal position to the reconstruction of the flagellum loss in fungi, one of the key morphological transitions associated with the colonization of land by the early fungi. We generated genome and transcriptome data from non-axenic zoospores of Olpidium bornovanus and used a metagenome approach to extract phylogenetically informative fungal markers. Our phylogenetic reconstruction strongly supported Olpidium as the closest zoosporic relative of the non-flagellated terrestrial fungi. Super-alignment analyses resolved Olpidium as sister to the non-flagellated terrestrial fungi, whereas a super-tree approach recovered different placements of Olpidium, but without strong support. Further investigations detected little conflicting signal among the sampled markers but revealed a potential polytomy in early fungal evolution associated with the branching order among Olpidium, Zoopagomycota and Mucoromycota. The branches defining the evolutionary relationships of these lineages were characterized by short branch lengths and low phylogenetic content and received equivocal support for alternative phylogenetic hypotheses from individual markers. These nodes were marked by important morphological innovations, including the transition to hyphal growth and the loss of flagellum, which enabled early fungi to explore new niches and resulted in rapid and temporally concurrent Precambrian diversifications of the ancestors of several phyla of fungi.

The recent accumulation of newly discovered fungal-bacterial mutualisms challenges the paradigm that fungi and bacteria are natural antagonists. To understand the mechanisms that govern the establishment and maintenance over evolutionary time of mutualisms between fungi and bacteria, we studied a symbiosis of the fungus Rhizopus microsporus (Mucoromycotina) and its Burkholderia endobacteria. We found that nonhost R. microsporus, as well as other mucoralean fungi, interact antagonistically with endobacteria derived from the host and are not invaded by them. Comparison of gene expression profiles of host and nonhost fungi during interaction with endobacteria revealed dramatic changes in expression of lipid metabolic genes in the host. Analysis of the host lipidome confirmed that symbiosis establishment was accompanied by specific changes in the fungal lipid profile. Diacylglycerol kinase (DGK) activity was important for these lipid metabolic changes, as its inhibition altered the fungal lipid profile and caused a shift in the host-bacterial interaction into an antagonism. We conclude that adjustments in host lipid metabolism during symbiosis establishment, mediated by DGKs, are required for the mutualistic outcome of the Rhizopus-Burkholderia symbiosis. In addition, the neutral and phospholipid profiles of R. microsporus provide important insights into lipid metabolism in an understudied group of oleaginous Mucoromycotina. Lastly, our study revealed that the DGKs involved in the symbiosis form a previously uncharacterized clade of DGK domain proteins.

© 2018, This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply. Aspergillus section Nigri comprises filamentous fungi relevant to biomedicine, bioenergy, health, and biotechnology. To learn more about what genetically sets these species apart, as well as about potential applications in biotechnology and biomedicine, we sequenced 23 genomes de novo, forming a full genome compendium for the section (26 species), as well as 6 Aspergillus niger isolates. This allowed us to quantify both inter- and intraspecies genomic variation. We further predicted 17,903 carbohydrate-active enzymes and 2,717 secondary metabolite gene clusters, which we condensed into 455 distinct families corresponding to compound classes, 49% of which are only found in single species. We performed metabolomics and genetic engineering to correlate genotypes to phenotypes, as demonstrated for the metabolite aurasperone, and by heterologous transfer of citrate production to Aspergillus nidulans. Experimental and computational analyses showed that both secondary metabolism and regulation are key factors that are significant in the delineation of Aspergillus species.

International audience; Ustilaginomycotina is home to a broad array of fungi including important plant pathogens collectively called smut fungi. Smuts are biotrophs that produce characteristic perennating propagules called teliospores, one of which, Ustilago maydis, is a model genetic organism. Broad exploration of smut biology has been hampered by limited phylogenetic resolution of Ustilaginiomycotina as well as an overall lack of genomic data for members of this subphylum. In this study, we sequenced eight Ustilaginomycotina genomes from previously unrepresented lineages, deciphered ordinal-level phylogenetic relationships for the subphylum, and performed comparative analyses. Unlike other Basidiomycota subphyla, all sampled Ustilaginomycotina genomes are relatively small and compact. Ancestral state reconstruction analyses indicate that teliospore formation was present at the origin of the subphylum. Divergence time estimation dates the divergence of most extant smut fungi after that of grasses (Poaceae). However, we found limited conservation of well-characterized genes related to smut pathogenesis from U. maydis, indicating dissimilar pathogenic mechanisms exist across other smut lineages. The genomes of Malasseziomycetes are highly diverged from the other sampled Ustilaginomycotina, likely due to their unique history as mammal-associated lipophilic yeasts. Despite extensive genomic data, the phylogenetic placement of this class remains ambiguous. Although the sampled Ustilaginomycotina members lack many core enzymes for plant cell wall decomposition and starch catabolism, we identified several novel carbohydrate active enzymes potentially related to pectin breakdown. Finally, ∼50% of Ustilaginomycotina species-specific genes are present in previously undersampled and rare lineages, highlighting the importance of exploring fungal diversity as a resource for novel gene discovery.

Endogonales (Mucoromycotina), composed of Endogonaceae and Densosporaceae, is the only known non-Dikarya order with ectomycorrhizal members. They also form mycorrhizal-like association with some nonspermatophyte plants. It has been recently proposed that Endogonales were among the earliest mycorrhizal partners with land plants. It remains unknown whether Endogonales possess genomes with mycorrhizal-lifestyle signatures and whether Endogonales originated around the same time as land plants did. We sampled sporocarp tissue from four Endogonaceae collections and performed shotgun genome sequencing. After binning the metagenome data, we assembled and annotated the Endogonaceae genomes. We performed comparative analysis on plant-cell-wall-degrading enzymes (PCWDEs) and small secreted proteins (SSPs). We inferred phylogenetic placement of Endogonaceae and estimated the ages of Endogonaceae and Endogonales with expanded taxon sampling. Endogonaceae have large genomes with high repeat content, low diversity of PCWDEs, but without elevated SSP/secretome ratios. Dating analysis estimated that Endogonaceae originated in the Permian-Triassic boundary and Endogonales originated in the mid-late Silurian. Mycoplasma-related endobacterium sequences were identified in three Endogonaceae genomes. Endogonaceae genomes possess typical signatures of mycorrhizal lifestyle. The early origin of Endogonales suggests that the mycorrhizal association between Endogonales and plants might have played an important role during the colonization of land by plants.

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells., Author summary Genetically identical cells growing in the same environment can vary in their cellular state and behavior. Such heterogeneity may explain why some cells in an isogenic population can survive sudden severe environmental stress whereas other cells succumb. Cell-to-cell variation in gene expression has been linked to variable stress survival, but how and why transcript levels vary across the transcriptome in single cells is only beginning to emerge. Here, we used single-cell RNA sequencing (scRNA-seq) to measure cell-to-cell heterogeneity in the transcriptome of budding yeast (Saccharomyces cerevisiae). We find surprising patterns of variation across known sets of transcription factor targets, indicating that cells vary in their transcriptome profile both before and after stress exposure. scRNA-seq analysis combined with live-cell imaging of transcription factor activation dynamics revealed some cells in which the stress response was coordinately activated and other cells in which the traditional response was decoupled, suggesting unrecognized regulatory nuances that expand our understanding of stress response and survival.

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in singleS. cerevisiaecells treated with and without salt stress, to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells decouple expression of ribosomal protein genes from the environmental stress response, in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results reveal that cells maintain a cellular filter against decoupled bursts of transcription-factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells.

Setaria viridis is a rapid-life-cycle model panicoid grass. To identify genes that may contribute to inflorescence architecture and thus have the potential to influence grain yield in related crops such as maize, we conducted an N-nitroso-N-methylurea (NMU) mutagenesis of S. viridis and screened for visible inflorescence mutant phenotypes. Of the approximately 2,700 M2 families screened, we identified four recessive sparse panicle mutants (spp1-spp4) characterized by reduced and uneven branching of the inflorescence. To identify the gene underlying the sparse panicle1 (spp1) phenotype, we performed bulked segregant analysis and deep sequencing to fine map it to an approximately 1 Mb interval. Within this interval, we identified disruptive mutations in two genes. Complementation tests between spp1 and spp3 revealed they were allelic, and deep sequencing of spp3 identified an independent disruptive mutation in SvAUX1 (AUXIN1), one of the two genes in the ∼1 Mb interval and the only gene disruption shared between spp1 and spp3. SvAUX1 was found to affect both inflorescence development and root gravitropism in S. viridis. A search for orthologous mutant alleles in maize confirmed a very similar role of ZmAUX1 in maize, which highlights the utility of S. viridis in accelerating functional genomic studies in maize.

Summary Mixia osmundae (Basidiomycota, Pucciniomycotina) represents a monotypic class containing an unusual fern pathogen with incompletely understood biology. We sequenced and analyzed the genome of M. osmundae, focusing on genes that may provide some insight into its mode of pathogenicity and reproductive biology. Mixia osmundae has the smallest plant pathogenic basidiomycete genome sequenced to date, at 13.6 Mb, with very few repeats, high gene density, and relatively few significant gene family gains. The genome shows that the yeast state of M. osmundae is haploid and the lack of segregation of mating genes suggests that the spores produced on Osmunda spp. fronds are probably asexual. However, our finding of a complete complement of mating and meiosis genes suggests the capacity to undergo sexual reproduction. Analyses of carbohydrate active enzymes suggest that this fungus is a biotroph with the ability to break down several plant cell wall components. Analyses of publicly available sequence data show that other Mixia members may exist on other plant hosts and with a broader distribution than previously known.