Distinguishing different modes of growth using single-cell data

Mechanistic basis for non-exponential bacterial growth

Bacterial populations typically exhibit exponential growth under resource-rich conditions, yet individual cells often deviate from this pattern. Recent work has shown that the elongation rates of Escherichia coli and Caulobacter crescentus increase throughout the cell cycle (super-exponential growth), while Bacillus subtilis displays a mid-cycle minimum (convex growth), and Mycobacterium tuberculosis grows linearly. Here, we develop a single-cell model linking gene expression, proteome allocation, and mass growth to explain these diverse growth trajectories. By calibrating model parameters with experimental data, we show that DNA-proportional mRNA transcription produces near-exponential growth, whereas deviations from this proportionality yield the observed non-exponential growth patterns. Analysis of gene expression perturbations reveals that ribosome expression primarily controls dry mass growth rate, whereas envelope expression more strongly affects cell elongation rate. Fitting our model to single-cell experimental data reproduces convex, super-exponential, and linear modes of growth, demonstrating how envelope and ribosome expression schedules drive cell-cycle-specific behaviors. These findings provide a mechanistic basis for non-exponential single-cell growth and offer insights into how bacterial cells dynamically regulate elongation rates within each generation.

From noisy cell size control to population growth: When variability can be beneficial

Single-cell experiments revealed substantial variability in generation times, growth rates, but also in birth and division sizes between genetically identical cells. Understanding how these fluctuations determine the fitness of the population, i.e., its growth rate, is necessary in any quantitative theory of evolution. Here, we develop a biologically relevant model which accounts for the stochasticity in single-cell growth rates, birth sizes, and division sizes. We derive expressions for the population growth rate and mean birth size in the population in terms of single-cell fluctuations. Allowing division sizes to fluctuate reveals how the mechanism of cell size control (timer, sizer, and adder) influences population growth. Surprisingly, we find that fluctuations in single-cell growth rates can be beneficial for population growth when slow-growing cells tend to divide at smaller sizes than fast-growing cells. Our framework is not limited to exponentially growing cells like Escherichia coli, and we derive similar expressions for cells with linear and bilinear growth laws, such as Mycobacterium tuberculosis and fission yeast Schizosaccharomyces pombe, respectively.

Accurate phenotype-to-genotype mapping of high-diversity yeast libraries by heat-shock-electroporation (HEEL)

High-throughput DNA transformation techniques are invaluable when generating high-diversity mutant libraries, a cornerstone of successful protein engineering. However, transformation efficiencies have a direct correlation with the probability of introducing multiple DNA molecules into each cell, although reliable library screenings require cells that contain a single unique genotype. Thus, transformation methods that yield a high multiplicity of transformations are unsuitable for high-diversity library screenings. Here, we describe an innovative yeast library transformation method that is both simple and highly efficient. Our dual heat-shock and electroporation approach (HEEL) creates high-quality DNA libraries by increasing the fraction of mono-transformed yeast cells from 20% to over 70% of all transformed cells, thus allowing for near-perfect phenotype-to-genotype associations. HEEL also allows more than 107 yeast cells per reaction to be transformed with a circular plasmid molecule, which corresponds to an almost 100-fold improvement compared with current yeast transformation methods. To further refine our library screening approach, we integrated an automated yeast genotyping workflow with a dual-barcode design that employs both a single nucleotide polymorphism and a high-diversity region. This design allows for robust identification and quantification of unique genotypes within a heterogeneous population using standard Sanger sequencing. Our findings demonstrate that the longstanding trade-off between the size and quality of transformed yeast libraries can be overcome. By employing the HEEL method, large DNA libraries can be transformed into yeast with high-efficiency, while maintaining high library quality, essential for successful mutant screenings. This advancement holds significant promise for the fields of molecular biology and protein engineering.IMPORTANCEWith the recent expansion of artificial intelligence in the field of synthetic biology, there has never been a greater need for high-quality data and reliable measurements of phenotype-to-genotype relationships. However, one major obstacle to creating accurate computer-based models is the current abundance of low-quality phenotypic measurements originating from numerous high-throughput but low-resolution assays. Rather than increasing the quantity of measurements, new studies should aim to generate as accurate measurements as possible. The HEEL methodology presented here aims to address this issue by minimizing the problem of multi-plasmid uptake during high-throughput yeast DNA transformations, which leads to the creation of heterogeneous cellular genotypes. HEEL should enable highly accurate phenotype-to-genotype measurements going forward, which could be used to construct better computer-based models.

Genome concentration limits cell growth and modulates proteome composition in Escherichia coli

Defining the cellular factors that drive growth rate and proteome composition is essential for understanding and manipulating cellular systems. In bacteria, ribosome concentration is known to be a constraining factor of cell growth rate, while gene concentration is usually assumed not to be limiting. Here, using single-molecule tracking, quantitative single-cell microscopy, and modeling, we show that genome dilution in Escherichia coli cells arrested for DNA replication limits total RNA polymerase activity within physiological cell sizes across tested nutrient conditions. This rapid-onset limitation on bulk transcription results in sub-linear scaling of total active ribosomes with cell size and sub-exponential growth. Such downstream effects on bulk translation and cell growth are near-immediately detectable in a nutrient-rich medium, but delayed in nutrient-poor conditions, presumably due to cellular buffering activities. RNA sequencing and tandem-mass-tag mass spectrometry experiments further reveal that genome dilution remodels the relative abundance of mRNAs and proteins with cell size at a global level. Altogether, our findings indicate that chromosome concentration is a limiting factor of transcription and a global modulator of the transcriptome and proteome composition in E. coli. Experiments in Caulobacter crescentus and comparison with eukaryotic cell studies identify broadly conserved DNA concentration-dependent scaling principles of gene expression.

Single-cell imaging of the Mycobacterium tuberculosis cell cycle reveals linear and heterogenous growth

Difficulties in antibiotic treatment of Mycobacterium tuberculosis (Mtb) are partly thought to be due to heterogeneity in growth. Although the ability of bacterial pathogens to regulate growth is crucial to control homeostasis, virulence and drug responses, single-cell growth and cell cycle behaviours of Mtb are poorly characterized. Here we use time-lapse, single-cell imaging of Mtb coupled with mathematical modelling to observe asymmetric growth and heterogeneity in cell size, interdivision time and elongation speed. We find that, contrary to Mycobacterium smegmatis, Mtb initiates cell growth not only from the old pole but also from new poles or both poles. Whereas most organisms grow exponentially at the single-cell level, Mtb has a linear growth mode. Our data show that the growth behaviour of Mtb diverges from that of model bacteria, provide details into how Mtb grows and creates heterogeneity and suggest that growth regulation may also diverge from that in other bacteria.

Determining the rate-limiting processes for cell division in Escherichia coli

A critical cell cycle checkpoint for most bacteria is the onset of constriction when the septal peptidoglycan synthesis starts. According to the current understanding, the arrival of FtsN to midcell triggers this checkpoint in Escherichia coli. Recent structural and in vitro data suggests that recruitment of FtsN to the Z-ring leads to a conformational switch in actin-like FtsA, which links FtsZ protofilaments to the cell membrane and acts as a hub for the late divisome proteins. Here, we investigate this putative pathway using in vivo measurements and stochastic cell cycle modeling at moderately fast growth rates. Quantitatively upregulating protein concentrations and determining the resulting division timings shows that FtsN and FtsA numbers are not rate-limiting for the division in E. coli. However, at higher overexpression levels, they affect divisions: FtsN by accelerating and FtsA by inhibiting them. At the same time, we find that the FtsZ numbers in the cell are one of the rate-limiting factors for cell divisions in E. coli. Altogether, these findings suggest that instead of FtsN, accumulation of FtsZ in the Z-ring is one of the main drivers of the onset of constriction in E. coli at faster growth rates.

Mechanics and Morphology of Proliferating Cell Collectives with Self-Inhibiting Growth

We study the dynamics of proliferating cell collectives whose microscopic constituents' growth is inhibited by macroscopic growth-induced stress. Discrete particle simulations of a growing collective show the emergence of concentric-ring patterns in cell size whose spatiotemporal structure is closely tied to the individual cell's stress response. Motivated by these observations, we derive a multiscale continuum theory whose parameters map directly to the discrete model. Analytical solutions of this theory show the concentric patterns arise from anisotropically accumulated resistance to growth over many cell cycles. This Letter shows how purely mechanical processes can affect the internal patterning and morphology of cell collectives, and provides a concise theoretical framework for connecting the micro- to macroscopic dynamics of proliferating matter.

Differential growth regulates asymmetric size partitioning in Caulobacter crescentus

Cell size regulation has been extensively studied in symmetrically dividing cells, but the mechanisms underlying the control of size asymmetry in asymmetrically dividing bacteria remain elusive. Here, we examine the control of asymmetric division in Caulobacter crescentus, a bacterium that produces daughter cells with distinct fates and morphologies upon division. Through comprehensive analysis of multi-generational growth and shape data, we uncover a tightly regulated cell size partitioning mechanism. We find that errors in division site positioning are promptly corrected early in the division cycle through differential growth. Our analysis reveals a negative feedback between the size of daughter cell compartments and their growth rates, wherein the larger compartment grows slower to achieve a homeostatic size partitioning ratio at division. To explain these observations, we propose a mechanistic model of differential growth, in which equal amounts of growth regulators are partitioned into daughter cell compartments of unequal sizes and maintained over time via size-independent synthesis.

Effects of Molecular Noise on Cell Size Control

Cells employ control strategies to maintain a stable size. Dividing at a target size (the "sizer" strategy) is thought to produce the tightest size distribution. However, this result follows from phenomenological models that ignore the molecular mechanisms required to implement the strategy. Here we investigate a simple mechanistic model for exponentially growing cells whose division is triggered at a molecular abundance threshold. We find that size noise inherits the molecular noise and is consequently minimized not by the sizer but by the "adder" strategy, where a cell divides after adding a target amount to its birth size. We derive a lower bound on size noise that agrees with publicly available data from six microfluidic studies on Escherichia coli bacteria.

Gene expression in growing cells: A biophysical primer

Cell growth and gene expression, two essential elements of all living systems, have long been the focus of biophysical interrogation. Advances in experimental single-cell methods have invigorated theoretical studies into these processes. However, until recently, there was little dialog between the two areas of study. In particular, most theoretical models for gene regulation assumed gene activity to be oblivious to the progression of the cell cycle between birth and division. But, in fact, there are numerous ways in which the periodic character of all cellular observables can modulate gene expression. The molecular factors required for transcription and translation-RNA polymerase, transcription factors, ribosomes-increase in number during the cell cycle, but are also diluted due to the continuous increase in cell volume. The replication of the genome changes the dosage of those same cellular players but also provides competing targets for regulatory binding. Finally, cell division reduces their number again, and so forth. Stochasticity is inherent to all these biological processes, manifested in fluctuations in the synthesis and degradation of new cellular components as well as the random partitioning of molecules at each cell division event. The notion of gene expression as stationary is thus hard to justify. In this review, we survey the emerging paradigm of cell-cycle regulated gene expression, with an emphasis on the global expression patterns rather than gene-specific regulation. We discuss recent experimental reports where cell growth and gene expression were simultaneously measured in individual cells, providing first glimpses into the coupling between the two, and motivating several questions. How do the levels of gene expression products - mRNA and protein - scale with the cell volume and cell-cycle progression? What are the molecular origins of the observed scaling laws, and when do they break down to yield non-canonical behavior? What are the consequences of cell-cycle dependence for the heterogeneity ("noise") in gene expression within a cell population? While the experimental findings, not surprisingly, differ among genes, organisms, and environmental conditions, several theoretical models have emerged that attempt to reconcile these differences and form a unifying framework for understanding gene expression in growing cells.

Super-exponential growth and stochastic size dynamics in rod-like bacteria

Proliferating bacterial cells exhibit stochastic growth and size dynamics, but the regulation of noise in bacterial growth and morphogenesis remains poorly understood. A quantitative understanding of morphogenetic noise control, and how it changes under different growth conditions, would provide better insights into cell-to-cell variability and intergenerational fluctuations in cell physiology. Using multigenerational growth and width data of single Escherichia coli and Caulobacter crescentus cells, we deduce the equations governing growth and size dynamics of rod-like bacterial cells. Interestingly, we find that both E. coli and C. crescentus cells deviate from exponential growth within the cell cycle. In particular, the exponential growth rate increases during the cell cycle irrespective of nutrient or temperature conditions. We propose a mechanistic model that explains the emergence of super-exponential growth from autocatalytic production of ribosomes coupled to the rate of cell elongation and surface area synthesis. Using this new model and statistical inference on large datasets, we construct the Langevin equations governing cell growth and size dynamics of E. coli cells in different nutrient conditions. The single-cell level model predicts how noise in intragenerational and intergenerational processes regulate variability in cell morphology and generation times, revealing quantitative strategies for cellular resource allocation and morphogenetic noise control in different growth conditions.

Using conditional independence tests to elucidate causal links in cell cycle regulation in Escherichia coli

How cells regulate their cell cycles is a central question for cell biology. Models of cell size homeostasis have been proposed for bacteria, archaea, yeast, plant, and mammalian cells. New experiments bring forth high volumes of data suitable for testing existing models of cell size regulation and proposing new mechanisms. In this paper, we use conditional independence tests in conjunction with data of cell size at key cell cycle events (birth, initiation of DNA replication, and constriction) in the model bacterium Escherichia coli to select between the competing cell cycle models. We find that in all growth conditions that we study, the division event is controlled by the onset of constriction at midcell. In slow growth, we corroborate a model where replication-related processes control the onset of constriction at midcell. In faster growth, we find that the onset of constriction is affected by additional cues beyond DNA replication. Finally, we also find evidence for the presence of additional cues triggering initiations of DNA replication apart from the conventional notion where the mother cells solely determine the initiation event in the daughter cells via an adder per origin model. The use of conditional independence tests is a different approach in the context of understanding cell cycle regulation and it can be used in future studies to further explore the causal links between cell events.

Correlation between speed and turning naturally arises for sparsely sampled cell movements

Mechanisms regulating cell movement are not fully understood. One feature of cell movement that determines how far cells displace from an initial position is persistence, the ability to perform movements in a direction similar to the previous movement direction. Persistence is thus determined by turning angles (TA) between two sequential displacements and can be characterized by an average TA or persistence time. Recent studies documenting T cell movement in zebrafish found that a cell's average speed and average TA are negatively correlated, suggesting a fundamental cell-intrinsic program whereby cells with a lower TA (and larger persistence time) are intrinsically faster (or faster cells turn less). In this paper we confirm the existence of the correlation between turning and speed for six different datasets on 3D movement of CD8 T cells in murine lymph nodes or liver. Interestingly, the negative correlation between TA and speed was observed in experiments in which liver-localized CD8 T cells rapidly displace due to floating with the blood flow, suggesting that other mechanisms besides cell-intrinsic program may be at play. By simulating correlated random walks using two different frameworks (one based on the von Mises-Fisher (vMF) distribution and another based on the Ornstein-Uhlenbeck (OU) process) we show that the negative correlation between speed and turning naturally arises when cell trajectories are sub-sampled, i.e. when the frequency of sampling is lower than frequency at which cells typically make movements. This effect is strongest when the sampling frequency is of the order of magnitude of the inverse of persistence time of cells and when cells vary in persistence time. The effect arises in part due to the sensitivity of estimated cell speeds to the frequency of imaging whereby less frequent imaging results in slower speeds. Interestingly, by using estimated persistence times for cells in two of our datasets and simulating cell movements using the OU process, we could partially reproduce the experimentally observed correlation between TA and speed without a cell-intrinsic program linking the two processes. Our results thus suggest that sub-sampling may contribute to (and perhaps fully explains) the observed correlation between speed and turning at least for some cell trajectory data and emphasize the role of sampling frequency in the inference of critical cellular parameters of cell motility such as speeds.

Analytical cell size distribution: lineage-population bias and parameter inference

We derive analytical steady-state cell size distributions for size-controlled cells in single-lineage experiments, such as the mother machine, which are fundamentally different from batch cultures where populations of cells grow freely. For exponential single-cell growth, characterizing most bacteria, the lineage-population bias is obtained explicitly. In addition, if volume is evenly split between the daughter cells at division, we show that cells are on average smaller in populations than in lineages. For more general power-law growth rates and deterministic volume partitioning, both symmetric and asymmetric, we derive the exact lineage distribution. This solution is in good agreement with Escherichia coli mother machine data and can be used to infer cell cycle parameters such as the strength of the size control and the asymmetry of the division. When introducing stochastic volume partitioning, we derive the large-size and small-size tails of the lineage distribution and show that the lineage-population bias only depends on the single-cell growth rate. These asymptotic behaviours are extended to the adder model of cell size control. When considering noisy single-cell growth rate, we derive the large-size lineage and population distributions. Finally, we show that introducing noise, either on the volume partitioning or on the single-cell growth rate, can cancel the lineage-population bias.