The rise of single-cell transcriptomics in yeast.

The field of single-cell omics has transformed our understanding of biological processes and is constantly advancing both experimentally and computationally. One of the most significant developments is the ability to measure the transcriptome of individual cells by single-cell RNA-seq (scRNA-seq), which was pioneered in higher eukaryotes. While yeast has served as a powerful model organism in which to test and develop transcriptomic technologies, the implementation of scRNA-seq has been significantly delayed in this organism, mainly because of technical constraints associated with its intrinsic characteristics, namely the presence of a cell wall, a small cell size and little amounts of RNA. In this review, we examine the current technologies for scRNA-seq in yeast and highlight their strengths and weaknesses. Additionally, we explore opportunities for developing novel technologies and the potential outcomes of implementing single-cell transcriptomics and extension to other modalities. Undoubtedly, scRNA-seq will be invaluable for both basic and applied yeast research, providing unique insights into fundamental biological processes.

Positive feedback induces switch between distributive and processive phosphorylation of Hog1.

Cellular decision making often builds on ultrasensitive MAPK pathways. The phosphorylation mechanism of MAP kinase has so far been described as either distributive or processive, with distributive mechanisms generating ultrasensitivity in theoretical analyses. However, the in vivo mechanism of MAP kinase phosphorylation and its activation dynamics remain unclear. Here, we characterize the regulation of the MAP kinase Hog1 in Saccharomyces cerevisiae via topologically different ODE models, parameterized on multimodal activation data. Interestingly, our best fitting model switches between distributive and processive phosphorylation behavior regulated via a positive feedback loop composed of an affinity and a catalytic component targeting the MAP kinase-kinase Pbs2. Indeed, we show that Hog1 directly phosphorylates Pbs2 on serine 248 (S248), that cells expressing a non-phosphorylatable (S248A) or phosphomimetic (S248E) mutant show behavior that is consistent with simulations of disrupted or constitutively active affinity feedback and that Pbs2-S248E shows significantly increased affinity to Hog1 in vitro. Simulations further suggest that this mixed Hog1 activation mechanism is required for full sensitivity to stimuli and to ensure robustness to different perturbations.

Structural disruption of BAF chromatin remodeller impairs neuroblastoma metastasis by reverting an invasiveness epigenomic program.

BACKGROUND: Epigenetic programming during development is essential for determining cell lineages, and alterations in this programming contribute to the initiation of embryonal tumour development. In neuroblastoma, neural crest progenitors block their course of natural differentiation into sympathoadrenergic cells, leading to the development of aggressive and metastatic paediatric cancer. Research of the epigenetic regulators responsible for oncogenic epigenomic networks is crucial for developing new epigenetic-based therapies against these tumours. Mammalian switch/sucrose non-fermenting (mSWI/SNF) ATP-dependent chromatin remodelling complexes act genome-wide translating epigenetic signals into open chromatin states. The present study aimed to understand the contribution of mSWI/SNF to the oncogenic epigenomes of neuroblastoma and its potential as a therapeutic target. METHODS: Functional characterisation of the mSWI/SNF complexes was performed in neuroblastoma cells using proteomic approaches, loss-of-function experiments, transcriptome and chromatin accessibility analyses, and in vitro and in vivo assays. RESULTS: Neuroblastoma cells contain three main mSWI/SNF subtypes, but only BRG1-associated factor (BAF) complex disruption through silencing of its key structural subunits, ARID1A and ARID1B, impairs cell proliferation by promoting cell cycle blockade. Genome-wide chromatin remodelling and transcriptomic analyses revealed that BAF disruption results in the epigenetic repression of an extensive invasiveness-related expression program involving integrins, cadherins, and key mesenchymal regulators, thereby reducing adhesion to the extracellular matrix and the subsequent invasion in vitro and drastically inhibiting the initiation and growth of neuroblastoma metastasis in vivo. CONCLUSIONS: We report a novel ATPase-independent role for the BAF complex in maintaining an epigenomic program that allows neuroblastoma invasiveness and metastasis, urging for the development of new BAF pharmacological structural disruptors for therapeutic exploitation in metastatic neuroblastoma.

Data-driven identification of inherent features of eukaryotic stress-responsive genes.

Living organisms are continuously challenged by changes in their environment that can propagate to stresses at the cellular level, such as rapid changes in osmolarity or oxygen tension. To survive these sudden changes, cells have developed stress-responsive mechanisms that tune cellular processes. The response of Saccharomyces cerevisiae to osmostress includes a massive reprogramming of gene expression. Identifying the inherent features of stress-responsive genes is of significant interest for understanding the basic principles underlying the rewiring of gene expression upon stress. Here, we generated a comprehensive catalog of osmostress-responsive genes from 5 independent RNA-seq experiments. We explored 30 features of yeast genes and found that 25 (83%) were distinct in osmostress-responsive genes. We then identified 13 non-redundant minimal osmostress gene traits and used statistical modeling to rank the most stress-predictive features. Intriguingly, the most relevant features of osmostress-responsive genes are the number of transcription factors targeting them and gene conservation. Using data on HeLa samples, we showed that the same features that define yeast osmostress-responsive genes can predict osmostress-responsive genes in humans, but with changes in the rank-ordering of feature-importance. Our study provides a holistic understanding of the basic principles of the regulation of stress-responsive gene expression across eukaryotes.

Sensitive high-throughput single-cell RNA-seq reveals within-clonal transcript correlations in yeast populations.

Single-cell RNA sequencing has revealed extensive cellular heterogeneity within many organisms, but few methods have been developed for microbial clonal populations. The yeast genome displays unusually dense transcript spacing, with interleaved and overlapping transcription from both strands, resulting in a minuscule but complex pool of RNA that is protected by a resilient cell wall. Here, we have developed a sensitive, scalable and inexpensive yeast single-cell RNA-seq (yscRNA-seq) method that digitally counts transcript start sites in a strand- and isoform-specific manner. YscRNA-seq detects the expression of low-abundance, noncoding RNAs and at least half of the protein-coding genome in each cell. In clonal cells, we observed a negative correlation for the expression of sense-antisense pairs, whereas paralogs and divergent transcripts co-expressed. By combining yscRNA-seq with index sorting, we uncovered a linear relationship between cell size and RNA content. Although we detected an average of ~3.5 molecules per gene, the number of expressed isoforms is restricted at the single-cell level. Remarkably, the expression of metabolic genes is highly variable, whereas their stochastic expression primes cells for increased fitness towards the corresponding environmental challenge. These findings suggest that functional transcript diversity acts as a mechanism that provides a selective advantage to individual cells within otherwise transcriptionally heterogeneous populations.

Yeast Single-cell RNA-seq, Cell by Cell and Step by Step.

Single-cell RNA-seq (scRNA-seq) has become an established method for uncovering the intrinsic complexity within populations. Even within seemingly homogenous populations of isogenic yeast cells, there is a high degree of heterogeneity that originates from a compact and pervasively transcribed genome. Research with microorganisms such as yeast represents a major challenge for single-cell transcriptomics, due to their small size, rigid cell wall, and low RNA content per cell. Because of these technical challenges, yeast-specific scRNA-seq methodologies have recently started to appear, each one of them relying on different cell-isolation and library-preparation methods. Consequently, each approach harbors unique strengths and weaknesses that need to be considered. We have recently developed a yeast single-cell RNA-seq protocol (yscRNA-seq), which is inexpensive, high-throughput and easy-to-implement, tailored to the unique needs of yeast. yscRNA-seq provides a unique platform that combines single-cell phenotyping via index sorting with the incorporation of unique molecule identifiers on transcripts that allows to digitally count the number of molecules in a strand- and isoform-specific manner. Here, we provide a detailed, step-by-step description of the experimental and computational steps of yscRNA-seq protocol. This protocol will ease the implementation of yscRNA-seq in other laboratories and provide guidelines for the development of novel technologies.

H3K4 monomethylation dictates nucleosome dynamics and chromatin remodeling at stress-responsive genes.

Chromatin remodeling is essential for proper adaptation to extracellular stimuli. The p38-related Hog1 SAPK is an important regulator of transcription that mediates chromatin remodeling upon stress. Hog1 targets the RSC chromatin remodeling complex to stress-responsive genes and rsc deficient cells display reduced induction of gene expression. Here we show that the absence of H3K4 methylation, either achieved by deletion of the SET1 methyltransferase or by amino acid substitution of H3K4, bypasses the requirement of RSC for stress-responsive gene expression. Monomethylation of H3K4 is specifically inhibiting RSC-independent chromatin remodeling and thus, it prevents osmostress-induced gene expression. The absence of H3K4 monomethylation permits that the association of alternative remodelers with stress-responsive genes and the Swr1 complex (SWR-C) is instrumental in the induction of gene expression upon stress. Accordingly, the absence of SWR-C or histone H2A.Z results in compromised chromatin remodeling and impaired gene expression in the absence of RSC and H3K4 methylation. These results indicate that expression of stress-responsive genes is controlled by two remodeling mechanisms: RSC in the presence of monomethylated H3K4, and SWR-C in the absence of H3K4 monomethylation. Our findings point to a novel role for H3K4 monomethylation in dictating the specificity of chromatin remodeling, adding an extra layer of regulation to the transcriptional stress response.

Hog1 targets Whi5 and Msa1 transcription factors to downregulate cyclin expression upon stress.

Yeast cells have developed complex mechanisms to cope with extracellular insults. An increase in external osmolarity leads to activation of the stress-activated protein kinase Hog1, which is the main regulator of adaptive responses, such as gene expression and cell cycle progression, that are essential for cellular survival. Upon osmostress, the G1-to-S transition is regulated by Hog1 through stabilization of the cyclin-dependent kinase inhibitor Sic1 and the downregulation of G1 cyclin expression by an unclear mechanism. Here, we show that Hog1 interacts with and phosphorylates components of the core cell cycle transcriptional machinery such as Whi5 and the coregulator Msa1. Phosphorylation of these two transcriptional regulators by Hog1 is essential for inhibition of G1 cyclin expression, for control of cell morphogenesis, and for maximal cell survival upon stress. The control of both Whi5 and Msa1 by Hog1 also revealed the necessity for proper coordination of budding and DNA replication. Thus, Hog1 regulates G1 cyclin transcription upon osmostress to ensure coherent passage through Start.

A novel role for lncRNAs in cell cycle control during stress adaptation.

Eukaryotic cells have developed sophisticated systems to constantly monitor changes in the extracellular environment and to orchestrate a proper cellular response. To maximize survival, cells delay cell-cycle progression in response to environmental changes. In response to extracellular insults, stress-activated protein kinases (SAPKs) modulate cell-cycle progression and gene expression. In yeast, osmostress induces activation of the p38-related SAPK Hog1, which plays a key role in reprogramming gene expression upon osmostress. Genomic analysis has revealed the existence of a large number of long non-coding RNAs (lncRNAs) with different functions in a variety of organisms, including yeast. Upon osmostress, hundreds of lncRNAs are induced by the SAPK p38/Hog1. One gene that expresses Hog1-dependent lncRNA in an antisense orientation is the CDC28 gene, which encodes CDK1 kinase that controls the cell cycle in yeast. Cdc28 lncRNA mediates the induction of CDC28 expression and this increase in the level of Cdc28 results in more efficient re-entry of the cells into the cell cycle after stress. Thus, the control of lncRNA expression as a new mechanism for the regulation of cell-cycle progression opens new avenues to understand how stress adaptation can be accomplished in response to changing environments.

Control of Cdc28 CDK1 by a stress-induced lncRNA.

Genomic analysis has revealed the existence of a large number of long noncoding RNAs (lncRNAs) with different functions in a variety of organisms, including yeast. Cells display dramatic changes of gene expression upon environmental changes. Upon osmostress, hundreds of stress-responsive genes are induced by the stress-activated protein kinase (SAPK) p38/Hog1. Using whole-genome tiling arrays, we found that Hog1 induces a set of lncRNAs upon stress. One of the genes expressing a Hog1-dependent lncRNA in antisense orientation is CDC28, the cyclin-dependent kinase 1 (CDK1) that controls the cell cycle in yeast. Cdc28 lncRNA mediates the establishment of gene looping and the relocalization of Hog1 and RSC from the 3' UTR to the +1 nucleosome to induce CDC28 expression. The increase in the levels of Cdc28 results in cells able to reenter the cell cycle more efficiently after stress. This may represent a general mechanism to prime expression of genes needed after stresses are alleviated.

Initiation of the transcriptional response to hyperosmotic shock correlates with the potential for volume recovery.

The control of activity and localization of transcription factors is critical for appropriate transcriptional responses. In eukaryotes, signal transduction components such as mitogen-activated protein kinase (MAPK) shuttle into the nucleus to activate transcription. It is not known in detail how different amounts of nuclear MAPK over time affect the transcriptional response. In the present study, we aimed to address this issue by studying the high osmolarity glycerol (HOG) system in Saccharomyces cerevisiae. We employed a conditional osmotic system, which changes the period of the MAPK Hog1 signal independent of the initial stress level. We determined the dynamics of the Hog1 nuclear localization and cell volume by single-cell analysis in well-controlled microfluidics systems and compared the responses with the global transcriptional output of cell populations. We discovered that the onset of the initial transcriptional response correlates with the potential of cells for rapid adaptation; cells that are capable of recovering quickly initiate the transcriptional responses immediately, whereas cells that require longer time to adapt also respond later. This is reflected by Hog1 nuclear localization, Hog1 promoter association and the transcriptional response, but not Hog1 phosphorylation, suggesting that a presently uncharacterized rapid adaptive mechanism precedes the Hog1 nuclear response. Furthermore, we found that the period of Hog1 nuclear residence affects the amplitude of the transcriptional response rather than the spectrum of responsive genes.

Hog1 bypasses stress-mediated down-regulation of transcription by RNA polymerase II redistribution and chromatin remodeling.

BACKGROUND: Cells are subjected to dramatic changes of gene expression upon environmental changes. Stress causes a general down-regulation of gene expression together with the induction of a set of stress-responsive genes. The p38-related stress-activated protein kinase Hog1 is an important regulator of transcription upon osmostress in yeast. RESULTS: Genome-wide localization studies of RNA polymerase II (RNA Pol II) and Hog1 showed that stress induced major changes in RNA Pol II localization, with a shift toward stress-responsive genes relative to housekeeping genes. RNA Pol II relocalization required Hog1, which was also localized to stress-responsive loci. In addition to RNA Pol II-bound genes, Hog1 also localized to RNA polymerase III-bound genes, pointing to a wider role for Hog1 in transcriptional control than initially expected. Interestingly, an increasing association of Hog1 with stress-responsive genes was strongly correlated with chromatin remodeling and increased gene expression. Remarkably, MNase-Seq analysis showed that although chromatin structure was not significantly altered at a genome-wide level in response to stress, there was pronounced chromatin remodeling for those genes that displayed Hog1 association. CONCLUSION: Hog1 serves to bypass the general down-regulation of gene expression that occurs in response to osmostress, and does so both by targeting RNA Pol II machinery and by inducing chromatin remodeling at stress-responsive loci.

Time-dependent quantitative multicomponent control of the G1-S network by the stress-activated protein kinase Hog1 upon osmostress.

Control of cell cycle progression by stress-activated protein kinases (SAPKs) is essential for cell adaptation to extracellular stimuli. Exposure of yeast to hyperosmotic stress activates the SAPK Hog1, which delays cell cycle progression through G1 by direct phosphorylation of the cyclin-dependent kinase (CDK) inhibitor Sic1 and by inhibition of the transcription of the genes encoding the G1 cyclins Cln1 and 2. Additional targets of Hog1 may also play a role in this response. We used mathematical modeling and quantitative in vivo experiments to define the contributions of individual components of the G1-S network downstream of Hog1 to this stress-induced delay in the cell cycle. The length of the arrest depended on the degree of stress and the temporal proximity of the onset of the stress to the commitment to cell division, called "Start." Hog1-induced inhibition of the transcription of the gene encoding cyclin Clb5, rather than that of the gene encoding Cln2, prevented entry into S phase upon osmostress. By controlling the accumulation of specific cyclins, Hog1 delayed bud morphogenesis (through Clns) and delayed DNA replication (through Clb5). Hog1-mediated phosphorylation and degradation of Sic1 at Start prevented residual activity of the cyclin/CDK complex Clb5/Cdc28 from initiating DNA replication before adaptation to the stress. Thus, our work defines distinct temporal roles for the actions of Hog1 on Sic1 and cyclins in mediating G1 arrest upon hyperosmotic stress.

Control of Ubp3 ubiquitin protease activity by the Hog1 SAPK modulates transcription upon osmostress.

Protein ubiquitylation is a key process in the regulation of many cellular processes. The balance between the activity of ubiquitin ligases and that of proteases controls the level of ubiquitylation. In response to extracellular stimuli, stress-activated protein kinases (SAPK) modulate gene expression to maximize cell survival. In yeast, the Hog1 SAPK has a key role in reprogramming the gene expression pattern required for cell survival upon osmostress. Here, we show that the Ubp3 ubiquitin protease is a target for the Hog1 SAPK to modulate gene expression. ubp3 mutant cells are defective in expression of osmoresponsive genes. Hog1 interacts with and phosphorylates Ubp3 at serine 695, which is essential to determine the extent of transcriptional activation in response to osmostress. Furthermore, Ubp3 is recruited to osmoresponsive genes to modulate transcriptional initiation as well as elongation. Therefore, Ubp3 activity responds to external stimuli and is required for transcriptional activation upon osmostress.

Transient activation of the HOG MAPK pathway regulates bimodal gene expression.

Mitogen-activated protein kinase (MAPK) cascades are conserved signaling modules that control many cellular processes by integrating intra- and extracellular cues. The p38/Hog1 MAPK is transiently activated in response to osmotic stress, leading to rapid translocation into the nucleus and induction of a specific transcriptional program. When investigating the dynamic interplay between Hog1 activation and Hog1-driven gene expression, we found that Hog1 activation increases linearly with stimulus, whereas the transcriptional output is bimodal. Modeling predictions, corroborated by single-cell experiments, established that a slow stochastic transition from a repressed to an activated transcriptional state in conjunction with transient Hog1 activation generates this behavior. Together, these findings provide a molecular mechanism by which a cell can impose a transcriptional threshold in response to a linear signaling behavior.