Random partitioning of molecules at cell division

On Information Extraction and Decoding Mechanisms Improved by Noisy Amplification in Signaling Pathways

The cells need to process information about extracellular stimuli. They encode, transmit and decode the information to elicit an appropriate response. Studies aimed at understanding how such information is decoded in the signaling pathways to generate a specific cellular response have become essential. Eukaryotic cells decode information through two different mechanisms: the feed-forward loop and the promoter affinity. Here, we investigate how these two mechanisms improve information transmission. A detailed comparison is made between the stochastic model of the MAPK/ERK pathway and a stochastic minimal decoding model. The maximal amount of transmittable information was computed. The results suggest that the decoding mechanism of the MAPK/ERK pathway improve the channel capacity because it behaves as a noisy amplifier. We show a positive dependence between the noisy amplification and the amount of information extracted. Additionally, we show that the extrinsic noise can be tuned to improve information transmission. This investigation has revealed that the feed-forward loop and the promoter affinity motifs extract information thanks to processes of amplification and noise addition. Moreover, the channel capacity is enhanced when both decoding mechanisms are coupled. Altogether, these findings suggest novel characteristics in how decoding mechanisms improve information transmission.

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.

Cell population heterogeneity driven by stochastic partition and growth optimality

A fundamental question in biology is how cell populations evolve into different subtypes based on homogeneous processes at the single cell level. Here we show that population bimodality can emerge even when biological processes are homogenous at the cell level and the environment is kept constant. Our model is based on the stochastic partitioning of a cell component with an optimal copy number. We show that the existence of unimodal or bimodal distributions depends on the variance of partition errors and the growth rate tolerance around the optimal copy number. In particular, our theory provides a consistent explanation for the maintenance of aneuploid states in a population. The proposed model can also be relevant for other cell components such as mitochondria and plasmids, whose abundances affect the growth rate and are subject to stochastic partition at cell division.

Quantitative single-cell live imaging links HES5 dynamics with cell-state and fate in murine neurogenesis

During embryogenesis cells make fate decisions within complex tissue environments. The levels and dynamics of transcription factor expression regulate these decisions. Here, we use single cell live imaging of an endogenous HES5 reporter and absolute protein quantification to gain a dynamic view of neurogenesis in the embryonic mammalian spinal cord. We report that dividing neural progenitors show both aperiodic and periodic HES5 protein fluctuations. Mathematical modelling suggests that in progenitor cells the HES5 oscillator operates close to its bifurcation boundary where stochastic conversions between dynamics are possible. HES5 expression becomes more frequently periodic as cells transition to differentiation which, coupled with an overall decline in HES5 expression, creates a transient period of oscillations with higher fold expression change. This increases the decoding capacity of HES5 oscillations and correlates with interneuron versus motor neuron cell fate. Thus, HES5 undergoes complex changes in gene expression dynamics as cells differentiate.

Distribution of Initiation Times Reveals Mechanisms of Transcriptional Regulation in Single Cells

Transcription is the dominant point of control of gene expression. Biochemical studies have revealed key molecular components of transcription and their interactions, but the dynamics of transcription initiation in cells is still poorly understood. This state of affairs is being remedied with experiments that observe transcriptional dynamics in single cells using fluorescent reporters. Quantitative information about transcription initiation dynamics can also be extracted from experiments that use electron micrographs of RNA polymerases caught in the act of transcribing a gene (Miller spreads). Inspired by these data, we analyze a general stochastic model of transcription initiation and elongation and compute the distribution of transcription initiation times. We show that different mechanisms of initiation leave distinct signatures in the distribution of initiation times that can be compared to experiments. We analyze published data from micrographs of RNA polymerases transcribing ribosomal RNA genes in Escherichia coli and compare the observed distributions of interpolymerase distances with the predictions from previously hypothesized mechanisms for the regulation of these genes. Our analysis demonstrates the potential of measuring the distribution of time intervals between initiation events as a probe for dissecting mechanisms of transcription initiation in live cells.

Is the cell really a machine?

It has become customary to conceptualize the living cell as an intricate piece of machinery, different to a man-made machine only in terms of its superior complexity. This familiar understanding grounds the conviction that a cell's organization can be explained reductionistically, as well as the idea that its molecular pathways can be construed as deterministic circuits. The machine conception of the cell owes a great deal of its success to the methods traditionally used in molecular biology. However, the recent introduction of novel experimental techniques capable of tracking individual molecules within cells in real time is leading to the rapid accumulation of data that are inconsistent with an engineering view of the cell. This paper examines four major domains of current research in which the challenges to the machine conception of the cell are particularly pronounced: cellular architecture, protein complexes, intracellular transport, and cellular behaviour. It argues that a new theoretical understanding of the cell is emerging from the study of these phenomena which emphasizes the dynamic, self-organizing nature of its constitution, the fluidity and plasticity of its components, and the stochasticity and non-linearity of its underlying processes.

Mathematical modelling reveals unexpected inheritance and variability patterns of cell cycle parameters in mammalian cells

The cell cycle is the fundamental process of cell populations, it is regulated by environmental cues and by intracellular checkpoints. Cell cycle variability in clonal cell population is caused by stochastic processes such as random partitioning of cellular components to progeny cells at division and random interactions among biomolecules in cells. One of the important biological questions is how the dynamics at the cell cycle scale, which is related to family dependencies between the cell and its descendants, affects cell population behavior in the long-run. We address this question using a “mechanistic” model, built based on observations of single cells over several cell generations, and then extrapolated in time. We used cell pedigree observations of NIH 3T3 cells including FUCCI markers, to determine patterns of inheritance of cell-cycle phase durations and single-cell protein dynamics. Based on that information we developed a hybrid mathematical model, involving bifurcating autoregression to describe stochasticity of partitioning and inheritance of cell-cycle-phase times, and an ordinary differential equation system to capture single-cell protein dynamics. Long-term simulations, concordant with in vitro experiments, demonstrated the model reproduced the main features of our data and had homeostatic properties. Moreover, heterogeneity of cell cycle may have important consequences during population development. We discovered an effect similar to genetic drift, amplified by family relationships among cells. In consequence, the progeny of a single cell with a short cell cycle time had a high probability of eventually dominating the population, due to the heritability of cell-cycle phases. Patterns of epigenetic heritability in proliferating cells are important for understanding long-term trends of cell populations which are either required to provide the influx of maturing cells (such as hematopoietic stem cells) or which started proliferating uncontrollably (such as cancer cells).

All cells in multicellular organisms obey orchestrated sequences of signals to ensure developmental and homeostatic fitness under a variety of external stimuli. However, there also exist self-perpetuating stem-cell populations, the function of which is to provide a steady supply of differentiated progenitors that in turn ensure persistence of organism functions. This “cell production engine” is an important element of biological homeostasis. A similar process, albeit distorted in many respects, plays a major role in cancer development; here the robustness of homeostasis contributes to difficulty in eradication of malignancy. An important role in homeostasis seems to be played by generation of heterogeneity among cell phenotypes, which then can be shaped by selection and other genetic forces. In the present paper, we present a model of a cultured cell population, which factors in relationships among related cells and the dynamics of cell growth and important proteins regulating cell division. We find that the model not only maintains homeostasis, but that it also responds to perturbations in a manner that is similar to that exhibited by the Wright-Fisher model of population genetics. The model-cell population can become dominated by the progeny of the fittest individuals, without invoking advantageous mutations. If confirmed, this may provide an alternative mode of evolution of cell populations.

Quantifying Cancer Epithelial-Mesenchymal Plasticity and its Association with Stemness and Immune Response

Cancer cells can acquire a spectrum of stable hybrid epithelial/mesenchymal (E/M) states during epithelial-mesenchymal transition (EMT). Cells in these hybrid E/M phenotypes often combine epithelial and mesenchymal features and tend to migrate collectively commonly as small clusters. Such collectively migrating cancer cells play a pivotal role in seeding metastases and their presence in cancer patients indicates an adverse prognostic factor. Moreover, cancer cells in hybrid E/M phenotypes tend to be more associated with stemness which endows them with tumor-initiation ability and therapy resistance. Most recently, cells undergoing EMT have been shown to promote immune suppression for better survival. A systematic understanding of the emergence of hybrid E/M phenotypes and the connection of EMT with stemness and immune suppression would contribute to more effective therapeutic strategies. In this review, we first discuss recent efforts combining theoretical and experimental approaches to elucidate mechanisms underlying EMT multi-stability (i.e., the existence of multiple stable phenotypes during EMT) and the properties of hybrid E/M phenotypes. Following we discuss non-cell-autonomous regulation of EMT by cell cooperation and extracellular matrix. Afterwards, we discuss various metrics that can be used to quantify EMT spectrum. We further describe possible mechanisms underlying the formation of clusters of circulating tumor cells. Last but not least, we summarize recent systems biology analysis of the role of EMT in the acquisition of stemness and immune suppression.

Asymmetric Inheritance of Cell Fate Determinants: Focus on RNA

During the last decade, and mainly primed by major developments in high-throughput sequencing technologies, the catalogue of RNA molecules harbouring regulatory functions has increased at a steady pace. Current evidence indicates that hundreds of mammalian RNAs have regulatory roles at several levels, including transcription, translation/post-translation, chromatin structure, and nuclear architecture, thus suggesting that RNA molecules are indeed mighty controllers in the flow of biological information. Therefore, it is logical to suggest that there must exist a series of molecular systems that safeguard the faithful inheritance of RNA content throughout cell division and that those mechanisms must be tightly controlled to ensure the successful segregation of key molecules to the progeny. Interestingly, whilst a handful of integral components of mammalian cells seem to follow a general pattern of asymmetric inheritance throughout division, the fate of RNA molecules largely remains a mystery. Herein, we will discuss current concepts of asymmetric inheritance in a wide range of systems, including prions, proteins, and finally RNA molecules, to assess overall the biological impact of RNA inheritance in cellular plasticity and evolutionary fitness.

Phenotypic Heterogeneity in Bacterial Quorum Sensing Systems

Quorum sensing is usually thought of as a collective behavior in which all members of a population partake. However, over the last decade, several reports of phenotypic heterogeneity in quorum sensing-related gene expression have been put forward, thus challenging this view. In the respective systems, cells of isogenic populations did not contribute equally to autoinducer production or target gene activation, and in some cases, the fraction of contributing cells was modulated by environmental factors. Here, we look into potential origins of these incidences and into how initial cell-to-cell variations might be amplified to establish distinct phenotypic heterogeneity. We furthermore discuss potential functions heterogeneity in bacterial quorum sensing systems could serve: as a preparation for environmental fluctuations (bet hedging), as a more cost-effective way of producing public goods (division of labor), as a loophole for genotypic cooperators when faced with non-contributing mutants (cheat protection), or simply as a means to fine-tune the output of the population as a whole (output modulation). We illustrate certain aspects of these recent developments with the model organisms Sinorhizobium meliloti, Sinorhizobium fredii and Bacillus subtilis, which possess quorum sensing systems of different complexity, but all show phenotypic heterogeneity therein.

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.

Exactly solvable model of gene expression in a proliferating bacterial cell population with stochastic protein bursts and protein partitioning

Many of the existing stochastic models of gene expression contain the first-order decay reaction term that may describe active protein degradation or dilution. If the model variable is interpreted as the molecule number, and not concentration, the decay term may also approximate the loss of protein molecules due to cell division as a continuous degradation process. The seminal model of that kind leads to gamma distributions of protein levels, whose parameters are defined by the mean frequency of protein bursts and mean burst size. However, such models (whether interpreted in terms of molecule numbers or concentrations) do not correctly account for the noise due to protein partitioning between daughter cells. We present an exactly solvable stochastic model of gene expression in cells dividing at random times, where we assume description in terms of molecule numbers with a constant mean protein burst size. The model is based on a population balance equation supplemented with protein production in random bursts. If protein molecules are partitioned equally between daughter cells, we obtain at steady state the analytical expressions for probability distributions similar in shape to gamma distributions, yet with quite different values of mean burst size and mean burst frequency than would result from fitting of the classical continuous-decay model to these distributions. For random partitioning of protein molecules between daughter cells, we obtain the moment equations for the protein number distribution and thus the analytical formulas for the squared coefficient of variation.

Maximum entropy and population heterogeneity in continuous cell cultures

Continuous cultures of mammalian cells are complex systems displaying hallmark phenomena of nonlinear dynamics, such as multi-stability, hysteresis, as well as sharp transitions between different metabolic states. In this context mathematical models may suggest control strategies to steer the system towards desired states. Although even clonal populations are known to exhibit cell-to-cell variability, most of the currently studied models assume that the population is homogeneous. To overcome this limitation, we use the maximum entropy principle to model the phenotypic distribution of cells in a chemostat as a function of the dilution rate. We consider the coupling between cell metabolism and extracellular variables describing the state of the bioreactor and take into account the impact of toxic byproduct accumulation on cell viability. We present a formal solution for the stationary state of the chemostat and show how to apply it in two examples. First, a simplified model of cell metabolism where the exact solution is tractable, and then a genome-scale metabolic network of the Chinese hamster ovary (CHO) cell line. Along the way we discuss several consequences of heterogeneity, such as: qualitative changes in the dynamical landscape of the system, increasing concentrations of byproducts that vanish in the homogeneous case, and larger population sizes.

Mitochondrial Heterogeneity

Cell-to-cell heterogeneity drives a range of (patho)physiologically important phenomena, such as cell fate and chemotherapeutic resistance. The role of metabolism, and particularly of mitochondria, is increasingly being recognized as an important explanatory factor in cell-to-cell heterogeneity. Most eukaryotic cells possess a population of mitochondria, in the sense that mitochondrial DNA (mtDNA) is held in multiple copies per cell, where the sequence of each molecule can vary. Hence, intra-cellular mitochondrial heterogeneity is possible, which can induce inter-cellular mitochondrial heterogeneity, and may drive aspects of cellular noise. In this review, we discuss sources of mitochondrial heterogeneity (variations between mitochondria in the same cell, and mitochondrial variations between supposedly identical cells) from both genetic and non-genetic perspectives, and mitochondrial genotype-phenotype links. We discuss the apparent homeostasis of mtDNA copy number, the observation of pervasive intra-cellular mtDNA mutation (which is termed "microheteroplasmy"), and developments in the understanding of inter-cellular mtDNA mutation ("macroheteroplasmy"). We point to the relationship between mitochondrial supercomplexes, cristal structure, pH, and cardiolipin as a potential amplifier of the mitochondrial genotype-phenotype link. We also discuss mitochondrial membrane potential and networks as sources of mitochondrial heterogeneity, and their influence upon the mitochondrial genome. Finally, we revisit the idea of mitochondrial complementation as a means of dampening mitochondrial genotype-phenotype links in light of recent experimental developments. The diverse sources of mitochondrial heterogeneity, as well as their increasingly recognized role in contributing to cellular heterogeneity, highlights the need for future single-cell mitochondrial measurements in the context of cellular noise studies.

Intrinsic and extrinsic noise of gene expression in lineage trees

Cell-to-cell heterogeneity is driven by stochasticity in intracellular reactions and the population dynamics. While these sources are usually studied separately, we develop an agent-based framework that accounts for both factors while tracking every single cell of a growing population. Apart from the common intrinsic variability, the framework also predicts extrinsic noise without the need to introduce fluctuating rate constants. Instead, extrinsic fluctuations are explained by cell cycle fluctuations and differences in cell age. We provide explicit formulas to quantify mean molecule numbers, intrinsic and extrinsic noise statistics in two-colour experiments. We find that these statistics differ significantly depending on the experimental setup used to observe the cells. We illustrate this fact using (i) averages over an isolated cell lineage tracked over many generations as observed in the mother machine, (ii) population snapshots with known cell ages as recorded in time-lapse microscopy, and (iii) snapshots with unknown cell ages as measured from static images or flow cytometry. Applying the method to models of stochastic gene expression and feedback regulation elucidates that isolated lineages, as compared to snapshot data, can significantly overestimate the mean number of molecules, overestimate extrinsic noise but underestimate intrinsic noise and have qualitatively different sensitivities to cell cycle fluctuations.

Association of mitochondria with microtubules inhibits mitochondrial fission by precluding assembly of the fission protein Dnm1

Mitochondria are organized as tubular networks in the cell and undergo fission and fusion. Although several of the molecular players involved in mediating mitochondrial dynamics have been identified, the precise cellular cues that initiate mitochondrial fission or fusion remain largely unknown. In fission yeast (Schizosaccharomyces pombe), mitochondria are organized along microtubule bundles. Here, we employed deletions of kinesin-like proteins to perturb microtubule dynamics and used high-resolution and time-lapse fluorescence microscopy, revealing that mitochondrial lengths mimic microtubule lengths. Furthermore, we determined that compared with WT cells, mutant cells with long microtubules exhibit fewer mitochondria, and mutant cells with short microtubules have an increased number of mitochondria because of reduced mitochondrial fission in the former and elevated fission in the latter. Correspondingly, upon onset of closed mitosis in fission yeast, wherein interphase microtubules assemble to form the spindle within the nucleus, we observed increased mitochondrial fission. We found that the consequent rise in the mitochondrial copy number is necessary to reduce partitioning errors during independent segregation of mitochondria between daughter cells. We also discovered that the association of mitochondria with microtubules physically impedes the assembly of the fission protein Dnm1 around mitochondria, resulting in inhibition of mitochondrial fission. Taken together, we demonstrate a mechanism for the regulation of mitochondrial fission that is dictated by the interaction between mitochondria and the microtubule cytoskeleton.

Noise in bacterial gene expression

The expression level of a gene can fluctuate significantly between individuals within a population of genetically identical cells. The resultant phenotypic heterogeneity could be exploited by bacteria to adapt to changing environmental conditions. Noise is hence a genome-wide phenomenon that arises from the stochastic nature of the biochemical reactions that take place during gene expression and the relatively low abundance of the molecules involved. The production of mRNA and proteins therefore occurs in bursts, with alternating episodes of high and low activity during transcription and translation. Single-cell and single-molecule studies demonstrated that noise within gene expression is influenced by a combination of both intrinsic and extrinsic factors. However, our mechanistic understanding of this process at the molecular level is still rather limited. Further investigation is necessary that takes into account the detailed knowledge of gene regulation gained from biochemical studies.

Massive mitochondrial DNA content in diplonemid and kinetoplastid protists

The mitochondrial DNA of diplonemid and kinetoplastid protists is known for its suite of bizarre features, including the presence of concatenated circular molecules, extensive trans‐splicing and various forms of RNA editing. Here we report on the existence of another remarkable characteristic: hyper‐inflated DNA content. We estimated the total amount of mitochondrial DNA in four kinetoplastid species (Trypanosoma brucei, Trypanoplasma borreli, Cryptobia helicis, and Perkinsela sp.) and the diplonemid Diplonema papillatum. Staining with 4′,6‐diamidino‐2‐phenylindole and RedDot1 followed by color deconvolution and quantification revealed massive inflation in the total amount of DNA in their organelles. This was further confirmed by electron microscopy. The most extreme case is the ∼260 Mbp of DNA in the mitochondrion of Diplonema, which greatly exceeds that in its nucleus; this is, to our knowledge, the largest amount of DNA described in any organelle. Perkinsela sp. has a total mitochondrial DNA content ~6.6× greater than its nuclear genome. This mass of DNA occupies most of the volume of the Perkinsela cell, despite the fact that it contains only six protein‐coding genes. Why so much DNA? We propose that these bloated mitochondrial DNAs accumulated by a ratchet‐like process. Despite their excessive nature, the synthesis and maintenance of these mtDNAs must incur a relatively low cost, considering that diplonemids are one of the most ubiquitous and speciose protist groups in the ocean. © 2018 The Authors. IUBMB Life published by Wiley Periodicals, Inc. on behalf of International Union of Biochemistry and Molecular Biology., 70(12):1267–1274, 2018

Making sense of snapshot data: ergodic principle for clonal cell populations

Population growth is often ignored when quantifying gene expression levels across clonal cell populations. We develop a framework for obtaining the molecule number distributions in an exponentially growing cell population taking into account its age structure. In the presence of generation time variability, the average acquired across a population snapshot does not obey the average of a dividing cell over time, apparently contradicting ergodicity between single cells and the population. Instead, we show that the variation observed across snapshots with known cell age is captured by cell histories, a single-cell measure obtained from tracking an arbitrary cell of the population back to the ancestor from which it originated. The correspondence between cells of known age in a population with their histories represents an ergodic principle that provides a new interpretation of population snapshot data. We illustrate the principle using analytical solutions of stochastic gene expression models in cell populations with arbitrary generation time distributions. We further elucidate that the principle breaks down for biochemical reactions that are under selection, such as the expression of genes conveying antibiotic resistance, which gives rise to an experimental criterion with which to probe selection on gene expression fluctuations.

A Symmetric Molecule Produced by Mycobacteria Generates Cell-Length Asymmetry during Cell-Division and Thereby Cell-Length Heterogeneity

Diadenosine polyphosphates, Ap_(2-7)A, which contain two adenosines in a 5',5' linkage through phosphodiester bonds involving 2-7 phosphates, regulate diverse cellular functions in all organisms, from bacteria to humans, under normal and stress conditions. We had earlier reported consistent occurrence of asymmetric constriction during division (ACD) in 20-30% of dividing mycobacterial cells in culture, irrespective of different growth media, implying exogenous action of some factor of mycobacterial origin. Consistent with this premise, concentrated culture supernatant (CCS), but not the equivalent volume-wise concentrated unused medium, dramatically enhanced the ACD proportion to 70-90%. Mass spectrometry and biochemical analyses of the bioactive fraction from CCS revealed the ACD-effecting factor to be Ap₆A. Synthetic Ap₆A showed a mass spectrometry profile, biochemical characteristics, and bioactivity identical to native Ap₆A in the CCS. Thus, the present work reveals a novel role for Ap₆A in generating cell-length asymmetry during mycobacterial cell-division and thereby cell-length heterogeneity in the population.

Division rate, cell size and proteome allocation: impact on gene expression noise and implications for the dynamics of genetic circuits

The cell division rate, size and gene expression programmes change in response to external conditions. These global changes impact on average concentrations of biomolecule and their variability or noise. Gene expression is inherently stochastic, and noise levels of individual proteins depend on synthesis and degradation rates as well as on cell-cycle dynamics. We have modelled stochastic gene expression inside growing and dividing cells to study the effect of division rates on noise in mRNA and protein expression. We use assumptions and parameters relevant to Escherichia coli, for which abundant quantitative data are available. We find that coupling of transcription, but not translation rates to the rate of cell division can result in protein concentration and noise homeostasis across conditions. Interestingly, we find that the increased cell size at fast division rates, observed in E. coli and other unicellular organisms, buffers noise levels even for proteins with decreased expression at faster growth. We then investigate the functional importance of these regulations using gene regulatory networks that exhibit bi-stability and oscillations. We find that network topology affects robustness to changes in division rate in complex and unexpected ways. In particular, a simple model of persistence, based on global physiological feedback, predicts increased proportion of persister cells at slow division rates. Altogether, our study reveals how cell size regulation in response to cell division rate could help controlling gene expression noise. It also highlights that understanding circuits' robustness across growth conditions is key for the effective design of synthetic biological systems.

Division rate, cell size and proteome allocation: impact on gene expression noise and implications for the dynamics of genetic circuits

The cell division rate, size and gene expression programmes change in response to external conditions. These global changes impact on average concentrations of biomolecule and their variability or noise. Gene expression is inherently stochastic, and noise levels of individual proteins depend on synthesis and degradation rates as well as on cell-cycle dynamics. We have modelled stochastic gene expression inside growing and dividing cells to study the effect of division rates on noise in mRNA and protein expression. We use assumptions and parameters relevant to Escherichia coli, for which abundant quantitative data are available. We find that coupling of transcription, but not translation rates to the rate of cell division can result in protein concentration and noise homeostasis across conditions. Interestingly, we find that the increased cell size at fast division rates, observed in E. coli and other unicellular organisms, buffers noise levels even for proteins with decreased expression at faster growth. We then investigate the functional importance of these regulations using gene regulatory networks that exhibit bi-stability and oscillations. We find that network topology affects robustness to changes in division rate in complex and unexpected ways. In particular, a simple model of persistence, based on global physiological feedback, predicts increased proportion of persister cells at slow division rates. Altogether, our study reveals how cell size regulation in response to cell division rate could help controlling gene expression noise. It also highlights that understanding circuits' robustness across growth conditions is key for the effective design of synthetic biological systems.

Statistics and simulation of growth of single bacterial cells: illustrations with B. subtilis and E. coli

The inherent stochasticity of molecular reactions prevents us from predicting the exact state of single-cells in a population. However, when a population grows at steady-state, the probability to observe a cell with particular combinations of properties is fixed. Here we validate and exploit existing theory on the statistics of single-cell growth in order to predict the probability of phenotypic characteristics such as cell-cycle times, volumes, accuracy of division and cell-age distributions, using real-time imaging data for Bacillus subtilis and Escherichia coli. Our results show that single-cell growth-statistics can accurately be predicted from a few basic measurements. These equations relate different phenotypic characteristics, and can therefore be used in consistency tests of experimental single-cell growth data and prediction of single-cell statistics. We also exploit these statistical relations in the development of a fast stochastic-simulation algorithm of single-cell growth and protein expression. This algorithm greatly reduces computational burden, by recovering the statistics of growing cell-populations from the simulation of only one of its lineages. Our approach is validated by comparison of simulations and experimental data. This work illustrates a methodology for the prediction, analysis and tests of consistency of single-cell growth and protein expression data from a few basic statistical principles.

Cultivation and quantitative single-cell analysis of Saccharomyces cerevisiae on a multifunctional microfluidic device

Here, we present a multifunctional microfluidic device whose integrative design enables to combine cell culture studies and quantitative single cell biomolecule analysis. The platform consists of 32 analysis units providing two key features; first, a micrometer-sized trap for hydrodynamic capture of a single Saccharomyces cerevisiae (S. cerevisiae) yeast cell; second, a convenient double-valve configuration surrounding the trap. Actuating of the outer valve with integrated opening results in a partial isolation in a volume of 11.8 nL, i.e. the cell surrounding fluid can be exchanged diffusion-based without causing shear stress or cell loss. Actuation of the inner ring-shaped valve isolates the trapped cell completely in a small analysis volume of 230 pL. The device was used to determine the growth rate of yeast cells (S. cerevisiae) under under optimum and oxidative stress conditions. In addition, we successfully quantified the cofactor beta-nicotinamide adenine dinucleotide phosphate (NAD(P)H) in single and few cells exposed to the different microenvironments. In conclusion, the microdevice enables to analyze the influence of an external stress factor on the cellular fitness in a fast and more comprehensive way as cell growth and intracellular biomolecule levels can be investigated.

Effects of growth rate and promoter activity on single-cell protein expression

Protein expression in a single cell depends on its global physiological state. Moreover, genetically-identical cells exhibit variability (noise) in protein expression, arising from the stochastic nature of biochemical processes, cell growth and division. While it is well understood how cellular growth rate influences mean protein expression, little is known about the relationship between growth rate and noise in protein expression. Here we quantify this relationship in Bacillus subtilis by a novel combination of experiments and theory. We measure the effects of promoter activity and growth rate on the expression of a fluorescent protein in single cells. We disentangle the observed protein expression noise into protein-specific and systemic contributions, using theory and variance decomposition. We find that noise in protein expression depends solely on mean expression levels, regardless of whether expression is set by promoter activity or growth rate, and that noise increases linearly with growth rate. Our results can aid studies of (synthetic) gene circuits of single cells and their condition dependence.

Taking chances and making mistakes: non-genetic phenotypic heterogeneity and its consequences for surviving in dynamic environments

Natural selection has shaped the strategies for survival and growth of microorganisms. The success of microorganisms depends not only on slow evolutionary tuning but also on the ability to adapt to unpredictable changes in their environment. In principle, adaptive strategies range from purely deterministic mechanisms to those that exploit the randomness intrinsic to many cellular and molecular processes. Depending on the environment and selective pressures, particular strategies can lie somewhere along this continuum. In recent years, non-genetic cell-to-cell differences have received a lot of attention, not least because of their potential impact on the ability of microbial populations to survive in dynamic environments. Using several examples, we describe the origins of spontaneous and induced mechanisms of phenotypic adaptation. We identify some of the commonalities of these examples and consider the potential role of chance and constraints in microbial phenotypic adaptation.

Quantifying ribosome dynamics in Escherichia coli using fluorescence

Ribosomes are a crucial component of the physiological state of a cell. Therefore, we aimed to monitor ribosome dynamics using a fast and easy fluorescence readout. Using fluorescent-labeled ribosomal proteins, the dynamics of ribosomes during batch cultivation and during nutritional shift conditions was investigated. The fluorescence readout was compared to the cellular rRNA content determined by capillary gel electrophoresis with laser-induced fluorescence detection during exponentially accelerating and decelerating growth. We found a linear correlation between the observed fluorescence and the extracted rRNA content throughout cultivation, demonstrating the applicability of this method. Moreover, the results show that ribosome dynamics, as a result of slowing growth, are accompanied by the passive effect of dilution of preexisting ribosomes, de novo ribosome synthesis and ribosome degradation. In light of the challenging task of deciphering ribosome regulatory mechanisms, our approach of using fluorescence to follow ribosome dynamics will allow more comprehensive studies of biological systems.

Organelles - understanding noise and heterogeneity in cell biology at an intermediate scale

Many studies over the years have shown that non-genetic mechanisms for producing cell-to-cell variation can lead to highly variable behaviors across genetically identical populations of cells. Most work to date has focused on gene expression noise as the primary source of phenotypic heterogeneity, yet other sources may also contribute. In this Commentary, we explore organelle-level heterogeneity as a potential secondary source of cellular 'noise' that contributes to phenotypic heterogeneity. We explore mechanisms for generating organelle heterogeneity and present evidence of functional links between organelle morphology and cellular behavior. Given the many instances in which molecular-level heterogeneity has been linked to phenotypic heterogeneity, we posit that organelle heterogeneity may similarly contribute to overall phenotypic heterogeneity and underline the importance of studying organelle heterogeneity to develop a more comprehensive understanding of phenotypic heterogeneity. Finally, we conclude with a discussion of the medical challenges associated with phenotypic heterogeneity and outline how improved methods for characterizing and controlling this heterogeneity may lead to improved therapeutic strategies and outcomes for patients.

An integrated modelling framework from cells to organism based on a cohort of digital embryos

We conducted a quantitative comparison of developing sea urchin embryos based on the analysis of five digital specimens obtained by automatic processing of in toto 3D+ time image data. These measurements served the reconstruction of a prototypical cell lineage tree able to predict the spatiotemporal cellular organisation of a normal sea urchin blastula. The reconstruction was achieved by designing and tuning a multi-level probabilistic model that reproduced embryo-level dynamics from a small number of statistical parameters characterising cell proliferation, cell surface area and cell volume evolution along the cell lineage. Our resulting artificial prototype was embedded in 3D space by biomechanical agent-based modelling and simulation, which allowed a systematic exploration and optimisation of free parameters to fit the experimental data and test biological hypotheses. The spherical monolayered blastula and the spatial arrangement of its different cell types appeared tightly constrained by cell stiffness, cell-adhesion parameters and blastocoel turgor pressure.

Effects of cell-cycle-dependent expression on random fluctuations in protein levels

Expression of many genes varies as a cell transitions through different cell-cycle stages. How coupling between stochastic expression and cell cycle impacts cell-to-cell variability (noise) in the level of protein is not well understood. We analyse a model where a stable protein is synthesized in random bursts, and the frequency with which bursts occur varies within the cell cycle. Formulae quantifying the extent of fluctuations in the protein copy number are derived and decomposed into components arising from the cell cycle and stochastic processes. The latter stochastic component represents contributions from bursty expression and errors incurred during partitioning of molecules between daughter cells. These formulae reveal an interesting trade-off: cell-cycle dependencies that amplify the noise contribution from bursty expression also attenuate the contribution from partitioning errors. We investigate the existence of optimum strategies for coupling expression to the cell cycle that minimize the stochastic component. Intriguingly, results show that a zero production rate throughout the cell cycle, with expression only occurring just before cell division, minimizes noise from bursty expression for a fixed mean protein level. By contrast, the optimal strategy in the case of partitioning errors is to make the protein just after cell division. We provide examples of regulatory proteins that are expressed only towards the end of the cell cycle, and argue that such strategies enhance robustness of cell-cycle decisions to the intrinsic stochasticity of gene expression.

Temperature-Dependent Model of Multi-step Transcription Initiation in Escherichia coli Based on Live Single-Cell Measurements

Transcription kinetics is limited by its initiation steps, which differ between promoters and with intra- and extracellular conditions. Regulation of these steps allows tuning both the rate and stochasticity of RNA production. We used time-lapse, single-RNA microscopy measurements in live Escherichia coli to study how the rate-limiting steps in initiation of the P_lac/ara-1 promoter change with temperature and induction scheme. For this, we compared detailed stochastic models fit to the empirical data in maximum likelihood sense using statistical methods. Using this analysis, we found that temperature affects the rate limiting steps unequally, as nonlinear changes in the closed complex formation suffice to explain the differences in transcription dynamics between conditions. Meanwhile, a similar analysis of the P_tetA promoter revealed that it has a different rate limiting step configuration, with temperature regulating different steps. Finally, we used the derived models to explore a possible cause for why the identified steps are preferred as the main cause for behavior modifications with temperature: we find that transcription dynamics is either insensitive or responds reciprocally to changes in the other steps. Our results suggests that different promoters employ different rate limiting step patterns that control not only their rate and variability, but also their sensitivity to environmental changes.

Temperature affects the behavior of cells, such as their growth rate. However, it is not well understood how these changes result from the changes at the single molecule level. We observed the production of individual RNA molecules in live cells under a wide range of temperatures. This allowed us to determine not only how fast they are produced, but also how much variability there is in this process. Next, we fit a stochastic model to the data to identify which rate-limiting steps during RNA production are responsible for the observed differences between conditions. We found that genes differ in how their RNA production is limited by different steps and in how these are affected by the temperature, which explains why different genes respond differently to temperature fluctuations.

Stochastic Kinetics of Nascent RNA

The stochastic kinetics of transcription is typically inferred from the distribution of RNA numbers in individual cells. However, cellular RNA reflects additional processes downstream of transcription, hampering this analysis. In contrast, nascent (actively transcribed) RNA closely reflects the kinetics of transcription. We present a theoretical model for the stochastic kinetics of nascent RNA, which we solve to obtain the probability distribution of nascent RNA per gene. The model allows us to evaluate the kinetic parameters of transcription from single-cell measurements of nascent RNA. The model also predicts surprising discontinuities in the distribution of nascent RNA, a feature which we verify experimentally.

Intercellular Variability in Protein Levels from Stochastic Expression and Noisy Cell Cycle Processes

Inside individual cells, expression of genes is inherently stochastic and manifests as cell-to-cell variability or noise in protein copy numbers. Since proteins half-lives can be comparable to the cell-cycle length, randomness in cell-division times generates additional intercellular variability in protein levels. Moreover, as many mRNA/protein species are expressed at low-copy numbers, errors incurred in partitioning of molecules between two daughter cells are significant. We derive analytical formulas for the total noise in protein levels when the cell-cycle duration follows a general class of probability distributions. Using a novel hybrid approach the total noise is decomposed into components arising from i) stochastic expression; ii) partitioning errors at the time of cell division and iii) random cell-division events. These formulas reveal that random cell-division times not only generate additional extrinsic noise, but also critically affect the mean protein copy numbers and intrinsic noise components. Counter intuitively, in some parameter regimes, noise in protein levels can decrease as cell-division times become more stochastic. Computations are extended to consider genome duplication, where transcription rate is increased at a random point in the cell cycle. We systematically investigate how the timing of genome duplication influences different protein noise components. Intriguingly, results show that noise contribution from stochastic expression is minimized at an optimal genome-duplication time. Our theoretical results motivate new experimental methods for decomposing protein noise levels from synchronized and asynchronized single-cell expression data. Characterizing the contributions of individual noise mechanisms will lead to precise estimates of gene expression parameters and techniques for altering stochasticity to change phenotype of individual cells.

Inside individual cells, gene products often occur at low molecular counts and are subject to considerable stochastic fluctuations (noise) in copy numbers over time. An important consequence of noisy expression is that the level of a protein can vary considerably even among genetically identical cells exposed to the same environment. Such non-genetic phenotypic heterogeneity is physiologically relevant and critically influences diverse cellular processes. In addition to noise sources inherent in gene product synthesis, recent experimental studies have uncovered additional noise mechanisms that critically effect expression. For example, the time within the cell cycle when a gene duplicates, and the time taken to complete cell cycle are governed by random processes. The key contribution of this work is development of novel mathematical results quantifying how cell cycle-related noise sources combine with stochastic expression to drive intercellular variability in protein molecular counts. Derived formulas lead to many counterintuitive results, such as increasing randomness in the timing of cell division can lower noise in the level of a protein. Finally, these results inform experimental strategies to systematically dissect the contributions of different noise sources in the expression of a gene of interest.

Tissue homogeneity requires inhibition of unequal gene silencing during development

Multicellular organisms can generate and maintain homogenous populations of cells that make up individual tissues. However, cellular processes that can disrupt homogeneity and how organisms overcome such disruption are unknown. We found that ∼100-fold differences in expression from a repetitive DNA transgene can occur between intestinal cells in Caenorhabditis elegans These differences are caused by gene silencing in some cells and are actively suppressed by parental and zygotic factors such as the conserved exonuclease ERI-1. If unsuppressed, silencing can spread between some cells in embryos but can be repeat specific and independent of other homologous loci within each cell. Silencing can persist through DNA replication and nuclear divisions, disrupting uniform gene expression in developed animals. Analysis at single-cell resolution suggests that differences between cells arise during early cell divisions upon unequal segregation of an initiator of silencing. Our results suggest that organisms with high repetitive DNA content, which include humans, could use similar developmental mechanisms to achieve and maintain tissue homogeneity.

In silico analysis of division times of Escherichia coli populations as a function of the partitioning scheme of non-functional proteins

Recent evidence suggests that cells employ functionally asymmetric partitioning schemes in division to cope with aging. We explore various schemes in silico, with a stochastic model of Escherichia coli that includes gene expression, non-functional proteins generation, aggregation and polar retention, and molecule partitioning in division. The model is implemented in SGNS2, which allows stochastic, multi-delayed reactions within hierarchical, transient, interlinked compartments. After setting parameter values of non-functional proteins’ generation and effects that reproduce realistic intracellular and population dynamics, we investigate how the spatial organization of non-functional proteins affects mean division times of cell populations in lineages and, thus, mean cell numbers over time. We find that division times decrease for increasingly asymmetric partitioning. Also, increasing the clustering of non-functional proteins decreases division times. Increasing the bias in polar segregation further decreases division times, particularly if the bias favors the older pole and aggregates’ polar retention is robust. Finally, we show that the non-energy consuming retention of inherited non-functional proteins at the older pole via nucleoid occlusion is a source of functional asymmetries and, thus, is advantageous. Our results suggest that the mechanisms of intracellular organization of non-functional proteins, including clustering and polar retention, affect the vitality of E. coli populations.

Cell Geometry: How Cells Count and Measure Size

The cell represents a highly organized state of living matter in which numerous geometrical parameters are under dynamic regulation in order to match the form of a cell with its function. Cells appear capable of regulating not only the total quantity of their internal organelles, but also the size and number of those organelles. The regulation of three parameters, size, number, and total quantity, can in principle be accomplished by regulating the production or growth of organelles, their degradation or disassembly, and their partitioning among daughter cells during division. Any or all of these steps could in principle be under regulation. But if organelle assembly or disassembly is regulated by number or size, how would the cell know how many copies of an organelle it has, or how big they are?

The importance of controlling mRNA turnover during cell proliferation

Microbial gene expression depends not only on specific regulatory mechanisms, but also on cellular growth because important global parameters, such as abundance of mRNAs and ribosomes, could be growth rate dependent. Understanding these global effects is necessary to quantitatively judge gene regulation. In the last few years, transcriptomic works in budding yeast have shown that a large fraction of its genes is coordinately regulated with growth rate. As mRNA levels depend simultaneously on synthesis and degradation rates, those studies were unable to discriminate the respective roles of both arms of the equilibrium process. We recently analyzed 80 different genomic experiments and found a positive and parallel correlation between both RNA polymerase II transcription and mRNA degradation with growth rates. Thus, the total mRNA concentration remains roughly constant. Some gene groups, however, regulate their mRNA concentration by uncoupling mRNA stability from the transcription rate. Ribosome-related genes modulate their transcription rates to increase mRNA levels under fast growth. In contrast, mitochondria-related and stress-induced genes lower mRNA levels by reducing mRNA stability or the transcription rate, respectively. We critically review here these results and analyze them in relation to their possible extrapolation to other organisms and in relation to the new questions they open.

Cathepsin B launches an apoptotic exit effort upon cell death-associated disruption of lysosomes

The release of cathepsin proteases from disrupted lysosomes results in lethal cellular autodigestion. Lysosomal disruption-related cell death is highly variable, showing both apoptotic and necrotic outcomes. As the substrate spectrum of lysosomal proteases encompasses the apoptosis-regulating proteins of the Bcl-2 family, their degradation could influence the cell death outcome upon lysosomal disruption. We used Förster resonance energy transfer (FRET)-based biosensors to image the real-time degradation of the Bcl-2-family members, Bcl-xl, Bax and Bid, in living cells undergoing lysosomal lysis and identified an early chain of proteolytic events, initiated by the release of cathepsin B, which directs cells toward apoptosis. In this apoptotic exit strategy, cathepsin B’s proteolytic activity results in apoptosis-inducing Bid and removes apoptosis-preventing Bcl-xl. Cathepsin B furthermore appears to degrade a cystein protease that would otherwise have eliminated apoptosis-supporting Bax, indirectly keeping cellular levels of the Bax protein up. The concerted effort of these three early events shifts the balance of cell fate away from necrosis and toward apoptosis.

Single Cell Analysis of a Bacterial Sender-Receiver System

Monitoring gene expression dynamics on the single cell level provides important information on cellular heterogeneity and stochasticity, and potentially allows for more accurate quantitation of gene expression processes. We here study bacterial senders and receivers genetically engineered with components of the quorum sensing system derived from Aliivibrio fischeri on the single cell level using microfluidics-based bacterial chemostats and fluorescence video microscopy. We track large numbers of bacteria over extended periods of time, which allows us to determine bacterial lineages and filter out subpopulations within a heterogeneous population. We quantitatively determine the dynamic gene expression response of receiver bacteria to varying amounts of the quorum sensing inducer N-3-oxo-C6-homoserine lactone (AHL). From this we construct AHL response curves and characterize gene expression dynamics of whole bacterial populations by investigating the statistical distribution of gene expression activity over time. The bacteria are found to display heterogeneous induction behavior within the population. We therefore also characterize gene expression in a homogeneous bacterial subpopulation by focusing on single cell trajectories derived only from bacteria with similar induction behavior. The response at the single cell level is found to be more cooperative than that obtained for the heterogeneous total population. For the analysis of systems containing both AHL senders and receiver cells, we utilize the receiver cells as ‘bacterial sensors’ for AHL. Based on a simple gene expression model and the response curves obtained in receiver-only experiments, the effective AHL concentration established by the senders and their ‘sending power’ is determined.

Accurate concentration control of mitochondria and nucleoids

All cellular materials are partitioned between daughters at cell division, but by various mechanisms and with different accuracy. In the yeast Schizosaccharomyces pombe, the mitochondria are pushed to the cell poles by the spindle. We found that mitochondria spatially reequilibrate just before division, and that the mitochondrial volume and DNA-containing nucleoids instead segregate in proportion to the cytoplasm inherited by each daughter. However, nucleoid partitioning errors are suppressed by control at two levels: Mitochondrial volume is actively distributed throughout a cell, and nucleoids are spaced out in semiregular arrays within mitochondria. During the cell cycle, both mitochondria and nucleoids appear to be produced without feedback, creating a net control of fluctuations that is just accurate enough to avoid substantial growth defects.

Robustness of the Process of Nucleoid Exclusion of Protein Aggregates in Escherichia coli

Escherichia coli segregates protein aggregates to the poles by nucleoid exclusion. Combined with cell divisions, this generates heterogeneous aggregate distributions in subsequent cell generations. We studied the robustness of this process with differing medium richness and antibiotics stress, which affect nucleoid size, using multimodal, time-lapse microscopy of live cells expressing both a fluorescently tagged chaperone (IbpA), which identifies in vivo the location of aggregates, and HupA-mCherry, a fluorescent variant of a nucleoid-associated protein. We find that the relative sizes of the nucleoid's major and minor axes change widely, in a positively correlated fashion, with medium richness and antibiotic stress. The aggregate's distribution along the major cell axis also changes between conditions and in agreement with the nucleoid exclusion phenomenon. Consequently, the fraction of aggregates at the midcell region prior to cell division differs between conditions, which will affect the degree of asymmetries in the partitioning of aggregates between cells of future generations. Finally, from the location of the peak of anisotropy in the aggregate displacement distribution, the nucleoid relative size, and the spatiotemporal aggregate distribution, we find that the exclusion of detectable aggregates from midcell is most pronounced in cells with mid-sized nucleoids, which are most common under optimal conditions. We conclude that the aggregate management mechanisms of E. coli are significantly robust but are not immune to stresses due to the tangible effect that these have on nucleoid size.

IMPORTANCE

Escherichia coli segregates protein aggregates to the poles by nucleoid exclusion. From live single-cell microscopy studies of the robustness of this process to various stresses known to affect nucleoid size, we find that nucleoid size and aggregate preferential locations change concordantly between conditions. Also, the degree of influence of the nucleoid on aggregate positioning differs between conditions, causing aggregate numbers at midcell to differ in cell division events, which will affect the degree of asymmetries in the partitioning of aggregates between cells of future generations. Finally, we find that aggregate segregation to the cell poles is most pronounced in cells with mid-sized nucleoids. We conclude that the energy-free process of the midcell exclusion of aggregates partially loses effectiveness under stressful conditions.

Stochasticity in the miR-9/Hes1 oscillatory network can account for clonal heterogeneity in the timing of differentiation

Recent studies suggest that cells make stochastic choices with respect to differentiation or division. However, the molecular mechanism underlying such stochasticity is unknown. We previously proposed that the timing of vertebrate neuronal differentiation is regulated by molecular oscillations of a transcriptional repressor, HES1, tuned by a post-transcriptional repressor, miR-9. Here, we computationally model the effects of intrinsic noise on the Hes1/miR-9 oscillator as a consequence of low molecular numbers of interacting species, determined experimentally. We report that increased stochasticity spreads the timing of differentiation in a population, such that initially equivalent cells differentiate over a period of time. Surprisingly, inherent stochasticity also increases the robustness of the progenitor state and lessens the impact of unequal, random distribution of molecules at cell division on the temporal spread of differentiation at the population level. This advantageous use of biological noise contrasts with the view that noise needs to be counteracted.