Systematic Analysis of the Determinants of Gene Expression Noise in Embryonic Stem Cells

Multimodal screen identifies noise-regulatory proteins

Gene-expression noise can influence cell-fate choices across pathology and physiology. However, a crucial question persists: do regulatory proteins or pathways exist that control noise independently of mean expression levels? Our integrative approach, combining single-cell RNA sequencing with proteomics and regulator enrichment analysis, identifies 32 putative noise regulators. SON, a nuclear speckle-associated protein, alters transcriptional noise without changing mean expression levels. Furthermore, SON's noise control can propagate to the protein level. Long-read and total RNA sequencing shows that SON's noise control does not significantly change isoform usage or splicing efficiency. Moreover, SON depletion reduces state switching in pluripotent mouse embryonic stem cells and impacts their fate choice during differentiation. Collectively, we demonstrate a class of proteins that control noise orthogonally to mean expression levels. This work serves as a proof of concept that can identify other functional noise regulators throughout development and disease progression.

EpiCHAOS: a metric to quantify epigenomic heterogeneity in single-cell data

Epigenetic heterogeneity is a fundamental property of biological systems and is recognized as a potential driver of tumor plasticity and therapy resistance. Single-cell epigenomics technologies have been widely employed to study epigenetic variation between-but not within-cellular clusters. We introduce epiCHAOS: a quantitative metric of cell-to-cell heterogeneity, applicable to any single-cell epigenomics data type. After validation in synthetic datasets, we apply epiCHAOS to investigate global and region-specific patterns of epigenetic heterogeneity across diverse biological systems. EpiCHAOS provides an excellent approximation of stemness and plasticity in development and malignancy, making it a valuable addition to single-cell cancer epigenomics analyses.

Preexisting cell state rather than stochastic noise confers high or low infection susceptibility of human lung epithelial cells to adenovirus

Viruses display large variability across all stages of their life cycle, including entry, gene expression, replication, assembly, and egress. We previously reported that the immediate early adenovirus (AdV) E1A transcripts accumulate in human lung epithelial A549 cancer cells with high variability, mostly independent of the number of incoming viral genomes, but somewhat correlated to the cell cycle state at the time of inoculation. Here, we leveraged the classical Luria-Delbrück fluctuation analysis to address whether infection variability primarily arises from the cell state or stochastic noise. The E1A expression was measured by the expression of green fluorescent protein (GFP) from the endogenous E1A promoter in AdV-C5_E1A-FS2A-GFP and found to be highly correlated with the viral plaque formation, indicating reliability of the reporter virus. As an ensemble, randomly picked clonal A549 cell isolates displayed significantly higher coefficients of variation in the E1A expression than technical noise, indicating a phenotypic variability larger than noise. The underlying cell state determining infection variability was maintained for at least 9 weeks of cell cultivation. Our results indicate that preexisting cell states tune adenovirus infection in favor of the cell or the virus. These findings have implications for antiviral strategies and gene therapy applications.IMPORTANCEViral infections are known for their variability. Underlying mechanisms are still incompletely understood but have been associated with particular cell states, for example, the eukaryotic cell division cycle in DNA virus infections. A cell state is the collective of biochemical, morphological, and contextual features owing to particular conditions or at random. It affects how intrinsic or extrinsic cues trigger a response, such as cell division or anti-viral state. Here, we provide evidence that cell states with a built-in memory confer high or low susceptibility of clonal human epithelial cells to adenovirus infection. Results are reminiscent of the Luria-Delbrück fluctuation test with bacteriophage infections back in 1943, which demonstrated that mutations, in the absence of selective pressure prior to infection, cause infection resistance rather than being a consequence of infection. Our findings of dynamic cell states conferring adenovirus infection susceptibility uncover new challenges for the prediction and treatment of viral infections.

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Self-organization underlies developmental robustness in plants

Development is a self-organized process that builds on cells and their interactions. Cells are heterogeneous in gene expression, growth, and division; yet how development is robust despite such heterogeneity is a fascinating question. Here, we review recent progress on this topic, highlighting how developmental robustness is achieved through self-organization. We will first discuss sources of heterogeneity, including stochastic gene expression, heterogeneity in growth rate and direction, and heterogeneity in division rate and precision. We then discuss cellular mechanisms that buffer against such noise, including Paf1C- and miRNA-mediated denoising, spatiotemporal growth averaging and compensation, mechanisms to improve cell division precision, and coordination of growth rate and developmental timing between different parts of an organ. We also discuss cases where such heterogeneity is not buffered but utilized for development. Finally, we highlight potential directions for future studies of noise and developmental robustness.

Effect of genomic and cellular environments on gene expression noise

BACKGROUND

Individual cells from isogenic populations often display large cell-to-cell differences in gene expression. This "noise" in expression derives from several sources, including the genomic and cellular environment in which a gene resides. Large-scale maps of genomic environments have revealed the effects of epigenetic modifications and transcription factor occupancy on mean expression levels, but leveraging such maps to explain expression noise will require new methods to assay how expression noise changes at locations across the genome.

RESULTS

To address this gap, we present Single-cell Analysis of Reporter Gene Expression Noise and Transcriptome (SARGENT), a method that simultaneously measures the noisiness of reporter genes integrated throughout the genome and the global mRNA profiles of individual reporter-gene-containing cells. Using SARGENT, we perform the first comprehensive genome-wide survey of how genomic locations impact gene expression noise. We find that the mean and noise of expression correlate with different histone modifications. We quantify the intrinsic and extrinsic components of reporter gene noise and, using the associated mRNA profiles, assign the extrinsic component to differences between the CD24+ "stem-like" substate and the more "differentiated" substate. SARGENT also reveals the effects of transgene integrations on endogenous gene expression, which will help guide the search for "safe-harbor" loci.

CONCLUSIONS

Taken together, we show that SARGENT is a powerful tool to measure both the mean and noise of gene expression at locations across the genome and that the data generatd by SARGENT reveals important insights into the regulation of gene expression noise genome-wide.

BASiCS workflow: a step-by-step analysis of expression variability using single cell RNA sequencing data

Cell-to-cell gene expression variability is an inherent feature of complex biological systems, such as immunity and development. Single-cell RNA sequencing is a powerful tool to quantify this heterogeneity, but it is prone to strong technical noise. In this article, we describe a step-by-step computational workflow that uses the BASiCS Bioconductor package to robustly quantify expression variability within and between known groups of cells (such as experimental conditions or cell types). BASiCS uses an integrated framework for data normalisation, technical noise quantification and downstream analyses, propagating statistical uncertainty across these steps. Within a single seemingly homogeneous cell population, BASiCS can identify highly variable genes that exhibit strong heterogeneity as well as lowly variable genes with stable expression. BASiCS also uses a probabilistic decision rule to identify changes in expression variability between cell populations, whilst avoiding confounding effects related to differences in technical noise or in overall abundance. Using a publicly available dataset, we guide users through a complete pipeline that includes preliminary steps for quality control, as well as data exploration using the scater and scran Bioconductor packages. The workflow is accompanied by a Docker image that ensures the reproducibility of our results.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Increased heterogeneity in expression of genes associated with cancer progression and drug resistance

Fluctuations in the number of regulatory molecules and differences in timings of molecular events can generate variation in gene expression among genetically identical cells in the same environmental condition. This variation, termed as expression noise, can create differences in metabolic state and cellular functions, leading to phenotypic heterogeneity. Expression noise and phenotypic heterogeneity have been recognized as important contributors to intra-tumor heterogeneity, and have been associated with cancer growth, progression, and therapy resistance. However, how expression noise changes with cancer progression in actual cancer patients has remained poorly explored. Such an analysis, through identification of genes with increasing expression noise, can provide valuable insights into generation of intra-tumor heterogeneity, and could have important implications for understanding immune-suppression, drug tolerance and therapy resistance. In this work, we performed a genome-wide identification of changes in gene expression noise with cancer progression using single-cell RNA-seq data of lung adenocarcinoma patients at different stages of cancer. We identified 37 genes in epithelial cells that showed an increasing noise trend with cancer progression, many of which were also associated with cancer growth, EMT and therapy resistance. We found that expression of several of these genes was positively associated with expression of mitochondrial genes, suggesting an important role of mitochondria in generation of heterogeneity. In addition, we uncovered substantial differences in sample-specific noise profiles which could have implications for personalized prognosis and treatment.

A concerted neuron-astrocyte program declines in ageing and schizophrenia

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

Concerted neuron-astrocyte gene expression declines in aging and schizophrenia

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a striking relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA-seq to analyze the prefrontal cortex of 191 human donors ages 22-97 years, including healthy individuals and persons with schizophrenia. Latent-factor analysis of these data revealed that in persons whose cortical neurons more strongly expressed genes for synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the Synaptic Neuron-and-Astrocyte Program (SNAP). In schizophrenia and aging - two conditions that involve declines in cognitive flexibility and plasticity 1,2 - cells had divested from SNAP: astrocytes, glutamatergic (excitatory) neurons, and GABAergic (inhibitory) neurons all reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy persons of similar age, may underlie many aspects of normal human interindividual differences and be an important point of convergence for multiple kinds of pathophysiology.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.

Transcription factor expression is the main determinant of variability in gene co-activity

Many genes are co-expressed and form genomic domains of coordinated gene activity. However, the regulatory determinants of domain co-activity remain unclear. Here, we leverage human individual variation in gene expression to characterize the co-regulatory processes underlying domain co-activity and systematically quantify their effect sizes. We employ transcriptional decomposition to extract from RNA expression data an expression component related to co-activity revealed by genomic positioning. This strategy reveals close to 1,500 co-activity domains, covering most expressed genes, of which the large majority are invariable across individuals. Focusing specifically on domains with high variability in co-activity reveals that contained genes have a higher sharing of eQTLs, a higher variability in enhancer interactions, and an enrichment of binding by variably expressed transcription factors, compared to genes within non-variable domains. Through careful quantification of the relative contributions of regulatory processes underlying co-activity, we find transcription factor expression levels to be the main determinant of gene co-activity. Our results indicate that distal trans effects contribute more than local genetic variation to individual variation in co-activity domains.

Quantifying stimulus-response specificity to probe the functional state of macrophages

Immune sentinel macrophages initiate responses to pathogens via hundreds of immune response genes. Each immune threat demands a tailored response, suggesting that the capacity for stimulus-specific gene expression is a key functional hallmark of healthy macrophages. To quantify this property, termed "stimulus-response specificity" (SRS), we developed a single-cell experimental workflow and analytical approaches based on information theory and machine learning. We found that the response specificity of macrophages is driven by combinations of specific immune genes that show low cell-to-cell heterogeneity and are targets of separate signaling pathways. The "response specificity profile," a systematic comparison of multiple stimulus-response distributions, was distinctly altered by polarizing cytokines, and it enabled an assessment of the functional state of macrophages. Indeed, the response specificity profile of peritoneal macrophages from old and obese mice showed characteristic differences, suggesting that SRS may be a basis for measuring the functional state of innate immune cells. A record of this paper's transparent peer review process is included in the supplemental information.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Stochastic variation in gene products ("noise") is an inescapable by-product of gene expression. Noise must be minimized to allow for the reliable execution of cellular functions. However, noise cannot be suppressed beyond an intrinsic lower limit. For constitutively expressed genes, this limit is believed to be Poissonian, meaning that the variance in mRNA numbers cannot be lower than their mean. Here, we show that several cell division genes in fission yeast have mRNA variances significantly below this limit, which cannot be explained by the classical gene expression model for low-noise genes. Our analysis reveals that multiple steps in both transcription and mRNA degradation are essential to explain this sub-Poissonian variance. The sub-Poissonian regime differs qualitatively from previously characterized noise regimes, a hallmark being that cytoplasmic noise is reduced when the mRNA export rate increases. Our study re-defines the lower limit of eukaryotic gene expression noise and identifies molecular requirements for ultra-low noise which are expected to support essential cell functions.

Genome-wide inference reveals that feedback regulations constrain promoter-dependent transcriptional burst kinetics

Gene expression in mammalian cells is highly variable and episodic, resulting in a series of discontinuous bursts of mRNAs. A challenge is to understand how static promoter architecture and dynamic feedback regulations dictate bursting on a genome-wide scale. Although single-cell RNA sequencing (scRNA-seq) provides an opportunity to address this challenge, effective analytical methods are scarce. We developed an interpretable and scalable inference framework, which combined experimental data with a mechanistic model to infer transcriptional burst kinetics (sizes and frequencies) and feedback regulations. Applying this framework to scRNA-seq data generated from embryonic mouse fibroblast cells, we found Simpson's paradoxes, i.e. genome-wide burst kinetics exhibit different characteristics in two cases without and with distinguishing feedback regulations. We also showed that feedbacks differently modulate burst frequencies and sizes and conceal the effects of transcription start site distributions on burst kinetics. Notably, only in the presence of positive feedback, TATA genes are expressed with high burst frequencies and enhancer-promoter interactions mainly modulate burst frequencies. The developed inference method provided a flexible and efficient way to investigate transcriptional burst kinetics and the obtained results would be helpful for understanding cell development and fate decision.

Transcription factor binding process is the primary driver of noise in gene expression

Noise in expression of individual genes gives rise to variations in activity of cellular pathways and generates heterogeneity in cellular phenotypes. Phenotypic heterogeneity has important implications for antibiotic persistence, mutation penetrance, cancer growth and therapy resistance. Specific molecular features such as the presence of the TATA box sequence and the promoter nucleosome occupancy have been associated with noise. However, the relative importance of these features in noise regulation is unclear and how well these features can predict noise has not yet been assessed. Here through an integrated statistical model of gene expression noise in yeast we found that the number of regulating transcription factors (TFs) of a gene was a key predictor of noise, whereas presence of the TATA box and the promoter nucleosome occupancy had poor predictive power. With an increase in the number of regulatory TFs, there was a rise in the number of cooperatively binding TFs. In addition, an increased number of regulatory TFs meant more overlaps in TF binding sites, resulting in competition between TFs for binding to the same region of the promoter. Through modeling of TF binding to promoter and application of stochastic simulations, we demonstrated that competition and cooperation among TFs could increase noise. Thus, our work uncovers a process of noise regulation that arises out of the dynamics of gene regulation and is not dependent on any specific transcription factor or specific promoter sequence.

Promoter sequence and architecture determine expression variability and confer robustness to genetic variants

Genetic and environmental exposures cause variability in gene expression. Although most genes are affected in a population, their effect sizes vary greatly, indicating the existence of regulatory mechanisms that could amplify or attenuate expression variability. Here, we investigate the relationship between the sequence and transcription start site architectures of promoters and their expression variability across human individuals. We find that expression variability can be largely explained by a promoter's DNA sequence and its binding sites for specific transcription factors. We show that promoter expression variability reflects the biological process of a gene, demonstrating a selective trade-off between stability for metabolic genes and plasticity for responsive genes and those involved in signaling. Promoters with a rigid transcription start site architecture are more prone to have variable expression and to be associated with genetic variants with large effect sizes, while a flexible usage of transcription start sites within a promoter attenuates expression variability and limits genotypic effects. Our work provides insights into the variable nature of responsive genes and reveals a novel mechanism for supplying transcriptional and mutational robustness to essential genes through multiple transcription start site regions within a promoter.

Comparative studies of genomic and epigenetic factors influencing transcriptional variation in two insect species

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.

Enhanced transcriptional heterogeneity mediated by NF-κB super-enhancers

The transcription factor NF-κB, which plays an important role in cell fate determination, is involved in the activation of super-enhancers (SEs). However, the biological functions of the NF-κB SEs in gene control are not fully elucidated. We investigated the characteristics of NF-κB-mediated SE activity using fluorescence imaging of RelA, single-cell transcriptome and chromatin accessibility analyses in anti-IgM-stimulated B cells. The formation of cell stimulation-induced nuclear RelA foci was abolished in the presence of hexanediol, suggesting an underlying process of liquid-liquid phase separation. The gained SEs induced a switch-like expression and enhanced cell-to-cell variability in transcriptional response. These properties were correlated with the number of gained cis-regulatory interactions, while switch-like gene induction was associated with the number of NF-κB binding sites in SE. Our study suggests that NF-κB SEs have an important role in the transcriptional regulation of B cells possibly through liquid condensate formation consisting of macromolecular interactions.

NF-κB produces an all-or-none activation response upon the activation of B cell receptors. These dynamics modulate the amplitude and frequency of target mRNA induction in cell populations. In this research, we performed single-cell assessment of chromatin accessibility and RNA expression, coupled with fluorescence imaging to characterize the nuclear dynamics of NF-κB proteins in B cell upon receptor stimulation. We found that upon cellular activation, NF-κB-mediated long-range activation of enhancers cooperatively evoked RNA production. In addition, predicted DNA contacts brought by open chromatin led to the high heterogeneity of RNA levels in cell populations. Stimuli-dependent NF-κB foci formation was further inhibited by 1,6-hexanediol (liquid-liquid phase separation inhibitor) and JQ1 (coactivator protein BRD4 inhibitor). We thus propose that nuclear NF-κB plays an important role in the transcriptional regulation of B cell development possibly through the formation of liquid condensates.

Hematopoietic differentiation is characterized by a transient peak of entropy at a single-cell level

Background

Mature blood cells arise from hematopoietic stem cells in the bone marrow by a process of differentiation along one of several different lineage trajectories. This is often represented as a series of discrete steps of increasing progenitor cell commitment to a given lineage, but as for differentiation in general, whether the process is instructive or stochastic remains controversial. Here, we examine this question by analyzing single-cell transcriptomic data from human bone marrow cells, assessing cell-to-cell variability along the trajectories of hematopoietic differentiation into four different types of mature blood cells. The instructive model predicts that cells will be following the same sequence of instructions and that there will be minimal variability of gene expression between them throughout the process, while the stochastic model predicts a role for cell-to-cell variability when lineage commitments are being made.

Results

Applying Shannon entropy to measure cell-to-cell variability among human hematopoietic bone marrow cells at the same stage of differentiation, we observed a transient peak of gene expression variability occurring at characteristic points in all hematopoietic differentiation pathways. Strikingly, the genes whose cell-to-cell variation of expression fluctuated the most over the course of a given differentiation trajectory are pathway-specific genes, whereas genes which showed the greatest variation of mean expression are common to all pathways. Finally, we showed that the level of cell-to-cell variation is increased in the most immature compartment of hematopoiesis in myelodysplastic syndromes.

Conclusions

These data suggest that human hematopoietic differentiation could be better conceptualized as a dynamical stochastic process with a transient stage of cellular indetermination, and strongly support the stochastic view of differentiation. They also highlight the need to consider the role of stochastic gene expression in complex physiological processes and pathologies such as cancers, paving the way for possible noise-based therapies through epigenetic regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01264-9.

Comparison and evaluation of statistical error models for scRNA-seq

BACKGROUND

Heterogeneity in single-cell RNA-seq (scRNA-seq) data is driven by multiple sources, including biological variation in cellular state as well as technical variation introduced during experimental processing. Deconvolving these effects is a key challenge for preprocessing workflows. Recent work has demonstrated the importance and utility of count models for scRNA-seq analysis, but there is a lack of consensus on which statistical distributions and parameter settings are appropriate.

RESULTS

Here, we analyze 59 scRNA-seq datasets that span a wide range of technologies, systems, and sequencing depths in order to evaluate the performance of different error models. We find that while a Poisson error model appears appropriate for sparse datasets, we observe clear evidence of overdispersion for genes with sufficient sequencing depth in all biological systems, necessitating the use of a negative binomial model. Moreover, we find that the degree of overdispersion varies widely across datasets, systems, and gene abundances, and argues for a data-driven approach for parameter estimation.

CONCLUSIONS

Based on these analyses, we provide a set of recommendations for modeling variation in scRNA-seq data, particularly when using generalized linear models or likelihood-based approaches for preprocessing and downstream analysis.

Single-cell Transcriptomes Reveal Characteristics of MicroRNA in Gene Expression Noise Reduction

Isogenic cells growing in identical environments show cell-to-cell variations because of the stochasticity in gene expression. High levels of variation or noise can disrupt robust gene expression and result in tremendous consequences for cell behaviors. In this work, we showed evidence from single-cell RNA sequencing data analysis that microRNAs (miRNAs) can reduce gene expression noise at the mRNA level in mouse cells. We identified that the miRNA expression level, number of targets, target pool abundance, and miRNA–target interaction strength are the key features contributing to noise repression. miRNAs tend to work together in cooperative subnetworks to repress target noise synergistically in a cell type-specific manner. By building a physical model of post-transcriptional regulation and observing in synthetic gene circuits, we demonstrated that accelerated degradation with elevated transcriptional activation of the miRNA target provides resistance to extrinsic fluctuations. Together, through the integrated analysis of single-cell RNA and miRNA expression profiles, we demonstrated that miRNAs are important post-transcriptional regulators for reducing gene expression noise and conferring robustness to biological processes.

Rapid Epigenetic Adaptation in Animals and Its Role in Invasiveness

Invasive species represent a serious ecological threat for many ecosystems worldwide and provide a unique opportunity to investigate rapid adaptation and evolution. Genetic variation allows populations of organisms to be both robust and adaptable to different environmental conditions over evolutionary timeframes. In contrast, invasive animals can rapidly adapt to new environments, with minimal genetic diversity. Thus, the extent to which environmental effects can trigger epigenetic responses is particularly interesting for understanding the role of epigenetics in rapid adaptation. In this review, we provide a brief overview of the different epigenetic mechanisms that control gene expression, and emphasize the importance of epigenetics for environmental adaptation. We also discuss recent publications that provide important examples for the role of epigenetic mechanisms in environmental adaptation. Furthermore, we present an overview of the current knowledge about epigenetic modulation as an adaptive strategy for invasive species. A particularly interesting example is provided by the marbled crayfish, a novel, monoclonal freshwater crayfish species that has colonized diverse habitats within a few years. Finally, we address important limitations of current approaches and highlight the potential importance of less well-known mechanisms for non-genetic organismal adaptation.

Imaging of native transcription and transcriptional dynamics in vivo using a tagged Argonaute protein

A flexible method to image unmodified transcripts and transcription in vivo would be a valuable tool to understand the regulation and dynamics of transcription. Here, we present a novel approach to follow native transcription, with fluorescence microscopy, in live C. elegans. By using the fluorescently tagged Argonaute protein NRDE-3, programmed by exposure to defined dsRNA to bind to nascent transcripts of the gene of interest, we demonstrate transcript labelling of multiple genes, at the transcription site and in the cytoplasm. This flexible approach does not require genetic manipulation, and can be easily scaled up by relying on whole-genome dsRNA libraries. We apply this method to image the transcriptional dynamics of the heat-shock inducible gene hsp-4 (a member of the hsp70 family), as well as two transcription factors: ttx-3 (a LHX2/9 orthologue) in embryos, and hlh-1 (a MyoD orthologue) in larvae, respectively involved in neuronal and muscle development.

Interplay between genetic, epigenetic, and gene expression variability: Considering complexity in evolvability

Genetic variability, epigenetic variability, and gene expression variability (noise) are generally considered independently in their relationship with phenotypic variation. However, they appear to be intrinsically interconnected and influence it in combination. The study of the interplay between genetic and epigenetic variability has the longest history. This article rather considers the introduction of gene expression variability in its relationships with the two others and reviews for the first time experimental evidences over the four relationships connected to gene expression noise. They show how introducing this third source of variability complicates the way of thinking evolvability and the emergence of biological novelty. Finally, cancer cells are proposed to be an ideal model to decipher the dynamic interplay between genetic, epigenetic, and gene expression variability when one of them is either experimentally increased or therapeutically targeted. This interplay is also discussed in an evolutionary perspective in the context of cancer cell drug resistance.

Pluripotent stem cell-derived models of neurological diseases reveal early transcriptional heterogeneity

BACKGROUND

Many neurodegenerative diseases develop only later in life, when cells in the nervous system lose their structure or function. In many forms of neurodegenerative diseases, this late-onset phenomenon remains largely unexplained.

RESULTS

Analyzing single-cell RNA sequencing from Alzheimer's disease (AD) and Huntington's disease (HD) patients, we find increased transcriptional heterogeneity in disease-state neurons. We hypothesize that transcriptional heterogeneity precedes neurodegenerative disease pathologies. To test this idea experimentally, we use juvenile forms (72Q; 180Q) of HD iPSCs, differentiate them into committed neuronal progenitors, and obtain single-cell expression profiles. We show a global increase in gene expression variability in HD. Autophagy genes become more stable, while energy and actin-related genes become more variable in the mutant cells. Knocking down several differentially variable genes results in increased aggregate formation, a pathology associated with HD. We further validate the increased transcriptional heterogeneity in CHD8+/- cells, a model for autism spectrum disorder.

CONCLUSIONS

Overall, our results suggest that although neurodegenerative diseases develop over time, transcriptional regulation imbalance is present already at very early developmental stages. Therefore, an intervention aimed at this early phenotype may be of high diagnostic value.

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue-specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late-stage Drosophila embryos to analyze the properties of promoter types. Using tissue-specific Pol II ChIP-seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC-seq data and have different expression characteristics in single-cell RNA-seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue-specific genes have evolved to use two different strategies for their differential expression across tissues.

Identifying chromatin features that regulate gene expression distribution

Gene expression variability, differences in the number of mRNA per cell across a population of cells, is ubiquitous across diverse organisms with broad impacts on cellular phenotypes. The role of chromatin in regulating average gene expression has been extensively studied. However, what aspects of the chromatin contribute to gene expression variability is still underexplored. Here we addressed this problem by leveraging chromatin diversity and using a systematic investigation of randomly integrated expression reporters to identify what aspects of chromatin microenvironment contribute to gene expression variability. Using DNA barcoding and split-pool decoding, we created a large library of isogenic reporter clones and identified reporter integration sites in a massive and parallel manner. By mapping our measurements of reporter expression at different genomic loci with multiple epigenetic profiles including the enrichment of transcription factors and the distance to different chromatin states, we identified new factors that impact the regulation of gene expression distributions.

Buffering noise: KAT2A modular contributions to stabilization of transcription and cell identity in cancer and development

KAT2A is a histone acetyltransferase recently identified as a vulnerability in at least some forms of Acute Myeloid Leukemia (AML). Its loss or inhibition prompts leukemia stem cells out of self-renewal and into differentiation with ultimate exhaustion of the leukemia pool. We have recently linked the Kat2a requirement in AML to control of transcriptional noise, reflecting an evolutionary-conserved role of Kat2a in promoting burst-like promoter activity and stabilizing gene expression. We suggest that through this role, Kat2a contributes to preservation of cell identity. KAT2A exerts its acetyltransferase activity in the context of two macromolecular complexes, Spt-Ada-Gcn5-Acetyltransferase (SAGA) and Ada-Two-A-Containing (ATAC), but the specific contribution of each complex to stabilization of gene expression is currently unknown. By reviewing specific gene targets and requirements of the two complexes in cancer and development, we suggest that SAGA regulates lineage-specific programs, and ATAC maintains biosynthetic activity through control of ribosomal protein and translation-associated genes, on which cells may be differentially dependent. While our data suggest that KAT2A-mediated regulation of transcriptional noise in AML may be exerted through ATAC, we discuss potential caveats and probe general vs. complex-specific contributions of KAT2A to transcriptional stability, with implications for control and perturbation of cell identity.

Inter-embryo gene expression variability recapitulates the hourglass pattern of evo-devo

BACKGROUND

The evolution of embryological development has long been characterized by deep conservation. In animal development, the phylotypic stage in mid-embryogenesis is more conserved than either early or late stages among species within the same phylum. Hypotheses to explain this hourglass pattern have focused on purifying the selection of gene regulation. Here, we propose an alternative-genes are regulated in different ways at different stages and have different intrinsic capacities to respond to perturbations on gene expression.

RESULTS

To eliminate the influence of natural selection, we quantified the expression variability of isogenetic single embryo transcriptomes throughout fly Drosophila melanogaster embryogenesis. We found that the expression variability is lower at the phylotypic stage, supporting that the underlying regulatory architecture in this stage is more robust to stochastic variation on gene expression. We present evidence that the phylotypic stage is also robust to genetic variations on gene expression. Moreover, chromatin regulation appears to play a key role in the variation and evolution of gene expression.

CONCLUSIONS

We suggest that a phylum-level pattern of embryonic conservation can be explained by the intrinsic difference of gene regulatory mechanisms in different stages.

Predictive features of gene expression variation reveal mechanistic link with differential expression

For most biological processes, organisms must respond to extrinsic cues, while maintaining essential gene expression programmes. Although studied extensively in single cells, it is still unclear how variation is controlled in multicellular organisms. Here, we used a machine-learning approach to identify genomic features that are predictive of genes with high versus low variation in their expression across individuals, using bulk data to remove stochastic cell-to-cell variation. Using embryonic gene expression across 75 Drosophila isogenic lines, we identify features predictive of expression variation (controlling for expression level), many of which are promoter-related. Genes with low variation fall into two classes reflecting different mechanisms to maintain robust expression, while genes with high variation seem to lack both types of stabilizing mechanisms. Applying this framework to humans revealed similar predictive features, indicating that promoter architecture is an ancient mechanism to control expression variation. Remarkably, expression variation features could also partially predict differential expression after diverse perturbations in both Drosophila and humans. Differential gene expression signatures may therefore be partially explained by genetically encoded gene-specific features, unrelated to the studied treatment.

Allele-specific single-cell RNA sequencing reveals different architectures of intrinsic and extrinsic gene expression noises

Gene expression noise refers to the variation of the expression level of a gene among isogenic cells in the same environment, and has two sources: extrinsic noise arising from the disparity of the cell state and intrinsic noise arising from the stochastic process of gene expression in the same cell state. Due to the low throughput of the existing method for measuring the two noise components, the architectures of intrinsic and extrinsic expression noises remain elusive. Using allele-specific single-cell RNA sequencing, we here estimate the two noise components of 3975 genes in mouse fibroblast cells. Our analyses verify predicted influences of several factors such as the TATA-box and microRNA targeting on intrinsic or extrinsic noises and reveal gene function-associated noise trends implicating the action of natural selection. These findings unravel differential regulations, optimizations, and biological consequences of intrinsic and extrinsic noises and can aid the construction of desired synthetic circuits.

Quantitative genetic analysis deciphers the impact of cis and trans regulation on cell-to-cell variability in protein expression levels

Identifying the factors that shape protein expression variability in complex multi-cellular organisms has primarily focused on promoter architecture and regulation of single-cell expression in cis. However, this targeted approach has to date been unable to identify major regulators of cell-to-cell gene expression variability in humans. To address this, we have combined single-cell protein expression measurements in the human immune system using flow cytometry with a quantitative genetics analysis. For the majority of proteins whose variability in expression has a heritable component, we find that genetic variants act in trans, with notably fewer variants acting in cis. Furthermore, we highlight using Mendelian Randomization that these variability-Quantitative Trait Loci might be driven by the cis regulation of upstream genes. This indicates that natural selection may balance the impact of gene regulation in cis with downstream impacts on expression variability in trans.

Loss of Kat2a enhances transcriptional noise and depletes acute myeloid leukemia stem-like cells

Acute Myeloid Leukemia (AML) is an aggressive hematological malignancy with abnormal progenitor self-renewal and defective white blood cell differentiation. Its pathogenesis comprises subversion of transcriptional regulation, through mutation and by hijacking normal chromatin regulation. Kat2a is a histone acetyltransferase central to promoter activity, that we recently associated with stability of pluripotency networks, and identified as a genetic vulnerability in AML. Through combined chromatin profiling and single-cell transcriptomics of a conditional knockout mouse, we demonstrate that Kat2a contributes to leukemia propagation through preservation of leukemia stem-like cells. Kat2a loss impacts transcription factor binding and reduces transcriptional burst frequency in a subset of gene promoters, generating enhanced variability of transcript levels. Destabilization of target programs shifts leukemia cell fate out of self-renewal into differentiation. We propose that control of transcriptional variability is central to leukemia stem-like cell propagation, and establish a paradigm exploitable in different tumors and distinct stages of cancer evolution.

Single-Cell Expression Variability Implies Cell Function

As single-cell RNA sequencing (scRNA-seq) data becomes widely available, cell-to-cell variability in gene expression, or single-cell expression variability (scEV), has been increasingly appreciated. However, it remains unclear whether this variability is functionally important and, if so, what are its implications for multi-cellular organisms. Here, we analyzed multiple scRNA-seq data sets from lymphoblastoid cell lines (LCLs), lung airway epithelial cells (LAECs), and dermal fibroblasts (DFs) and, for each cell type, selected a group of homogenous cells with highly similar expression profiles. We estimated the scEV levels for genes after correcting the mean-variance dependency in that data and identified 465, 466, and 364 highly variable genes (HVGs) in LCLs, LAECs, and DFs, respectively. Functions of these HVGs were found to be enriched with those biological processes precisely relevant to the corresponding cell type’s function, from which the scRNA-seq data used to identify HVGs were generated—e.g., cytokine signaling pathways were enriched in HVGs identified in LCLs, collagen formation in LAECs, and keratinization in DFs. We repeated the same analysis with scRNA-seq data from induced pluripotent stem cells (iPSCs) and identified only 79 HVGs with no statistically significant enriched functions; the overall scEV in iPSCs was of negligible magnitude. Our results support the “variation is function” hypothesis, arguing that scEV is required for cell type-specific, higher-level system function. Thus, quantifying and characterizing scEV are of importance for our understating of normal and pathological cellular processes.

Multivariable regulation of gene expression plasticity in metazoans

An important capacity of genes is the rapid change of expression levels to cope with the environment, known as expression responsiveness or plasticity. Elucidating the genomic mechanisms determining expression plasticity is critical for understanding the molecular basis of phenotypic plasticity, fitness and adaptation. In this study, we systematically quantified gene expression plasticity in four metazoan species by integrating changes of expression levels under a large number of genetic and environmental conditions. From this, we demonstrated that expression plasticity measures a distinct feature of gene expression that is orthogonal to other well-studied features, including gene expression level and tissue specificity/broadness. Expression plasticity is conserved across species with important physiological implications. The magnitude of expression plasticity is highly correlated with gene function and genes with high plasticity are implicated in disease susceptibility. Genome-wide analysis identified many conserved promoter cis-elements, trans-acting factors (such as CTCF), and gene body histone modifications (H3K36me3, H3K79me2 and H4K20me1) that are significantly associated with expression plasticity. Analysis of expression changes in perturbation experiments further validated a causal role of specific transcription factors and histone modifications. Collectively, this work reveals the general properties, physiological implications and multivariable regulation of gene expression plasticity in metazoans, extending the mechanistic understanding of gene regulation.

Challenges in measuring and understanding biological noise

Biochemical reactions are intrinsically stochastic, leading to variation in the production of mRNAs and proteins within cells. In the scientific literature, this source of variation is typically referred to as 'noise'. The observed variability in molecular phenotypes arises from a combination of processes that amplify and attenuate noise. Our ability to quantify cell-to-cell variability in numerous biological contexts has been revolutionized by recent advances in single-cell technology, from imaging approaches through to 'omics' strategies. However, defining, accurately measuring and disentangling the stochastic and deterministic components of cell-to-cell variability is challenging. In this Review, we discuss the sources, impact and function of molecular phenotypic variability and highlight future directions to understand its role.

Identification of genes under dynamic post-transcriptional regulation from time-series epigenomic data

Aim: Prediction of genes under dynamic post-transcriptional regulation from epigenomic data. Materials & methods: We used time-series profiles of chromatin immunoprecipitation-seq data of histone modifications from differentiation of mesenchymal progenitor cells toward adipocytes and osteoblasts to predict gene expression levels at five time points in both lineages and estimated the deviation of those predictions from the RNA-seq measured expression levels using linear regression. Results & conclusion: The genes with biggest changes in their estimated stability across the time series are enriched for noncoding RNAs and lineage-specific biological processes. Clustering mRNAs according to their stability dynamics allows identification of post-transcriptionally coregulated mRNAs and their shared regulators through sequence enrichment analysis. We identify miR-204 as an early induced adipogenic microRNA targeting Akr1c14 and Il1rl1.

Gene expression variability: the other dimension in transcriptome analysis

Transcriptome sequencing is a powerful technique to study molecular changes that underlie the differences in physiological conditions and disease progression. A typical question that is posed in such studies is finding genes with significant changes between sample groups. In this respect expression variability is regarded as a nuisance factor that is primarily of technical origin and complicates the data analysis. However, it is becoming apparent that the biological variation in gene expression might be an important molecular phenotype that can affect physiological parameters. In this review we explore the recent literature on technical and biological variability in gene expression, sources of expression variability, (epi-)genetic hallmarks, and evolutionary constraints in genes with robust and variable gene expression. We provide an overview of recent findings on effects of external cues, such as diet and aging, on expression variability and on other biological phenomena that can be linked to it. We discuss metrics and tools that were developed for quantification of expression variability and highlight the importance of future studies in this direction. To assist the adoption of expression variability analysis, we also provide a detailed description and computer code, which can easily be utilized by other researchers. We also provide a reanalysis of recently published data to highlight the value of the analysis method.

Proteotyping pluripotency with mass spectrometry

INTRODUCTION

Pluripotency emerges transiently during embryogenesis in two main forms with different developmental potential, termed naïve and primed states. Importantly, these pluripotent states can be recapitulated in vitro under specific culture conditions, representing a unique model to study the regulatory principles of development and cellular plasticity. Areas covered: A complex network of signaling pathways that senses intrinsic and extrinsic cues controls the fine balance between self-renewal and differentiation. Much of our knowledge on this tight regulation originates from epigenetic and transcriptomic approaches. However, the presence of post-transcriptional and post-translational mechanisms demands a direct assessment of the proteome in its multiple facets. Mass spectrometry-based proteomics is now a mature technique and has started to deliver new insights in the stem cell field. Expert opinion: Here, we review our current understanding on the mechanisms that dictate the spectrum of pluripotency levels. We put special emphasis on the emerging proteomic studies that focused on the molecular properties behind the naïve and primed states. In addition, we hypothesize on the impact that future developments in proteomic technologies can have to improve our view of pluripotency.

Ribosome Stoichiometry: From Form to Function

The existence of eukaryotic ribosomes with distinct ribosomal protein (RP) stoichiometry and regulatory roles in protein synthesis has been speculated for over 60 years. Recent advances in mass spectrometry (MS) and high-throughput analysis have begun to identify and characterize distinct ribosome stoichiometry in yeast and mammalian systems. In addition to RP stoichiometry, ribosomes host a vast array of protein modifications, effectively expanding the number of human RPs from 80 to many thousands of distinct proteoforms. Is it possible that these proteoforms combine to function as a 'ribosome code' to tune protein synthesis? We outline the specific benefits that translational regulation by specialized ribosomes can offer and discuss the means and methodologies available to correlate and characterize RP stoichiometry with function. We highlight previous research with a focus on formulating hypotheses that can guide future experiments and crack the ribosome code.

Gene expression variability across cells and species shapes innate immunity

As the first line of defence against pathogens, cells mount an innate immune response, which varies widely from cell to cell. The response must be potent but carefully controlled to avoid self-damage. How these constraints have shaped the evolution of innate immunity remains poorly understood. Here we characterize the innate immune response's transcriptional divergence between species and variability in expression among cells. Using bulk and single-cell transcriptomics in fibroblasts and mononuclear phagocytes from different species, challenged with immune stimuli, we map the architecture of the innate immune response. Transcriptionally diverging genes, including those that encode cytokines and chemokines, vary across cells and have distinct promoter structures. Conversely, genes that are involved in the regulation of this response, such as those that encode transcription factors and kinases, are conserved between species and display low cell-to-cell variability in expression. We suggest that this expression pattern, which is observed across species and conditions, has evolved as a mechanism for fine-tuned regulation to achieve an effective but balanced response.

The methylome of the marbled crayfish links gene body methylation to stable expression of poorly accessible genes

Background

The parthenogenetic marbled crayfish (Procambarus virginalis) is a novel species that has rapidly invaded and colonized various different habitats. Adaptation to different environments appears to be independent of the selection of genetic variants, but epigenetic programming of the marbled crayfish genome remains to be understood.

Results

Here, we provide a comprehensive analysis of DNA methylation in marbled crayfish. Whole-genome bisulfite sequencing of multiple replicates and different tissues revealed a methylation pattern that is characterized by gene body methylation of housekeeping genes. Interestingly, this pattern was largely tissue invariant, suggesting a function that is unrelated to cell fate specification. Indeed, integrative analysis of DNA methylation, chromatin accessibility and mRNA expression patterns revealed that gene body methylation correlated with limited chromatin accessibility and stable gene expression, while low-methylated genes often resided in chromatin with higher accessibility and showed increased expression variation. Interestingly, marbled crayfish also showed reduced gene body methylation and higher gene expression variability when compared with their noninvasive mother species, Procambarus fallax.

Conclusions

Our results provide novel insights into invertebrate gene body methylation and its potential role in adaptive gene regulation.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0229-6) contains supplementary material, which is available to authorized users.

Epigenetic Variability Confounds Transcriptome but Not Proteome Profiling for Coexpression-based Gene Function Prediction

Genes are often coexpressed with their genomic neighbors, even if these are functionally unrelated. For small expression changes driven by genetic variation within the same cell type, non-functional mRNA coexpression is not propagated to the protein level. However, it is unclear if protein levels are also buffered against any non-functional mRNA coexpression accompanying large, regulated changes in the gene expression program, such as those occurring during cell differentiation. Here, we address this question by analyzing mRNA and protein expression changes for housekeeping genes across 20 mouse tissues. We find that a large proportion of mRNA coexpression is indeed non-functional and does not lead to coexpressed proteins. Chromosomal proximity of genes explains a proportion of this nonfunctional mRNA coexpression. However, the main driver of non-functional mRNA coexpression across mouse tissues is epigenetic similarity. Both factors together provide an explanation for why monitoring protein coexpression outperforms mRNA coexpression data in gene function prediction. Furthermore, this suggests that housekeeping genes translocating during evolution within genomic subcompartments might maintain their broad expression pattern.

CpG island composition differences are a source of gene expression noise indicative of promoter responsiveness

BACKGROUND

Population phenotypic variation can arise from genetic differences between individuals, or from cellular heterogeneity in an isogenic group of cells or organisms. The emergence of gene expression differences between genetically identical cells is referred to as gene expression noise, the sources of which are not well understood.

RESULTS

In this work, by studying gene expression noise between multiple cell lineages and mammalian species, we find consistent evidence of a role for CpG islands as sources of gene expression noise. Variation in noise among CpG island promoters can be partially attributed to differences in island size, in which short islands have noisier gene expression. Building on these findings, we investigate the potential for short CpG islands to act as fast response elements to environmental stimuli. Specifically, we find that these islands are enriched amongst primary response genes in SWI/SNF-independent stimuli, suggesting that expression noise is an indicator of promoter responsiveness.

CONCLUSIONS

Thus, through the integration of single-cell RNA expression profiling, chromatin landscape and temporal gene expression dynamics, we have uncovered a role for short CpG island promoters as fast response elements.