Living with noisy genes: how cells function reliably with inherent variability in gene expression

An Ocular Gene-Set Expression Library for Heritability Partition and Cell Line Enrichment Analyses

Purpose

Use of genome-wide association studies (GWASs) in combination with transcriptomic arrays of different tissues or cell lines can reveal relevant cellular profiles and significantly improve understanding of the mechanisms of diseases. However, due to difficulty of access, few ocular transcriptomics datasets are available. This work aimed to create and make available an expression library platform that can be used with popular and versatile tools such as the linkage disequilibrium score (LDSC) regression techniques to identify specific cell lines enriched in ocular diseases.

Methods

We used transcriptome information from six publicly available single-cell and single-nucleus RNA sequence datasets to obtain matrices of normalized gene expression and estimated enrichment of the 10% most expressed transcripts in each cell line. We tested for type 1 error using simulated GWAS datasets and then used LDSC regression analyses to study the enrichment in different eye diseases.

Results

Gene expression databases for over 197 ocular cell lines were generated. Simulations of 1000 random GWASs of varying sample sizes showed no genomic inflation. Cell line-specific analyses of GWAS results found that genes near single nucleotide polymorphisms (SNPs) associated with refractive error were significantly enriched in cones (P = 0.008), rods (P = 0.002) and peripheral retina Müller cells (P = 0.002), juxtacanalicular cribriform cells (P = 0.0017), stromal keratocytes (P = 0.0018), and one beam-cell subpopulation (P = 0.0028) in primary open-angle glaucoma, emphasizing the importance of intraocular pressure in its etiology.

Conclusions

We have built a structured ocular transcriptomics library to estimate cell line enrichment among association signals from genome-wide association studies that may be extended by incorporating other datasets. We identified cells that appear important in the genetics of common eye diseases.

Dynamical model of antibiotic responses linking expression of resistance genes to metabolism explains emergence of heterogeneity during drug exposures

Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.

Exact solution of a three-stage model of stochastic gene expression including cell-cycle dynamics

The classical three-stage model of stochastic gene expression predicts the statistics of single cell mRNA and protein number fluctuations as a function of the rates of promoter switching, transcription, translation, degradation and dilution. While this model is easily simulated, its analytical solution remains an unsolved problem. Here we modify this model to explicitly include cell-cycle dynamics and then derive an exact solution for the time-dependent joint distribution of mRNA and protein numbers. We show large differences between this model and the classical model which captures cell-cycle effects implicitly via effective first-order dilution reactions. In particular we find that the Fano factor of protein numbers calculated from a population snapshot measurement are underestimated by the classical model whereas the correlation between mRNA and protein can be either over- or underestimated, depending on the timescales of mRNA degradation and promoter switching relative to the mean cell-cycle duration time.

The interplay between metabolic stochasticity and cAMP-CRP regulation in single E. coli cells

The inherent stochasticity of metabolism raises a critical question for understanding homeostasis: are cellular processes regulated in response to internal fluctuations? Here, we show that, in E. coli cells under constant external conditions, catabolic enzyme expression continuously responds to metabolic fluctuations. The underlying regulatory feedback is enabled by the cyclic AMP (cAMP) and cAMP receptor protein (CRP) system, which controls catabolic enzyme expression based on metabolite concentrations. Using single-cell microscopy, genetic constructs in which this feedback is disabled, and mathematical modeling, we show how fluctuations circulate through the metabolic and genetic network at sub-cell-cycle timescales. Modeling identifies four noise propagation modes, including one specific to CRP regulation. Together, these modes correctly predict noise circulation at perturbed cAMP levels. The cAMP-CRP system may thus have evolved to control internal metabolic fluctuations in addition to external growth conditions. We conjecture that second messengers may more broadly function to achieve cellular homeostasis.

Dynamical model of antibiotic responses linking expression of resistance to metabolism explains emergence of heterogeneity during drug exposures

Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistance tet operon in E. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.

Quantifying Parameter Interdependence in Stochastic Discrete Models of Biochemical Systems

Stochastic modeling of biochemical processes at the cellular level has been the subject of intense research in recent years. The Chemical Master Equation is a broadly utilized stochastic discrete model of such processes. Numerous important biochemical systems consist of many species subject to many reactions. As a result, their mathematical models depend on many parameters. In applications, some of the model parameters may be unknown, so their values need to be estimated from the experimental data. However, the problem of parameter value inference can be quite challenging, especially in the stochastic setting. To estimate accurately the values of a subset of parameters, the system should be sensitive with respect to variations in each of these parameters and they should not be correlated. In this paper, we propose a technique for detecting collinearity among models' parameters and we apply this method for selecting subsets of parameters that can be estimated from the available data. The analysis relies on finite-difference sensitivity estimations and the singular value decomposition of the sensitivity matrix. We illustrated the advantages of the proposed method by successfully testing it on several models of biochemical systems of practical interest.

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.

Inferencing Bulk Tumor and Single-Cell Multi-Omics Regulatory Networks for Discovery of Biomarkers and Therapeutic Targets

There are insufficient accurate biomarkers and effective therapeutic targets in current cancer treatment. Multi-omics regulatory networks in patient bulk tumors and single cells can shed light on molecular disease mechanisms. Integration of multi-omics data with large-scale patient electronic medical records (EMRs) can lead to the discovery of biomarkers and therapeutic targets. In this review, multi-omics data harmonization methods were introduced, and common approaches to molecular network inference were summarized. Our Prediction Logic Boolean Implication Networks (PLBINs) have advantages over other methods in constructing genome-scale multi-omics networks in bulk tumors and single cells in terms of computational efficiency, scalability, and accuracy. Based on the constructed multi-modal regulatory networks, graph theory network centrality metrics can be used in the prioritization of candidates for discovering biomarkers and therapeutic targets. Our approach to integrating multi-omics profiles in a patient cohort with large-scale patient EMRs such as the SEER-Medicare cancer registry combined with extensive external validation can identify potential biomarkers applicable in large patient populations. These methodologies form a conceptually innovative framework to analyze various available information from research laboratories and healthcare systems, accelerating the discovery of biomarkers and therapeutic targets to ultimately improve cancer patient survival outcomes.

Stochastic Fluctuations Drive Non-genetic Evolution of Proliferation in Clonal Cancer Cell Populations

Evolutionary dynamics allows us to understand many changes happening in a broad variety of biological systems, ranging from individuals to complete ecosystems. It is also behind a number of remarkable organizational changes that happen during the natural history of cancers. These reflect tumour heterogeneity, which is present at all cellular levels, including the genome, proteome and phenome, shaping its development and interrelation with its environment. An intriguing observation in different cohorts of oncological patients is that tumours exhibit an increased proliferation as the disease progresses, while the timescales involved are apparently too short for the fixation of sufficient driver mutations to promote explosive growth. Here, we discuss how phenotypic plasticity, emerging from a single genotype, may play a key role and provide a ground for a continuous acceleration of the proliferation rate of clonal populations with time. We address this question by combining the analysis of real-time growth of non-small-cell lung carcinoma cells (N-H460) together with stochastic and deterministic mathematical models that capture proliferation trait heterogeneity in clonal populations to elucidate the contribution of phenotypic transitions on tumour growth dynamics.

Approximating solutions of the Chemical Master equation using neural networks

The Chemical Master Equation (CME) provides an accurate description of stochastic biochemical reaction networks in well-mixed conditions, but it cannot be solved analytically for most systems of practical interest. Although Monte Carlo methods provide a principled means to probe system dynamics, the large number of simulations typically required can render the estimation of molecule number distributions and other quantities infeasible. In this article, we aim to leverage the representational power of neural networks to approximate the solutions of the CME and propose a framework for the Neural Estimation of Stochastic Simulations for Inference and Exploration (Nessie). Our approach is based on training neural networks to learn the distributions predicted by the CME from relatively few stochastic simulations. We show on biologically relevant examples that simple neural networks with one hidden layer can capture highly complex distributions across parameter space, thereby accelerating computationally intensive tasks such as parameter exploration and inference.

Cytoplasmic LMO2-LDB1 Complex Activates STAT3 Signaling through Interaction with gp130-JAK in Glioma Stem Cells

The oncogenic role of nuclear LIM domain only 2 (LMO2) as a transcriptional regulator is well established, but its function in the cytoplasm is largely unknown. Here, we identified LMO2 as a cytoplasmic activator for signal transducer and activator of transcription 3 (STAT3) signaling in glioma stem cells (GSCs) through biochemical and bioinformatics analyses. LMO2 increases STAT3 phosphorylation by interacting with glycoprotein 130 (gp130) and Janus kinases (JAKs). LMO2-driven activation of STAT3 signaling requires the LDB1 protein and leads to increased expression of an inhibitor of differentiation 1 (ID1), a master regulator of cancer stemness. Our findings indicate that the cytoplasmic LMO2-LDB1 complex plays a crucial role in the activation of the GSC signaling cascade via interaction with gp130 and JAK1/2. Thus, LMO2-LDB1 is a bona fide oncogenic protein complex that activates either the JAK-STAT signaling cascade in the cytoplasm or direct transcriptional regulation in the nucleus.

CellDynaMo–stochastic reaction-diffusion-dynamics model: Application to search-and-capture process of mitotic spindle assembly

We introduce a Stochastic Reaction-Diffusion-Dynamics Model (SRDDM) for simulations of cellular mechanochemical processes with high spatial and temporal resolution. The SRDDM is mapped into the CellDynaMo package, which couples the spatially inhomogeneous reaction-diffusion master equation to account for biochemical reactions and molecular transport within the Langevin Dynamics (LD) framework to describe dynamic mechanical processes. This computational infrastructure allows the simulation of hours of molecular machine dynamics in reasonable wall-clock time. We apply SRDDM to test performance of the Search-and-Capture of mitotic spindle assembly by simulating, in three spatial dimensions, dynamic instability of elastic microtubules anchored in two centrosomes, movement and deformations of geometrically realistic centromeres with flexible kinetochores and chromosome arms. Furthermore, the SRDDM describes the mechanics and kinetics of Ndc80 linkers mediating transient attachments of microtubules to the chromosomal kinetochores. The rates of these attachments and detachments depend upon phosphorylation states of the Ndc80 linkers, which are regulated in the model by explicitly accounting for the reactions of Aurora A and B kinase enzymes undergoing restricted diffusion. We find that there is an optimal rate of microtubule-kinetochore detachments which maximizes the accuracy of the chromosome connections, that adding chromosome arms to kinetochores improve the accuracy by slowing down chromosome movements, that Aurora A and kinetochore deformations have a small positive effect on the attachment accuracy, and that thermal fluctuations of the microtubules increase the rates of kinetochore capture and also improve the accuracy of spindle assembly.

The CellDynaMo package models, in 3D, any cellular subsystem where sufficient detail of the macromolecular players and the kinetics of relevant reactions are available. The package is based on the Stochastic Reaction-Diffusion-Dynamics model that combines the stochastic description of chemical kinetics, Brownian diffusion-based description of molecular transport, and Langevin dynamics-based representation of mechanical processes most pertinent to the system. We apply the model to test the Search-and-Capture mechanism of mitotic spindle assembly. We find that there is an optimal rate of microtubule-kinetochore detachments which maximizes the accuracy of chromosome connections, that chromosome arms improve the attachment accuracy by slowing down chromosome movements, that Aurora A kinase and kinetochore deformations have small positive effects on the accuracy, and that thermal fluctuations of the microtubules increase the rates of kinetochore capture and also improve the accuracy.

Environment-driven shifts in inter-individual variation and phenotypic integration within subnetworks of the mussel transcriptome and proteome

The environment can alter the magnitude of phenotypic variation among individuals, potentially influencing evolutionary trajectories. However, environmental influences on variation are complex and remain understudied. Populations in heterogeneous environments might exhibit more variation, the amount of variation could differ between benign and stressful conditions, and/or variation might manifest in different ways among stages of the gene‐to‐protein expression cascade or among physiological functions. Here, we explore these three issues by quantifying patterns of inter‐individual variation in both transcript and protein expression levels among California mussels, Mytilus californianus Conrad. Mussels were exposed to five ecologically relevant treatments that varied in the mean and interindividual heterogeneity of body temperature. To target a diverse set of physiological functions, we assessed variation within 19 expression subnetworks, including canonical stress‐response pathways and empirically derived coexpression clusters that represent a diffuse set of cellular processes. Variation in expression was particularly pronounced in the treatments with high mean and heterogeneous body temperatures. However, with few exceptions, environment‐dependent shifts of variation in the transcriptome were not reflected in the proteome. A metric of phenotypic integration provided evidence for a greater degree of constraint on relative expression levels (i.e., stronger correlation) within expression subnetworks in benign, homogeneous environments. Our results suggest that environments that are more stressful on average – and which also tend to be more heterogeneous – can relax these expression constraints and reduce phenotypic integration within biochemical subnetworks. Context‐dependent “unmasking” of functional variation may contribute to interindividual differences in physiological phenotype and performance in stressful environments.

Mapping the dynamic transfer functions of eukaryotic gene regulation

Biological information can be encoded within the dynamics of signaling components, which has been implicated in a broad range of physiological processes including stress response, oncogenesis, and stem cell differentiation. To study the complexity of information transfer across the eukaryotic promoter, we screened 119 dynamic conditions-modulating the pulse frequency, amplitude, and pulse width of light-regulating the binding of an epigenome editor to a fluorescent reporter. This system revealed tunable gene expression and filtering behaviors and provided a quantification of the limit to the amount of information that can be reliably transferred across a single promoter as ∼1.7 bits. Using a library of over 100 orthogonal chromatin regulators, we further determined that chromatin state could be used to tune mutual information and expression levels, as well as completely alter the input-output transfer function of the promoter. This system unlocks the information-rich content of eukaryotic gene regulation.

Improving the performance of machine learning models for biotechnology: The quest for deus ex machina

Machine learning is becoming an integral part of the Design-Build-Test-Learn cycle in biotechnology. Machine learning models learn from collected datasets such as omics data and predict a defined outcome, which has led to both production improvements and predictive tools in the field. Robust prediction of the behavior of microbial cell factories and production processes not only greatly increases our understanding of the function of such systems, but also provides significant savings of development time. However, many pitfalls when modeling biological data - bad fit, noisy data, model instability, low data quantity and imbalances in the data - cause models to suffer in their performance. Here we provide an accessible, in-depth analysis on the problems created by these pitfalls, as well as means of their detection and mediation, with a focus on supervised learning. Assessing the state of the art, we show that, currently, in-depth analyses of model performance are often absent and must be improved. This review provides a toolbox for the analysis of model robustness and performance, and simultaneously proposes a standard for the community to facilitate future work. It is further accompanied by an interactive online tutorial on the discussed issues.

Temporal changes in the gene expression heterogeneity during brain development and aging

Cells in largely non-mitotic tissues such as the brain are prone to stochastic (epi-)genetic alterations that may cause increased variability between cells and individuals over time. Although increased inter-individual heterogeneity in gene expression was previously reported, whether this process starts during development or if it is restricted to the aging period has not yet been studied. The regulatory dynamics and functional significance of putative aging-related heterogeneity are also unknown. Here we address these by a meta-analysis of 19 transcriptome datasets from three independent studies, covering diverse human brain regions. We observed a significant increase in inter-individual heterogeneity during aging (20 + years) compared to postnatal development (0 to 20 years). Increased heterogeneity during aging was consistent among different brain regions at the gene level and associated with lifespan regulation and neuronal functions. Overall, our results show that increased expression heterogeneity is a characteristic of aging human brain, and may influence aging-related changes in brain functions.

Analytical distributions for detailed models of stochastic gene expression in eukaryotic cells

The stochasticity of gene expression presents significant challenges to the modeling of genetic networks. A two-state model describing promoter switching, transcription, and messenger RNA (mRNA) decay is the standard model of stochastic mRNA dynamics in eukaryotic cells. Here, we extend this model to include mRNA maturation, cell division, gene replication, dosage compensation, and growth-dependent transcription. We derive expressions for the time-dependent distributions of nascent mRNA and mature mRNA numbers, provided two assumptions hold: 1) nascent mRNA dynamics are much faster than those of mature mRNA; and 2) gene-inactivation events occur far more frequently than gene-activation events. We confirm that thousands of eukaryotic genes satisfy these assumptions by using data from yeast, mouse, and human cells. We use the expressions to perform a sensitivity analysis of the coefficient of variation of mRNA fluctuations averaged over the cell cycle, for a large number of genes in mouse embryonic stem cells, identifying degradation and gene-activation rates as the most sensitive parameters. Furthermore, it is shown that, despite the model's complexity, the time-dependent distributions predicted by our model are generally well approximated by the negative binomial distribution. Finally, we extend our model to include translation, protein decay, and auto-regulatory feedback, and derive expressions for the approximate time-dependent protein-number distributions, assuming slow protein decay. Our expressions enable us to study how complex biological processes contribute to the fluctuations of gene products in eukaryotic cells, as well as allowing a detailed quantitative comparison with experimental data via maximum-likelihood methods.

mRNA with Mammalian Codon Bias Accumulates in Yeast Mutants with Constitutive Stress Granules

Stress granules and P bodies are cytoplasmic structures assembled in response to various stress factors and represent sites of temporary storage or decay of mRNAs. Depending on the source of stress, the formation of these structures may be driven by distinct mechanisms, but several stresses have been shown to stabilize mRNAs via inhibition of deadenylation. A recent study identified yeast gene deletion mutants with constitutive stress granules and elevated P bodies; however, the mechanisms which trigger its formation remain poorly understood. Here, we investigate the possibility of accumulating mRNA with mammalian codon bias, which we termed the model RNA, in these mutants. We found that the model RNA accumulates in dcp2 and xrn1 mutants and in four mutants with constitutive stress granules overlapping with P bodies. However, in eight other mutants with constitutive stress granules, the model RNA is downregulated, or its steady state levels vary. We further suggest that the accumulation of the model RNA is linked to its protection from the main mRNA surveillance path. However, there is no obvious targeting of the model RNA to stress granules or P bodies. Thus, accumulation of the model RNA and formation of constitutive stress granules occur independently and only some paths inducing formation of constitutive stress granules will stabilize mRNA as well.

Quantitation and Comparison of Phenotypic Heterogeneity Among Single Cells of Monoclonal Microbial Populations

Phenotypic heterogeneity within microbial populations arises even when the cells are exposed to putatively constant and homogeneous conditions. The outcome of this phenomenon can affect the whole function of the population, resulting in, for example, new "adapted" metabolic strategies and impacting its fitness at given environmental conditions. Accounting for phenotypic heterogeneity becomes thus necessary, due to its relevance in medical and applied microbiology as well as in environmental processes. Still, a comprehensive evaluation of this phenomenon requires a common and unique method of quantitation, which allows for the comparison between different studies carried out with different approaches. Consequently, in this study, two widely applicable indices for quantitation of heterogeneity were developed. The heterogeneity coefficient (HC) is valid when the population follows unimodal activity, while the differentiation tendency index (DTI) accounts for heterogeneity implying outbreak of subpopulations and multimodal activity. We demonstrated the applicability of HC and DTI for heterogeneity quantitation on stable isotope probing with nanoscale secondary ion mass spectrometry (SIP-nanoSIMS), flow cytometry, and optical microscopy datasets. The HC was found to provide a more accurate and precise measure of heterogeneity, being at the same time consistent with the coefficient of variation (CV) applied so far. The DTI is able to describe the differentiation in single-cell activity within monoclonal populations resolving subpopulations with low cell abundance, individual cells with similar phenotypic features (e.g., isotopic content close to natural abundance, as detected with nanoSIMS). The developed quantitation approach allows for a better understanding on the impact and the implications of phenotypic heterogeneity in environmental, medical and applied microbiology, microbial ecology, cell biology, and biotechnology.

Modulation of first-passage time for gene expression via asymmetric cell division

How to balance the size of exponentially growing cells has always been a focus of biologists. Recent experiments have uncovered that the cell is divided into two daughter cells only when the level of time-keeper protein reaches a fixed threshold and cell division in prokaryote is not completely symmetric. The timing of cell division is essentially random because gene expression is stochastic, but cells seen to manage to have precise timing of cell division events. Although the inter-cellular variability of gene expression has attracted much attention, the randomness of event timing has been rarely studied. In our analysis, the timing of cell division is formulated as the first-passage time (denoted by FPT) for time-keeper protein’s level to cross a critical threshold firstly, we derive exact analytical formulae for the mean and noise of FPT based on stochastic gene expression model with asymmetric cell division. The results of numerical simulation show that the regulatory factors (division rate, newborn cell size, exponential growth rate and threshold) have significant influence on the mean and noise of FPT. We also show that both the increase of division rate and newborn cell size could reduce the mean of FPT and increase the noise of FPT, the larger the exponential growth rate is, the smaller the mean and noise of FPT will be; and the larger the threshold value is, the higher the mean of FPT is and the lower the noise is. In addition, compared with symmetric division, asymmetric division can reduce the mean of FPT and improve the noise of FPT. In summary, our results provide insight into the relationship between regulatory factors and FPT and reveal that asymmetric division is an effective mechanism to shorten the mean of FPT.

Operating regimes in a single enzymatic cascade at ensemble-level

Single enzymatic cascade, ubiquitously found in cellular signaling networks, is a phosphorylation-dephosphorylation reaction cycle causing a transition between inactive and active states of a protein catalysed by kinase and phosphatase, respectively. Steady-state information processing ability of such a cycle (e.g., MAPK cascade) has been classified into four qualitatively different operating regimes, viz., hyperbolic (H), signal-transducing (ST), threshold-hyperbolic (TH) and ultrasensitive (U). These four regimes represent qualitatively different dose-response curves, that is, relationship between concentrations of input kinase (e.g., pMEK) and response activated protein (e.g., pERK). Regimes were identified using a deterministic model accounting for population-averaged behavior only. Operating regimes can be strongly influenced by the inherently present cell-to-cell variability in an ensemble of cells which is captured in the form of pMEK and pERK distributions using reporter-based single-cell experimentation. In this study, we show that such experimentally acquired snapshot pMEK and pERK distribution data of a single MAPK cascade can be directly used to infer the underlying operating regime even in the absence of a dose-response curve. This deduction is possible primarily due to the presence of a monotonic relationship between experimental observables R_IQR, ratio of the inter-quartile range of the pERK and pMEK distribution pairs and R_M, ratio of the medians of the distribution pair. We demonstrate this relationship by systematic analysis of a quasi-steady state approximated model superimposed with an input gamma distribution constrained by the stimulus strength specific pMEK distribution measured on Jurkat-T cells stimulated with PMA. As a first, we show that introduction of cell-to-cell variability only in the upstream kinase achieved by superimposition of an appropriate input pMEK distribution on the dose-response curve can predict bimodal response pERK distribution in ST regime. Implementation of the proposed method on the input-response distribution pair obtained in stimulated Jurkat-T cells revealed that while low-dosage PMA stimulation preserves the H regime observed in resting cells, high-dosage causes H to ST regime transition.

Chemically active nanodroplets in a multi-component fluid

We introduce a model of chemically active particles of a multi-component fluid that can change their interactions with other particles depending on their state. Since such switching of interactions can only be maintained by the input of chemical energy, the system is inherently non-equilibrium. Focusing on a scenario where the equilibrium interactions would lead to condensation into a liquid droplet, and despite the relative simplicity of the interaction rules, these systems display a wealth of interesting and novel behaviors such as oscillations of droplet size and molecular sorting, and raise the possibility of spatio-temporal control of chemical reactions on the nanoscale.

The nonlinear dynamics and fluctuations of mRNA levels in cell cycle coupled transcription

Gene transcription is a noisy process, and cell division cycle is an important source of gene transcription noise. In this work, we develop a mathematical approach by coupling transcription kinetics with cell division cycles to delineate how they are combined to regulate transcription output and noise. In view of gene dosage, a cell cycle is divided into an early stage S1 and a late stage S2. The analytical forms for the mean and the noise of mRNA numbers are given in each stage. The analysis based on these formulas predicts precisely the fold change r* of mRNA numbers from S1 to S2 measured in a mouse embryonic stem cell line. When transcription follows similar kinetics in both stages, r* buffers against DNA dosage variation and r* ∈ (1, 2). Numerical simulations suggest that increasing cell cycle durations up-regulates transcription with less noise, whereas rapid stage transitions induce highly noisy transcription. A minimization of the transcription noise is observed when transcription homeostasis is attained by varying a single kinetic rate. When the transcription level scales with cellular volume, either by reducing the transcription burst frequency or by increasing the burst size in S2, the noise shows only a minor variation over a wide range of cell cycle stage durations. The reduction level in the burst frequency is nearly a constant, whereas the increase in the burst size is conceivably sensitive, when responding to a large random variation of the cell cycle durations and the gene duplication time.

Gene transcription in single cells is inherently a stochastic process, resulting in a large variability in the number of transcripts and constituting the phenotypic heterogeneity in cell population. Cell division cycle has global effects on transcriptional outputs, and is thought to be an additional source of transcription noise. In this work, we develop a hybrid model to delineate the combined contribution of transcription activities and cell divisions in the variability of transcript counts. By working with the analytical forms of the mean and the noise of mRNA numbers, we show that if the transcription kinetic rates do not change considerably, then the average mRNA level is increased about 1 to 2 folds from earlier to later cell cycle stages. When transcription homeostasis is attained by varying a single kinetic rate between the two cell cycle stages, we find no significant changes in the transcription noise, and the homeostasis nearly minimizes the noise. In our continuous study on the transcript concentration homeostasis that the transcription level scales with the cellular volume, we find only minor variations of the noise if the homeostasis is maintained either by reducing the transcription burst frequency or by increasing the burst size in late cell cycle phase, in the face of a large cell cycle stage duration variation. The reduction in the burst frequency is relative robust, while the increase in the burst size is conceivably sensitive, to the large random variation of the cell cycle durations and the gene duplication time.

Autoregulation of bacterial gene expression: lessons from the MazEF toxin-antitoxin system

Autoregulation is the direct modulation of gene expression by the product of the corresponding gene. Autoregulation of bacterial gene expression has been mostly studied at the transcriptional level, when a protein acts as the cognate transcriptional repressor. A recent study investigating dynamics of the bacterial toxin-antitoxin MazEF system has shown how autoregulation at both the transcriptional and post-transcriptional levels affects the heterogeneity of Escherichia coli populations. Toxin-antitoxin systems hold a crucial but still elusive part in bacterial response to stress. This perspective highlights how these modules can also serve as a great model system for investigating basic concepts in gene regulation. However, as the genomic background and environmental conditions substantially influence toxin activation, it is important to study (auto)regulation of toxin-antitoxin systems in well-defined setups as well as in conditions that resemble the environmental niche.

On the duration of the microbial lag phase

When faced with environmental changes, microbes enter a lag phase during which cell growth is arrested, allowing cells to adapt to the new situation. The discovery of the lag phase started the field of gene regulation and led to the unraveling of underlying mechanisms. However, the factors determining the exact duration and dynamics of the lag phase remain largely elusive. Naively, one would expect that cells adapt as quickly as possible, so they can resume growth and compete with other organisms. However, recent studies show that the lag phase can last from several hours up to several days. Moreover, some cells within the same population take much longer than others, despite being genetically identical. In addition, the lag phase duration is also influenced by the past, with recent exposure to a given environment leading to a quicker adaptation when that environment returns. Genome-wide screens in Saccharomyces cerevisiae on carbon source shifts now suggest that the length of the lag phase, the heterogeneity in lag times of individual cells, and the history-dependent behavior are not determined by the time it takes to induce a few specific genes related to uptake and metabolism of a new carbon source. Instead, a major shift in general metabolism, and in particular a switch between fermentation and respiration, is the major bottleneck that determines lag duration. This suggests that there may be a fitness trade-off between complete adaptation of a cell’s metabolism to a given environment, and a short lag phase when the environment changes.

Synthetic STARR-seq reveals how DNA shape and sequence modulate transcriptional output and noise

The binding of transcription factors to short recognition sequences plays a pivotal role in controlling the expression of genes. The sequence and shape characteristics of binding sites influence DNA binding specificity and have also been implicated in modulating the activity of transcription factors downstream of binding. To quantitatively assess the transcriptional activity of tens of thousands of designed synthetic sites in parallel, we developed a synthetic version of STARR-seq (synSTARR-seq). We used the approach to systematically analyze how variations in the recognition sequence of the glucocorticoid receptor (GR) affect transcriptional regulation. Our approach resulted in the identification of a novel highly active functional GR binding sequence and revealed that sequence variation both within and flanking GR’s core binding site can modulate GR activity without apparent changes in DNA binding affinity. Notably, we found that the sequence composition of variants with similar activity profiles was highly diverse. In contrast, groups of variants with similar activity profiles showed specific DNA shape characteristics indicating that DNA shape may be a better predictor of activity than DNA sequence. Finally, using single cell experiments with individual enhancer variants, we obtained clues indicating that the architecture of the response element can independently tune expression mean and cell-to cell variability in gene expression (noise). Together, our studies establish synSTARR as a powerful method to systematically study how DNA sequence and shape modulate transcriptional output and noise.

The expression level of genes is controlled by transcription factors, which are proteins that bind to genomic response elements that contain their recognition DNA sequence. Importantly, genes are not simply turned on but need to be expressed at the right level. This is, at least in part, assured by the sequence composition of genomic response elements. Here, we studied how the recognition DNA sequence influences gene regulation by a transcription factor called the glucocorticoid receptor. Specifically, we developed a method to test the activity of variants in a highly parallelized setting where everything is kept identical except for the sequence of the binding site. The systematic analysis of tens of thousands of sequence variants facilitated the identification of a previously unknown sequence variant with high activity. Moreover, we report how sequence variation of the response element influences cell-to-cell variability in expression levels. Finally, we observe similar activity profiles for distinct sequence variants that share similar three-dimensional DNA shape characteristics arguing that the three-dimensional perception of DNA by the glucocorticoid receptor, modulates its activity towards individual target genes.

Noise⁻Seeded Developmental Pattern Formation in Filamentous Cyanobacteria

Under nitrogen-poor conditions, multicellular cyanobacteria such as Anabaena sp. PCC 7120 undergo a process of differentiation, forming nearly regular, developmental patterns of individual nitrogen-fixing cells, called heterocysts, interspersed between intervals of vegetative cells that carry out photosynthesis. Developmental pattern formation is mediated by morphogen species that can act as activators and inhibitors, some of which can diffuse along filaments. We survey the limitations of the classical, deterministic Turing mechanism that has been often invoked to explain pattern formation in these systems, and then, focusing on a simpler system governed by birth-death processes, we illustrate pedagogically a recently proposed paradigm that provides a much more robust description of pattern formation: stochastic Turing patterns. We emphasize the essential role that cell-to-cell differences in molecular numbers—caused by inevitable fluctuations in gene expression—play, the so called demographic noise, in seeding the formation of stochastic Turing patterns over a much larger region of parameter space, compared to their deterministic counterparts.

Hollow Nanoneedle-Electroporation System To Extract Intracellular Protein Repetitively and Nondestructively

Techniques used to understand the dynamic expression of intracellular proteins are critical in both fundamental biological research and biomedical engineering. Various methods for analyzing proteins have been developed, but these methods require the extraction of intracellular proteins from the cells resulting in cell lysis and subsequent protein purifications from the lysate, which limits the potential of repetitive extraction from the same set of viable cells to track dynamic intracellular protein expression. Therefore, it is crucial to develop novel methods that enable nondestructive and repeated extraction of intracellular proteins. This work reports a hollow nanoneedle-electroporation system for the repeated extraction of intracellular proteins from living cells. Hollow nanoneedles with ∼450 nm diameter were fabricated by a material deposition and etching process, followed by integration with a microfluidic device. Long-lasting electrical pulses were coupled with the nanoneedles to permeate the cell membrane, allowing intracellular contents to diffuse into the microfluidic channels located below the cells via hollow nanoneedles. Using lactate dehydrogenase B (LDHB) as the model intracellular protein, the nanoneedle-electroporation system effectively and repeatedly extracted LDHB from the same set of cells at different time points, followed by quantitative analysis of LDHB via standard enzyme-linked immunosorbent assay. Our work demonstrated an efficient method to nondestructively probe intracellular protein levels and monitor the dynamic protein expression, with great potential to help understanding cell behaviors and functions.

Environmental properties of cells improve machine learning-based phenotype recognition accuracy

To answer major questions of cell biology, it is often essential to understand the complex phenotypic composition of cellular systems precisely. Modern automated microscopes produce vast amounts of images routinely, making manual analysis nearly impossible. Due to their efficiency, machine learning-based analysis software have become essential tools to perform single-cell-level phenotypic analysis of large imaging datasets. However, an important limitation of such methods is that they do not use the information gained from the cellular micro- and macroenvironment: the algorithmic decision is based solely on the local properties of the cell of interest. Here, we present how various features from the surrounding environment contribute to identifying a cell and how such additional information can improve single-cell-level phenotypic image analysis. The proposed methodology was tested for different sizes of Euclidean and nearest neighbour-based cellular environments both on tissue sections and cell cultures. Our experimental data verify that the surrounding area of a cell largely determines its entity. This effect was found to be especially strong for established tissues, while it was somewhat weaker in the case of cell cultures. Our analysis shows that combining local cellular features with the properties of the cell’s neighbourhood significantly improves the accuracy of machine learning-based phenotyping.

On the networked architecture of genotype spaces and its critical effects on molecular evolution

Evolutionary dynamics is often viewed as a subtle process of change accumulation that causes a divergence among organisms and their genomes. However, this interpretation is an inheritance of a gradualistic view that has been challenged at the macroevolutionary, ecological and molecular level. Actually, when the complex architecture of genotype spaces is taken into account, the evolutionary dynamics of molecular populations becomes intrinsically non-uniform, sharing deep qualitative and quantitative similarities with slowly driven physical systems: nonlinear responses analogous to critical transitions, sudden state changes or hysteresis, among others. Furthermore, the phenotypic plasticity inherent to genotypes transforms classical fitness landscapes into multiscapes where adaptation in response to an environmental change may be very fast. The quantitative nature of adaptive molecular processes is deeply dependent on a network-of-networks multilayered structure of the map from genotype to function that we begin to unveil.

Microfluidics and single-cell microscopy to study stochastic processes in bacteria

Bacteria have molecules present in low and fluctuating numbers that randomize cell behaviors. Understanding these stochastic processes and their impact on cells has, until recently, been limited by the lack of single-cell measurement methods. Here, we review recent developments in microfluidics that enable following individual cells over long periods of time under precisely controlled conditions, and counting individual fluorescent molecules in many cells. We showcase discoveries that were made possible using these devices in various aspects of microbiology, such as antibiotic tolerance/persistence, cell-size control, cell-fate determination, DNA damage response, and synthetic biology.

Minimum Action Path Theory Reveals the Details of Stochastic Transitions Out of Oscillatory States

Cell state determination is the outcome of intrinsically stochastic biochemical reactions. Transitions between such states are studied as noise-driven escape problems in the chemical species space. Escape can occur via multiple possible multidimensional paths, with probabilities depending nonlocally on the noise. Here we characterize the escape from an oscillatory biochemical state by minimizing the Freidlin-Wentzell action, deriving from it the stochastic spiral exit path from the limit cycle. We also use the minimized action to infer the escape time probability density function.

Influence of multiplicative stochastic variation on translational elongation rates

Experimental data indicate that stochastic effects exerted at the level of translation contribute substantially to the variation in abundance of proteins expressed at moderate to high levels. This study analyzes the theoretical consequences of fluctuations in residue-specific elongation rates during translation. A simple analytical framework shows that rate variation during elongation gives rise to protein production rates that consist of sums of products of random variables. Simulations show that because the contribution to total variation of products of random variables greatly exceeds that of sums of random variables, the overall distribution exhibits approximately log-normal behavior. Empirical fits of the data can be satisfied by either sums of log-normal distributions, or sums of log-normal and log-logistic distributions. Elongation rate stochastic variation offers an accounting for a major component of biological variation. The analysis provided here highlights a probability distribution that is a natural extension of the Poisson and has broad applicability to many types of multiplicative noise processes.

Likelihood for transcriptions in a genetic regulatory system under asymmetric stable Lévy noise

This work is devoted to investigating the evolution of concentration in a genetic regulation system, when the synthesis reaction rate is under additive and multiplicative asymmetric stable Lévy fluctuations. By focusing on the impact of skewness (i.e., non-symmetry) in the probability distributions of noise, we find that via examining the mean first exit time (MFET) and the first escape probability (FEP), the asymmetric fluctuations, interacting with nonlinearity in the system, lead to peculiar likelihood for transcription. This includes, in the additive noise case, realizing higher likelihood of transcription for larger positive skewness (i.e., asymmetry) index β, causing a stochastic bifurcation at the non-Gaussianity index value α = 1 (i.e., it is a separating point or line for the likelihood for transcription), and achieving a turning point at the threshold value β≈-0.5 (i.e., beyond which the likelihood for transcription suddenly reversed for α values). The stochastic bifurcation and turning point phenomena do not occur in the symmetric noise case (β = 0). While in the multiplicative noise case, non-Gaussianity index value α = 1 is a separating point or line for both the MFET and the FEP. We also investigate the noise enhanced stability phenomenon. Additionally, we are able to specify the regions in the whole parameter space for the asymmetric noise, in which we attain desired likelihood for transcription. We have conducted a series of numerical experiments in "regulating" the likelihood of gene transcription by tuning asymmetric stable Lévy noise indexes. This work offers insights for possible ways of achieving gene regulation in experimental research.

Stochastic gene expression in Arabidopsis thaliana

Although plant development is highly reproducible, some stochasticity exists. This developmental stochasticity may be caused by noisy gene expression. Here we analyze the fluctuation of protein expression in Arabidopsis thaliana. Using the photoconvertible KikGR marker, we show that the protein expressions of individual cells fluctuate over time. A dual reporter system was used to study extrinsic and intrinsic noise of marker gene expression. We report that extrinsic noise is higher than intrinsic noise and that extrinsic noise in stomata is clearly lower in comparison to several other tissues/cell types. Finally, we show that cells are coupled with respect to stochastic protein expression in young leaves, hypocotyls and roots but not in mature leaves. Our data indicate that stochasticity of gene expression can vary between tissues/cell types and that it can be coupled in a non-cell-autonomous manner.

Noisy gene expression can cause stochasticity in the expression of plant traits. Here, Araújo et al. use a dual reporter system of protein expression in Arabidopsis to show that expression noise is lowest in stomata relative to other tissues and that leaf cells are coupled with respect to noise.

Translation initiation events on structured eukaryotic mRNAs generate gene expression noise

Gene expression stochasticity plays a major role in biology, creating non-genetic cellular individuality and influencing multiple processes, including differentiation and stress responses. We have addressed the lack of knowledge about posttranscriptional contributions to noise by determining cell-to-cell variations in the abundance of mRNA and reporter protein in yeast. Two types of structural element, a stem–loop and a poly(G) motif, not only inhibit translation initiation when inserted into an mRNA 5΄ untranslated region, but also generate noise. The noise-enhancing effect of the stem–loop structure also remains operational when combined with an upstream open reading frame. This has broad significance, since these elements are known to modulate the expression of a diversity of eukaryotic genes. Our findings suggest a mechanism for posttranscriptional noise generation that will contribute to understanding of the generally poor correlation between protein-level stochasticity and transcriptional bursting. We propose that posttranscriptional stochasticity can be linked to cycles of folding/unfolding of a stem–loop structure, or to interconversion between higher-order structural conformations of a G-rich motif, and have created a correspondingly configured computational model that generates fits to the experimental data. Stochastic events occurring during the ribosomal scanning process can therefore feature alongside transcriptional bursting as a source of noise.

Lack of 14-3-3 proteins in Saccharomyces cerevisiae results in cell-to-cell heterogeneity in the expression of Pho4-regulated genes SPL2 and PHO84

Background

Ion homeostasis is an essential property of living organisms. The yeast Saccharomyces cerevisiae is an ideal model organism to investigate ion homeostasis at all levels. In this yeast genes involved in high-affinity phosphate uptake (PHO genes) are strongly induced during both phosphate and potassium starvation, indicating a link between phosphate and potassium homeostasis. However, the signal transduction processes involved are not completely understood. As 14-3-3 proteins are key regulators of signal transduction processes, we investigated the effect of deletion of the 14-3-3 genes BMH1 or BMH2 on gene expression during potassium starvation and focused especially on the expression of genes involved in phosphate uptake.

Results

Genome-wide analysis of the effect of disruption of either BMH1 or BMH2 revealed that the mRNA levels of the PHO genes PHO84 and SPL2 are greatly reduced in the mutant strains compared to the levels in wild type strains. This was especially apparent at standard potassium and phosphate concentrations. Furthermore the promoter of these genes is less active after deletion of BMH1. Microscopic and flow cytometric analysis of cells with GFP-tagged SPL2 showed that disruption of BMH1 resulted in two populations of genetically identical cells, cells expressing the protein and the majority of cells with no detectible expression. Heterogeneity was also observed for the expression of GFP under control of the PHO84 promoter. Upon deletion of PHO80 encoding a regulator of the transcription factor Pho4, the effect of the BMH1 deletion on SPL2 and PHO84 promoter was lost, suggesting that the BMH1 deletion mainly influences processes upstream of the Pho4 transcription factor.

Conclusion

Our data indicate that that yeast cells can be in either of two states, expressing or not expressing genes required for high-affinity phosphate uptake and that 14-3-3 proteins are involved in the process(es) that establish the activation state of the PHO regulon.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4105-8) contains supplementary material, which is available to authorized users.

Optimal information transfer in enzymatic networks: A field theoretic formulation

Signaling in enzymatic networks is typically triggered by environmental fluctuations, resulting in a series of stochastic chemical reactions, leading to corruption of the signal by noise. For example, information flow is initiated by binding of extracellular ligands to receptors, which is transmitted through a cascade involving kinase-phosphatase stochastic chemical reactions. For a class of such networks, we develop a general field-theoretic approach to calculate the error in signal transmission as a function of an appropriate control variable. Application of the theory to a simple push-pull network, a module in the kinase-phosphatase cascade, recovers the exact results for error in signal transmission previously obtained using umbral calculus [Hinczewski and Thirumalai, Phys. Rev. X 4, 041017 (2014)2160-330810.1103/PhysRevX.4.041017]. We illustrate the generality of the theory by studying the minimal errors in noise reduction in a reaction cascade with two connected push-pull modules. Such a cascade behaves as an effective three-species network with a pseudointermediate. In this case, optimal information transfer, resulting in the smallest square of the error between the input and output, occurs with a time delay, which is given by the inverse of the decay rate of the pseudointermediate. Surprisingly, in these examples the minimum error computed using simulations that take nonlinearities and discrete nature of molecules into account coincides with the predictions of a linear theory. In contrast, there are substantial deviations between simulations and predictions of the linear theory in error in signal propagation in an enzymatic push-pull network for a certain range of parameters. Inclusion of second-order perturbative corrections shows that differences between simulations and theoretical predictions are minimized. Our study establishes that a field theoretic formulation of stochastic biological signaling offers a systematic way to understand error propagation in networks of arbitrary complexity.

Independent regulation of gene expression level and noise by histone modifications

The inherent stochasticity generates substantial gene expression variation among isogenic cells under identical conditions, which is frequently referred to as gene expression noise or cell-to-cell expression variability. Similar to (average) expression level, expression noise is also subject to natural selection. Yet it has been observed that noise is negatively correlated with expression level, which manifests as a potential constraint for simultaneous optimization of both. Here, we studied expression noise in human embryonic cells with computational analysis on single-cell RNA-seq data and in yeast with flow cytometry experiments. We showed that this coupling is overcome, to a certain degree, by a histone modification strategy in multiple embryonic developmental stages in human, as well as in yeast. Importantly, this epigenetic strategy could fit into a burst-like gene expression model: promoter-localized histone modifications (such as H3K4 methylation) are associated with both burst size and burst frequency, which together influence expression level, while gene-body-localized ones (such as H3K79 methylation) are more associated with burst frequency, which influences both expression level and noise. We further knocked out the only “writer” of H3K79 methylation in yeast, and observed that expression noise is indeed increased. Consistently, dosage sensitive genes, such as genes in the Wnt signaling pathway, tend to be marked with gene-body-localized histone modifications, while stress responding genes, such as genes regulating autophagy, tend to be marked with promoter-localized ones. Our findings elucidate that the “division of labor” among histone modifications facilitates the independent regulation of expression level and noise, extend the “histone code” hypothesis to include expression noise, and shed light on the optimization of transcriptome in evolution.

Gene expression noise, or cell-to-cell expression variability, has been a topic of intense interest for more than a decade. The prevailing model of “burst-like transcription” mediated by the promoter transitions between on and off states explains the formation of noise in eukaryotes. Albeit widely accepted, the cis- elements that determine the burst frequency and burst size remain largely unknown. Here we systematically examined the relationship between transcriptional burst frequency/size and all major histone modifications in various cell types, including human embryonic cells, mouse embryonic stem cells, and yeast, and found that histone markers can be divided into two groups based on their associations with burst frequency/ size. Coincidently, promoter-localized histone markers are associated with both burst size and burst frequency whereas gene-body-localized ones are more associated with burst frequency. We further knocked out a gene that is responsible for “writing” a gene-body histone mark in yeast, and found that burst frequency is indeed reduced. Our findings reveal a new mechanism of transcriptional burst regulation and shed light on the simultaneous optimization of gene expression level and noise in evolution.

Identification of potential new treatment response markers and therapeutic targets using a Gaussian process-based method in lapatinib insensitive breast cancer models

Molecularly targeted therapeutics hold promise of revolutionizing treatments of advanced malignancies. However, a large number of patients do not respond to these treatments. Here, we take a systems biology approach to understand the molecular mechanisms that prevent breast cancer (BC) cells from responding to lapatinib, a dual kinase inhibitor that targets human epidermal growth factor receptor 2 (HER2) and epidermal growth factor receptor (EGFR). To this end, we analysed temporal gene expression profiles of four BC cell lines, two of which respond and the remaining two do not respond to lapatinib. For this analysis, we developed a Gaussian process based algorithm which can accurately find differentially expressed genes by analysing time course gene expression profiles at a fraction of the computational cost of other state-of-the-art algorithms. Our analysis identified 519 potential genes which are characteristic of lapatinib non-responsiveness in the tested cell lines. Data from the Genomics of Drug Sensitivity in Cancer (GDSC) database suggested that the basal expressions 120 of the above genes correlate with the response of BC cells to HER2 and/or EGFR targeted therapies. We selected 27 genes from the larger panel of 519 genes for experimental verification and 16 of these were successfully validated. Further bioinformatics analysis identified vitamin D receptor (VDR) as a potential target of interest for lapatinib non-responsive BC cells. Experimentally, calcitriol, a commonly used reagent for VDR targeted therapy, in combination with lapatinib additively inhibited proliferation in two HER2 positive cell lines, lapatinib insensitive MDA-MB-453 and lapatinib resistant HCC 1954-L cells.

Single-Cell and Single-Molecule Analysis of Gene Expression Regulation

Recent advancements in single-cell and single-molecule imaging technologies have resolved biological processes in time and space that are fundamental to understanding the regulation of gene expression. Observations of single-molecule events in their cellular context have revealed highly dynamic aspects of transcriptional and post-transcriptional control in eukaryotic cells. This approach can relate transcription with mRNA abundance and lifetimes. Another key aspect of single-cell analysis is the cell-to-cell variability among populations of cells. Definition of heterogeneity has revealed stochastic processes, determined characteristics of under-represented cell types or transitional states, and integrated cellular behaviors in the context of multicellular organisms. In this review, we discuss novel aspects of gene expression of eukaryotic cells and multicellular organisms revealed by the latest advances in single-cell and single-molecule imaging technology.

Nondestructive nanostraw intracellular sampling for longitudinal cell monitoring

Here, we report a method for time-resolved, longitudinal extraction and quantitative measurement of intracellular proteins and mRNA from a variety of cell types. Cytosolic contents were repeatedly sampled from the same cell or population of cells for more than 5 d through a cell-culture substrate, incorporating hollow 150-nm-diameter nanostraws (NS) within a defined sampling region. Once extracted, the cellular contents were analyzed with conventional methods, including fluorescence, enzymatic assays (ELISA), and quantitative real-time PCR. This process was nondestructive with >95% cell viability after sampling, enabling long-term analysis. It is important to note that the measured quantities from the cell extract were found to constitute a statistically significant representation of the actual contents within the cells. Of 48 mRNA sequences analyzed from a population of cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs), 41 were accurately quantified. The NS platform samples from a select subpopulation of cells within a larger culture, allowing native cell-to-cell contact and communication even during vigorous activity such as cardiomyocyte beating. This platform was applied both to cell lines and to primary cells, including CHO cells, hiPSC-CMs, and human astrocytes derived in 3D cortical spheroids. By tracking the same cell or group of cells over time, this method offers an avenue to understand dynamic cell behavior, including processes such as induced pluripotency and differentiation.

Numerical Simulations Reveal Randomness of Cu(II) Induced Aβ Peptide Dimerization under Conditions Present in Glutamatergic Synapses

The interactions between the Aβ1-40 molecules species and the copper ions (Cu(II)) were intensively investigated due to their potential role in the development of the Alzheimer Disease (AD). The rate and the mechanism of the Cu(II)-Aβ complexes formation determines the aggregation pathway of the Aβ species, starting from smaller but more cytotoxic oligomers and ending up in large Aβ plaques, being the main hallmark of the AD. In our study we exploit the existing knowledge on the Cu(II)-Aβ interactions and create the theoretical model of the initial phase of the copper- driven Aβ aggregation mechanism. The model is based on the direct solution of the Chemical Master Equations, which capture the inherent stochastics of the considered system. In our work we argue that due to a strong Cu(II) affinity to Aβ and temporal accessibility of the Cu(II) ions during normal synaptic activity the aggregation driven by Cu(II) dominates the pure Aβ aggregation. We also demonstrate the dependence of the formation of different Cu(II)-Aβ complexes on the concentrations of reagents and the synaptic activity. Our findings correspond to recent experimental results and give a sound hypothesis on the AD development mechanisms.

Fluctuation sensitivity of a transcriptional signaling cascade

The internal biochemical state of a cell is regulated by a vast transcriptional network that kinetically correlates the concentrations of numerous proteins. Fluctuations in protein concentration that encode crucial information about this changing state must compete with fluctuations caused by the noisy cellular environment in order to successfully transmit information across the network. Oftentimes, one protein must regulate another through a sequence of intermediaries, and conventional wisdom, derived from the data processing inequality of information theory, leads us to expect that longer sequences should lose more information to noise. Using the metric of mutual information to characterize the fluctuation sensitivity of transcriptional signaling cascades, we find, counter to this expectation, that longer chains of regulatory interactions can instead lead to enhanced informational efficiency. We derive an analytic expression for the mutual information from a generalized chemical kinetics model that we reduce to simple, mass-action kinetics by linearizing for small fluctuations about the basal biological steady state, and we find that at long times this expression depends only on a simple ratio of protein production to destruction rates and the length of the cascade. We place bounds on the values of these parameters by requiring that the mutual information be at least one bit-otherwise, any received signal would be indistinguishable from noise-and we find not only that nature has devised a way to circumvent the data processing inequality, but that it must be circumvented to attain this one-bit threshold. We demonstrate how this result places informational and biochemical efficiency at odds with one another by correlating high transcription factor binding affinities with low informational output, and we conclude with an analysis of the validity of our assumptions and propose how they might be tested experimentally.

Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion

Quantitative mechanistic models are valuable tools for disentangling biochemical pathways and for achieving a comprehensive understanding of biological systems. However, to be quantitative the parameters of these models have to be estimated from experimental data. In the presence of significant stochastic fluctuations this is a challenging task as stochastic simulations are usually too time-consuming and a macroscopic description using reaction rate equations (RREs) is no longer accurate. In this manuscript, we therefore consider moment-closure approximation (MA) and the system size expansion (SSE), which approximate the statistical moments of stochastic processes and tend to be more precise than macroscopic descriptions. We introduce gradient-based parameter optimization methods and uncertainty analysis methods for MA and SSE. Efficiency and reliability of the methods are assessed using simulation examples as well as by an application to data for Epo-induced JAK/STAT signaling. The application revealed that even if merely population-average data are available, MA and SSE improve parameter identifiability in comparison to RRE. Furthermore, the simulation examples revealed that the resulting estimates are more reliable for an intermediate volume regime. In this regime the estimation error is reduced and we propose methods to determine the regime boundaries. These results illustrate that inference using MA and SSE is feasible and possesses a high sensitivity.

In this manuscript, we introduce efficient methods for parameter estimation for stochastic processes. The stochasticity of chemical reactions can influence the average behavior of the considered system. For some biological systems, a microscopic, stochastic description is computationally intractable but a macroscopic, deterministic description too inaccurate. This inaccuracy manifests itself in an error in parameter estimates, which impede the predictive power of the proposed model. Until now, no rigorous analysis on the magnitude of the estimation error exists. We show by means of two simulation examples that using mesoscopic descriptions based on the system size expansions and moment-closure approximations can reduce this estimation error compared to inference using a macroscopic description. This reduction is most pronounced in an intermediate volume regime where the influence of stochasticity on the average behavior is moderately strong. For the JAK/STAT pathway where experimental data is available, we show that one parameter that was not structurally identifiable when using a macroscopic description becomes structurally identifiable when using a mesoscopic description for parameter estimation.

Entrainment and Control of Bacterial Populations: An in Silico Study over a Spatially Extended Agent Based Model

We extend a spatially explicit agent based model (ABM) developed previously to investigate entrainment and control of the emergent behavior of a population of synchronized oscillating cells in a microfluidic chamber. Unlike most of the work in models of control of cellular systems which focus on temporal changes, we model individual cells with spatial dependencies which may contribute to certain behavioral responses. We use the model to investigate the response of both open loop and closed loop strategies, such as proportional control (P-control), proportional-integral control (PI-control) and proportional-integral-derivative control (PID-control), to heterogeinities and growth in the cell population, variations of the control parameters and spatial effects such as diffusion in the spatially explicit setting of a microfluidic chamber setup. We show that, as expected from the theory of phase locking in dynamical systems, open loop control can only entrain the cell population in a subset of forcing periods, with a wide variety of dynamical behaviors obtained outside these regions of entrainment. Closed-loop control is shown instead to guarantee entrainment in a much wider region of control parameter space although presenting limitations when the population size increases over a certain threshold. In silico tracking experiments are also performed to validate the ability of classical control approaches to achieve other reference behaviors such as a desired constant output or a linearly varying one. All simulations are carried out in BSim, an advanced agent-based simulator of microbial population which is here extended ad hoc to include the effects of control strategies acting onto the population.