Mammalian gene expression variability is explained by underlying cell state

Cell-cycle dependence of bursty gene expression: insights from fitting mechanistic models to single-cell RNA-seq data

Bursty gene expression is characterized by two intuitive parameters, burst frequency and burst size, the cell-cycle dependence of which has not been extensively profiled at the transcriptome level. In this study, we estimate the burst parameters per allele in the G1 and G2/M cell-cycle phases for thousands of mouse genes by fitting mechanistic models of gene expression to messenger RNA count data, obtained by sequencing of single cells whose cell-cycle position has been inferred using a deep-learning method. We find that upon DNA replication, the median burst frequency approximately halves, while the burst size remains mostly unchanged. Genome-wide distributions of the burst parameter ratios between the G2/M and G1 phases are broad, indicating substantial heterogeneity in transcriptional regulation. We also observe a significant negative correlation between the burst frequency and size ratios, suggesting that regulatory processes do not independently control the burst parameters. We show that to accurately estimate the burst parameter ratios, mechanistic models must explicitly account for gene copy number variation and extrinsic noise due to the coupling of transcription to cell age across the cell cycle, but corrections for technical noise due to imperfect capture of RNA molecules in sequencing experiments are less critical.

Gradient matching accelerates mixed-effects inference for biochemical networks

MOTIVATION

Single-cell time series data often exhibit significant variability within an isogenic cell population. When modeling intracellular processes, it is therefore more appropriate to infer parameter distributions that reflect this variability, rather than fitting the population average to obtain a single point estimate. The Global Two-Stage (GTS) approach for nonlinear mixed-effects (NLME) models is a simple and modular method commonly used for this purpose. However, this method is computationally intensive due to its repeated use of nonconvex optimization and numerical integration of the underlying system.

RESULTS

We propose the Gradient Matching GTS (GMGTS) method as an efficient alternative to GTS. Gradient matching offers an integration-free approach to parameter estimation that is particularly powerful for systems that are linear in the unknown parameters, such as biochemical networks modeled by mass action kinetics. By incorporating gradient matching into the GTS framework, we expand its capabilities through uncertainty propagation calculations and an iterative estimation scheme for partially observed systems. Comparisons between GMGTS and GTS across various inference setups show that our method significantly reduces computational demands, facilitating the application of complex NLME models in systems biology.

AVAILABILITY AND IMPLEMENTATION

A Matlab implementation of GMGTS is provided at https://github.com/yulanvanoppen/GMGTS (DOI: http://doi.org/10.5281/zenodo.14884457).

Transient power-law behaviour following induction distinguishes between competing models of stochastic gene expression

What features of transcription can be learnt by fitting mathematical models of gene expression to mRNA count data? Given a suite of models, fitting to data selects an optimal one, thus identifying a probable transcriptional mechanism. Whilst attractive, the utility of this methodology remains unclear. Here, we sample steady-state, single-cell mRNA count distributions from parameters in the physiological range, and show they cannot be used to confidently estimate the number of inactive gene states, i.e. the number of rate-limiting steps in transcriptional initiation. Distributions from over 99% of the parameter space generated using models with 2, 3, or 4 inactive states can be well fit by one with a single inactive state. However, we show that for many minutes following induction, eukaryotic cells show an increase in the mean mRNA count that obeys a power law whose exponent equals the sum of the number of states visited from the initial inactive to the active state and the number of rate-limiting post-transcriptional processing steps. Our study shows that estimation of the exponent from eukaryotic data can be sufficient to determine a lower bound on the total number of regulatory steps in transcription initiation, splicing, and nuclear export.

Dynamic and Biphasic Regulation of Cell Migration by Ras

Ras has traditionally been regarded as a positive regulator and therapeutic target due to its role in cell proliferation, but recent findings indicate a more nuanced role in cell migration, where suppressed Ras activity can unexpectedly promote migration. To clarify this complexity, we systematically modulate Ras activity using various RasGEF and RasGAP proteins and assess their effects on migration dynamics. Leveraging optogenetics, we assess the immediate, non-transcriptional effects of Ras signaling on migration. Local RasGEF recruitment to the plasma membrane induces protrusions and new fronts to effectively guide migration, even in the absence of GPCR/G-protein signaling whereas global recruitment causes immediate cell spreading halting cell migration. Local RasGAP recruitment suppresses protrusions, generates new backs, and repels cells whereas global relocation either eliminates all protrusions to inhibit migration or preserves a single protrusion to maintain polarity. Consistent local and global increases or decreases in signal transduction and cytoskeletal activities accompany these morphological changes. Additionally, we performed cortical tension measurements and found that RasGEFs generally increase cortical tension while RasGAPs decrease it. Our results reveal a biphasic relationship between Ras activity and cellular dynamics, reinforcing our previous findings that optimal Ras activity and cortical tension are critical for efficient migration.

Significance

This study challenges the traditional view of Ras as solely a positive regulator of cell functions by controlling of gene expression. Using optogenetics to rapidly modulate Ras activity in Dictyostelium , we demonstrate a biphasic relationship between Ras activity and migration: both excessive and insufficient Ras activity impair cell movement. Importantly, these effects occur rapidly, independent of transcriptional changes, revealing the mechanism by which Ras controls cell migration. The findings suggest that optimal Ras activity and cortical tension are crucial for efficient migration, and that targeting Ras in cancer therapy should consider the cell's initial state, aiming to push Ras activity outside the optimal range for migration. This nuanced understanding of the role of Ras in migration has significant implications for developing more effective cancer treatments, as simply inhibiting Ras might inadvertently promote metastasis in certain contexts.

Priestia megaterium cells are primed for surviving lethal doses of antibiotics and chemical stress

Antibiotic resistant infections kill millions worldwide yearly. However, a key factor in recurrent infections is antibiotic persisters. Persisters are not inherently antibiotic-resistant but can withstand antibiotic exposure by entering a non-dividing state. This tolerance often results in prolonged antibiotic usage, increasing the likelihood of developing resistant strains. Here, we show the existence of "primed cells" in the Gram-positive bacterium Priestia megaterium, formerly known as Bacillus megaterium. These cells are pre-adapted to become persisters prior to lethal antibiotic stress. Remarkably, this prepared state is passed down through multiple generations via epigenetic memory, enhancing survival against antibiotics and other chemical stress. Previously, two distinct types of persisters were proposed: Type I and Type II, formed during stationary and log phases, respectively. However, our findings reveal that primed cells contribute to an increase in persisters during transition and stationary phases, with no evidence supporting distinct phenotypes between Type I and Type II persisters.

Normalization of gene counts affects principal components-based exploratory analysis of RNA-sequencing data

Normalization of gene expression count data is an essential step of in the analysis of RNA-sequencing data. Its statistical analysis has been mostly addressed in the context of differential expression analysis, that is in the univariate setting. However, relationships among genes and samples are better explored and quantified using multivariate exploratory data analysis tools like Principal Component Analysis (PCA). In this study we investigate how normalization impacts the PCA model and its interpretation, considering twelve different widely used normalization methods that were applied on simulated and experimental data. Correlation patterns in the normalized data were explored using both summary statistics and Covariance Simultaneous Component Analysis. The impact of normalization on the PCA solution was assessed by exploring the model complexity, the quality of sample clustering in the low-dimensional PCA space and gene ranking in the model fit to normalized data. PCA models upon normalization were interpreted in the context gene enrichment pathway analysis. We found that although PCA score plots are often similar independently form the normalization used, biological interpretation of the models can depend heavily on the normalization method applied.

Robustness and variability in Caenorhabditis elegans dauer gene expression

Both plasticity and robustness are pervasive features of developmental programs. The dauer in Caenorhabditis elegans is an arrested, hypometabolic alternative to the third larval stage of the nematode. Dauers undergo dramatic tissue remodeling and extensive physiological, metabolic, behavioral, and gene expression changes compared to conspecifics that continue development and can be induced by several adverse environments or genetic mutations that act as independent and parallel inputs into the larval developmental program. Therefore, dauer induction is an example of phenotypic plasticity. However, whether gene expression in dauer larvae induced to arrest development by different genetic or environmental triggers is invariant or varies depending on their route into dauer has not been examined. By using RNA-sequencing to characterize gene expression in different types of dauer larvae and computing the variance and concordance within Gene Ontologies (GO) and gene expression networks, we find that the expression patterns within most pathways are strongly correlated between dauer larvae, suggestive of transcriptional robustness. However, gene expression within specific defense pathways, pathways regulating some morphological traits, and several metabolic pathways differ between the dauer larvae. We speculate that the transcriptional robustness of core dauer pathways allows for the buffering of variation in the expression of genes involved in adaptation, allowing the dauers induced by different stimuli to survive in and exploit different niches.

Evaluation of Integrity of Allogeneic Bone Processed with High Hydrostatic Pressure: A Pilot Animal Study

Processing of bone allografts with strong acids and γ-sterilization results in decreased biomechanical properties and reduction in osteogenecity and osteoconductivity. High hydrostatic pressure (HHP) treatment could be a gentle alternative to processing techniques usually applied. HHP is known to induce devitalization of cancellous bone while preserving biomechanical stability and molecules that induce cell differentiation. Here, a specific HHP protocol for devitalization of cancellous bone was applied to rabbit femoral bone. Allogeneic bone cylinders were subsequently implanted into a defect in the lateral condyles of rabbit femora and were compared to autologous bone grafts. Analysis of bone integration 4 and 12 weeks postoperatively revealed no differences between autografts and HHP-treated allografts regarding the expression of genes characteristic for bone remodeling, showing expression niveous comparable to original bone cylinder. Furthermore, biomechanical properties were evaluated 12 weeks postoperatively. Autografts and HHP-treated allografts both showed a yield strength ranging between 2 and 2.5 MPa and an average bone mass density of 250 mg/cm2. Furthermore, histological analysis of the region of interest revealed a rate of 5 to 10% BPM-2 and approximately 40% osteocalcin-positive staining, with no marked differences between allografts and autografts demonstrating comparable matrix deposition in the graft region. A suitable graft integrity was pointed out by μCT imaging in both groups, supporting the biomechanical data. In summary, the integrity of HHP-treated cancellous bone allografts showed similar results to untreated autografts. Hence, HHP treatment may represent a gentle and effective alternative to existing processing techniques for bone allografts.

Multiplexed spatial transcriptomics methods and the application of expansion microscopy

While spatial transcriptomics has undeniably revolutionized our ability to study cellular organization, it has driven the development of a great number of innovative transcriptomics methods, which can be classified into in situ sequencing (ISS) methods, in situ hybridization (ISH) techniques, and next-generation sequencing (NGS)-based sequencing with region capture. These technologies not only refine our understanding of cellular processes, but also open up new possibilities for breakthroughs in various research domains. One challenge of spatial transcriptomics experiments is the limitation of RNA detection due to optical crowding of RNA in the cells. Expansion microscopy (ExM), characterized by the controlled enlargement of biological specimens, offers a means to achieve super-resolution imaging, overcoming the diffraction limit inherent in conventional microscopy and enabling precise visualization of RNA in spatial transcriptomics methods. In this review, we elaborate on ISS, ISH and NGS-based spatial transcriptomic protocols and on how performance of these techniques can be extended by the combination of these protocols with ExM. Moving beyond the techniques and procedures, we highlight the broader implications of transcriptomics in biology and medicine. These include valuable insight into the spatial organization of gene expression in cells within tissues, aid in the identification and the distinction of cell types and subpopulations and understanding of molecular mechanisms and intercellular changes driving disease development.

Analysis of a detailed multi-stage model of stochastic gene expression using queueing theory and model reduction

We introduce a biologically detailed, stochastic model of gene expression describing the multiple rate-limiting steps of transcription, nuclear pre-mRNA processing, nuclear mRNA export, cytoplasmic mRNA degradation and translation of mRNA into protein. The processes in sub-cellular compartments are described by an arbitrary number of processing stages, thus accounting for a significantly finer molecular description of gene expression than conventional models such as the telegraph, two-stage and three-stage models of gene expression. We use two distinct tools, queueing theory and model reduction using the slow-scale linear-noise approximation, to derive exact or approximate analytic expressions for the moments or distributions of nuclear mRNA, cytoplasmic mRNA and protein fluctuations, as well as lower bounds for their Fano factors in steady-state conditions. We use these to study the phase diagram of the stochastic model; in particular we derive parametric conditions determining three types of transitions in the properties of mRNA fluctuations: from sub-Poissonian to super-Poissonian noise, from high noise in the nucleus to high noise in the cytoplasm, and from a monotonic increase to a monotonic decrease of the Fano factor with the number of processing stages. In contrast, protein fluctuations are always super-Poissonian and show weak dependence on the number of mRNA processing stages. Our results delineate the region of parameter space where conventional models give qualitatively incorrect results and provide insight into how the number of processing stages, e.g. the number of rate-limiting steps in initiation, splicing and mRNA degradation, shape stochastic gene expression by modulation of molecular memory.

Gram-positive bacteria are primed for surviving lethal doses of antibiotics and chemical stress

Antibiotic resistance kills millions worldwide yearly. However, a major contributor to recurrent infections lies in a small fraction of bacterial cells, known as persisters. These cells are not inherently antibiotic-resistant, yet they lead to increased antibiotic usage, raising the risk of developing resistant progenies. In a bacterial population, individual cells exhibit considerable fluctuations in their gene expression levels despite being cultivated under identical, stable conditions. This variability in cell-to-cell characteristics (phenotypic diversity) within an isogenic population enables persister cells to withstand antibiotic exposure by entering a non-dividing state. We recently showed the existence of "primed cells" in E. coli. Primed cells are dividing cells prepared for antibiotic stress before encountering it and are more prone to form persisters. They also pass their "prepared state" down for several generations through epigenetic memory. Here, we show that primed cells are common among distant bacterial lineages, allowing for survival against antibiotics and other chemical stress, and form in different growth phases. They are also responsible for increased persister levels in transition and stationary phases compared to the log phase. We tested and showed that the Gram-positive bacterium Bacillus megaterium, evolutionarily very distant from E. coli, forms primed cells and has a transient epigenetic memory that is maintained for 7 generations or more. We showed this using ciprofloxacin and the non-antibiotic chemical stress fluoride. It is well established that persister levels are higher in the stationary phase than in the log phase, and B. megaterium persisters levels are nearly identical from the early to late-log phase but are ~2-fold and ~4-fold higher in the transition and stationary phase, respectively. It was previously proposed that there are two distinct types of persisters: Type II forms in the log phase, while Type I forms in the stationary phase. However, we show that primed cells lead to increased persisters in the transition and stationary phase and found no evidence of Type I or II persisters with distant phenotypes. Overall, we have provided substantial evidence of the importance of primed cells and their transitory epigenetic memories to surviving stress.

Complementing Cell Taxonomies with a Multicellular Analysis of Tissues

The application of single-cell molecular profiling coupled with spatial technologies has enabled charting of cellular heterogeneity in reference tissues and in disease. This new wave of molecular data has highlighted the expected diversity of single-cell dynamics upon shared external queues and spatial organizations. However, little is known about the relationship between single-cell heterogeneity and the emergence and maintenance of robust multicellular processes in developed tissues and its role in (patho)physiology. Here, we present emerging computational modeling strategies that use increasingly available large-scale cross-condition single-cell and spatial datasets to study multicellular organization in tissues and complement cell taxonomies. This perspective should enable us to better understand how cells within tissues collectively process information and adapt synchronized responses in disease contexts and to bridge the gap between structural changes and functions in tissues.

Effects of mRNA conformational switching on translational noise in gene circuits

Intragenic translational heterogeneity describes the variation in translation at the level of transcripts for an individual gene. A factor that contributes to this source of variation is the mRNA structure. Both the composition of the thermodynamic ensemble, i.e., the stationary distribution of mRNA structures, and the switching dynamics between those play a role. The effect of the switching dynamics on intragenic translational heterogeneity remains poorly understood. We present a stochastic translation model that accounts for mRNA structure switching and is derived from a Markov model via approximate stochastic filtering. We assess the approximation on various timescales and provide a method to quantify how mRNA structure dynamics contributes to translational heterogeneity. With our approach, we allow quantitative information on mRNA switching from biophysical experiments or coarse-grain molecular dynamics simulations of mRNA structures to be included in gene regulatory chemical reaction network models without an increase in the number of species. Thereby, our model bridges a gap between mRNA structure kinetics and gene expression models, which we hope will further improve our understanding of gene regulatory networks and facilitate genetic circuit design.

Resource allocation in mammalian systems

Cells execute biological functions to support phenotypes such as growth, migration, and secretion. Complementarily, each function of a cell has resource costs that constrain phenotype. Resource allocation by a cell allows it to manage these costs and optimize their phenotypes. In fact, the management of resource constraints (e.g., nutrient availability, bioenergetic capacity, and macromolecular machinery production) shape activity and ultimately impact phenotype. In mammalian systems, quantification of resource allocation provides important insights into higher-order multicellular functions; it shapes intercellular interactions and relays environmental cues for tissues to coordinate individual cells to overcome resource constraints and achieve population-level behavior. Furthermore, these constraints, objectives, and phenotypes are context-dependent, with cells adapting their behavior according to their microenvironment, resulting in distinct steady-states. This review will highlight the biological insights gained from probing resource allocation in mammalian cells and tissues.

Integrating single-cell transcriptomics with cellular phenotypes: cell morphology, Ca2+ imaging and electrophysiology

I review recent technological advancements in coupling single-cell transcriptomics with cellular phenotypes including morphology, calcium signaling, and electrophysiology. Single-cell RNA sequencing (scRNAseq) has revolutionized cell type classifications by capturing the transcriptional diversity of cells. A new wave of methods to integrate scRNAseq and biophysical measurements is facilitating the linkage of transcriptomic data to cellular function, which provides physiological insight into cellular states. I briefly discuss critical factors of these phenotypical characterizations such as timescales, information content, and analytical tools. Dedicated sections focus on the integration with cell morphology, calcium imaging, and electrophysiology (patch-seq), emphasizing their complementary roles. I discuss their application in elucidating cellular states, refining cell type classifications, and uncovering functional differences in cell subtypes. To illustrate the practical applications and benefits of these methods, I highlight their use in tissues with excitable cell-types such as the brain, pancreatic islets, and the retina. The potential of combining functional phenotyping with spatial transcriptomics for a detailed mapping of cell phenotypes in situ is explored. Finally, I discuss open questions and future perspectives, emphasizing the need for a shift towards broader accessibility through increased throughput.

Quantifying and correcting bias in transcriptional parameter inference from single-cell data

The snapshot distribution of mRNA counts per cell can be measured using single-molecule fluorescence in situ hybridization or single-cell RNA sequencing. These distributions are often fit to the steady-state distribution of the two-state telegraph model to estimate the three transcriptional parameters for a gene of interest: mRNA synthesis rate, the switching on rate (the on state being the active transcriptional state), and the switching off rate. This model assumes no extrinsic noise, i.e., parameters do not vary between cells, and thus estimated parameters are to be understood as approximating the average values in a population. The accuracy of this approximation is currently unclear. Here, we develop a theory that explains the size and sign of estimation bias when inferring parameters from single-cell data using the standard telegraph model. We find specific bias signatures depending on the source of extrinsic noise (which parameter is most variable across cells) and the mode of transcriptional activity. If gene expression is not bursty then the population averages of all three parameters are overestimated if extrinsic noise is in the synthesis rate; underestimation occurs if extrinsic noise is in the switching on rate; both underestimation and overestimation can occur if extrinsic noise is in the switching off rate. We find that some estimated parameters tend to infinity as the size of extrinsic noise approaches a critical threshold. In contrast when gene expression is bursty, we find that in all cases the mean burst size (ratio of the synthesis rate to the switching off rate) is overestimated while the mean burst frequency (the switching on rate) is underestimated. We estimate the size of extrinsic noise from the covariance matrix of sequencing data and use this together with our theory to correct published estimates of transcriptional parameters for mammalian genes.

A guide to ERK dynamics, part 2: downstream decoding

Signaling by the extracellular signal-regulated kinase (ERK) pathway controls many cellular processes, including cell division, death, and differentiation. In this second installment of a two-part review, we address the question of how the ERK pathway exerts distinct and context-specific effects on multiple processes. We discuss how the dynamics of ERK activity induce selective changes in gene expression programs, with insights from both experiments and computational models. With a focus on single-cell biosensor-based studies, we summarize four major functional modes for ERK signaling in tissues: adjusting the size of cell populations, gradient-based patterning, wave propagation of morphological changes, and diversification of cellular gene expression states. These modes of operation are disrupted in cancer and other related diseases and represent potential targets for therapeutic intervention. By understanding the dynamic mechanisms involved in ERK signaling, there is potential for pharmacological strategies that not only simply inhibit ERK, but also restore functional activity patterns and improve disease outcomes.

A guide to ERK dynamics, part 1: mechanisms and models

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

Gene transcription regulation by ER at the single cell and allele level

In this short review we discuss the current view of how the estrogen receptor (ER), a pivotal member of the nuclear receptor superfamily of transcription factors, regulates gene transcription at the single cell and allele level, focusing on in vitro cell line models. We discuss central topics and new trends in molecular biology including phenotypic heterogeneity, single cell sequencing, nuclear phase separated condensates, single cell imaging, and image analysis methods, with particular focus on the methodologies and results that have been reported in the last few years using microscopy-based techniques. These observations augment the results from biochemical assays that lead to a much more complex and dynamic view of how ER, and arguably most transcription factors, act to regulate gene transcription.

Light inducible protein degradation in E. coli with the LOVdeg tag

Molecular tools for optogenetic control allow for spatial and temporal regulation of cell behavior. In particular, light controlled protein degradation is a valuable mechanism of regulation because it can be highly modular, used in tandem with other control mechanisms, and maintain functionality throughout growth phases. Here, we engineered LOVdeg, a tag that can be appended to a protein of interest for inducible degradation in Escherichia coli using blue light. We demonstrate the modularity of LOVdeg by using it to tag a range of proteins, including the LacI repressor, CRISPRa activator, and the AcrB efflux pump. Additionally, we demonstrate the utility of pairing the LOVdeg tag with existing optogenetic tools to enhance performance by developing a combined EL222 and LOVdeg system. Finally, we use the LOVdeg tag in a metabolic engineering application to demonstrate post-translational control of metabolism. Together, our results highlight the modularity and functionality of the LOVdeg tag system, and introduce a powerful new tool for bacterial optogenetics.

Quantitative determination of fluorescence labeling implemented in cell cultures

BACKGROUND

Labeling efficiency is a crucial parameter in fluorescence applications, especially when studying biomolecular interactions. Current approaches for estimating the yield of fluorescent labeling have critical drawbacks that usually lead them to be inaccurate or not quantitative.

RESULTS

We present a method to quantify fluorescent-labeling efficiency that addresses the critical issues marring existing approaches. The method operates in the same conditions of the target experiments by exploiting a ratiometric evaluation with two fluorophores used in sequential reactions. We show the ability of the protocol to extract reliable quantification for different fluorescent probes, reagents concentrations, and reaction timing and to optimize labeling performance. As paradigm, we consider the labeling of the membrane-receptor TrkA through 4'-phosphopantetheinyl transferase Sfp in living cells, visualizing the results by TIRF microscopy. This investigation allows us to find conditions for demanding single and multi-color single-molecule studies requiring high degrees of labeling.

CONCLUSIONS

The developed method allows the quantitative determination and the optimization of staining efficiency in any labeling strategy based on stable reactions.

Automatic design of gene regulatory mechanisms for spatial pattern formation

Synthetic developmental biology aims to engineer gene regulatory mechanisms (GRMs) for understanding and producing desired multicellular patterns and shapes. However, designing GRMs for spatial patterns is a current challenge due to the nonlinear interactions and feedback loops in genetic circuits. Here we present a methodology to automatically design GRMs that can produce any given spatial pattern. The proposed approach uses two orthogonal morphogen gradients acting as positional information signals in a multicellular tissue area or culture, which constitutes a continuous field of engineered cells implementing the same designed GRM. To efficiently design both the circuit network and the interaction mechanisms-including the number of genes necessary for the formation of the target pattern-we developed an automated algorithm based on high-performance evolutionary computation. The tolerance of the algorithm can be configured to design GRMs that are either simple to produce approximate patterns or complex to produce precise patterns. We demonstrate the approach by automatically designing GRMs that can produce a diverse set of synthetic spatial expression patterns by interpreting just two orthogonal morphogen gradients. The proposed framework offers a versatile approach to systematically design and discover pattern-producing genetic circuits.

Evaluating single-cell variability in proteasomal decay

Gene expression is a stochastic process that leads to variability in mRNA and protein abundances even within an isogenic population of cells grown in the same environment. This variation, often called gene-expression noise, has typically been attributed to transcriptional and translational processes while ignoring the contributions of protein decay variability across cells. Here we estimate the single-cell protein decay rates of two degron GFPs in Saccharomyces cerevisiae using time-lapse microscopy. We find substantial cell-to-cell variability in the decay rates of the degron GFPs. We evaluate cellular features that explain the variability in the proteasomal decay and find that the amount of 20s catalytic beta subunit of the proteasome marginally explains the observed variability in the degron GFP half-lives. We propose alternate hypotheses that might explain the observed variability in the decay of the two degron GFPs. Overall, our study highlights the importance of studying the kinetics of the decay process at single-cell resolution and that decay rates vary at the single-cell level, and that the decay process is stochastic. A complex model of decay dynamics must be included when modeling stochastic gene expression to estimate gene expression noise.

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.

SLAM-seq reveals early transcriptomic response mechanisms upon glutamine deprivation in Chinese hamster ovary cells

Mammalian cells frequently encounter subtle perturbations during recombinant protein production. Identifying the genetic factors that govern the cellular stress response can facilitate targeted genetic engineering to obtain production cell lines that demonstrate a higher stress tolerance. To simulate nutrient stress, Chinese hamster ovary (CHO) cells were transferred into a glutamine(Q)-free medium and transcriptional dynamics using thiol(SH)-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) along with standard RNA-seq of stressed and unstressed cells were investigated. The SLAM-seq method allows differentiation between actively transcribed, nascent mRNA, and total (previously present) mRNA in the sample, adding an additional, time-resolved layer to classic RNA-sequencing. The cells tackle amino acid (AA) limitation by inducing the integrated stress response (ISR) signaling pathway, reflected in Atf4 overexpression in the early hours post Q deprivation, leading to subsequent activation of its targets, Asns, Atf3, Ddit3, Eif4ebp1, Gpt2, Herpud1, Slc7a1, Slc7a11, Slc38a2, Trib3, and Vegfa. The GCN2-eIF2α-ATF4 pathway is confirmed by a significant halt in transcription of translation-related genes at 24 h post Q deprivation. The downregulation of lipid synthesis indicates the inhibition of the mTOR pathway, further confirmed by overexpression of Sesn2. Furthermore, SLAM-seq detects short-lived transcription factors, such as Egr1, that would have been missed in standard experimental designs with RNA-seq. Our results describe the successful establishment of SLAM-seq in CHO cells and therefore facilitate its future use in other scenarios where dynamic transcriptome profiling in CHO cells is essential.

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.

Making use of noise in biological systems

Disorder and noise are inherent in biological systems. They are required to provide systems with the advantages required for proper functioning. Noise is a part of the flexibility and plasticity of biological systems. It provides systems with increased routes, improves information transfer, and assists in response triggers. This paper reviews recent studies on noise at the genome, cellular, and whole organ levels. We focus on the need to use noise in system engineering. We present some of the challenges faced in studying noise. Optimizing the efficiency of complex systems requires a degree of variability in their functions within certain limits. Constrained noise can be considered a method for improving system robustness by regulating noise levels in continuously dynamic settings. The digital pill-based artificial intelligence (AI)-based platform is the first to implement second-generation AI comprising variability-based signatures. This platform enhances the efficacy of the therapeutic regimens. Systems requiring variability and mechanisms regulating noise are mandatory for understanding biological functions.

Probing transient memory of cellular states using single-cell lineages

The inherent stochasticity in the gene product levels can drive single cells within an isoclonal population to different phenotypic states. The dynamic nature of this intercellular variation, where individual cells can transition between different states over time, makes it a particularly hard phenomenon to characterize. We reviewed recent progress in leveraging the classical Luria-Delbrück experiment to infer the transient heritability of the cellular states. Similar to the original experiment, individual cells were first grown into cell colonies, and then, the fraction of cells residing in different states was assayed for each colony. We discuss modeling approaches for capturing dynamic state transitions in a growing cell population and highlight formulas that identify the kinetics of state switching from the extent of colony-to-colony fluctuations. The utility of this method in identifying multi-generational memory of the both expression and phenotypic states is illustrated across diverse biological systems from cancer drug resistance, reactivation of human viruses, and cellular immune responses. In summary, this fluctuation-based methodology provides a powerful approach for elucidating cell-state transitions from a single time point measurement, which is particularly relevant in situations where measurements lead to cell death (as in single-cell RNA-seq or drug treatment) or cause an irreversible change in cell physiology.

Gene expression model inference from snapshot RNA data using Bayesian non-parametrics

Gene expression models, which are key towards understanding cellular regulatory response, underlie observations of single-cell transcriptional dynamics. Although RNA expression data encode information on gene expression models, existing computational frameworks do not perform simultaneous Bayesian inference of gene expression models and parameters from such data. Rather, gene expression models-composed of gene states, their connectivities and associated parameters-are currently deduced by pre-specifying gene state numbers and connectivity before learning associated rate parameters. Here we propose a method to learn full distributions over gene states, state connectivities and associated rate parameters, simultaneously and self-consistently from single-molecule RNA counts. We propagate noise from fluctuating RNA counts over models by treating models themselves as random variables. We achieve this within a Bayesian non-parametric paradigm. We demonstrate our method on the Escherichia coli lacZ pathway and the Saccharomyces cerevisiae STL1 pathway, and verify its robustness on synthetic data.

Coupling gene expression dynamics to cell size dynamics and cell cycle events: Exact and approximate solutions of the extended telegraph model

The standard model describing the fluctuations of mRNA numbers in single cells is the telegraph model which includes synthesis and degradation of mRNA, and switching of the gene between active and inactive states. While commonly used, this model does not describe how fluctuations are influenced by the cell cycle phase, cellular growth and division, and other crucial aspects of cellular biology. Here, we derive the analytical time-dependent solution of an extended telegraph model that explicitly considers the doubling of gene copy numbers upon DNA replication, dependence of the mRNA synthesis rate on cellular volume, gene dosage compensation, partitioning of molecules during cell division, cell-cycle duration variability, and cell-size control strategies. Based on the time-dependent solution, we obtain the analytical distributions of transcript numbers for lineage and population measurements in steady-state growth and also find a linear relation between the Fano factor of mRNA fluctuations and cell volume fluctuations. We show that generally the lineage and population distributions in steady-state growth cannot be accurately approximated by the steady-state solution of extrinsic noise models, i.e. a telegraph model with parameters drawn from probability distributions. This is because the mRNA lifetime is often not small enough compared to the cell cycle duration to erase the memory of division and replication. Accurate approximations are possible when this memory is weak, e.g. for genes with bursty expression and for which there is sufficient gene dosage compensation when replication occurs.

A guidebook of spatial transcriptomic technologies, data resources and analysis approaches

Advances in transcriptomic technologies have deepened our understanding of the cellular gene expression programs of multicellular organisms and provided a theoretical basis for disease diagnosis and therapy. However, both bulk and single-cell RNA sequencing approaches lose the spatial context of cells within the tissue microenvironment, and the development of spatial transcriptomics has made overall bias-free access to both transcriptional information and spatial information possible. Here, we elaborate development of spatial transcriptomic technologies to help researchers select the best-suited technology for their goals and integrate the vast amounts of data to facilitate data accessibility and availability. Then, we marshal various computational approaches to analyze spatial transcriptomic data for various purposes and describe the spatial multimodal omics and its potential for application in tumor tissue. Finally, we provide a detailed discussion and outlook of the spatial transcriptomic technologies, data resources and analysis approaches to guide current and future research on spatial transcriptomics.

Quantifying how post-transcriptional noise and gene copy number variation bias transcriptional parameter inference from mRNA distributions

Transcriptional rates are often estimated by fitting the distribution of mature mRNA numbers measured using smFISH (single molecule fluorescence in situ hybridization) with the distribution predicted by the telegraph model of gene expression, which defines two promoter states of activity and inactivity. However, fluctuations in mature mRNA numbers are strongly affected by processes downstream of transcription. In addition, the telegraph model assumes one gene copy but in experiments, cells may have two gene copies as cells replicate their genome during the cell cycle. While it is often presumed that post-transcriptional noise and gene copy number variation affect transcriptional parameter estimation, the size of the error introduced remains unclear. To address this issue, here we measure both mature and nascent mRNA distributions of GAL10 in yeast cells using smFISH and classify each cell according to its cell cycle phase. We infer transcriptional parameters from mature and nascent mRNA distributions, with and without accounting for cell cycle phase and compare the results to live-cell transcription measurements of the same gene. We find that: (i) correcting for cell cycle dynamics decreases the promoter switching rates and the initiation rate, and increases the fraction of time spent in the active state, as well as the burst size; (ii) additional correction for post-transcriptional noise leads to further increases in the burst size and to a large reduction in the errors in parameter estimation. Furthermore, we outline how to correctly adjust for measurement noise in smFISH due to uncertainty in transcription site localisation when introns cannot be labelled. Simulations with parameters estimated from nascent smFISH data, which is corrected for cell cycle phases and measurement noise, leads to autocorrelation functions that agree with those obtained from live-cell imaging.

Quantifying the phenotypic information in mRNA abundance

Quantifying the dependency between mRNA abundance and downstream cellular phenotypes is a fundamental open problem in biology. Advances in multimodal single-cell measurement technologies provide an opportunity to apply new computational frameworks to dissect the contribution of individual genes and gene combinations to a given phenotype. Using an information theory approach, we analyzed multimodal data of the expression of 83 genes in the Ca2+ signaling network and the dynamic Ca2+ response in the same cell. We found that the overall expression levels of these 83 genes explain approximately 60% of Ca2+ signal entropy. The average contribution of each single gene was 17%, revealing a large degree of redundancy between genes. Using different heuristics, we estimated the dependency between the size of a gene set and its information content, revealing that on average, a set of 53 genes contains 54% of the information about Ca2+ signaling. Our results provide the first direct quantification of information content about complex cellular phenotype that exists in mRNA abundance measurements.

The emerging landscape of spatial profiling technologies

Improved scale, multiplexing and resolution are establishing spatial nucleic acid and protein profiling methods as a major pillar for cellular atlas building of complex samples, from tissues to full organisms. Emerging methods yield omics measurements at resolutions covering the nano- to microscale, enabling the charting of cellular heterogeneity, complex tissue architectures and dynamic changes during development and disease. We present an overview of the developing landscape of in situ spatial genome, transcriptome and proteome technologies, exemplify their impact on cell biology and translational research, and discuss current challenges for their community-wide adoption. Among many transformative applications, we envision that spatial methods will map entire organs and enable next-generation pathology.

Cellular Environment and Phenotypic Heterogeneity: How Data-Driven Modeling Finds the Smoking Gun

The cellular environment modifies cellular phenotypes, in particular, the stress response phenotype, which easily exhibits high phenotypic heterogeneity due to the common characteristics of its regulatory networks. The aim of this work is to quantify and interpret the impact of collagen type I, a major component of the cellular environment, on the phenotypic heterogeneity of the cellular response. Our approach combines in an original way the monitoring of the response of a single cell and the mathematical modeling of the network. After a detailed statistical description of the phenotypic heterogeneity of the cellular response, the mathematical modeling explains how the observed changes can be explained by an induced increase in the average expression of a central protein of the regulatory network. The predictions of the data-driven model are fully consistent with the biochemical measurements performed. The framework presented here is also a new general methodology to study phenotypic heterogeneity, although we focus here on the response to proteotoxic stress in HeLa cells.

Mechanisms of cellular mRNA transcript homeostasis

For most genes, mRNA transcript abundance scales with cell size to ensure a constant concentration. Scaling of mRNA synthesis rates with cell size plays an important role, with regulation of the activity and abundance of RNA polymerase II (Pol II) now emerging as a key point of control. However, there is also considerable evidence for feedback mechanisms that kinetically couple the rates of mRNA synthesis, nuclear export, and degradation to allow cells to compensate for changes in one by adjusting the others. Researchers are beginning to integrate results from these different fields to reveal the mechanisms underlying transcript homeostasis. This will be crucial for moving beyond our current understanding of relative gene expression towards an appreciation of how absolute transcript levels are linked to other aspects of the cellular phenotype.

Epigenetic fluctuations underlie gene expression timescales and variability

Isogenic populations of mammalian cells exhibit significant gene expression variability. This variability can be separated into two sources, cis, or allele-specific sources, and trans and global processes. Furthermore, each source of variability has its own timescale. Fast timescales will result in rapid fluctuation of gene expression whereas slow timescales will result in longer persistence of gene expression levels over time. Here we investigated sources of gene expression that are intrinsic, i.e. coming from cis-regulatory factors and follow slow timescales. To do so, we developed a reporter system that isolates allele-specific variability and measures its persistence in imaging and long-term fluctuation analysis experiments. Our results identify a new source of gene expression variability that is allele-specific but that is fluctuating on timescales of days. We hypothesized that allele-specific fluctuations of epigenetic regulatory factors are responsible for the newly discovered allele-specific and slow source of gene expression variability. Using mathematical modeling we showed that the addition of this effect to the two-state model is sufficient to account for all empirical observation. Furthermore, using direct assays of chromatin markers we find fluctuation in H3K4me3 levels that match the observed changes in gene expression levels providing direct experimental support of our model. Collectively, our work shows that slow fluctuations of regulatory chromatin modifications contribute to the variability in gene expression.

miR-9a regulates levels of both rhomboid mRNA and protein in the early Drosophila melanogaster embryo

MicroRNAs can have subtle and combinatorial effects on the levels of the targets and pathways they act on. Studying the consequences of a single microRNA knockout often proves difficult as many such knockouts exhibit phenotypes only under stress conditions. This has often led to the hypothesis that microRNAs buffer the effects of intrinsic and environmental stochasticity on gene expression. Observing and understanding this buffering effect entails quantitative analysis of microRNA and target expression in single cells. To this end, we have employed single-molecule fluorescence in situ hybridization, immunofluorescence, and high-resolution confocal microscopy to investigate the effects of miR-9a loss on the expression of the serine-protease Rhomboid in Drosophila melanogaster early embryos. Our single-cell quantitative approach shows that spatially, the rhomboid mRNA pattern is identical in WT and miR-9a knockout embryos. However, we find that the number of mRNA molecules per cell is higher when miR-9a is absent, and the level and temporal accumulation of rhomboid protein shows a more dramatic increase in the miR-9a knockout. Specifically, we see accumulation of rhomboid protein in miR-9a mutants by stage 5, much earlier than in WT. The data, therefore, show that miR-9a functions in the regulation of rhomboid mRNA and protein levels. While further work is required to establish whether rhomboid is a direct target of miR-9 in Drosophila, our results further establish the miR-9 family microRNAs as conserved regulators of timing in neurogenic processes. This study shows the power of single-cell quantification as an experimental tool to study phenotypic consequences of microRNA mis-regulation.

Behavioral Neuroscience in the Era of Genomics: Tools and Lessons for Analyzing High-Dimensional Datasets

Behavioral neuroscience underwent a technology-driven revolution with the emergence of machine-vision and machine-learning technologies. These technological advances facilitated the generation of high-resolution, high-throughput capture and analysis of complex behaviors. Therefore, behavioral neuroscience is becoming a data-rich field. While behavioral researchers use advanced computational tools to analyze the resulting datasets, the search for robust and standardized analysis tools is still ongoing. At the same time, the field of genomics exploded with a plethora of technologies which enabled the generation of massive datasets. This growth of genomics data drove the emergence of powerful computational approaches to analyze these data. Here, we discuss the composition of a large behavioral dataset, and the differences and similarities between behavioral and genomics data. We then give examples of genomics-related tools that might be of use for behavioral analysis and discuss concepts that might emerge when considering the two fields together.

A single-cell resolved cell-cell communication model explains lineage commitment in hematopoiesis

Cells do not make fate decisions independently. Arguably, every cell-fate decision occurs in response to environmental signals. In many cases, cell-cell communication alters the dynamics of the internal gene regulatory network of a cell to initiate cell-fate transitions, yet models rarely take this into account. Here, we have developed a multiscale perspective to study the granulocyte-monocyte versus megakaryocyte-erythrocyte fate decisions. This transition is dictated by the GATA1-PU.1 network: a classical example of a bistable cell-fate system. We show that, for a wide range of cell communication topologies, even subtle changes in signaling can have pronounced effects on cell-fate decisions. We go on to show how cell-cell coupling through signaling can spontaneously break the symmetry of a homogenous cell population. Noise, both intrinsic and extrinsic, shapes the decision landscape profoundly, and affects the transcriptional dynamics underlying this important hematopoietic cell-fate decision-making system.

This article has an associated ‘The people behind the papers’ interview.

Summary: Through theory and computational modeling, cell-cell communication is revealed to be a crucial and under-appreciated determinant of cell-fate decision-making during hematopoiesis.

Reduction in gene expression noise by targeted increase in accessibility at gene loci

Many eukaryotic genes are expressed in randomly initiated bursts that are punctuated by periods of quiescence. Here, we show that the intermittent access of the promoters to transcription factors through relatively impervious chromatin contributes to this "noisy" transcription. We tethered a nuclease-deficient Cas9 fused to a histone acetyl transferase at the promoters of two endogenous genes in HeLa cells. An assay for transposase-accessible chromatin using sequencing showed that the activity of the histone acetyl transferase altered the chromatin architecture locally without introducing global changes in the nucleus and rendered the targeted promoters constitutively accessible. We measured the gene expression variability from the gene loci by performing single-molecule fluorescence in situ hybridization against mature messenger RNAs (mRNAs) and by imaging nascent mRNA molecules present at active gene loci in single cells. Because of the increased accessibility of the promoter to transcription factors, the transcription from two genes became less noisy, even when the average levels of expression did not change. In addition to providing evidence for chromatin accessibility as a determinant of the noise in gene expression, our study offers a mechanism for controlling gene expression noise which is otherwise unavoidable.

PRC1 drives Polycomb-mediated gene repression by controlling transcription initiation and burst frequency

The Polycomb repressive system plays a fundamental role in controlling gene expression during mammalian development. To achieve this, Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) bind target genes and use histone modification-dependent feedback mechanisms to form Polycomb chromatin domains and repress transcription. The inter-relatedness of PRC1 and PRC2 activity at these sites has made it difficult to discover the specific components of Polycomb chromatin domains that drive gene repression and to understand mechanistically how this is achieved. Here, by exploiting rapid degron-based approaches and time-resolved genomics, we kinetically dissect Polycomb-mediated repression and discover that PRC1 functions independently of PRC2 to counteract RNA polymerase II binding and transcription initiation. Using single-cell gene expression analysis, we reveal that PRC1 acts uniformly within the cell population and that repression is achieved by controlling transcriptional burst frequency. These important new discoveries provide a mechanistic and conceptual framework for Polycomb-dependent transcriptional control.

Stimulus-specific responses in innate immunity: Multilayered regulatory circuits

Immune sentinel cells initiate immune responses to pathogens and tissue injury and are capable of producing highly stimulus-specific responses. Insight into the mechanisms underlying such specificity has come from the identification of regulatory factors and biochemical pathways, as well as the definition of signaling circuits that enable combinatorial and temporal coding of information. Here, we review the multi-layered molecular mechanisms that underlie stimulus-specific gene expression in macrophages. We categorize components of inflammatory and anti-pathogenic signaling pathways into five layers of regulatory control and discuss unifying mechanisms determining signaling characteristics at each layer. In this context, we review mechanisms that enable combinatorial and temporal encoding of information, identify recurring regulatory motifs and principles, and present strategies for integrating experimental and computational approaches toward the understanding of signaling specificity in innate immunity.

Precision of Protein Thermometry

Temperature sensing is a ubiquitous cell behavior, but the fundamental limits to the precision of temperature sensing are poorly understood. Unlike in chemical concentration sensing, the precision of temperature sensing is not limited by extrinsic fluctuations in the temperature field itself. Instead, we find that precision is limited by the intrinsic copy number, turnover, and binding kinetics of temperature-sensitive proteins. Developing a model based on the canonical TlpA protein, we find that a cell can estimate temperature to within 2%. We compare this prediction with in vivo data on temperature sensing in bacteria.

Bridging scales: From cell biology to physiology using in situ single-cell technologies

Biological organization crosses multiple spatial scales: from molecular, cellular, to tissues and organs. The proliferation of molecular profiling technologies enables increasingly detailed cataloging of the components at each scale. However, the scarcity of spatial profiling has made it challenging to bridge across these scales. Emerging technologies based on highly multiplexed in situ profiling are paving the way to study the spatial organization of cells and tissues in greater detail. These new technologies provide the data needed to cross the scale from cell biology to physiology and identify the fundamental principles that govern tissue organization. Here, we provide an overview of these key technologies and discuss the current and future insights these powerful techniques enable.

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.

The circadian oscillator analysed at the single-transcript level

The circadian clock is an endogenous and self-sustained oscillator that anticipates daily environmental cycles. While rhythmic gene expression of circadian genes is well-described in populations of cells, the single-cell mRNA dynamics of multiple core clock genes remain largely unknown. Here we use single-molecule fluorescence in situ hybridisation (smFISH) at multiple time points to measure pairs of core clock transcripts, Rev-erbα (Nr1d1), Cry1 and Bmal1, in mouse fibroblasts. The mean mRNA level oscillates over 24 h for all three genes, but mRNA numbers show considerable spread between cells. We develop a probabilistic model for multivariate mRNA counts using mixtures of negative binomials, which accounts for transcriptional bursting, circadian time and cell-to-cell heterogeneity, notably in cell size. Decomposing the mRNA variability into distinct noise sources shows that clock time contributes a small fraction of the total variability in mRNA number between cells. Thus, our results highlight the intrinsic biological challenges in estimating circadian phase from single-cell mRNA counts and suggest that circadian phase in single cells is encoded post-transcriptionally.

Perspective: Controlling Epidermal Terminal Differentiation with Transcriptional Bursting and RNA Bodies

The outer layer of the skin, the epidermis, is the principal barrier to the external environment: post-mitotic cells terminally differentiate to form a tough outer cornified layer of enucleate and flattened cells that confer the majority of skin barrier function. Nuclear degradation is required for correct cornified envelope formation. This process requires mRNA translation during the process of nuclear destruction. In this review and perspective, we address the biology of transcriptional bursting and the formation of ribonuclear particles in model organisms including mammals, and then examine the evidence that these phenomena occur as part of epidermal terminal differentiation.

Using single-plant-omics in the field to link maize genes to functions and phenotypes

Most of our current knowledge on plant molecular biology is based on experiments in controlled laboratory environments. However, translating this knowledge from the laboratory to the field is often not straightforward, in part because field growth conditions are very different from laboratory conditions. Here, we test a new experimental design to unravel the molecular wiring of plants and study gene-phenotype relationships directly in the field. We molecularly profiled a set of individual maize plants of the same inbred background grown in the same field and used the resulting data to predict the phenotypes of individual plants and the function of maize genes. We show that the field transcriptomes of individual plants contain as much information on maize gene function as traditional laboratory-generated transcriptomes of pooled plant samples subject to controlled perturbations. Moreover, we show that field-generated transcriptome and metabolome data can be used to quantitatively predict individual plant phenotypes. Our results show that profiling individual plants in the field is a promising experimental design that could help narrow the lab-field gap.