Homeostasis of protein and mRNA concentrations in growing cells

A coarse-grained bacterial cell model for resource-aware analysis and design of synthetic gene circuits

Within a cell, synthetic and native genes compete for expression machinery, influencing cellular process dynamics through resource couplings. Models that simplify competitive resource binding kinetics can guide the design of strategies for countering these couplings. However, in bacteria resource availability and cell growth rate are interlinked, which complicates resource-aware biocircuit design. Capturing this interdependence requires coarse-grained bacterial cell models that balance accurate representation of metabolic regulation against simplicity and interpretability. We propose a coarse-grained E. coli cell model that combines the ease of simplified resource coupling analysis with appreciation of bacterial growth regulation mechanisms and the processes relevant for biocircuit design. Reliably capturing known growth phenomena, it provides a unifying explanation to disparate empirical relations between growth and synthetic gene expression. Considering a biomolecular controller that makes cell-wide ribosome availability robust to perturbations, we showcase our model's usefulness in numerically prototyping biocircuits and deriving analytical relations for design guidance.

Influence of cell volume on the gene transcription rate

Cells vary in volume throughout their life cycle and in many other circumstances, while their genome remains identical. Hence, the RNA production factory must adapt to changing needs, while maintaining the same production lines. This paradox is resolved by different mechanisms in distinct cells and circumstances. RNA polymerases have evolved to cope with the particular circumstances of each case and the different characteristics of the several RNA molecule types, especially their stabilities. Here we review current knowledge on these issues. We focus on the yeast Saccharomyces cerevisiae, where many of the studies have been performed, although we compare and discuss the results obtained in other eukaryotes and propose several ideas and questions to be tested and solved in the future. TAKE AWAY.

Resource allocation in mammalian systems

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

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

Transcription-replication interactions reveal bacterial genome regulation

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

Cell size homeostasis is tightly controlled throughout the cell cycle

To achieve a stable size distribution over multiple generations, proliferating cells require a means of counteracting stochastic noise in the rate of growth, the time spent in various phases of the cell cycle, and the imprecision in the placement of the plane of cell division. In the most widely accepted model, cell size is thought to be regulated at the G1/S transition, such that cells smaller than a critical size pause at the end of G1 phase until they have accumulated mass to a predetermined size threshold, at which point the cells proceed through the rest of the cell cycle. However, a model, based solely on a specific size checkpoint at G1/S, cannot readily explain why cells with deficient G1/S control mechanisms are still able to maintain a very stable cell size distribution. Furthermore, such a model would not easily account for stochastic variation in cell size during the subsequent phases of the cell cycle, which cannot be anticipated at G1/S. To address such questions, we applied computationally enhanced quantitative phase microscopy (ceQPM) to populations of cultured human cell lines, which enables highly accurate measurement of cell dry mass of individual cells throughout the cell cycle. From these measurements, we have evaluated the factors that contribute to maintaining cell mass homeostasis at any point in the cell cycle. Our findings reveal that cell mass homeostasis is accurately maintained, despite disruptions to the normal G1/S machinery or perturbations in the rate of cell growth. Control of cell mass is generally not confined to regulation of the G1 length. Instead mass homeostasis is imposed throughout the cell cycle. In the cell lines examined, we find that the coefficient of variation (CV) in dry mass of cells in the population begins to decline well before the G1/S transition and continues to decline throughout S and G2 phases. Among the different cell types tested, the detailed response of cell growth rate to cell mass differs. However, in general, when it falls below that for exponential growth, the natural increase in the CV of cell mass is effectively constrained. We find that both mass-dependent cell cycle regulation and mass-dependent growth rate modulation contribute to reducing cell mass variation within the population. Through the interplay and coordination of these 2 processes, accurate cell mass homeostasis emerges. Such findings reveal previously unappreciated and very general principles of cell size control in proliferating cells. These same regulatory processes might also be operative in terminally differentiated cells. Further quantitative dynamical studies should lead to a better understanding of the underlying molecular mechanisms of cell size control.

Hsf1 and the molecular chaperone Hsp90 support a 'rewiring stress response' leading to an adaptive cell size increase in chronic stress

Cells are exposed to a wide variety of internal and external stresses. Although many studies have focused on cellular responses to acute and severe stresses, little is known about how cellular systems adapt to sublethal chronic stresses. Using mammalian cells in culture, we discovered that they adapt to chronic mild stresses of up to two weeks, notably proteotoxic stresses such as heat, by increasing their size and translation, thereby scaling the amount of total protein. These adaptations render them more resilient to persistent and subsequent stresses. We demonstrate that Hsf1, well known for its role in acute stress responses, is required for the cell size increase, and that the molecular chaperone Hsp90 is essential for coupling the cell size increase to augmented translation. We term this translational reprogramming the 'rewiring stress response', and propose that this protective process of chronic stress adaptation contributes to the increase in size as cells get older, and that its failure promotes aging.

Cytosolic concentrations of actin binding proteins and the implications for in vivo F-actin turnover

Understanding how numerous actin-binding proteins (ABPs) work in concert to control the assembly, organization, and turnover of the actin cytoskeleton requires quantitative information about the levels of each component. Here, we measured the cellular concentrations of actin and the majority of the conserved ABPs in Saccharomyces cerevisiae, as well as the free (cytosolic) fractions of each ABP. The cellular concentration of actin is estimated to be 13.2 µM, with approximately two-thirds in the F-actin form and one-third in the G-actin form. Cellular concentrations of ABPs range from 12.4 to 0.85 µM (Tpm1> Pfy1> Cof1> Abp1> Srv2> Abp140> Tpm2> Aip1> Cap1/2> Crn1> Sac6> Twf1> Arp2/3> Scp1). The cytosolic fractions of all ABPs are unexpectedly high (0.6-0.9) and remain so throughout the cell cycle. Based on these numbers, we speculate that F-actin binding sites are limited in vivo, which leads to high cytosolic levels of ABPs, and in turn helps drive the rapid assembly and turnover of cellular F-actin structures.

RNA polymerase II dynamics and mRNA stability feedback scale mRNA amounts with cell size

A fundamental feature of cellular growth is that total protein and RNA amounts increase with cell size to keep concentrations approximately constant. A key component of this is that global transcription rates increase in larger cells. Here, we identify RNA polymerase II (RNAPII) as the limiting factor scaling mRNA transcription with cell size in budding yeast, as transcription is highly sensitive to the dosage of RNAPII but not to other components of the transcriptional machinery. Our experiments support a dynamic equilibrium model where global RNAPII transcription at a given size is set by the mass action recruitment kinetics of unengaged nucleoplasmic RNAPII to the genome. However, this only drives a sub-linear increase in transcription with size, which is then partially compensated for by a decrease in mRNA decay rates as cells enlarge. Thus, limiting RNAPII and feedback on mRNA stability work in concert to scale mRNA amounts with cell size.

Physiological stress drives the emergence of a Salmonella subpopulation through ribosomal RNA regulation

Bacteria undergo cycles of growth and starvation to which they must adapt swiftly. One important strategy for adjusting growth rates relies on ribosomal levels. Although high ribosomal levels are required for fast growth, their dynamics during starvation remain unclear. Here, we analyzed ribosomal RNA (rRNA) content of individual Salmonella cells by using fluorescence in situ hybridization (rRNA-FISH) and measured a dramatic decrease in rRNA numbers only in a subpopulation during nutrient limitation, resulting in a bimodal distribution of cells with high and low rRNA content. During nutritional upshifts, the two subpopulations were associated with distinct phenotypes. Using a transposon screen coupled with rRNA-FISH, we identified two mutants, DksA and RNase I, acting on rRNA transcription shutdown and degradation, which abolished the formation of the subpopulation with low rRNA content. Our work identifies a bacterial mechanism for regulation of ribosomal bimodality that may be beneficial for population survival during starvation.

Intermediate metabolites and molecular correlates of one‑carbon and nutrient metabolism differ in tissues from Holstein fetuses

Methionine and folate cycles along with transsulfuration comprise the one‑carbon metabolism (OCM) pathway. Amino acids and other nutrients feed into OCM, which is central to cellular function. mRNA abundance, proteins (Western blotting), and metabolites (GC-MC) associated with OCM were used to characterize these mechanisms in fetal tissues. Liver, whole intestine, and semitendinosus muscle were harvested from fetuses in 6 multiparous Holstein cows (37 kg milk/d, 100 d gestation). Data were analyzed using PROC MIXED (SAS 9.4). Protein abundance of BHMT was greatest (P < 0.01) in liver suggesting active remethylation of homocysteine to methionine. This idea was supported by the greater (P < 0.05) mRNA of CBS, BHMT, MTR, SHMT1, and MAT1A (encoding OCM enzymes) in liver. The antioxidant protein GPX3 had greatest (P < 0.05) abundance in liver, whereas the glutathione-transferase GSTM1 was 5-fold greater (P < 0.05) in intestine than liver and muscle. Greatest concentrations of glycine, serine, and taurine along with lower cysteine underscored the relevance of OCM in fetal liver. Phosphoethanolamine concentration was greatest (4-fold, P < 0.05) in intestine and along with the greatest (P < 0.05) mRNA of SLC44A1 (choline transporter), CHKA, and CEPT1 underscored the importance of the CDP-choline pathway. Greatest (P < 0.05) mRNA of PPARA, CPT1A, and HMGCS2 along with lower PCK1 in liver highlighted a potential reliance on fatty acid oxidation. In contrast, greater (P < 0.05) concentration of myo-inositol in muscle and intestine suggested both tissues rely on glucose as main source of energy. Future research should address how environmental inputs such as maternal nutrition alter these pathways in fetal tissues and their phenotypic outcomes.

Genome dilution by cell growth drives starvation-like proteome remodeling in mammalian and yeast cells

Cell size is tightly controlled in healthy tissues and single-celled organisms, but it remains unclear how size influences cell physiology. Increasing cell size was recently shown to remodel the proteomes of cultured human cells, demonstrating that large and small cells of the same type can be biochemically different. Here, we corroborate these results in mouse hepatocytes and extend our analysis using yeast. We find that size-dependent proteome changes are highly conserved and mostly independent of metabolic state. As eukaryotic cells grow larger, the dilution of the genome elicits a starvation-like proteome phenotype, suggesting that growth in large cells is limited by the genome in a manner analogous to a limiting nutrient. We also demonstrate that the proteomes of replicatively-aged yeast are primarily determined by their large size. Overall, our data suggest that genome concentration is a universal determinant of proteome content in growing cells.

Measuring prion propagation in single bacteria elucidates a mechanism of loss

Prions are self-propagating protein aggregates formed by specific proteins that can adopt alternative folds. Prions were discovered as the cause of the fatal transmissible spongiform encephalopathies in mammals, but prions can also constitute nontoxic protein-based elements of inheritance in fungi and other species. Prion propagation has recently been shown to occur in bacteria for more than a hundred cell divisions, yet a fraction of cells in these lineages lost the prion through an unknown mechanism. Here, we investigate prion propagation in single bacterial cells as they divide using microfluidics and fluorescence microscopy. We show that the propagation occurs in two distinct modes. In a fraction of the population, cells had multiple small visible aggregates and lost the prion through random partitioning of aggregates to one of the two daughter cells at division. In the other subpopulation, cells had a stable large aggregate localized to the pole; upon division the mother cell retained this polar aggregate and a daughter cell was generated that contained small aggregates. Extending our findings to prion domains from two orthologous proteins, we observe similar propagation and loss properties. Our findings also provide support for the suggestion that bacterial prions can form more than one self-propagating state. We implement a stochastic version of the molecular model of prion propagation from yeast and mammals that recapitulates all the observed single-cell properties. This model highlights challenges for prion propagation that are unique to prokaryotes and illustrates the conservation of fundamental characteristics of prion propagation.

Bacterial cell-size changes resulting from altering the relative expression of Min proteins

The timing of cell division, and thus cell size in bacteria, is determined in part by the accumulation dynamics of the protein FtsZ, which forms the septal ring. FtsZ localization depends on membrane-associated Min proteins, which inhibit FtsZ binding to the cell pole membrane. Changes in the relative concentrations of Min proteins can disrupt FtsZ binding to the membrane, which in turn can delay cell division until a certain cell size is reached, in which the dynamics of Min proteins frees the cell membrane long enough to allow FtsZ ring formation. Here, we study the effect of Min proteins relative expression on the dynamics of FtsZ ring formation and cell size in individual Escherichia coli bacteria. Upon inducing overexpression of minE, cell size increases gradually to a new steady-state value. Concurrently, the time required to initiate FtsZ ring formation grows as the size approaches the new steady-state, at which point the ring formation initiates as early as before induction. These results highlight the contribution of Min proteins to cell size control, which may be partially responsible for the size fluctuations observed in bacterial populations, and may clarify how the size difference acquired during asymmetric cell division is offset.

Optimizing Arsenic Therapy by Selectively Targeting Leukemia Cells

Arsenic, in the simple form of arsenic trioxide, is currently marketed for the treatment of acute promyelocytic leukemia. Due to the multifaceted mechanisms of action of arsenic, it has also shown promise in other types of leukemias but is hindered by its toxic effects toward normal cells. This research has aimed to determine whether tumor-homing peptide complexes of arsenic can be designed and developed to strategically target specific cancers. The end goal is to achieve dose reduction and decreased side effects of the resultant arsenic therapeutic agent. In this article, we present the synthesis, characterization, and stability studies of a new class of As-peptide complexes designed to target leukemia. In vitro biological studies of the most stable complex show 1000 times greater toxicity toward leukemia cells over human blood cells, indicating potential for progression to in vivo studies.

Out-of-equilibrium gene expression fluctuations in the presence of extrinsic noise

Cell-to-cell variability in protein concentrations is strongly affected by extrinsic noise, especially for highly expressed genes. Extrinsic noise can be due to fluctuations of several possible cellular factors connected to cell physiology and to the level of key enzymes in the expression process. However, how to identify the predominant sources of extrinsic noise in a biological system is still an open question. This work considers a general stochastic model of gene expression with extrinsic noise represented as fluctuations of the different model rates, and focuses on the out-of-equilibrium expression dynamics. Combining analytical calculations with stochastic simulations, we characterize how extrinsic noise shapes the protein variability during gene activation or inactivation, depending on the prevailing source of extrinsic variability, on its intensity and timescale. In particular, we show that qualitatively different noise profiles can be identified depending on which are the fluctuating parameters. This indicates an experimentally accessible way to pinpoint the dominant sources of extrinsic noise using time-coarse experiments.

Implications of differential size-scaling of cell-cycle regulators on cell size homeostasis

Accurate timing of division and size homeostasis is crucial for cells. A potential mechanism for cells to decide the timing of division is the differential scaling of regulatory protein copy numbers with cell size. However, it remains unclear whether such a mechanism can lead to robust growth and division, and how the scaling behaviors of regulatory proteins influence the cell size distribution. Here we study a mathematical model combining gene expression and cell growth, in which the cell-cycle activators scale superlinearly with cell size while the inhibitors scale sublinearly. The cell divides once the ratio of their concentrations reaches a threshold value. We find that the cell can robustly grow and divide within a finite range of the threshold value with the cell size proportional to the ploidy. In a stochastic version of the model, the cell size at division is uncorrelated with that at birth. Also, the more differential the cell-size scaling of the cell-cycle regulators is, the narrower the cell-size distribution is. Intriguingly, our model with multiple regulators rationalizes the observation that after the deletion of a single regulator, the coefficient of variation of cell size remains roughly the same though the average cell size changes significantly. Our work reveals that the differential scaling of cell-cycle regulators provides a robust mechanism of cell size control.

Scaling of nuclear numbers and their spatial arrangement in skeletal muscle cell size regulation

Many cells display considerable functional plasticity and depend on the regulation of numerous organelles and macromolecules for their maintenance. In large cells, organelles also need to be carefully distributed to supply the cell with essential resources and regulate intracellular activities. Having multiple copies of the largest eukaryotic organelle, the nucleus, epitomizes the importance of scaling gene products to large cytoplasmic volumes in skeletal muscle fibers. Scaling of intracellular constituents within mammalian muscle fibers is, however, poorly understood, but according to the myonuclear domain hypothesis, a single nucleus supports a finite amount of cytoplasm and is thus postulated to act autonomously, causing the nuclear number to be commensurate with fiber volume. In addition, the orderly peripheral distribution of myonuclei is a hallmark of normal cell physiology, as nuclear mispositioning is associated with impaired muscle function. Because underlying structures of complex cell behaviors are commonly formalized by scaling laws and thus emphasize emerging principles of size regulation, the work presented herein offers more of a unified conceptual platform based on principles from physics, chemistry, geometry, and biology to explore cell size-dependent correlations of the largest mammalian cell by means of scaling.

PyEcoLib: a python library for simulating stochastic cell size dynamics

Recently, there has been an increasing need for tools to simulate cell size regulation due to important applications in cell proliferation and gene expression. However, implementing the simulation usually presents some difficulties, as the division has a cycle-dependent occurrence rate. In this article, we gather a recent theoretical framework inPyEcoLib, a python-based library to simulate the stochastic dynamics of the size of bacterial cells. This library can simulate cell size trajectories with an arbitrarily small sampling period. In addition, this simulator can include stochastic variables, such as the cell size at the beginning of the experiment, the cycle duration timing, the growth rate, and the splitting position. Furthermore, from a population perspective, the user can choose between tracking a single lineage or all cells in a colony. They can also simulate the most common division strategies (adder, timer, and sizer) using the division rate formalism and numerical methods. As an example of PyecoLib applications, we explain how to couple size dynamics with gene expression predicting, from simulations, how the noise in protein levels increases by increasing the noise in division timing, the noise in growth rate and the noise in cell splitting position. The simplicity of this library and its transparency about the underlying theoretical framework yield the inclusion of cell size stochasticity in complex models of gene expression.

Mycobacterium tuberculosis grows linearly at the single-cell level with larger variability than model organisms

The ability of bacterial pathogens to regulate growth is crucial to control homeostasis, virulence, and drug response. Yet, we do not understand the growth and cell cycle behaviors of Mycobacterium tuberculosis (Mtb), a slow-growing pathogen, at the single-cell level. Here, we use time-lapse imaging and mathematical modeling to characterize these fundamental properties of Mtb. Whereas most organisms grow exponentially at the single-cell level, we find that Mtb exhibits a unique linear growth mode. Mtb growth characteristics are highly variable from cell-to-cell, notably in their growth speeds, cell cycle timing, and cell sizes. Together, our study demonstrates that growth behavior of Mtb diverges from what we have learned from model bacteria. Instead, Mtb generates a heterogeneous population while growing slowly and linearly. Our study provides a new level of detail into how Mtb grows and creates heterogeneity, and motivates more studies of growth behaviors in bacterial pathogens.

Engineering synthetic biomolecular condensates

The concept of phase-separation-mediated formation of biomolecular condensates provides a new framework to understand cellular organization and cooperativity-dependent cellular functions. With growing understanding of how biological systems drive phase separation and how cellular functions are encoded by biomolecular condensates, opportunities have emerged for cellular control through engineering of synthetic biomolecular condensates. In this Review, we discuss how to construct synthetic biomolecular condensates and how they can regulate cellular functions. We first describe the fundamental principles by which biomolecular components can drive phase separation. Next, we discuss the relationship between the properties of condensates and their cellular functions, which informs the design of components to create programmable synthetic condensates. Finally, we describe recent applications of synthetic biomolecular condensates for cellular control and discuss some of the design considerations and prospective applications.

Predictive model of transcriptional elongation control identifies trans regulatory factors from chromatin signatures

Promoter-proximal Polymerase II (Pol II) pausing is a key rate-limiting step for gene expression. DNA and RNA-binding trans-acting factors regulating the extent of pausing have been identified. However, we lack a quantitative model of how interactions of these factors determine pausing, therefore the relative importance of implicated factors is unknown. Moreover, previously unknown regulators might exist. Here we address this gap with a machine learning model that accurately predicts the extent of promoter-proximal Pol II pausing from large-scale genome and transcriptome binding maps and gene annotation and sequence composition features. We demonstrate high accuracy and generalizability of the model by validation on an independent cell line which reveals the model's cell line agnostic character. Model interpretation in light of prior knowledge about molecular functions of regulatory factors confirms the interconnection of pausing with other RNA processing steps. Harnessing underlying feature contributions, we assess the relative importance of each factor, quantify their predictive effects and systematically identify previously unknown regulators of pausing. We additionally identify 16 previously unknown 7SK ncRNA interacting RNA-binding proteins predictive of pausing. Our work provides a framework to further our understanding of the regulation of the critical early steps in transcriptional elongation.

The myonuclear domain in adult skeletal muscle fibres: past, present and future

Most cells in the body are mononuclear whereas skeletal muscle fibres are uniquely multinuclear. The nuclei of muscle fibres (myonuclei) are usually situated peripherally which complicates the equitable distribution of gene products. Myonuclear abundance can also change under conditions such as hypertrophy and atrophy. Specialised zones in muscle fibres have different functions and thus distinct synthetic demands from myonuclei. The complex structure and regulatory requirements of multinuclear muscle cells understandably led to the hypothesis that myonuclei govern defined 'domains' to maintain homeostasis and facilitate adaptation. The purpose of this review is to provide historical context for the myonuclear domain and evaluate its veracity with respect to mRNA and protein distribution resulting from myonuclear transcription. We synthesise insights from past and current in vitro and in vivo genetically modified models for studying the myonuclear domain under dynamic conditions. We also cover the most contemporary knowledge on mRNA and protein transport in muscle cells. Insights from emerging technologies such as single myonuclear RNA-sequencing further inform our discussion of the myonuclear domain. We broadly conclude: (1) the myonuclear domain can be flexible during muscle fibre growth and atrophy, (2) the mechanisms and role of myonuclear loss and motility deserve further consideration, (3) mRNA in muscle is actively transported via microtubules and locally restricted, but proteins may travel far from a myonucleus of origin and (4) myonuclear transcriptional specialisation extends beyond the classic neuromuscular and myotendinous populations. A deeper understanding of the myonuclear domain in muscle may promote effective therapies for ageing and disease.

Biology-inspired data-driven quality control for scientific discovery in single-cell transcriptomics

BACKGROUND

Quality control (QC) of cells, a critical first step in single-cell RNA sequencing data analysis, has largely relied on arbitrarily fixed data-agnostic thresholds applied to QC metrics such as gene complexity and fraction of reads mapping to mitochondrial genes. The few existing data-driven approaches perform QC at the level of samples or studies without accounting for biological variation.

RESULTS

We first demonstrate that QC metrics vary with both tissue and cell types across technologies, study conditions, and species. We then propose data-driven QC (ddqc), an unsupervised adaptive QC framework to perform flexible and data-driven QC at the level of cell types while retaining critical biological insights and improved power for downstream analysis. ddqc applies an adaptive threshold based on the median absolute deviation on four QC metrics (gene and UMI complexity, fraction of reads mapping to mitochondrial and ribosomal genes). ddqc retains over a third more cells when compared to conventional data-agnostic QC filters. Finally, we show that ddqc recovers biologically meaningful trends in gradation of gene complexity among cell types that can help answer questions of biological interest such as which cell types express the least and most number of transcripts overall, and ribosomal transcripts specifically.

CONCLUSIONS

ddqc retains cell types such as metabolically active parenchymal cells and specialized cells such as neutrophils which are often lost by conventional QC. Taken together, our work proposes a revised paradigm to quality filtering best practices-iterative QC, providing a data-driven QC framework compatible with observed biological diversity.

Principles of gene regulation quantitatively connect DNA to RNA and proteins in bacteria

Protein concentrations are set by a complex interplay between gene-specific regulatory processes and systemic factors, including cell volume and shared gene expression machineries. Elucidating this interplay is crucial for discerning and designing gene regulatory systems. We quantitatively characterized gene-specific and systemic factors that affect transcription and translation genome-wide for Escherichia coli across many conditions. The results revealed two design principles that make regulation of gene expression insulated from concentrations of shared machineries: RNA polymerase activity is fine-tuned to match translational output, and translational characteristics are similar across most messenger RNAs (mRNAs). Consequently, in bacteria, protein concentration is set primarily at the promoter level. A simple mathematical formula relates promoter activities and protein concentrations across growth conditions, enabling quantitative inference of gene regulation from omics data.

Robust replication initiation from coupled homeostatic mechanisms

The bacterium Escherichia coli initiates replication once per cell cycle at a precise volume per origin and adds an on average constant volume between successive initiation events, independent of the initiation size. Yet, a molecular model that can explain these observations has been lacking. Experiments indicate that E. coli controls replication initiation via titration and activation of the initiator protein DnaA. Here, we study by mathematical modelling how these two mechanisms interact to generate robust replication-initiation cycles. We first show that a mechanism solely based on titration generates stable replication cycles at low growth rates, but inevitably causes premature reinitiation events at higher growth rates. In this regime, the DnaA activation switch becomes essential for stable replication initiation. Conversely, while the activation switch alone yields robust rhythms at high growth rates, titration can strongly enhance the stability of the switch at low growth rates. Our analysis thus predicts that both mechanisms together drive robust replication cycles at all growth rates. In addition, it reveals how an origin-density sensor yields adder correlations.

Eukaryotic Cell Size Control and Its Relation to Biosynthesis and Senescence

The most fundamental feature of cellular form is size, which sets the scale of all cell biological processes. Growth, form, and function are all necessarily linked in cell biology, but we often do not understand the underlying molecular mechanisms nor their specific functions. Here, we review progress toward determining the molecular mechanisms that regulate cell size in yeast, animals, and plants, as well as progress toward understanding the function of cell size regulation. It has become increasingly clear that the mechanism of cell size regulation is deeply intertwined with basic mechanisms of biosynthesis, and how biosynthesis can be scaled (or not) in proportion to cell size. Finally, we highlight recent findings causally linking aberrant cell size regulation to cellular senescence and their implications for cancer therapies.

Concentration fluctuations in growing and dividing cells: Insights into the emergence of concentration homeostasis

Intracellular reaction rates depend on concentrations and hence their levels are often regulated. However classical models of stochastic gene expression lack a cell size description and cannot be used to predict noise in concentrations. Here, we construct a model of gene product dynamics that includes a description of cell growth, cell division, size-dependent gene expression, gene dosage compensation, and size control mechanisms that can vary with the cell cycle phase. We obtain expressions for the approximate distributions and power spectra of concentration fluctuations which lead to insight into the emergence of concentration homeostasis. We find that (i) the conditions necessary to suppress cell division-induced concentration oscillations are difficult to achieve; (ii) mRNA concentration and number distributions can have different number of modes; (iii) two-layer size control strategies such as sizer-timer or adder-timer are ideal because they maintain constant mean concentrations whilst minimising concentration noise; (iv) accurate concentration homeostasis requires a fine tuning of dosage compensation, replication timing, and size-dependent gene expression; (v) deviations from perfect concentration homeostasis show up as deviations of the concentration distribution from a gamma distribution. Some of these predictions are confirmed using data for E. coli, fission yeast, and budding yeast.

Increasing cell size remodels the proteome and promotes senescence

Cell size is tightly controlled in healthy tissues, but it is unclear how deviations in cell size affect cell physiology. To address this, we measured how the cell's proteome changes with increasing cell size. Size-dependent protein concentration changes are widespread and predicted by subcellular localization, size-dependent mRNA concentrations, and protein turnover. As proliferating cells grow larger, concentration changes typically associated with cellular senescence are increasingly pronounced, suggesting that large size may be a cause rather than just a consequence of cell senescence. Consistent with this hypothesis, larger cells are prone to replicative, DNA-damage-induced, and CDK4/6i-induced senescence. Size-dependent changes to the proteome, including those associated with senescence, are not observed when an increase in cell size is accompanied by an increase in ploidy. Together, our findings show how cell size could impact many aspects of cell physiology by remodeling the proteome and provide a rationale for cell size control and polyploidization.

Scaling of biosynthesis and metabolism with cell size

Cells adopt a size that is optimal for their function, and pushing them beyond this limit can cause cell aging and death by senescence or reduce proliferative potential. However, by increasing their genome copy number (ploidy), cells can increase their size dramatically and homeostatically maintain physiological properties such as biosynthesis rate. Recent studies investigating the relationship between cell size and rates of biosynthesis and metabolism under normal, polyploid, and pathological conditions are revealing new insights into how cells attain the best function or fitness for their size by tuning processes including transcription, translation, and mitochondrial respiration. A new frontier is to connect single-cell scaling relationships with tissue and whole-organism physiology, which promises to reveal molecular and evolutionary principles underlying the astonishing diversity of size observed across the tree of life.

Nucleo-cytoplasmic shuttling of RNA-binding factors: mRNA buffering and beyond

Gene expression is a highly regulated process that adapts RNAs and proteins content to the cellular context. Under steady-state conditions, mRNA homeostasis is robustly maintained by tight controls that act on both nuclear transcription and cytoplasmic mRNA stability. In recent years, it has been revealed that several RNA-binding proteins (RBPs) that perform functions in mRNA decay can move to the nucleus and regulate transcription. The RBPs involved in transcription can also travel to the cytoplasm and regulate mRNA degradation and/or translation. The multifaceted functions of these shuttling nucleo-cytoplasm RBPs have raised the possibility that they can act as mRNA metabolism coordinators. In addition, this indicates the existence of crosstalk mechanisms between the enzymatic machineries that drive the different mRNA life-cycle phases. The buffering of the mRNA concentration is the best known consequence of a transcription-degradation crosstalk counteraction, but alternative ways of RBP action can also imply enhanced gene regulation.

Mechanisms of cellular mRNA transcript homeostasis

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

Feedback from nuclear RNA on transcription promotes robust RNA concentration homeostasis in human cells

RNA concentration homeostasis involves coordinating RNA abundance and synthesis rates with cell size. Here, we study this in human cells by combining genome-wide perturbations with quantitative single-cell measurements. Despite relative ease in perturbing RNA synthesis, we find that RNA concentrations generally remain highly constant. Perturbations that would be expected to increase nuclear mRNA levels, including those targeting nuclear mRNA degradation or export, result in downregulation of RNA synthesis. This is associated with reduced abundance of transcription-associated proteins and protein states that are normally coordinated with RNA production in single cells, including RNA polymerase II (RNA Pol II) itself. Acute perturbations, elevation of nuclear mRNA levels, and mathematical modeling indicate that mammalian cells achieve robust mRNA concentration homeostasis by the mRNA-based negative feedback on transcriptional activity in the nucleus. This ultimately acts to coordinate RNA Pol II abundance with nuclear mRNA degradation and export rates and may underpin the scaling of mRNA abundance with cell size.

In vivo Evaluation of Non-viral NICD Plasmid-Loaded PLGA Nanoparticles in Developing Zebrafish to Improve Cardiac Functions

In the era of the advanced nanomaterials, use of nanoparticles has been highlighted in biomedical research. However, the demonstration of DNA plasmid delivery with nanoparticles for in vivo gene delivery experiments must be carefully tested due to many possible issues, including toxicity. The purpose of the current study was to deliver a Notch Intracellular Domain (NICD)-encoded plasmid via poly(lactic-co-glycolic acid) (PLGA) nanoparticles and to investigate the toxic environmental side effects for an in vivo experiment. In addition, we demonstrated the target delivery to the endothelium, including the endocardial layer, which is challenging to manipulate gene expression for cardiac functions due to the beating heart and rapid blood pumping. For this study, we used a zebrafish animal model and exposed it to nanoparticles at varying concentrations to observe for specific malformations over time for toxic effects of PLGA nanoparticles as a delivery vehicle. Our nanoparticles caused significantly less malformations than the positive control, ZnO nanoparticles. Additionally, the NICD plasmid was successfully delivered by PLGA nanoparticles and significantly increased Notch signaling related genes. Furthermore, our image based deep-learning analysis approach evaluated that the antibody conjugated nanoparticles were successfully bound to the endocardium to overexpress Notch related genes and improve cardiac function such as ejection fraction, fractional shortening, and cardiac output. This research demonstrates that PLGA nanoparticle-mediated target delivery to upregulate Notch related genes which can be a potential therapeutic approach with minimum toxic effects.

Double-edged role of resource competition in gene expression noise and control

Despite extensive investigation demonstrating that resource competition can significantly alter the deterministic behaviors of synthetic gene circuits, it remains unclear how resource competition contributes to the gene expression noise and how this noise can be controlled. Utilizing a two-gene circuit as a prototypical system, we uncover a surprising double-edged role of resource competition in gene expression noise: competition decreases noise through introducing a resource constraint but generates its own type of noise which we name as "resource competitive noise." Utilization of orthogonal resources enables retainment of the noise reduction conferred by resource constraint while removing the added resource competitive noise. The noise reduction effects are studied using three negative feedback types: negatively competitive regulation (NCR), local, and global controllers, each having four placement architectures in the protein biosynthesis pathway (mRNA or protein inhibition on transcription or translation). Our results show that both local and NCR controllers with mRNA-mediated inhibition are efficacious at reducing noise, with NCR controllers demonstrating a superior noise-reduction capability. We also find that combining feedback controllers with orthogonal resources can improve the local controllers. This work provides deep insights into the origin of stochasticity in gene circuits with resource competition and guidance for developing effective noise control strategies.

Protein concentration fluctuations in the high expression regime: Taylor's law and its mechanistic origin

Protein concentration in a living cell fluctuates over time due to noise in growth and division processes. In the high expression regime, variance of the protein concentration in a cell was found to scale with the square of the mean, which belongs to a general phenomenon called Taylor's law (TL). To understand the origin for these fluctuations, we measured protein concentration dynamics in single E. coli cells from a set of strains with a variable expression of fluorescent proteins. The protein expression is controlled by a set of constitutive promoters with different strength, which allows to change the expression level over 2 orders of magnitude without introducing noise from fluctuations in transcription regulators. Our data confirms the square TL, but the prefactor A has a cell-to-cell variation independent of the promoter strength. Distributions of the normalized protein concentration for different promoters are found to collapse onto the same curve. To explain these observations, we used a minimal mechanistic model to describe the stochastic growth and division processes in a single cell with a feedback mechanism for regulating cell division. In the high expression regime where extrinsic noise dominates, the model reproduces our experimental results quantitatively. By using a mean-field approximation in the minimal model, we showed that the stochastic dynamics of protein concentration is described by a Langevin equation with multiplicative noise. The Langevin equation has a scale invariance which is responsible for the square TL. By solving the Langevin equation, we obtained an analytical solution for the protein concentration distribution function that agrees with experiments. The solution shows explicitly how the prefactor A depends on strength of different noise sources, which explains its cell-to-cell variability. By using this approach to analyze our single-cell data, we found that the noise in production rate dominates the noise from cell division. The deviation from the square TL in the low expression regime can also be captured in our model by including intrinsic noise in the production rate.

Spatiotemporal localization of proteins in mycobacteria

Although prokaryotic organisms lack traditional organelles, they must still organize cellular structures in space and time, challenges that different species solve differently. To systematically define the subcellular architecture of mycobacteria, we perform high-throughput imaging of a library of fluorescently tagged proteins expressed in Mycobacterium smegmatis and develop a customized computational pipeline, MOMIA and GEMATRIA, to analyze these data. Our results establish a spatial organization network of over 700 conserved mycobacterial proteins and reveal a coherent localization pattern for many proteins of known function, including those in translation, energy metabolism, cell growth and division, as well as proteins of unknown function. Furthermore, our pipeline exploits morphologic proxies to enable a pseudo-temporal approximation of protein localization and identifies previously uncharacterized cell-cycle-dependent dynamics of essential mycobacterial proteins. Collectively, these data provide a systems perspective on the subcellular organization of mycobacteria and provide tools for the analysis of bacteria with non-standard growth characteristics.

Zhu et al. develop a two-stage image analysis pipeline, MOMIA and GEMATRIA, that efficiently models the spatial and temporal dynamics of over 700 conserved proteins in M. smegmatis. Through the analysis they report spatial constraints of mycobacterial ribosomes and membrane complexes and reconstruct temporal dynamics from still image data.

Cell cycle progression and transmitotic apoptosis resistance promote escape from extrinsic apoptosis

Extrinsic apoptosis relies on TNF-family receptor activation by immune cells or receptor-activating drugs. Here, we monitored cell cycle progression at a resolution of minutes to relate apoptosis kinetics and cell-to-cell heterogeneities in death decisions to cell cycle phases. Interestingly, we found that cells in S phase delay TRAIL receptor-induced death in favour of mitosis, thereby passing on an apoptosis-primed state to their offspring. This translates into two distinct fates, apoptosis execution post mitosis or cell survival from inefficient apoptosis. Transmitotic resistance is linked to Mcl-1 upregulation and its increased accumulation at mitochondria from mid-S phase onwards, which allows cells to pass through mitosis with activated caspase-8, and with cells escaping apoptosis after mitosis sustaining sublethal DNA damage. Antagonizing Mcl-1 suppresses cell cycle-dependent delays in apoptosis, prevents apoptosis-resistant progression through mitosis and averts unwanted survival after apoptosis induction. Cell cycle progression therefore modulates signal transduction during extrinsic apoptosis, with Mcl-1 governing decision making between death, proliferation and survival. Cell cycle progression thus is a crucial process from which cell-to-cell heterogeneities in fates and treatment outcomes emerge in isogenic cell populations during extrinsic apoptosis. This article has an associated First Person interview with the first author of the paper.

Distinguishing different modes of growth using single-cell data

Collection of high-throughput data has become prevalent in biology. Large datasets allow the use of statistical constructs such as binning and linear regression to quantify relationships between variables and hypothesize underlying biological mechanisms based on it. We discuss several such examples in relation to single-cell data and cellular growth. In particular, we show instances where what appears to be ordinary use of these statistical methods leads to incorrect conclusions such as growth being non-exponential as opposed to exponential and vice versa. We propose that the data analysis and its interpretation should be done in the context of a generative model, if possible. In this way, the statistical methods can be validated either analytically or against synthetic data generated via the use of the model, leading to a consistent method for inferring biological mechanisms from data. On applying the validated methods of data analysis to infer cellular growth on our experimental data, we find the growth of length in E. coli to be non-exponential. Our analysis shows that in the later stages of the cell cycle the growth rate is faster than exponential.

Heterogeneous recruitment abilities to RNA polymerases generate nonlinear scaling of gene expression with cell volume

While most genes' expression levels are proportional to cell volumes, some genes exhibit nonlinear scaling between their expression levels and cell volume. Therefore, their mRNA and protein concentrations change as the cell volume increases, which often have crucial biological functions such as cell-cycle regulation. However, the biophysical mechanism underlying the nonlinear scaling between gene expression and cell volume is still unclear. In this work, we show that the nonlinear scaling is a direct consequence of the heterogeneous recruitment abilities of promoters to RNA polymerases based on a gene expression model at the whole-cell level. Those genes with weaker (stronger) recruitment abilities than the average ability spontaneously exhibit superlinear (sublinear) scaling with cell volume. Analysis of the promoter sequences and the nonlinear scaling of Saccharomyces cerevisiae's mRNA levels shows that motifs associated with transcription regulation are indeed enriched in genes exhibiting nonlinear scaling, in concert with our model.

The Single-Cell Transcriptomic Analysis of Prefrontal Pyramidal Cells and Interneurons Reveals the Neuronal Expression of Genes Encoding Antimicrobial Peptides and Immune Proteins

The investigation of the molecular background of direct communication of neurons and immune cells in the brain is an important issue for understanding physiological and pathological processes in the nervous system. Direct contacts between brain-infiltrating immune cells and neurons, and the neuromodulatory effect of immune cell-derived regulatory peptides are well established. Several aspects of the role of immune and glial cells in the direct neuro-immune communication are also well known; however, there remain many questions regarding the molecular details of signaling from neurons to immune cells. Thus, we report here on the neuronal expression of genes encoding antimicrobial and immunomodulatory peptides, as well as proteins of immune cell-specific activation and communication mechanisms. In the present study, we analyzed the single-cell sequencing data of our previous transcriptomic work, obtained from electrophysiologically identified pyramidal cells and interneurons of the murine prefrontal cortex. We filtered out the genes that may be associated with the direct communication between immune cells and neurons and examined their expression pattern in the neuronal transcriptome. The expression of some of these genes by cortical neurons has not yet been reported. The vast majority of antimicrobial (~53%) and immune cell protein (~94%) transcripts was identified in the transcriptome of the 84 cells, owing to the high sensitivity of ultra-deep sequencing. Several of the antimicrobial and immune process-related protein transcripts showed cell type-specific or enriched expression. Individual neurons transcribed only a fraction of the investigated genes with low copy numbers probably due to the bursting kinetics of gene expression; however, the comparison of our data with available transcriptomic datasets from immune cells and neurons suggests the functional relevance of the reported findings. Accordingly, we propose further experimental and in silico studies on the neuronal expression of immune system-related genes and the potential role of the encoded proteins in neuroimmunological processes.

The total mRNA concentration buffering system in yeast is global rather than gene-specific

Gene expression in eukaryotes does not follow a linear process from transcription to translation and mRNA degradation. Instead it follows a circular process in which cytoplasmic mRNA decay crosstalks with nuclear transcription. In many instances, this crosstalk contributes to buffer mRNA at a roughly constant concentration. Whether the mRNA buffering concept operates on the total mRNA concentration or at the gene-specific level, and if the mechanism to do so is a global or a specific one, remain unknown. Here we assessed changes in mRNA concentrations and their synthesis rates along the transcriptome of aneuploid strains of the yeast Saccharomyces cerevisiae We also assessed mRNA concentrations and their synthesis rates in nonsense-mediated decay (NMD) targets in euploid strains. We found that the altered synthesis rates in the genes from the aneuploid chromosome and the changes in their mRNA stabilities were not counterbalanced. In addition, the stability of NMD targets was not specifically compensated by the changes in synthesis rate. We conclude that there is no genetic compensation of NMD mRNA targets in yeast, and total mRNA buffering uses mostly a global system rather than a gene-specific one.

Fundamental limits on the rate of bacterial growth and their influence on proteomic composition

Despite abundant measurements of bacterial growth rate, cell size, and protein content, we lack a rigorous understanding of what sets the scale of these quantities and when protein abundances should (or should not) depend on growth rate. Here, we estimate the basic requirements and physical constraints on steady-state growth by considering key processes in cellular physiology across a collection of Escherichia coli proteomic data covering ≈4,000 proteins and 36 growth rates. Our analysis suggests that cells are predominantly tuned for the task of cell doubling across a continuum of growth rates; specific processes do not limit growth rate or dictate cell size. We present a model of proteomic regulation as a function of nutrient supply that reconciles observed interdependences between protein synthesis, cell size, and growth rate and propose that a theoretical inability to parallelize ribosomal synthesis places a firm limit on the achievable growth rate. A record of this paper's transparent peer review process is included in the supplemental information.

A single-cell analytical approach to quantify activated caspase-3/7 during osteoblast proliferation, differentiation, and apoptosis

The protein heterogeneity at the single-cell level has been recognized to be vital for an understanding of various life processes during animal development. In addition, the knowledge of accurate quantity of relevant proteins at cellular level is essential for appropriate interpretation of diagnostic and therapeutic results. Some low-copy-number proteins are known to play a crucial role during cell proliferation, differentiation, and also in apoptosis. The fate decision is often based on the concentration of these proteins in the individual cells. This is likely to apply also for caspases, cysteine proteases traditionally associated with cell death via apoptosis but recently being discovered also as important factors in cell proliferation and differentiation. The hypothesis was tested in bone-related cells, where modulation of fate from apoptosis to proliferation/differentiation and vice versa is particularly challenging, e.g., towards anti-osteoporotic treatments and anti-cancer strategies. An ultrasensitive and highly selective method based on bioluminescence photon counting was used to quantify activated caspase-3/7 in order to demonstrate protein-level heterogeneity in individual cells within one population and to associate quantitative measurements with different cell fates (proliferation, differentiation, apoptosis). The results indicate a gradual increase of caspase-3/7 activation from the proliferative status to differentiation (more than three times) and towards apoptosis (more than six times). The findings clearly support one of the putative key mechanisms of non-apoptotic functions of pro-apoptotic caspases based on fine-tuning of their activation levels.

Modeling the impact of single-cell stochasticity and size control on the population growth rate in asymmetrically dividing cells

Microbial populations show striking diversity in cell growth morphology and lifecycle; however, our understanding of how these factors influence the growth rate of cell populations remains limited. We use theory and simulations to predict the impact of asymmetric cell division, cell size regulation and single-cell stochasticity on the population growth rate. Our model predicts that coarse-grained noise in the single-cell growth rate λ decreases the population growth rate, as previously seen for symmetrically dividing cells. However, for a given noise in λ we find that dividing asymmetrically can enhance the population growth rate for cells with strong size control (between a “sizer” and an “adder”). To reconcile this finding with the abundance of symmetrically dividing organisms in nature, we propose that additional constraints on cell growth and division must be present which are not included in our model, and we explore the effects of selected extensions thereof. Further, we find that within our model, epigenetically inherited generation times may arise due to size control in asymmetrically dividing cells, providing a possible explanation for recent experimental observations in budding yeast. Taken together, our findings provide insight into the complex effects generated by non-canonical growth morphologies.

How rapidly a population of single-celled organisms can grow will strongly impact their long-term success. Prior work has shown that many factors impact this population growth rate, including the rate at which single cells grow, random variability between cells, and whether cells regulate their own size. Here we show that cell division asymmetry can also have a strong impact on the population growth rate. We use theory and computer simulations to study the growth rate of cells that divide asymmetrically, producing one smaller cell and one larger cell with each cell division event. We show that variability in how fast single cells grow will still decrease the population growth rate, when asymmetry is moderate or size control is weak, but that cells with strong size control can diminish this decrease by dividing more asymmetrically. We also demonstrate that cell cycle lengths can be positively correlated for closely related cells when they both divide asymmetrically and regulate their size. This counter-intuitive result contrasts with previous findings based on cell size regulation in symmetrically dividing cells that if cells grow for “too long” in one cell cycle, this will be corrected for by reduced growth during a shorter, subsequent cell cycle.

Limits and Constraints on Mechanisms of Cell-Cycle Regulation Imposed by Cell Size-Homeostasis Measurements

High-throughput imaging has led to an explosion of observations about cell-size homeostasis across the kingdoms of life. Among bacteria, "adder" behavior-in which a constant size increment appears to be added during each cell cycle-is ubiquitous, while various eukaryotes show other size-homeostasis behaviors. Since interactions between cell-cycle progression and growth ultimately determine such behaviors, we developed a general model of cell-cycle regulation. Our analyses reveal a range of scenarios that are plausible but fail to regulate cell size, indicating that mechanisms of cell-cycle regulation are stringently limited by size-control requirements, and possibly why certain cell-cycle features are strongly conserved. Cell-cycle features can play unintuitive roles in altering size-homeostasis behaviors: noisy regulator production can enhance adder behavior, while Whi5-like inhibitor dilutors respond sensitively to perturbations to G2/M control and noisy G1/S checkpoints. Our model thus provides holistic insights into the mechanistic implications of size-homeostasis experimental measurements.

Coordination of gene expression noise with cell size: analytical results for agent-based models of growing cell populations

The chemical master equation and the Gillespie algorithm are widely used to model the reaction kinetics inside living cells. It is thereby assumed that cell growth and division can be modelled through effective dilution reactions and extrinsic noise sources. We here re-examine these paradigms through developing an analytical agent-based framework of growing and dividing cells accompanied by an exact simulation algorithm, which allows us to quantify the dynamics of virtually any intracellular reaction network affected by stochastic cell size control and division noise. We find that the solution of the chemical master equation-including static extrinsic noise-exactly agrees with the agent-based formulation when the network under study exhibits stochastic concentration homeostasis, a novel condition that generalizes concentration homeostasis in deterministic systems to higher order moments and distributions. We illustrate stochastic concentration homeostasis for a range of common gene expression networks. When this condition is not met, we demonstrate by extending the linear noise approximation to agent-based models that the dependence of gene expression noise on cell size can qualitatively deviate from the chemical master equation. Surprisingly, the total noise of the agent-based approach can still be well approximated by extrinsic noise models.

Exponential trajectories, cell size fluctuations, and the adder property in bacteria follow from simple chemical dynamics and division control

Experiments on steady-state bacterial cultures have uncovered several quantitative regularities at the system level. These include, first, the exponential growth of cell size with time and the balanced growth of intracellular chemicals between cell birth and division, which are puzzling given the nonlinear and decentralized chemical dynamics in the cell. We model a cell as a set of chemical populations undergoing nonlinear mass action kinetics in a container whose volume is a linear function of the chemical populations. This turns out to be a special class of dynamical systems that generically has attractors in which all populations grow exponentially with time at the same rate. This explains exponential balanced growth of bacterial cells without invoking any regulatory mechanisms and suggests that this could be a robust property of protocells as well. Second, we consider the hypothesis that cells commit themselves to division when a certain internal chemical population reaches a threshold of N molecules. We show that this hypothesis leads to a simple explanation of some of the variability observed across cells in a bacterial culture. In particular, it reproduces the adder property of cell size fluctuations observed recently in E. coli; the observed correlations among interdivision time, birth volume, and added volume in a generation; and the observed scale of the fluctuations (CV ≈ 10-30%) when N is between 10 and 100. Third, upon including a suitable regulatory mechanism that optimizes the growth rate of the cell, the model reproduces the observed bacterial growth laws including the dependence of the growth rate and ribosomal protein fraction on the medium. Thus, the models provide a framework for unifying diverse aspects of bacterial growth physiology under one roof. They also suggest new questions for experimental and theoretical enquiry.