Surface-to-volume scaling and aspect ratio preservation in rod-shaped bacteria
Abstract
Rod-shaped bacterial cells can readily adapt their lengths and widths in response to environmental changes. While many recent studies have focused on the mechanisms underlying bacterial cell size control, it remains largely unknown how the coupling between cell length and width results in robust control of rod-like bacterial shapes. In this study we uncover a conserved surface-to-volume scaling relation in Escherichia coli and other rod-shaped bacteria, resulting from the preservation of cell aspect ratio. To explain the mechanistic origin of aspect-ratio control, we propose a quantitative model for the coupling between bacterial cell elongation and the accumulation of an essential division protein, FtsZ. This model reveals a mechanism for why bacterial aspect ratio is independent of cell size and growth conditions, and predicts cell morphological changes in response to nutrient perturbations, antibiotics, MreB or FtsZ depletion, in quantitative agreement with experimental data.
https://doi.org/10.7554/eLife.47033.001Introduction
Cell morphology is an important adaptive trait that is crucial for bacterial growth, motility, nutrient uptake, and proliferation (Young, 2006). When rod-shaped bacteria grow in media with different nutrient availability, both cell length and width increase with growth rate (Schaechter et al., 1958; Si et al., 2017). At the single-cell level, control of cell volume in many rod-shaped cells is achieved via an adder mechanism, whereby cells elongate by a fixed length per division cycle (Amir, 2014; Campos et al., 2014; Taheri-Araghi et al., 2015; Deforet et al., 2015; Wallden et al., 2016; Banerjee et al., 2017). A recent study has linked the determination of cell size to a condition-dependent regulation of cell surface-to-volume ratio (Harris and Theriot, 2016). However, it remains largely unknown how cell length and width are coupled to regulate rod-like bacterial shapes in diverse growth conditions (Volkmer and Heinemann, 2011; Belgrave and Wolgemuth, 2013; Colavin et al., 2018; Shi et al., 2018).
Results
Here we investigated the relation between cell surface area () and cell volume () for E. coli cells grown under different nutrient conditions, challenged with antibiotics, protein overexpression or depletion, and single gene deletions (Nonejuie et al., 2013; Harris and Theriot, 2016; Si et al., 2017; Vadia et al., 2017; Campos et al., 2018; Gray et al., 2019). Collected surface and volume data span two orders of magnitude and exhibit a single power law in this regime: (Figure 1A). Specifically, during steady-state growth (Si et al., 2017), γ = 6.24 ± 0.04, suggesting an elegant geometric relation: . This surface-to-volume scaling with a constant prefactor, , is a consequence of tight control of cell aspect ratio (length/width) (Figure 1D), whose mechanistic origin has been puzzling for almost half a century (Zaritsky, 1975; Zaritsky, 2015). Specifically, for a sphero-cylindrical bacterium, implies . A constant thus defines a constant aspect ratio η = 4.14 ± 0.17 (Figure 1B-inset), with a coefficient of variation ∼14% (Figure 1B).
The surface-to-volume relation for steady-state growth, , results in a simple expression for cell surface-to-volume ratio: . Using the phenomenological nutrient growth law (Schaechter et al., 1958), where is the population growth rate, a negative correlation emerges between and :
with the cell volume at , and is the relative rate of increase in with (Figure 1C). In Equation (1) underlies an adaptive feedback response of the cell — at low nutrient concentrations, cells increase their surface-to-volume ratio to promote nutrient influx (Si et al., 2017; Harris and Theriot, 2018). Prediction from Equation (1) is in excellent agreement with the best fit to the experimental data. Furthermore, a constant aspect ratio of ≈4 implies and , where is the cell width, suggesting stronger geometric constraints than recently proposed (Harris and Theriot, 2018; Shi et al., 2018). Thus, knowing cell volume as a function of cell cycle parameters (Si et al., 2017) we can directly predict cell width and length under changes in growth media, in agreement with experimental data (Figure 1—figure supplement 1A–B). We further analysed cell shape data for 48 rod-shaped bacteria, one rod-shaped Archaea (H. vulcanii), two long spiral Spirochete, and one coccoid bacteria (Figure 1E). Collected data for all rod-shaped cells follow closely the relationship , while the long Spirochetes deviate from this curve (Figure 1D–E). Coccoid S. aureus also follows the universal scaling relation (with γ = 4.92), but maintains a much lower aspect ratio η = 1.38 ± 0.18 (Quach et al., 2016) when exposed to different antibiotics (Figure 1D–E). Therefore, aspect-ratio preservation likely emerges from a mechanism that is common to diverse rod-shaped and coccoid bacterial species.
To investigate how aspect ratio is regulated at the single cell level we analysed the morphologies of E. coli cells grown in the mother machine (Taheri-Araghi et al., 2015) (Figure 2A,B). For five different growth media, mean volume and surface area of newborn cells also follow the relationship , suggesting that a fixed aspect ratio is maintained on average. In the single-cell data, slight deviation from the 2/3 scaling is a consequence of large fluctuations in newborn cell lengths for a given cell width (Figure 2—figure supplement 1A–B). Importantly, the probability distribution of aspect ratio is independent of the growth media (Figure 2B), implying that cellular aspect ratio is independent of cell size as well as growth rate.
To explain the origin of aspect ratio homeostasis we developed a quantitative model for cell shape dynamics that accounts for the coupling between cell elongation and the accumulation of cell division proteins FtsZ (Figure 2C). Our model is thus only applicable to bacteria that divide using the FtsZ machinery. E. coli and other rod-like bacteria maintain a constant width during their cell cycle while elongating exponentially in length (Donachie et al., 1976; Taheri-Araghi et al., 2015): , with the elongation rate. Cell division is triggered when a constant length is added per division cycle — a mechanism that is captured by a model for threshold accumulation of division initiator proteins, produced at a rate proportional to cell size (Basan et al., 2015; Deforet et al., 2015; Ghusinga et al., 2016). While many molecular candidates have been suggested as initiators of division (Adams and Errington, 2009), a recent study (Si et al., 2019) has identified FtsZ as the key initiator protein that assembles a ring-like structure in the mid-cell region to trigger septation.
Dynamics of division protein accumulation can be described using a two-component model. First, a cytoplasmic component with abundance grows in proportion to cell size (), as ribosome content increases with cell size (Marguerat and Bähler, 2012). Second, a ring-bound component, , is assembled from the cytoplasmic pool at a constant rate. Dynamics of the cytoplasmic and ring-bound FtsZ are given by:
where is the constant production rate of cytoplasmic FtsZ, is the rate of binding of cytoplasmic FtsZ to the Z-ring, and is the rate of disassembly of Z-ring bound FtsZ. At the start of the cell cycle, we have (a constant) and . Cell divides when reaches a threshold amount, , required for the completion of ring assembly. A key ingredient of our model is that scales linearly with the cell circumference, , preserving the density of FtsZ in the ring. This is consistent with experimental findings that the total FtsZ scales with the cell width (Shi et al., 2017). Accumulation of FtsZ proteins, , follows the equation: , where is the production rate of division proteins, with at the start of the division cycle. We assume that , such that all the newly synthesized cytoplasmic proteins are recruited to the Z-ring at a rate much faster than growth rate (Söderström et al., 2018). As a result, cell division occurs when (Figure 2C). Upon division is reset to 0 for the two daughter cells. It is reasonable to assume that all the FtsZ proteins are in filamentous form at cell division, as the concentration of FtsZ in an average E. coli cell is in the range 4–10 μM, much higher than the critical concentration 1 μM (Erickson et al., 2010).
From the model it follows that during one division cycle cells grow by adding a length , which equals the homeostatic length of newborn cells. Furthermore, recent experiments suggest that the amount of FtsZ synthesised per unit cell length, , is constant (Si et al., 2019). This implies,
Aspect ratio homeostasis is thus achieved via a balance between the rates of cell elongation and division protein production, consistent with observations that FtsZ overexpression leads to minicells and FtsZ depletion induces elongated phenotypes (Potluri et al., 2012; Zheng et al., 2016). Indeed single cell E. coli data (Taheri-Araghi et al., 2015) show that is constant on average and independent of growth conditions (Figure 2D). Furthermore, added length correlates with cell width during one cell cycle implying that the cell width is a good predictor for added cell length (Figure 2—figure supplement 1C–D).
To predict cell-shape dynamics under perturbations to growth conditions we simulated our single-cell model (Figure 3, Materials and methods) with an additional equation for cell width that we derived from a recent model proposed by Harris and Theriot (2016): , where is the rate of surface area synthesis relative to volume and is a linearly increasing function of (Figure 3—figure supplement 1A). This model leads to an equation for the control of cell width for a sphero-cylinder shaped bacterium,
such that at steady-state. It then follows from Equation (4) that the added cell length . However, our model for division control is mechanistically different from Harris and Theriot (2016). In the latter, cells accumulate a threshold amount of excess surface area material to trigger septation, which does not lead to aspect ratio preservation. By contrast, we propose that cells divide when they accumulate a threshold amount of division proteins in the Z-ring, proportional to the cell diameter.
We simulated nutrient shift experiments using the coupled equations for cell length, width and division protein production (Materials and methods). When simulated cells are exposed to new nutrient conditions (Figure 3—figure supplement 1B–E), changes in cell width result in a transient increase in aspect ratio () during nutrient downshift, or a transient decrease in during nutrient upshift (Figure 3C). After nutrient shift, aspect ratio reaches its pre-stimulus homeostatic value over multiple generations. Typical timescale for transition to the new steady-state is controlled by the growth rate of the new medium (), such that the cell shape parameters reach a steady state faster in media with higher growth rate. This result is consistent with the experimental observation that newborn aspect ratio reaches equilibrium faster in fast growing media (Taheri-Araghi et al., 2015) (Figure 3—figure supplement 1F). In our model, cell shape changes are controlled by two parameters: the ratio that determines cell aspect ratio, and that controls cell width (Figure 4A). Nutrient upshift or downshift only changes the ratio while keeping the steady-state aspect ratio constant.
We further used our model to predict drastic shape changes leading to deviations from the homeostatic aspect ratio when cells are perturbed by FtsZ knockdown, MreB depletion, and antibiotic treatments that induce non steady state filamentation (Figure 4B). First, FtsZ depletion results in long cells while the width stays approximately constant, (Figure 4—figure supplement 1C), data from Zheng et al. (2016). We modelled FtsZ knockdown by decreasing and simulations quantitatively agree with experimental data. Second, MreB depletion increases the cell width and slightly decreases cell length while keeping growth rate constant (Zheng et al., 2016). We modelled MreB knockdown by decreasing as expected for disruption in cell wall synthesis machinery, while simultaneously increasing (Materials and methods). This increase in is consistent with a prior finding that in MreB mutant cells of various sizes, the total FtsZ scales with the cell width (Shi et al., 2017). Furthermore, cells treated with MreB inhibitor A22 induce envelope stress response system (Rcs) that in turn activates FtsZ overproduction (Carballès et al., 1999; Cho et al., 2014). Third, transient long filamentous cells result from exposure to high dosages of cell-wall targeting antibiotics that prevent cell division, or DNA-targeting antibiotics that induce filamentation via SOS response (Nonejuie et al., 2013). Cell-wall targeting antibiotics inhibit the activity of essential septum forming penicillin binding proteins, preventing cell septation. We modelled this response as an effective reduction in , while slightly decreasing surface synthesis rate (Materials and methods). For DNA targeting antibiotics, FtsZ is directly sequestered during SOS response resulting in delayed ring formation and septation (Chen et al., 2012). Surprisingly all filamentous cells have a similar aspect ratio of 11.0 ± 1.4, represented by a single curve in the - plane (Figure 4B).
Discussion
The conserved surface-to-volume scaling in diverse bacterial species, , is a direct consequence of aspect-ratio homeostasis at the single-cell level. We present a regulatory model (Figure 2C) where aspect-ratio control is the consequence of a constant ratio between the rate of cell elongation () and division protein accumulation (). Deviation from the homeostatic aspect ratio is a consequence of altered , as observed in filamentous cells, FtsZ or MreB depleted cells (Figure 4B). By contrast, drugs that target cell wall biogenesis, for example Fosfomycin, do not alter and maintain cellular aspect ratio (Figure 4—figure supplement 1C).
Our study suggests that cell width is an essential shape parameter for determining cell length in E. coli (Figure 2—figure supplement 1C–D). This is to be contrasted with B. subtilis, where cell width stays approximately constant across different media, while elongating in length (Sharpe et al., 1998). However, FtsZ recruitment in B. subtilis is additionally controlled by effector UgtP, which localises to the division site in a nutrient-dependent manner and prevents Z-ring assembly (Weart et al., 2007). This can be interpreted as a reduction in with increasing , within the framework of our model. As a result, B. subtilis aspect ratio () is predicted to increase with increasing growth rate.
Aspect ratio control may have several adaptive benefits. For instance, increasing cell surface-to-volume ratio under low nutrient conditions can result in an increased nutrient influx to promote cell growth (Figure 1C). Under translation inhibition by ribosome-targeting antibiotics, bacterial cells increase their volume while preserving aspect ratio (Harris and Theriot, 2016; Si et al., 2017). This leads to a reduction in surface-to-volume ratio to counter further antibiotic influx. Furthermore, recent studies have shown that the efficiency of swarming bacteria strongly depends on their aspect ratio (Ilkanaiv et al., 2017; Jeckel et al., 2019). The highest foraging speed has been observed for aspect ratios in the range 4–6 (Ilkanaiv et al., 2017), suggesting that the maintenance of an optimal aspect ratio may have evolutionary benefits for cell swarmers.
Materials and methods
Cell shape analysis
Request a detailed protocolBacterial cell surface area and volume are obtained directly from previous publications where these values were reported (Si et al., 2017; Campos et al., 2018; Gray et al., 2019), or they are calculated assuming a sphero-cylindrical cell geometry using reported values for population-averaged cell length and width (Vadia et al., 2017; Nonejuie et al., 2013; Zheng et al., 2016; Wright et al., 2015; Deforet et al., 2015; Desmarais et al., 2015; Ojkic et al., 2016; Carabetta et al., 2016; Lopez-Garrido et al., 2018). Single cell data are obtained from Suckjoon Jun lab (UCSD) (Taheri-Araghi et al., 2015). For number of cells analyzed per growth condition see Supplementary file 1. Intergeneration autocorrelation function (Figure 2—figure supplement 1D) of average cell width during one cell cycle is calculated using expression in Ojkic et al. (2014). For a spherocylinder of pole-to-pole length and width , the surface area is , and volume is given by . In the case of S. aureus, surface area and volume are computed assuming prolate spheroidal shape using reported population averaged values of cell major axis, , and minor axis (Quach et al., 2016). Surface area of a prolate spheroid is , and volume is .
Cell growth simulations
Request a detailed protocolWe simulated the single-cell model using the coupled equations for the dynamics of cell length , cell width , and division protein production (Figure 2C). In simulations, when reaches the threshold , the mother cell divides into two daughter cells whose lengths are 0.5 ± δ of the mother cell. Parameter is picked from Gaussian distribution (, ).
For nutrient shift simulations we simulated 105 asynchronous cells growing at a rate (Figure 3). In Equation 5, parameter is obtained from the fit to experimental data for vs (Figure 3—figure supplement 1A) (Si et al., 2017), giving , where is in units of h–1, and in h-1 μm-1. At t = 0 h we change corresponding to nutrient upshift ( 1.25, 2 h-1) or nutrient downshift ( 0.75, 0.25 h-1). We calculated population average of length and width (Figure 3B), and population average of aspect ratio of newborn cells (Figure 3C). Aspect ratio of newborn cells are binned in time and the bin average is calculated for a temporal bin size of 10 min. Examples of single cell traces during the nutrient shift are shown in Figure 3—figure supplement 1B–E.
FtsZ depletion experiment (Zheng et al., 2016) was simulated for while was reduced to 40% of its initial value. This is consistent with the reduction of relative mRNA to ∼40% corresponding to addition of 3 ng/ml of aTc to reduce ftsZ expression (Zheng et al., 2016). Our model predictions for the dependence of cell aspect ratio on is shown in Figure 4—figure supplement 1B.
Best fit for MreB depletion experiment (Zheng et al., 2016) was obtained for η ≈ 2.7, by simulating reduction in division protein production rate, , and by varying so that width spans range from 0.9 to 1.8 μm. The best fit for long filamentous cells (resulting from DNA or cell-wall targeting antibiotics) was obtained for η ≈ 11.0 . Filamentation was simulated by decreasing and so that spans the range from 0.9 to 1.4 μm as experimentally observed (Nonejuie et al., 2013).
Open Source Physics (www.compadre.org) Java was used for executing the simulations and Mathematica 11 for data analysis, model fitting, and data presentation.
Data availability
All data generated or analysed during this study are referenced in the manuscript and supporting files.
References
-
Bacterial cell division: assembly, maintenance and disassembly of the Z ringNature Reviews Microbiology 7:642–653.https://doi.org/10.1038/nrmicro2198
-
Cell size regulation in BacteriaPhysical Review Letters 112:208102.https://doi.org/10.1103/PhysRevLett.112.208102
-
Inflating bacterial cells by increased protein synthesisMolecular Systems Biology 11:836.https://doi.org/10.15252/msb.20156178
-
Evaluation of exposure to Brevundimonas diminuta and Pseudomonas aeruginosa during ShoweringJournal of Aerosol Science 114:77–93.https://doi.org/10.1016/j.jaerosci.2017.08.008
-
SulA inhibits assembly of FtsZ by a simple sequestration mechanismBiochemistry 51:3100–3109.https://doi.org/10.1021/bi201669d
-
Cell-Size homeostasis and the incremental rule in a bacterial pathogenBiophysical Journal 109:521–528.https://doi.org/10.1016/j.bpj.2015.07.002
-
High-throughput, highly sensitive analyses of bacterial morphogenesis using ultra performance liquid chromatographyJournal of Biological Chemistry 290:31090–31100.https://doi.org/10.1074/jbc.M115.661660
-
FtsZ in bacterial cytokinesis: cytoskeleton and force generator all in oneMicrobiology and Molecular Biology Reviews 74:504–528.https://doi.org/10.1128/MMBR.00021-10
-
Surface area to volume ratio: a natural variable for bacterial morphogenesisTrends in Microbiology 26:815–832.https://doi.org/10.1016/j.tim.2018.04.008
-
Effect of cell aspect ratio on swarming BacteriaPhysical Review Letters 118:.https://doi.org/10.1103/PhysRevLett.118.158002
-
Coordinating genome expression with cell sizeTrends in Genetics 28:560–565.https://doi.org/10.1016/j.tig.2012.07.003
-
Bistable forespore engulfment in Bacillus subtilis by a zipper mechanism in absence of the cell wallPLOS Computational Biology 10:e1003912.https://doi.org/10.1371/journal.pcbi.1003912
-
Helicobacter pylori in culture: an ultrastructural studyThe Hokkaido Journal of Medical Science 75:.
-
Bacillus subtilis cell cycle as studied by fluorescence microscopy: constancy of cell length at initiation of DNA replication and evidence for active nucleoid partitioningJournal of Bacteriology 180:547–555.
-
Mechanistic origin of Cell-Size control and homeostasis in BacteriaCurrent Biology 29:1760–1770.https://doi.org/10.1016/j.cub.2019.04.062
-
Cell-size control and homeostasis in BacteriaCurrent Biology 25:385–391.https://doi.org/10.1016/j.cub.2014.12.009
-
Shaping and moving a SpiroplasmaJournal of Molecular Microbiology and Biotechnology 7:78–87.https://doi.org/10.1159/000077872
-
The selective value of bacterial shapeMicrobiology and Molecular Biology Reviews 70:660–703.https://doi.org/10.1128/MMBR.00001-06
-
Dimensional rearrangement of Escherichia coli B/r cells during a nutritional shift-downJournal of Theoretical Biology 54:.https://doi.org/10.1099/00221287-139-11-2711
Article and author information
Author details
Funding
Royal Society (URF/R1/180187)
- Shiladitya Banerjee
Royal Society (RGF/EA/181044)
- Shiladitya Banerjee
Engineering and Physical Sciences Research Council (EP/R029822/1)
- Shiladitya Banerjee
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Suckjoon Jun lab (UCSD) for providing single cell shape data for E. coli, and Javier López-Garrido, Guillaume Charras, and Deb Sankar Banerjee for useful comments. SB gratefully acknowledges funding from EPSRC New Investigator Award EP/R029822/1, Royal Society Tata University Research Fellowship URF/R1/180187, and Royal Society grant RGF/EA/181044.
Copyright
© 2019, Ojkic et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 9,665
- views
-
- 673
- downloads
-
- 71
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
Planar cell polarity (PCP) – tissue-scale alignment of the direction of asymmetric localization of proteins at the cell-cell interface – is essential for embryonic development and physiological functions. Abnormalities in PCP can result in developmental imperfections, including neural tube closure defects and misaligned hair follicles. Decoding the mechanisms responsible for PCP establishment and maintenance remains a fundamental open question. While the roles of various molecules – broadly classified into “global” and “local” modules – have been well-studied, their necessity and sufficiency in explaining PCP and connecting their perturbations to experimentally observed patterns have not been examined. Here, we develop a minimal model that captures the proposed features of PCP establishment – a global tissue-level gradient and local asymmetric distribution of protein complexes. The proposed model suggests that while polarity can emerge without a gradient, the gradient not only acts as a global cue but also increases the robustness of PCP against stochastic perturbations. We also recapitulated and quantified the experimentally observed features of swirling patterns and domineering non-autonomy, using only three free model parameters - the rate of protein binding to membrane, the concentration of PCP proteins, and the gradient steepness. We explain how self-stabilizing asymmetric protein localizations in the presence of tissue-level gradient can lead to robust PCP patterns and reveal minimal design principles for a polarized system.
-
- Computational and Systems Biology
- Physics of Living Systems
Synthetic genetic oscillators can serve as internal clocks within engineered cells to program periodic expression. However, cell-to-cell variability introduces a dispersion in the characteristics of these clocks that drives the population to complete desynchronization. Here, we introduce the optorepressilator, an optically controllable genetic clock that combines the repressilator, a three-node synthetic network in E. coli, with an optogenetic module enabling to reset, delay, or advance its phase using optical inputs. We demonstrate that a population of optorepressilators can be synchronized by transient green light exposure or entrained to oscillate indefinitely by a train of short pulses, through a mechanism reminiscent of natural circadian clocks. Furthermore, we investigate the system’s response to detuned external stimuli observing multiple regimes of global synchronization. Integrating experiments and mathematical modeling, we show that the entrainment mechanism is robust and can be understood quantitatively from single cell to population level.