What Population Reveals about Individual Cell Identity: Single-Cell Parameter Estimation of Models of Gene Expression in Yeast

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).

Mathematical modelling, selection and hierarchical inference to determine the minimal dose in IFNα therapy against myeloproliferative neoplasms

Myeloproliferative neoplasms (MPN) are blood cancers that appear after acquiring a driver mutation in a hematopoietic stem cell. These hematological malignancies result in the overproduction of mature blood cells and, if not treated, induce a risk of cardiovascular events and thrombosis. Pegylated IFN$\alpha $ is commonly used to treat MPN, but no clear guidelines exist concerning the dose prescribed to patients. We applied a model selection procedure and ran a hierarchical Bayesian inference method to decipher how dose variations impact the response to the therapy. We inferred that IFN$\alpha $ acts on mutated stem cells by inducing their differentiation into progenitor cells; the higher the dose, the higher the effect. We found that the treatment can induce long-term remission when a sufficient (patient-dependent) dose is reached. We determined this minimal dose for individuals in a cohort of patients and estimated the most suitable starting dose to give to a new patient to increase the chances of being cured.

Exploring carbon source related localization and phosphorylation in the Snf1/Mig1 network using population and single cell-based approaches

The AMPK/SNF1 pathway governs energy balance in eukaryotic cells, notably influencing glucose de-repression. In S. cerevisiae, Snf1 is phosphorylated and hence activated upon glucose depletion. This activation is required but is not sufficient for mediating glucose de-repression, indicating further glucose-dependent regulation mechanisms. Employing fluorescence recovery after photobleaching (FRAP) in conjunction with non-linear mixed effects modelling, we explore the spatial dynamics of Snf1 as well as the relationship between Snf1 phosphorylation and its target Mig1 controlled by hexose sugars. Our results suggest that inactivation of Snf1 modulates Mig1 localization and that the kinetic of Snf1 localization to the nucleus is modulated by the presence of non-fermentable carbon sources. Our data offer insight into the true complexity of regulation of this central signaling pathway in orchestrating cellular responses to fluctuating environmental cues. These insights not only expand our understanding of glucose homeostasis but also pave the way for further studies evaluating the importance of Snf1 localization in relation to its phosphorylation state and regulation of downstream targets.

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.

Modeling early treatment response in AML from cell-free tumor DNA

Monitoring disease response after intensive chemotherapy for acute myeloid leukemia (AML) currently requires invasive bone marrow biopsies, imposing a significant burden on patients. In contrast, cell-free tumor DNA (ctDNA) in peripheral blood, carrying tumor-specific mutations, offers a less-invasive assessment of residual disease. However, the relationship between ctDNA levels and bone marrow blast kinetics remains unclear. We explored this in 10 AML patients with NPM1 and IDH2 mutations undergoing initial chemotherapy. Comparison of mathematical mixed-effect models showed that (1) inclusion of blast cell death in the bone marrow, (2) transition of ctDNA to peripheral blood, and (3) ctDNA decay in peripheral blood describes kinetics of blast cells and ctDNA best. The fitted model allows prediction of residual bone marrow blast content from ctDNA, and its scaling factor, representing clonal heterogeneity, correlates with relapse risk. Our study provides precise insights into blast and ctDNA kinetics, offering novel avenues for AML disease monitoring.

Advances in automated and reactive flow cytometry for synthetic biotechnology

Automated flow cytometry (FC) has been initially considered for bioprocess monitoring and optimization. More recently, new physical and software interfaces have been made available, facilitating the access to this technology for labs and industries. It also comes with new capabilities, such as being able to act on the cultivation conditions based on population data. This approach, known as reactive FC, extended the range of applications of automated FC to bioprocess control and the stabilization of cocultures, but also to the broad field of synthetic and systems biology for the characterization of gene circuits. However, several issues must be addressed before automated and reactive FC can be considered standard and modular technologies.

Optogenetic control of beta-carotene bioproduction in yeast across multiple lab-scales

Optogenetics arises as a valuable tool to precisely control genetic circuits in microbial cell factories. Light control holds the promise of optimizing bioproduction methods and maximizing yields, but its implementation at different steps of the strain development process and at different culture scales remains challenging. In this study, we aim to control beta-carotene bioproduction using optogenetics in Saccharomyces cerevisiae and investigate how its performance translates across culture scales. We built four lab-scale illumination devices, each handling different culture volumes, and each having specific illumination characteristics and cultivating conditions. We evaluated optogenetic activation and beta-carotene production across devices and optimized them both independently. Then, we combined optogenetic induction and beta-carotene production to make a light-inducible beta-carotene producer strain. This was achieved by placing the transcription of the bifunctional lycopene cyclase/phytoene synthase CrtYB under the control of the pC120 optogenetic promoter regulated by the EL222-VP16 light-activated transcription factor, while other carotenogenic enzymes (CrtI, CrtE, tHMG) were expressed constitutively. We show that illumination, culture volume and shaking impact differently optogenetic activation and beta-carotene production across devices. This enabled us to determine the best culture conditions to maximize light-induced beta-carotene production in each of the devices. Our study exemplifies the stakes of scaling up optogenetics in devices of different lab scales and sheds light on the interplays and potential conflicts between optogenetic control and metabolic pathway efficiency. As a general principle, we propose that it is important to first optimize both components of the system independently, before combining them into optogenetic producing strains to avoid extensive troubleshooting. We anticipate that our results can help designing both strains and devices that could eventually lead to larger scale systems in an effort to bring optogenetics to the industrial scale.

Modeling Cellular Signaling Variability Based on Single-Cell Data: The TGFβ-SMAD Signaling Pathway

Nongenetic heterogeneity is key to cellular decisions, as even genetically identical cells respond in very different ways to the same external stimulus, e.g., during cell differentiation or therapeutic treatment of disease. Strong heterogeneity is typically already observed at the level of signaling pathways that are the first sensors of external inputs and transmit information to the nucleus where decisions are made. Since heterogeneity arises from random fluctuations of cellular components, mathematical models are required to fully describe the phenomenon and to understand the dynamics of heterogeneous cell populations. Here, we review the experimental and theoretical literature on cellular signaling heterogeneity, with special focus on the TGFβ/SMAD signaling pathway.

Resolving Crosstalk Between Signaling Pathways Using Mathematical Modeling and Time-Resolved Single Cell Data

Crosstalk between signaling pathways can modulate the cellular response to stimuli and is therefore an important part of signal transduction. For a comprehensive understanding of cellular responses, identifying points of interaction between the underlying molecular networks is essential. Here, we present an approach that allows the systematic prediction of such interactions by perturbing one pathway and quantifying the concomitant alterations in the response of a second pathway. As the observed alterations contain information about the crosstalk, we use an ordinary differential equation-based model to extract this information by linking altered dynamics to individual processes. Consequently, we can predict the interaction points between two pathways. As an example, we employed our approach to investigate the crosstalk between the NF-κB and p53 signaling pathway. We monitored the response of p53 to genotoxic stress using time-resolved single cell data and perturbed NF-κB signaling by inhibiting the kinase IKK2. Employing a subpopulation-based modeling approach enabled us to identify multiple interaction points that are simultaneously affected by perturbation of NF-κB signaling. Hence, our approach can be used to analyze crosstalk between two signaling pathways in a systematic manner.

CyberSco.Py an open-source software for event-based, conditional microscopy

Timelapse fluorescence microscopy imaging is routinely used in quantitative cell biology. However, microscopes could become much more powerful investigation systems if they were endowed with simple unsupervised decision-making algorithms to transform them into fully responsive and automated measurement devices. Here, we report CyberSco.Py, Python software for advanced automated timelapse experiments. We provide proof-of-principle of a user-friendly framework that increases the tunability and flexibility when setting up and running fluorescence timelapse microscopy experiments. Importantly, CyberSco.Py combines real-time image analysis with automation capability, which allows users to create conditional, event-based experiments in which the imaging acquisition parameters and the status of various devices can be changed automatically based on the image analysis. We exemplify the relevance of CyberSco.Py to cell biology using several use case experiments with budding yeast. We anticipate that CyberSco.Py could be used to address the growing need for smart microscopy systems to implement more informative quantitative cell biology experiments.

Surface-enhanced Raman scattering: An emerging tool for sensing cellular function

Continuous long-term intracellular imaging and multiplexed monitoring of biomolecular changes associated with key cellular processes remains a challenge for the scientific community. Recently, surface-enhanced Raman scattering (SERS) has been demonstrated as a powerful spectroscopic tool in the field of biology owing to its significant advantages. Some of these include the ability to provide molecule-specific information with exquisite sensitivity, working with small volumes of precious samples, real-time monitoring, and optimal optical contrast. More importantly, the availability of a large number of novel Raman reporters with narrower full width at half maximum (FWHM) of spectral peaks/vibrational modes than conventional fluorophores has created a versatile palette of SERS-based probes that allow targeted multiplex sensing surpassing the detection sensitivity of even fluorescent probes. Due to its nondestructive nature, its applicability has been recognized for biological sensing, molecular imaging, and dynamic monitoring of complex intracellular processes. We critically discuss recent developments in this area with a focus on different applications where SERS has been used for obtaining information that remains elusive for conventional imaging methods. Current reports indicate that SERS has made significant inroads in the field of biology and has the potential to be used for in vivo human applications. This article is categorized under: Diagnostic Tools > In Vitro Nanoparticle-Based Sensing Nanotechnology Approaches to Biology > Nanoscale Systems in Biology Diagnostic Tools > Biosensing Diagnostic Tools > In Vivo Nanodiagnostics and Imaging.

One Cell at a Time: Advances in Single-Cell Methods and Instrumentation for Discovery in Aquatic Microbiology

Studying microbes from a single-cell perspective has become a major theme and interest within the field of aquatic microbiology. One emerging trend is the unfailing observation of heterogeneity in activity levels within microbial populations. Wherever researchers have looked, intra-population variability in biochemical composition, growth rates, and responses to varying environmental conditions has been evident and probably reflect coexisting genetically distinct strains of the same species. Such observations of heterogeneity require a shift away from bulk analytical approaches and development of new methods or adaptation of existing techniques, many of which were first pioneered in other, unrelated fields, e.g., material, physical, and biomedical sciences. Many co-opted approaches were initially optimized using model organisms. In a field with so few cultivable models, method development has been challenging but has also contributed tremendous insights, breakthroughs, and stimulated curiosity. In this perspective, we present a subset of methods that have been effectively applied to study aquatic microbes at the single-cell level. Opportunities and challenges for innovation are also discussed. We suggest future directions for aquatic microbiological research that will benefit from open access to sophisticated instruments and highly interdisciplinary collaborations.

Scalable and flexible inference framework for stochastic dynamic single-cell models

Understanding the inherited nature of how biological processes dynamically change over time and exhibit intra- and inter-individual variability, due to the different responses to environmental stimuli and when interacting with other processes, has been a major focus of systems biology. The rise of single-cell fluorescent microscopy has enabled the study of those phenomena. The analysis of single-cell data with mechanistic models offers an invaluable tool to describe dynamic cellular processes and to rationalise cell-to-cell variability within the population. However, extracting mechanistic information from single-cell data has proven difficult. This requires statistical methods to infer unknown model parameters from dynamic, multi-individual data accounting for heterogeneity caused by both intrinsic (e.g. variations in chemical reactions) and extrinsic (e.g. variability in protein concentrations) noise. Although several inference methods exist, the availability of efficient, general and accessible methods that facilitate modelling of single-cell data, remains lacking. Here we present a scalable and flexible framework for Bayesian inference in state-space mixed-effects single-cell models with stochastic dynamic. Our approach infers model parameters when intrinsic noise is modelled by either exact or approximate stochastic simulators, and when extrinsic noise is modelled by either time-varying, or time-constant parameters that vary between cells. We demonstrate the relevance of our approach by studying how cell-to-cell variation in carbon source utilisation affects heterogeneity in the budding yeast Saccharomyces cerevisiae SNF1 nutrient sensing pathway. We identify hexokinase activity as a source of extrinsic noise and deduce that sugar availability dictates cell-to-cell variability.

Parameter inference for stochastic biochemical models from perturbation experiments parallelised at the single cell level

Understanding and characterising biochemical processes inside single cells requires experimental platforms that allow one to perturb and observe the dynamics of such processes as well as computational methods to build and parameterise models from the collected data. Recent progress with experimental platforms and optogenetics has made it possible to expose each cell in an experiment to an individualised input and automatically record cellular responses over days with fine time resolution. However, methods to infer parameters of stochastic kinetic models from single-cell longitudinal data have generally been developed under the assumption that experimental data is sparse and that responses of cells to at most a few different input perturbations can be observed. Here, we investigate and compare different approaches for calculating parameter likelihoods of single-cell longitudinal data based on approximations of the chemical master equation (CME) with a particular focus on coupling the linear noise approximation (LNA) or moment closure methods to a Kalman filter. We show that, as long as cells are measured sufficiently frequently, coupling the LNA to a Kalman filter allows one to accurately approximate likelihoods and to infer model parameters from data even in cases where the LNA provides poor approximations of the CME. Furthermore, the computational cost of filtering-based iterative likelihood evaluation scales advantageously in the number of measurement times and different input perturbations and is thus ideally suited for data obtained from modern experimental platforms. To demonstrate the practical usefulness of these results, we perform an experiment in which single cells, equipped with an optogenetic gene expression system, are exposed to various different light-input sequences and measured at several hundred time points and use parameter inference based on iterative likelihood evaluation to parameterise a stochastic model of the system.

Two-level modeling approach to identify the regulatory dynamics capturing drug response heterogeneity in single-cells

Single-cell multimodal technologies reveal the scales of cellular heterogeneity impairing cancer treatment, yet cell response dynamics remain largely underused to decipher the mechanisms of drug resistance they take part in. As the phenotypic heterogeneity of a clonal cell population informs on the capacity of each single-cell to recapitulate the whole range of observed behaviors, we developed a modeling approach utilizing single-cell response data to identify regulatory reactions driving population heterogeneity in drug response. Dynamic data of hundreds of HeLa cells treated with TNF-related apoptosis-inducing ligand (TRAIL) were used to characterize the fate-determining kinetic parameters of an apoptosis receptor reaction model. Selected reactions sets were augmented to incorporate a mechanism that leads to the separation of the opposing response phenotypes. Using a positive feedback loop motif to identify the reaction set, we show that caspase-8 is able to encapsulate high levels of heterogeneity by introducing a response delay and amplifying the initial differences arising from natural protein expression variability. Our approach enables the identification of fate-determining reactions that drive the population response heterogeneity, providing regulatory targets to curb the cell dynamics of drug resistance.

Uncertainty propagation for deterministic models of biochemical networks using moment equations and the extended Kalman filter

Differential equation models of biochemical networks are frequently associated with a large degree of uncertainty in parameters and/or initial conditions. However, estimating the impact of this uncertainty on model predictions via Monte Carlo simulation is computationally demanding. A more efficient approach could be to track a system of low-order statistical moments of the state. Unfortunately, when the underlying model is nonlinear, the system of moment equations is infinite-dimensional and cannot be solved without a moment closure approximation which may introduce bias in the moment dynamics. Here, we present a new method to study the time evolution of the desired moments for nonlinear systems with polynomial rate laws. Our approach is based on solving a system of low-order moment equations by substituting the higher-order moments with Monte Carlo-based estimates from a small number of simulations, and using an extended Kalman filter to counteract Monte Carlo noise. Our algorithm provides more accurate and robust results compared to traditional Monte Carlo and moment closure techniques, and we expect that it will be widely useful for the quantification of uncertainty in biochemical model predictions.

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.

Stochastic models for single-cell data: Current challenges and the way forward

Although the quantity and quality of single-cell data have progressed rapidly, making quantitative predictions with single-cell stochastic models remains challenging. The stochastic nature of cellular processes leads to at least three challenges in building models with single-cell data: (a) because variability in single-cell data can be attributed to multiple different sources, it is difficult to rule out conflicting mechanistic models that explain the same data equally well; (b) the distinction between interesting biological variability and experimental variability is sometimes ambiguous; (c) the nonstandard distributions of single-cell data can lead to violations of the assumption of symmetric errors in least-squares fitting. In this review, we first discuss recent studies that overcome some of the challenges or set up a promising direction and then introduce some powerful statistical approaches utilized in these studies. We conclude that applying and developing statistical approaches could lead to further progress in building stochastic models for single-cell data.

Modeling and characterization of inter-individual variability in CD8 T cell responses in mice

To develop vaccines it is mandatory yet challenging to account for inter-individual variability during immune responses. Even in laboratory mice, T cell responses of single individuals exhibit a high heterogeneity that may come from genetic backgrounds, intra-specific processes (e.g. antigen-processing and presentation) and immunization protocols.

To account for inter-individual variability in CD8 T cell responses in mice, we propose a dynamical model coupled to a statistical, nonlinear mixed effects model. Average and individual dynamics during a CD8 T cell response are characterized in different immunization contexts (vaccinia virus and tumor). On one hand, we identify biological processes that generate inter-individual variability (activation rate of naive cells, the mortality rate of effector cells, and dynamics of the immunogen). On the other hand, introducing categorical covariates to analyze two different immunization regimens, we highlight the steps of the response impacted by immunogens (priming, differentiation of naive cells, expansion of effector cells and generation of memory cells). The robustness of the model is assessed by confrontation to new experimental data.

Our approach allows to investigate immune responses in various immunization contexts, when measurements are scarce or missing, and contributes to a better understanding of inter-individual variability in CD8 T cell immune responses.

Fine-Tuning of Energy Levels Regulates SUC2 via a SNF1-Dependent Feedback Loop

Nutrient sensing pathways are playing an important role in cellular response to different energy levels. In budding yeast, Saccharomyces cerevisiae, the sucrose non-fermenting protein kinase complex SNF1 is a master regulator of energy homeostasis. It is affected by multiple inputs, among which energy levels is the most prominent. Cells which are exposed to a switch in carbon source availability display a change in the gene expression machinery. It has been shown that the magnitude of the change varies from cell to cell. In a glucose rich environment Snf1/Mig1 pathway represses the expression of its downstream target, such as SUC2. However, upon glucose depletion SNF1 is activated which leads to an increase in SUC2 expression. Our single cell experiments indicate that upon starvation, gene expression pattern of SUC2 shows rapid increase followed by a decrease to initial state with high cell-to-cell variability. The mechanism behind this behavior is currently unknown. In this work we study the long-term behavior of the Snf1/Mig1 pathway upon glucose starvation with a microfluidics and non-linear mixed effect modeling approach. We show a negative feedback mechanism, involving Snf1 and Reg1, which reduces SUC2 expression after the initial strong activation. Snf1 kinase activity plays a key role in this feedback mechanism. Our systems biology approach proposes a negative feedback mechanism that works through the SNF1 complex and is controlled by energy levels. We further show that Reg1 likely is involved in the negative feedback mechanism.

A systematic approach to decipher crosstalk in the p53 signaling pathway using single cell dynamics

The transcription factors NF-κB and p53 are key regulators in the genotoxic stress response and are critical for tumor development. Although there is ample evidence for interactions between both networks, a comprehensive understanding of the crosstalk is lacking. Here, we developed a systematic approach to identify potential interactions between the pathways. We perturbed NF-κB signaling by inhibiting IKK2, a critical regulator of NF-κB activity, and monitored the altered response of p53 to genotoxic stress using single cell time lapse microscopy. Fitting subpopulation-specific computational p53 models to this time-resolved single cell data allowed to reproduce in a quantitative manner signaling dynamics and cellular heterogeneity for the unperturbed and perturbed conditions. The approach enabled us to untangle the integrated effects of IKK/ NF-κB perturbation on p53 dynamics and thereby derive potential interactions between both networks. Intriguingly, we find that a simultaneous perturbation of multiple processes is necessary to explain the observed changes in the p53 response. Specifically, we show interference with the activation and degradation of p53 as well as the degradation of Mdm2. Our results highlight the importance of the crosstalk and its potential implications in p53-dependent cellular functions.

Cells can respond to external and internal inputs by transducing information to the nucleus where transcription factors initiate corresponding cellular responses. Cellular signaling is mediated by several pathways; molecular networks that can interact with each other, which alters signal processing and modulates cellular responses. As deregulated signaling can lead to the development of tumors it is important to understand not only how signaling pathways function but also the contribution of their interaction on the signaling dynamics. Here, we analyzed the interplay of the IKK/ NF-κB and p53 pathway, which are both critical for the cellular response to DNA damage and have been implicated in tumor development. To systematically identify interaction points between both pathways, we perturbed IKK/ NF-κB signaling and tracked the changes in the response of p53 to DNA damage. Using computational methods, we show that several reactions in the p53 pathway are simultaneously affected by NF-κB signaling and that this combined action is necessary to explain altered behaviour of the p53 pathway. Hence, our results provide new insights into the interplay between the NF-κB and p53 pathway and help to gain a more comprehensive understanding of the crosstalk.

Random Parametric Perturbations of Gene Regulatory Circuit Uncover State Transitions in Cell Cycle

Many biological processes involve precise cellular state transitions controlled by complex gene regulation. Here, we use budding yeast cell cycle as a model system and explore how a gene regulatory circuit encodes essential information of state transitions. We present a generalized random circuit perturbation method for circuits containing heterogeneous regulation types and its usage to analyze both steady and oscillatory states from an ensemble of circuit models with random kinetic parameters. The stable steady states form robust clusters with a circular structure that are associated with cell cycle phases. This circular structure in the clusters is consistent with single-cell RNA sequencing data. The oscillatory states specify the irreversible state transitions along cell cycle progression. Furthermore, we identify possible mechanisms to understand the irreversible state transitions from the steady states. We expect this approach to be robust and generally applicable to unbiasedly predict dynamical transitions of a gene regulatory circuit.

Inheritance and variability of kinetic gene expression parameters in microbial cells: modeling and inference from lineage tree data

Motivation

Modern experimental technologies enable monitoring of gene expression dynamics in individual cells and quantification of its variability in isogenic microbial populations. Among the sources of this variability is the randomness that affects inheritance of gene expression factors at cell division. Known parental relationships among individually observed cells provide invaluable information for the characterization of this extrinsic source of gene expression noise. Despite this fact, most existing methods to infer stochastic gene expression models from single-cell data dedicate little attention to the reconstruction of mother–daughter inheritance dynamics.

Results

Starting from a transcription and translation model of gene expression, we propose a stochastic model for the evolution of gene expression dynamics in a population of dividing cells. Based on this model, we develop a method for the direct quantification of inheritance and variability of kinetic gene expression parameters from single-cell gene expression and lineage data. We demonstrate that our approach provides unbiased estimates of mother–daughter inheritance parameters, whereas indirect approaches using lineage information only in the post-processing of individual-cell parameters underestimate inheritance. Finally, we show on yeast osmotic shock response data that daughter cell parameters are largely determined by the mother, thus confirming the relevance of our method for the correct assessment of the onset of gene expression variability and the study of the transmission of regulatory factors.

Availability and implementation

Software code is available at https://github.com/almarguet/IdentificationWithARME. Lineage tree data is available upon request.

Supplementary information

Supplementary material is available at Bioinformatics online.

The Maximum Caliber Variational Principle for Nonequilibria

Ever since Clausius in 1865 and Boltzmann in 1877, the concepts of entropy and of its maximization have been the foundations for predicting how material equilibria derive from microscopic properties. But, despite much work, there has been no equally satisfactory general variational principle for nonequilibrium situations. However, in 1980, a new avenue was opened by E.T. Jaynes and by Shore and Johnson. We review here maximum caliber, which is a maximum-entropy-like principle that can infer distributions of flows over pathways, given dynamical constraints. This approach is providing new insights, particularly into few-particle complex systems, such as gene circuits, protein conformational reaction coordinates, network traffic, bird flocking, cell motility, and neuronal firing.

Maximum Entropy Framework for Predictive Inference of Cell Population Heterogeneity and Responses in Signaling Networks

Predictive models of signaling networks are essential for understanding cell population heterogeneity and designing rational interventions in disease. However, using computational models to predict heterogeneity of signaling dynamics is often challenging because of the extensive variability of biochemical parameters across cell populations. Here, we describe a maximum entropy-based framework for inference of heterogeneity in dynamics of signaling networks (MERIDIAN). MERIDIAN estimates the joint probability distribution over signaling network parameters that is consistent with experimentally measured cell-to-cell variability of biochemical species. We apply the developed approach to investigate the response heterogeneity in the EGFR/Akt signaling network. Our analysis demonstrates that a significant fraction of cells exhibits high phosphorylated Akt (pAkt) levels hours after EGF stimulation. Our findings also suggest that cells with high EGFR levels predominantly contribute to the subpopulation of cells with high pAkt activity. We also discuss how MERIDIAN can be extended to accommodate various experimental measurements.

Estimating numbers of intracellular molecules through analysing fluctuations in photobleaching

The impact of fluorescence microscopy has been limited by the difficulties of expressing measurements of fluorescent proteins in numbers of molecules. Absolute numbers enable the integration of results from different laboratories, empower mathematical modelling, and are the bedrock for a quantitative, predictive biology. Here we propose an estimator to infer numbers of molecules from fluctuations in the photobleaching of proteins tagged with Green Fluorescent Protein. Performing experiments in budding yeast, we show that our estimates of numbers agree, within an order of magnitude, with published biochemical measurements, for all six proteins tested. The experiments we require are straightforward and use only a wide-field fluorescence microscope. As such, our approach has the potential to become standard for those practising quantitative fluorescence microscopy.

Of Gene Expression and Cell Division Time: A Mathematical Framework for Advanced Differential Gene Expression and Data Analysis

Estimating fold changes of average mRNA and protein molecule counts per cell is the most common way to perform differential expression analysis. However, these gene expression data may be affected by cell division, an often-neglected phenomenon. Here, we develop a quantitative framework that links population-based mRNA and protein measurements to rates of gene expression in single cells undergoing cell division. The equations we derive are easy-to-use and widely robust against biological variability. They integrate multiple "omics" data into a coherent, quantitative description of single-cell gene expression and improve analysis when comparing systems or states with different cell division times. We explore these ideas in the context of resting versus activated B cells. Analyzing differences in protein synthesis rates enables to account for differences in cell division times. We demonstrate that this improves the resolution and hit rate of differential gene expression analysis when compared to analyzing population protein abundances alone.

Landscape perspectives of tumor, EMT, and development

A tumor is rarely fatal until becoming metastatic. Recent discoveries suggest that epithelial mesenchymal transition(EMT) is an important process which contributes to not only cancer metastasis but also increased stemness. Cancer cells with stem cell characteristics are called cancer stem cells (CSCs). We review recent efforts to quantify and delineate the relationship among EMT, CSC and tumor development. When the gene regulatory network is tightly regulated through the effectively fast regulatory binding, Cancer, Premalignant, Normal, CSC, stem cell (SC), Lesion and Hyperplasia states emerged. The corresponding landscape topography for all of these states can be quantified to a global way for uncovering the relationship among the tumor, metastasis, and development. On the other hand, phenotypic and functional heterogeneity is regarded as one of the greatest challenge in cancer treatment. Cancer and CSCs are heterogeneous and give rise to tumorigenic and non-tumorigenic cells during self-renewal, differentiation and epigenetic diversification. Further, if the gene regulatory network is weakly regulated through the effective slow regulatory binding (by DNA methylation or histone modification etc), multiple meta-stable states can emerge. This model can provide an epigenetic and physical rather than genetic and fixed origin of heterogeneity. Elucidating the origin of and dynamic nature of tumor cells will likely help better understand the cellular basis of therapeutic response, resistance, and relapse.

A microfluidic device for inferring metabolic landscapes in yeast monolayer colonies

Microbial colonies are fascinating structures in which growth and internal organization reflect complex morphogenetic processes. Here, we generated a microfluidics device with arrays of long monolayer yeast colonies to further global understanding of how intercellular metabolic interactions affect the internal structure of colonies within defined boundary conditions. We observed the emergence of stable glucose gradients using fluorescently labeled hexose transporters and quantified the spatial correlations with intra-colony growth rates and expression of other genes regulated by glucose availability. These landscapes depended on the external glucose concentration as well as secondary gradients, for example amino acid availability. This work demonstrates the regulatory genetic networks governing cellular physiological adaptation are the key to internal structuration of cellular assemblies. This approach could be used in the future to decipher the interplay between long-range metabolic interactions, cellular development and morphogenesis in more complex systems.

Hyperosmotic Stress Response Memory is Modulated by Gene Positioning in Yeast

Cellular memory is a critical ability that allows microorganisms to adapt to potentially detrimental environmental fluctuations. In the unicellular eukaryote Saccharomyces cerevisiae, cellular memory can take the form of faster or slower responses within the cell population to repeated stresses. Using microfluidics and fluorescence time-lapse microscopy, we studied how yeast responds to short, pulsed hyperosmotic stresses at the single-cell level by analyzing the dynamic behavior of the stress-responsive STL1 promoter (pSTL1) fused to a fluorescent reporter. We established that pSTL1 exhibits variable successive activation patterns following two repeated short stresses. Despite this variability, most cells exhibited a memory of the first stress as decreased pSTL1 activity in response to the second stress. Notably, we showed that genomic location is important for the memory effect, since displacement of the promoter to a pericentromeric chromatin domain decreased the transcriptional strength of pSTL1 and led to a loss of memory. This study provides a quantitative description of a cellular memory that includes single-cell variability and highlights the contribution of chromatin structure to stress memory.

Optimal parameter identification of synthetic gene networks using harmony search algorithm

Computational modeling of engineered gene circuits is an important while challenged task in systems biology. In order to describe and predict the response behaviors of genetic circuits using reliable model parameters, this paper applies an optimal experimental design(OED) method to obtain input signals. In order to obtain informative observations, this study focuses on maximizing Fisher information matrix(FIM)-based optimal criteria and to provide optimal inputs. Furthermore, this paper designs a two-stage optimization with the modified E-optimal criteria and applies harmony search(HS)-based OED algorithm to minimize estimation errors. The proposed optimal identification methodology involves estimation errors and the sample size to pursue a trade-off between estimation accuracy and measurement cost in modeling gene networks. The designed cost function takes two major factors into account, in which experimental costs are proportional to the number of time points. Experiments select two types of synthetic genetic networks to validate the effectiveness of the proposed HS-OED approach. Identification outcomes and analysis indicate the proposed HS-OED method outperforms two candidate OED approaches, with reduced computational effort.

A Simple and Flexible Computational Framework for Inferring Sources of Heterogeneity from Single-Cell Dynamics

Single-cell time-lapse data provide the means for disentangling sources of cell-to-cell and intra-cellular variability, a key step for understanding heterogeneity in cell populations. However, single-cell analysis with dynamic models is a challenging open problem: current inference methods address only single-gene expression or neglect parameter correlations. We report on a simple, flexible, and scalable method for estimating cell-specific and population-average parameters of non-linear mixed-effects models of cellular networks, demonstrating its accuracy with a published model and dataset. We also propose sensitivity analysis for identifying which biological sub-processes quantitatively and dynamically contribute to cell-to-cell variability. Our application to endocytosis in yeast demonstrates that dynamic models of realistic size can be developed for the analysis of single-cell data and that shifting the focus from single reactions or parameters to nuanced and time-dependent contributions of sub-processes helps biological interpretation. Generality and simplicity of the approach will facilitate customized extensions for analyzing single-cell dynamics.

Multi-experiment nonlinear mixed effect modeling of single-cell translation kinetics after transfection

Single-cell time-lapse studies have advanced the quantitative understanding of cellular pathways and their inherent cell-to-cell variability. However, parameters retrieved from individual experiments are model dependent and their estimation is limited, if based on solely one kind of experiment. Hence, methods to integrate data collected under different conditions are expected to improve model validation and information content. Here we present a multi-experiment nonlinear mixed effect modeling approach for mechanistic pathway models, which allows the integration of multiple single-cell perturbation experiments. We apply this approach to the translation of green fluorescent protein after transfection using a massively parallel read-out of micropatterned single-cell arrays. We demonstrate that the integration of data from perturbation experiments allows the robust reconstruction of cell-to-cell variability, i.e., parameter densities, while each individual experiment provides insufficient information. Indeed, we show that the integration of the datasets on the population level also improves the estimates for individual cells by breaking symmetries, although each of them is only measured in one experiment. Moreover, we confirmed that the suggested approach is robust with respect to batch effects across experimental replicates and can provide mechanistic insights into the nature of batch effects. We anticipate that the proposed multi-experiment nonlinear mixed effect modeling approach will serve as a basis for the analysis of cellular heterogeneity in single-cell dynamics.

Role of noise and parametric variation in the dynamics of gene regulatory circuits

Stochasticity in gene expression impacts the dynamics and functions of gene regulatory circuits. Intrinsic noises, including those that are caused by low copy number of molecules and transcriptional bursting, are usually studied by stochastic simulations. However, the role of extrinsic factors, such as cell-to-cell variability and heterogeneity in the microenvironment, is still elusive. To evaluate the effects of both the intrinsic and extrinsic noises, we develop a method, named sRACIPE, by integrating stochastic analysis with random circuit perturbation (RACIPE) method. RACIPE uniquely generates and analyzes an ensemble of models with random kinetic parameters. Previously, we have shown that the gene expression from random models form robust and functionally related clusters. In sRACIPE we further develop two stochastic simulation schemes, aiming to reduce the computational cost without sacrificing the convergence of statistics. One scheme uses constant noise to capture the basins of attraction, and the other one uses simulated annealing to detect the stability of states. By testing the methods on several synthetic gene regulatory circuits and an epithelial-mesenchymal transition network in squamous cell carcinoma, we demonstrate that sRACIPE can interpret the experimental observations from single-cell gene expression data. We observe that parametric variation (the spread of parameters around a median value) increases the spread of the gene expression clusters, whereas high noise merges the states. Our approach quantifies the robustness of a gene circuit in the presence of noise and sheds light on a new mechanism of noise-induced hybrid states. We have implemented sRACIPE as an R package.

Fitness effects of altering gene expression noise in Saccharomyces cerevisiae

Gene expression noise is an evolvable property of biological systems that describes differences in expression among genetically identical cells in the same environment. Prior work has shown that expression noise is heritable and can be shaped by selection, but the impact of variation in expression noise on organismal fitness has proven difficult to measure. Here, we quantify the fitness effects of altering expression noise for the TDH3 gene in Saccharomyces cerevisiae. We show that increases in expression noise can be deleterious or beneficial depending on the difference between the average expression level of a genotype and the expression level maximizing fitness. We also show that a simple model relating single-cell expression levels to population growth produces patterns consistent with our empirical data. We use this model to explore a broad range of average expression levels and expression noise, providing additional insight into the fitness effects of variation in expression noise.

RACIPE: a computational tool for modeling gene regulatory circuits using randomization

Background

One of the major challenges in traditional mathematical modeling of gene regulatory circuits is the insufficient knowledge of kinetic parameters. These parameters are often inferred from existing experimental data and/or educated guesses, which can be time-consuming and error-prone, especially for large networks.

Results

We present a user-friendly computational tool for the community to use our newly developed method named random circuit perturbation (RACIPE), to explore the robust dynamical features of gene regulatory circuits without the requirement of detailed kinetic parameters. Taking the network topology as the only input, RACIPE generates an ensemble of circuit models with distinct randomized parameters and uniquely identifies robust dynamical properties by statistical analysis. Here, we discuss the implementation of the software and the statistical analysis methods of RACIPE-generated data to identify robust gene expression patterns and the functions of genes and regulatory links. Finally, we apply the tool on coupled toggle-switch circuits and a published circuit of B-lymphopoiesis.

Conclusions

We expect our new computational tool to contribute to a more comprehensive and unbiased understanding of mechanisms underlying gene regulatory networks. RACIPE is a free open source software distributed under (Apache 2.0) license and can be downloaded from GitHub (https://github.com/simonhb1990/RACIPE-1.0).

Electronic supplementary material

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

Genomics of cellular proliferation in periodic environmental fluctuations

Living systems control cell growth dynamically by processing information from their environment. Although responses to a single environmental change have been intensively studied, little is known about how cells react to fluctuating conditions. Here, we address this question at the genomic scale by measuring the relative proliferation rate (fitness) of 3,568 yeast gene deletion mutants in out-of-equilibrium conditions: periodic oscillations between two environmental conditions. In periodic salt stress, fitness and its genetic variance largely depended on the oscillating period. Surprisingly, dozens of mutants displayed pronounced hyperproliferation under short stress periods, revealing unexpected controllers of growth under fast dynamics. We validated the implication of the high-affinity cAMP phosphodiesterase and of a regulator of protein translocation to mitochondria in this group. Periodic oscillations of extracellular methionine, a factor unrelated to salinity, also altered fitness but to a lesser extent and for different genes. The results illustrate how natural selection acts on mutations in a dynamic environment, highlighting unsuspected genetic vulnerabilities to periodic stress in molecular processes that are conserved across all eukaryotes.

Assigning function to natural allelic variation via dynamic modeling of gene network induction

More and more natural DNA variants are being linked to physiological traits. Yet, understanding what differences they make on molecular regulations remains challenging. Important properties of gene regulatory networks can be captured by computational models. If model parameters can be “personalized” according to the genotype, their variation may then reveal how DNA variants operate in the network. Here, we combined experiments and computations to visualize natural alleles of the yeast GAL3 gene in a space of model parameters describing the galactose response network. Alleles altering the activation of Gal3p by galactose were discriminated from those affecting its activity (production/degradation or efficiency of the activated protein). The approach allowed us to correctly predict that a non‐synonymous SNP would change the binding affinity of Gal3p with the Gal80p transcriptional repressor. Our results illustrate how personalizing gene regulatory models can be used for the mechanistic interpretation of genetic variants.

Single-cell systems biology: probing the basic unit of information flow

Gene expression varies across cells in a population or a tissue. This heterogeneity has come into sharp focus in recent years through developments in new imaging and sequencing technologies. However, our ability to measure variation has outpaced our ability to interpret it. Much of the variability may arise from random effects occurring in the processes of gene expression (transcription, RNA processing and decay, translation). The molecular basis of these effects is largely unknown. Likewise, a functional role of this variability in growth, differentiation and disease has only been elucidated in a few cases. In this review, we highlight recent experimental and theoretical advances for measuring and analyzing stochastic variation.

Taking control over microbial populations: Current approaches for exploiting biological noise in bioprocesses

Phenotypic plasticity of microbial cells has attracted much attention and several research efforts have been dedicated to the description of methods aiming at characterizing phenotypic heterogeneity and its impact on microbial populations. However, different approaches have also been suggested in order to take benefit from noise in a bioprocess perspective, e.g. by increasing the robustness or productivity of a microbial population. This review is dedicated to outline these controlling methods. A common issue, that has still to be addressed, is the experimental identification and the mathematical expression of noise. Indeed, the effective interfacing of microbial physiology with external parameters that can be used for controlling physiology depends on the acquisition of reliable signals. Latest technologies, like single cell microfluidics and advanced flow cytometric approaches, enable linking physiology, noise, heterogeneity in productive microbes with environmental cues and hence allow correctly mapping and predicting biological behavior via mathematical representations. However, like in the field of electronics, signals are perpetually subjected to noise. If appropriately interpreted, this noise can give an additional insight into the behavior of the individual cells within a microbial population of interest. This review focuses on recent progress made at describing, treating and exploiting biological noise in the context of microbial populations used in various bioprocess applications.

Parameter inference for stochastic single-cell dynamics from lineage tree data

Background

With the advance of experimental techniques such as time-lapse fluorescence microscopy, the availability of single-cell trajectory data has vastly increased, and so has the demand for computational methods suitable for parameter inference with this type of data. Most of currently available methods treat single-cell trajectories independently, ignoring the mother-daughter relationships and the information provided by the population structure. However, this information is essential if a process of interest happens at cell division, or if it evolves slowly compared to the duration of the cell cycle.

Results

In this work, we propose a Bayesian framework for parameter inference on single-cell time-lapse data from lineage trees. Our method relies on a combination of Sequential Monte Carlo for approximating the parameter likelihood function and Markov Chain Monte Carlo for parameter exploration. We demonstrate our inference framework on two simple examples in which the lineage tree information is crucial: one in which the cell phenotype can only switch at cell division and another where the cell state fluctuates slowly over timescales that extend well beyond the cell-cycle duration.

Conclusion

There exist several examples of biological processes, such as stem cell fate decisions or epigenetically controlled phase variation in bacteria, where the cell ancestry is expected to contain important information about the underlying system dynamics. Parameter inference methods that discard this information are expected to perform poorly for such type of processes. Our method provides a simple and computationally efficient way to take into account single-cell lineage tree data for the purpose of parameter inference and serves as a starting point for the development of more sophisticated and powerful approaches in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0425-1) contains supplementary material, which is available to authorized users.

Interrogating the topological robustness of gene regulatory circuits by randomization

One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.

Cells are able to robustly carry out their essential biological functions, possibly because of multiple layers of tight regulation via complex, yet well-designed, gene regulatory networks involving a substantial number of genes. State-of-the-art genomics technology has enabled the mapping of these large gene networks, yet it remains a tremendous challenge to elucidate their design principles and the regulatory mechanisms underlying their biological functions such as signal processing and decision-making. One of the key barriers is the absence of accurate kinetics for the regulatory interactions, especially from in vivo experiments. To this end, we have developed a new computational modeling method, Random Circuit Perturbation (RACIPE), to explore the dynamic behaviors of gene regulatory circuits without the requirement of detailed kinetic parameters. RACIPE takes a network topology as the input, and generates an unbiased ensemble of models with varying kinetic parameters. Each model is subjected to simulation, followed by statistical analysis for the ensemble. We tested RACIPE on several gene circuits, and found that the predicted gene expression patterns from all of the models converge to experimentally observed gene state clusters. We expect RACIPE to be a powerful method to identify the role of network topology in determining network operating principles.

Sensitivity estimation for stochastic models of biochemical reaction networks in the presence of extrinsic variability

Determining the sensitivity of certain system states or outputs to variations in parameters facilitates our understanding of the inner working of that system and is an essential design tool for the de novo construction of robust systems. In cell biology, the output of interest is often the response of a certain reaction network to some input (e.g., stressors or nutrients) and one aims to quantify the sensitivity of this response in the presence of parameter heterogeneity. We argue that for such applications, parametric sensitivities in their standard form do not paint a complete picture of a system's robustness since one assumes that all cells in the population have the same parameters and are perturbed in the same way. Here, we consider stochasticreaction networks in which the parameters are randomly distributed over the population and propose a new sensitivity index that captures the robustness of system outputs upon changes in the characteristics of the parameter distribution, rather than the parameters themselves. Subsequently, we make use of Girsanov's likelihood ratio method to construct a Monte Carlo estimator of this sensitivity index. However, it turns out that this estimator has an exceedingly large variance. To overcome this problem, we propose a novel estimation algorithm that makes use of a marginalization of the path distribution of stochasticreaction networks and leads to Rao-Blackwellized estimators with reduced variance.

Stochastic modeling and numerical simulation of gene regulatory networks with protein bursting

Gene expression is inherently stochastic. Advanced single-cell microscopy techniques together with mathematical models for single gene expression led to important insights in elucidating the sources of intrinsic noise in prokaryotic and eukaryotic cells. In addition to the finite size effects due to low copy numbers, translational bursting is a dominant source of stochasticity in cell scenarios involving few short lived mRNA transcripts with high translational efficiency (as is typically the case for prokaryotes), causing protein synthesis to occur in random bursts. In the context of gene regulation cascades, the Chemical Master Equation (CME) governing gene expression has in general no closed form solution, and the accurate stochastic simulation of the dynamics of complex gene regulatory networks is a major computational challenge. The CME associated to a single gene self regulatory motif has been previously approximated by a one dimensional time dependent partial integral differential equation (PIDE). However, to the best of our knowledge, multidimensional versions for such PIDE have not been developed yet. Here we propose a multidimensional PIDE model for regulatory networks involving multiple genes with self and cross regulations (in which genes can be regulated by different transcription factors) derived as the continuous counterpart of a CME with jump process. The model offers a reliable description of systems with translational bursting. In order to provide an efficient numerical solution, we develop a semilagrangian method to discretize the differential part of the PIDE, combined with a composed trapezoidal quadrature formula to approximate the integral term. We apply the model and numerical method to study sustained stochastic oscillations and the development of competence, a particular case of transient differentiation attained by certain bacterial cells under stress conditions. We found that the resulting probability distributions are distinguishable from those characteristic of other transient differentiation processes. In this way, they can be employed as markers or signatures that identify such phenomena from bacterial population experimental data, for instance. The computational efficiency of the semilagrangian method makes it suitable for purposes like model identification and parameter estimation from experimental data or, in combination with optimization routines, the design of gene regulatory networks under molecular noise.

Long-term tracking of budding yeast cells in brightfield microscopy: CellStar and the Evaluation Platform

With the continuous expansion of single cell biology, the observation of the behaviour of individual cells over extended durations and with high accuracy has become a problem of central importance. Surprisingly, even for yeast cells that have relatively regular shapes, no solution has been proposed that reaches the high quality required for long-term experiments for segmentation and tracking (S&T) based on brightfield images. Here, we present CellStar, a tool chain designed to achieve good performance in long-term experiments. The key features are the use of a new variant of parametrized active rays for segmentation, a neighbourhood-preserving criterion for tracking, and the use of an iterative approach that incrementally improves S&T quality. A graphical user interface enables manual corrections of S&T errors and their use for the automated correction of other, related errors and for parameter learning. We created a benchmark dataset with manually analysed images and compared CellStar with six other tools, showing its high performance, notably in long-term tracking. As a community effort, we set up a website, the Yeast Image Toolkit, with the benchmark and the Evaluation Platform to gather this and additional information provided by others.