Preassembled Cas9 Ribonucleoprotein-Mediated Gene Deletion Identifies the Carbon Catabolite Repressor and Its Target Genes in Coprinopsis cinerea.

Cre1 is an important transcription factor that regulates carbon catabolite repression (CCR) and is widely conserved across fungi. The cre1 gene has been extensively studied in several Ascomycota species, whereas its role in gene expression regulation in the Basidiomycota species remains poorly understood. Here, we identified and investigated the role of cre1 in Coprinopsis cinerea, a basidiomycete model mushroom that can efficiently degrade lignocellulosic plant wastes. We used a rapid and efficient gene deletion approach based on PCR-amplified split-marker DNA cassettes together with in vitro assembled Cas9-guide RNA ribonucleoproteins (Cas9 RNPs) to generate C. cinerea cre1 gene deletion strains. Gene expression profiling of two independent C. cinerea cre1 mutants showed significant deregulation of carbohydrate metabolism, plant cell wall degrading enzymes (PCWDEs), plasma membrane transporter-related and several transcription factor-encoding genes, among others. Our results support the notion that, like reports in the ascomycetes, Cre1 of C. cinerea orchestrates CCR through a combined regulation of diverse genes, including PCWDEs, transcription factors that positively regulate PCWDEs, and membrane transporters which could import simple sugars that can induce the expression of PWCDEs. Somewhat paradoxically, though in accordance with other Agaricomycetes, genes related to lignin degradation were mostly downregulated in cre1 mutants, indicating they fall under different regulation than other PCWDEs. The gene deletion approach and the data presented here will expand our knowledge of CCR in the Basidiomycota and provide functional hypotheses on genes related to plant biomass degradation. IMPORTANCE Mushroom-forming fungi include some of the most efficient lignocellulosic plant biomass degraders. They degrade dead plant materials by a battery of lignin-, cellulose-, hemicellulose-, and pectin-degrading enzymes, the encoding genes of which are under tight transcriptional control. One of the highest-level regulations of these metabolic enzymes is known as carbon catabolite repression, which is orchestrated by the transcription factor Cre1, and ensures that costly lignocellulose-degrading enzyme genes are expressed only when simple carbon sources (e.g., glucose) are not available. Here, we identified the Cre1 ortholog in a litter decomposer Agaricomycete, Coprinopsis cinerea, knocked it out, and characterized transcriptional changes in the mutants. We identified several dozen lignocellulolytic enzyme genes as well as membrane transporters and other transcription factors as putative target genes of C. cinerea cre1. These results extend knowledge on carbon catabolite repression to litter decomposer Basidiomycota.

Novel phylogenetic methods are needed for understanding gene function in the era of mega-scale genome sequencing.

Ongoing large-scale genome sequencing projects are forecasting a data deluge that will almost certainly overwhelm current analytical capabilities of evolutionary genomics. In contrast to population genomics, there are no standardized methods in evolutionary genomics for extracting evolutionary and functional (e.g. gene-trait association) signal from genomic data. Here, we examine how current practices of multi-species comparative genomics perform in this aspect and point out that many genomic datasets are under-utilized due to the lack of powerful methodologies. As a result, many current analyses emphasize gene families for which some functional data is already available, resulting in a growing gap between functionally well-characterized genes/organisms and the universe of unknowns. This leaves unknown genes on the 'dark side' of genomes, a problem that will not be mitigated by sequencing more and more genomes, unless we develop tools to infer functional hypotheses for unknown genes in a systematic manner. We provide an inventory of recently developed methods capable of predicting gene-gene and gene-trait associations based on comparative data, then argue that realizing the full potential of whole genome datasets requires the integration of phylogenetic comparative methods into genomics, a rich but underutilized toolbox for looking into the past.

Transcriptomic atlas of mushroom development reveals conserved genes behind complex multicellularity in fungi.

The evolution of complex multicellularity has been one of the major transitions in the history of life. In contrast to simple multicellular aggregates of cells, it has evolved only in a handful of lineages, including animals, embryophytes, red and brown algae, and fungi. Despite being a key step toward the evolution of complex organisms, the evolutionary origins and the genetic underpinnings of complex multicellularity are incompletely known. The development of fungal fruiting bodies from a hyphal thallus represents a transition from simple to complex multicellularity that is inducible under laboratory conditions. We constructed a reference atlas of mushroom formation based on developmental transcriptome data of six species and comparisons of >200 whole genomes, to elucidate the core genetic program of complex multicellularity and fruiting body development in mushroom-forming fungi (Agaricomycetes). Nearly 300 conserved gene families and >70 functional groups contained developmentally regulated genes from five to six species, covering functions related to fungal cell wall remodeling, targeted protein degradation, signal transduction, adhesion, and small secreted proteins (including effector-like orphan genes). Several of these families, including F-box proteins, expansin-like proteins, protein kinases, and transcription factors, showed expansions in Agaricomycetes, many of which convergently expanded in multicellular plants and/or animals too, reflecting convergent solutions to genetic hurdles imposed by complex multicellularity among independently evolved lineages. This study provides an entry point to studying mushroom development and complex multicellularity in one of the largest clades of complex eukaryotic organisms.

Unmatched Level of Molecular Convergence among Deeply Divergent Complex Multicellular Fungi.

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

Starvation- and xenobiotic-related transcriptomic responses of the sulfanilic acid-degrading bacterium, Novosphingobium resinovorum SA1.

Novosphingobium resinovorum SA1 was the first single isolate capable of degrading sulfanilic acid, a widely used representative of sulfonated aromatic compounds. The genome of the strain was recently sequenced, and here, we present whole-cell transcriptome analyses of cells exposed to sulfanilic acid as compared to cells grown on glucose. The comparison of the transcript profiles suggested that the primary impact of sulfanilic acid on the cell transcriptome was a starvation-like effect. The genes of the peripheral, central, and common pathways of sulfanilic acid biodegradation had distinct transcript profiles. The peripheral genes located on a plasmid had very high basal expressions which were hardly upregulated by sulfanilic acid. The genomic context and the codon usage preference of these genes suggested that they were acquired by horizontal gene transfer. The genes of the central pathways were remarkably inducible by sulfanilic acid indicating the presence of a substrate-specific regulatory system in the cells. Surprisingly, the genes of the common part of the metabolic pathway had low and sulfanilic acid-independent transcript levels. The approach applied resulted in the identification of the genes of proteins involved in auxiliary processes such as electron transfer, substrate and iron transports, sulfite oxidases, and sulfite transporters. The whole transcriptome analysis revealed that the cells exposed to xenobiotics had multiple responses including general starvation-like, substrate-specific, and substrate-related effects. From the results, we propose that the genes of the peripheral, central, and common parts of the pathway have been evolved independently.

Metabolic responses of Rhodococcus erythropolis PR4 grown on diesel oil and various hydrocarbons.

Rhodococcus erythropolis PR4 is able to degrade diesel oil, normal-, iso- and cycloparaffins and aromatic compounds. The complete DNA content of the strain was previously sequenced and numerous oxygenase genes were identified. In order to identify the key elements participating in biodegradation of various hydrocarbons, we performed a comparative whole transcriptome analysis of cells grown on hexadecane, diesel oil and acetate. The transcriptomic data for the most prominent genes were validated by RT-qPCR. The expression of two genes coding for alkane-1-monooxygenase enzymes was highly upregulated in the presence of hydrocarbon substrates. The transcription of eight phylogenetically diverse cytochrome P450 (cyp) genes was upregulated in the presence of diesel oil. The transcript levels of various oxygenase genes were determined in cells grown in an artificial mixture, containing hexadecane, cycloparaffin and aromatic compounds and six cyp genes were induced by this hydrocarbon mixture. Five of them were not upregulated by linear and branched hydrocarbons. The expression of fatty acid synthase I genes was downregulated by hydrocarbon substrates, indicating the utilization of external alkanes for fatty acid synthesis. Moreover, the transcription of genes involved in siderophore synthesis, iron transport and exopolysaccharide biosynthesis was also upregulated, indicating their important role in hydrocarbon metabolism. Based on the results, complex metabolic response profiles were established for cells grown on various hydrocarbons. Our results represent a functional annotation of a rhodococcal genome, provide deeper insight into molecular events in diesel/hydrocarbon utilization and suggest novel target genes for environmental monitoring projects.

ContScout: sensitive detection and removal of contamination from annotated genomes.

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

Snowball: a novel gene family required for developmental patterning of fruiting bodies of mushroom-forming fungi (Agaricomycetes).

The morphogenesis of sexual fruiting bodies of fungi is a complex process determined by a genetically encoded program. Fruiting bodies reached the highest complexity levels in the Agaricomycetes; yet, the underlying genetics is currently poorly known. In this work, we functionally characterized a highly conserved gene termed snb1, whose expression level increases rapidly during fruiting body initiation. According to phylogenetic analyses, orthologs of snb1 are present in almost all agaricomycetes and may represent a novel conserved gene family that plays a substantial role in fruiting body development. We disrupted snb1 using CRISPR/Cas9 in the agaricomycete model organism Coprinopsis cinerea. snb1 deletion mutants formed unique, snowball-shaped, rudimentary fruiting bodies that could not differentiate caps, stipes, and lamellae. We took advantage of this phenotype to study fruiting body differentiation using RNA-Seq analyses. This revealed differentially regulated genes and gene families that, based on wild-type RNA-Seq data, were upregulated early during development and showed tissue-specific expression, suggesting a potential role in differentiation. Taken together, the novel gene family of snb1 and the differentially expressed genes in the snb1 mutants provide valuable insights into the complex mechanisms underlying developmental patterning in the Agaricomycetes. IMPORTANCE: Fruiting bodies of mushroom-forming fungi (Agaricomycetes) are complex multicellular structures, with a spatially and temporally integrated developmental program that is, however, currently poorly known. In this study, we present a novel, conserved gene family, Snowball (snb), termed after the unique, differentiation-less fruiting body morphology of snb1 knockout strains in the model mushroom Coprinopsis cinerea. snb is a gene of unknown function that is highly conserved among agaricomycetes and encodes a protein of unknown function. A comparative transcriptomic analysis of the early developmental stages of differentiated wild-type and non-differentiated mutant fruiting bodies revealed conserved differentially expressed genes which may be related to tissue differentiation and developmental patterning fruiting body development.

Extreme overall mushroom genome expansion in Mycena s.s. irrespective of plant hosts or substrate specializations.

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.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom-forming fungi (Agaricomycetes).

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a 'developmental hourglass,' act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Comparative genomics reveals the origin of fungal hyphae and multicellularity.

Hyphae represent a hallmark structure of multicellular fungi. The evolutionary origins of hyphae and of the underlying genes are, however, hardly known. By systematically analyzing 72 complete genomes, we here show that hyphae evolved early in fungal evolution probably via diverse genetic changes, including co-option and exaptation of ancient eukaryotic (e.g. phagocytosis-related) genes, the origin of new gene families, gene duplications and alterations of gene structure, among others. Contrary to most multicellular lineages, the origin of filamentous fungi did not correlate with expansions of kinases, receptors or adhesive proteins. Co-option was probably the dominant mechanism for recruiting genes for hypha morphogenesis, while gene duplication was apparently less prevalent, except in transcriptional regulators and cell wall - related genes. We identified 414 novel gene families that show correlated evolution with hyphae and that may have contributed to its evolution. Our results suggest that hyphae represent a unique multicellular organization that evolved by limited fungal-specific innovations and gene duplication but pervasive co-option and modification of ancient eukaryotic functions.

Complete genome sequence of Novosphingobium resinovorum SA1, a versatile xenobiotic-degrading bacterium capable of utilizing sulfanilic acid.

Sulfanilic acid (4-aminobenzenesulfonic acid) is a sulfonated aromatic amine widely used in chemical industries for synthesis of various organic dyes and sulfa drugs. There are quite a few microbial co-cultures or single isolates capable of completely degrading this compound. Novosphingobium resinovorum SA1 was the first single bacterium which could utilize sulfanilic acid as its sole carbon, nitrogen and sulfur source. The strain has versatile catabolic routes for the bioconversion of numerous other aromatic compounds. Here, the complete genome sequence of the N. resinovorum SA1 strain is reported. The genome consists of a circular chromosome of 3.8 Mbp and four extrachromosomal elements between 67 and 1 759.8 kbp in size. Three alternative 3-ketoadipate pathways were identified on the plasmids. Sulfanilic acid is decomposed via a modified 3-ketoadipate pathway and the oxygenases involved form a phylogenetically separate branch on the tree. Sequence analysis of these elements might provide a genetic background for deeper insight into the versatile catabolic metabolism of various aromatic xenobiotics, including sulfanilic acid and its derivatives. Moreover, this is also a good model strain for understanding the role and evolution of multiple genetic elements within a single strain.