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

Here we investigated the relation between cell surface area ($S$) and cell volume ($V$) 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: $S=\gamma V^{2/3}$ (Figure 1A). Specifically, during steady-state growth (Si et al., 2017), γ = 6.24 ± 0.04, suggesting an elegant geometric relation: $S\approx 2\pi V^{2/3}$. This surface-to-volume scaling with a constant prefactor, $\gamma$, is a consequence of tight control of cell aspect ratio $\eta$ (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, $S=\gamma V^{2/3}$ implies $\gamma =\eta \pi {\left(\frac{\eta \pi }{4}-\frac{\pi }{12}\right)}^{-2/3}$. A constant $\gamma$ thus defines a constant aspect ratio η = 4.14 ± 0.17 (Figure 1B-inset), with a coefficient of variation ∼14% (Figure 1B).

Figure 1 with 1 supplement see all

Download asset Open asset

The surface-to-volume relation for steady-state growth, $S\approx 2\pi V^{2/3}$, results in a simple expression for cell surface-to-volume ratio: $S/V\approx 2\pi V^{-1/3}$. Using the phenomenological nutrient growth law $V=V_0{e}^{\alpha \kappa }$ (Schaechter et al., 1958), where $\kappa$ is the population growth rate, a negative correlation emerges between $S/V$ and $\kappa$:

with $V_0$ the cell volume at $\kappa =0$, and $\alpha$ is the relative rate of increase in $V$ with $\kappa$ (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 $V\approx \sqrt{8}{w}^3$ and $S\approx 4\pi w^2$, where $w$ 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 $S\approx 2\pi V^{2/3}$, while the long Spirochetes deviate from this curve (Figure 1D–E). Coccoid S. aureus also follows the universal scaling relation $S=\gamma V^{2/3}$ (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 $S=2\pi V^{2/3}$, 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.

Figure 2 with 1 supplement see all

Download asset Open asset

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 $L$ (Donachie et al., 1976; Taheri-Araghi et al., 2015): $\mathrm{d}L/\mathrm{d}t=kL$, with $k$ 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 $P_{\mathrm{c}}$ grows in proportion to cell size ($\propto L$), as ribosome content increases with cell size (Marguerat and Bähler, 2012). Second, a ring-bound component, $P_{\mathrm{r}}$, is assembled from the cytoplasmic pool at a constant rate. Dynamics of the cytoplasmic and ring-bound FtsZ are given by:

where $k_P$ is the constant production rate of cytoplasmic FtsZ, $k_{\mathrm{b}}$ is the rate of binding of cytoplasmic FtsZ to the Z-ring, and $k_{\mathrm{d}}$ is the rate of disassembly of Z-ring bound FtsZ. At the start of the cell cycle, we have $P_{\mathrm{c}}=P^{*}$ (a constant) and $P_r=0$. Cell divides when $P_{\mathrm{r}}$ reaches a threshold amount, $P_0$, required for the completion of ring assembly. A key ingredient of our model is that $P_0$ scales linearly with the cell circumference, $P_0=\rho \pi w$, preserving the density $\rho$ 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, $P=P_{\mathrm{c}}+P_{\mathrm{r}}-P^{*}$, follows the equation: $\mathrm{d}P/\mathrm{d}t=k_{P}L$, where $k_P$ is the production rate of division proteins, with $P=0$ at the start of the division cycle. We assume that $k_{\mathrm{b}}\gg k_{\mathrm{d}}$, 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 $P=P_0$ (Figure 2C). Upon division $P$ 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 $\mathrm{\Delta }L=P_{0}k/k_P$, which equals the homeostatic length of newborn cells. Furthermore, recent experiments suggest that the amount of FtsZ synthesised per unit cell length, $\mathrm{d}P/\mathrm{d}L$, 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 $\mathrm{\Delta }L/w$ 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): $\mathrm{d}S/\mathrm{d}t=\beta V$, where $\beta$ is the rate of surface area synthesis relative to volume and is a linearly increasing function of $k$ (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 $w=4k/\beta$ at steady-state. It then follows from Equation (4) that the added cell length $\mathrm{\Delta }L\propto k^2/\beta k_{\mathrm{P}}$. 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.

Figure 3 with 1 supplement see all

Download asset Open asset

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 ($\eta =L/w$) during nutrient downshift, or a transient decrease in $\eta$ 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 ($\propto k^{-1}$), 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 $k/k_P$ that determines cell aspect ratio, and $k/\beta$ that controls cell width (Figure 4A). Nutrient upshift or downshift only changes the ratio $k/\beta$ while keeping the steady-state aspect ratio $(\propto k/k_P)$ constant.

Figure 4 with 1 supplement see all

Download asset Open asset

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, $S\propto V^{0.95}$ (Figure 4—figure supplement 1C), data from Zheng et al. (2016). We modelled FtsZ knockdown by decreasing $k_P$ 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 $\beta$ as expected for disruption in cell wall synthesis machinery, while simultaneously increasing $k_P$ (Materials and methods). This increase in $k_P$ 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 $k_P$, while slightly decreasing surface synthesis rate $\beta$ (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 $S$-$V$ plane (Figure 4B).

The conserved surface-to-volume scaling in diverse bacterial species, $S\sim V^{2/3}$, 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 ($k$) and division protein accumulation ($k_P$). Deviation from the homeostatic aspect ratio is a consequence of altered $k/k_P$, 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 $k/k_P$ 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 $k_P$ with increasing $k$, within the framework of our model. As a result, B. subtilis aspect ratio ($\propto k/k_P$) 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.

All data generated or analysed during this study are referenced in the manuscript and supporting files.

1. 1. Adams DW
   2. Errington J

   (2009) Bacterial cell division: assembly, maintenance and disassembly of the Z ring

   Nature Reviews Microbiology 7:642–653.

   https://doi.org/10.1038/nrmicro2198

   • PubMed
   • Google Scholar
2. 1. Amir A

   (2014) Cell size regulation in Bacteria

   Physical Review Letters 112:208102.

   https://doi.org/10.1103/PhysRevLett.112.208102

   • Google Scholar
3. 1. Banerjee S
   2. Lo K
   3. Daddysman MK
   4. Selewa A
   5. Kuntz T
   6. Dinner AR
   7. Scherer NF

   (2017) Biphasic growth dynamics control cell division in Caulobacter crescentus

   Nature Microbiology 2:.

   https://doi.org/10.1038/nmicrobiol.2017.116

   • PubMed
   • Google Scholar
4. 1. Basan M
   2. Zhu M
   3. Dai X
   4. Warren M
   5. Sévin D
   6. Wang YP
   7. Hwa T

   (2015) Inflating bacterial cells by increased protein synthesis

   Molecular Systems Biology 11:836.

   https://doi.org/10.15252/msb.20156178

   • PubMed
   • Google Scholar
5. 1. Belgrave AM
   2. Wolgemuth CW

   (2013) Elasticity and biochemistry of growth relate replication rate to cell length and cross-link density in rod-shaped Bacteria

   Biophysical Journal 104:2607–2611.

   https://doi.org/10.1016/j.bpj.2013.04.028

   • PubMed
   • Google Scholar
6. 1. Campos M
   2. Surovtsev IV
   3. Kato S
   4. Paintdakhi A
   5. Beltran B
   6. Ebmeier SE
   7. Jacobs-Wagner C

   (2014) A constant size extension drives bacterial cell size homeostasis

   Cell 159:1433–1446.

   https://doi.org/10.1016/j.cell.2014.11.022

   • PubMed
   • Google Scholar
7. 1. Campos M
   2. Govers SK
   3. Irnov I
   4. Dobihal GS
   5. Cornet F
   6. Jacobs-Wagner C

   (2018) Genomewide phenotypic analysis of growth, cell morphogenesis, and cell cycle events in Escherichia coli

   Molecular Systems Biology 14:e7573.

   https://doi.org/10.15252/msb.20177573

   • PubMed
   • Google Scholar
8. 1. Carabetta VJ
   2. Greco TM
   3. Tanner AW
   4. Cristea IM
   5. Dubnau D

   (2016) Temporal regulation of the Bacillus subtilis Acetylome and Evidence for a Role of MreB Acetylation in Cell Wall Growth

   mSystems 1:e00005.

   https://doi.org/10.1128/mSystems.00005-16

   • PubMed
   • Google Scholar
9. 1. Carballès F
   2. Bertrand C
   3. Bouché JP
   4. Cam K

   (1999) Regulation of osmC Gene Expression by the Two-Component System rcsB-rcsC in Escherichia coli

   Molecular Microbiology 34:.

   https://doi.org/10.1128/JB.183.20.5870-5876.2001

   • Google Scholar
10. 1. Chattopadhyay S
    2. Perkins SD
    3. Shaw M
    4. Nichols TL

    (2017) Evaluation of exposure to Brevundimonas diminuta and Pseudomonas aeruginosa during Showering

    Journal of Aerosol Science 114:77–93.

    https://doi.org/10.1016/j.jaerosci.2017.08.008

    • PubMed
    • Google Scholar
11. 1. Chen Y
    2. Milam SL
    3. Erickson HP

    (2012) SulA inhibits assembly of FtsZ by a simple sequestration mechanism

    Biochemistry 51:3100–3109.

    https://doi.org/10.1021/bi201669d

    • PubMed
    • Google Scholar
12. 1. Cho S-H
    2. Szewczyk J
    3. Pesavento C
    4. Zietek M
    5. Banzhaf M
    6. Roszczenko P
    7. Asmar A
    8. Laloux G
    9. Hov A-K
    10. Leverrier P
    11. Van der Henst C
    12. Vertommen D
    13. Typas A
    14. Collet J-F

    (2014) Detecting envelope stress by monitoring β-Barrel assembly

    Cell 159:1652–1664.

    https://doi.org/10.1016/j.cell.2014.11.045

    • Google Scholar
13. 1. Colavin A
    2. Shi H
    3. Huang KC

    (2018) RodZ modulates geometric localization of the bacterial actin MreB to regulate cell shape

    Nature Communications 9:.

    https://doi.org/10.1038/s41467-018-03633-x

    • PubMed
    • Google Scholar
14. 1. Deforet M
    2. van Ditmarsch D
    3. Xavier JB

    (2015) Cell-Size homeostasis and the incremental rule in a bacterial pathogen

    Biophysical Journal 109:521–528.

    https://doi.org/10.1016/j.bpj.2015.07.002

    • PubMed
    • Google Scholar
15. 1. Desmarais SM
    2. Tropini C
    3. Miguel A
    4. Cava F
    5. Monds RD
    6. de Pedro MA
    7. Huang KC

    (2015) High-throughput, highly sensitive analyses of bacterial morphogenesis using ultra performance liquid chromatography

    Journal of Biological Chemistry 290:31090–31100.

    https://doi.org/10.1074/jbc.M115.661660

    • PubMed
    • Google Scholar
16. 1. Donachie WD
    2. Begg KJ
    3. Vicente M

    (1976) Cell length, cell growth and cell division

    Nature 264:328–333.

    https://doi.org/10.1038/264328a0

    • PubMed
    • Google Scholar
17. 1. Erickson HP
    2. Anderson DE
    3. Osawa M

    (2010) FtsZ in bacterial cytokinesis: cytoskeleton and force generator all in one

    Microbiology and Molecular Biology Reviews 74:504–528.

    https://doi.org/10.1128/MMBR.00021-10

    • PubMed
    • Google Scholar
18. 1. Ghusinga KR
    2. Vargas-Garcia CA
    3. Singh A

    (2016) A mechanistic stochastic framework for regulating bacterial cell division

    Scientific Reports 6:.

    https://doi.org/10.1038/srep30229

    • PubMed
    • Google Scholar
19. 1. Gray WT
    2. Govers SK
    3. Xiang Y
    4. Parry BR
    5. Campos M
    6. Kim S
    7. Jacobs-Wagner C

    (2019) Nucleoid size scaling and intracellular organization of translation across Bacteria

    Cell 177:1632–1648.

    https://doi.org/10.1016/j.cell.2019.05.017

    • PubMed
    • Google Scholar
20. 1. Harris LK
    2. Theriot JA

    (2016) Relative rates of surface and volume synthesis set bacterial cell size

    Cell 165:1479–1492.

    https://doi.org/10.1016/j.cell.2016.05.045

    • PubMed
    • Google Scholar
21. 1. Harris LK
    2. Theriot JA

    (2018) Surface area to volume ratio: a natural variable for bacterial morphogenesis

    Trends in Microbiology 26:815–832.

    https://doi.org/10.1016/j.tim.2018.04.008

    • PubMed
    • Google Scholar
22. 1. Ilkanaiv B
    2. Kearns DB
    3. Ariel G
    4. Be'er A

    (2017) Effect of cell aspect ratio on swarming Bacteria

    Physical Review Letters 118:.

    https://doi.org/10.1103/PhysRevLett.118.158002

    • PubMed
    • Google Scholar
23. 1. Jeckel H
    2. Jelli E
    3. Hartmann R
    4. Singh PK
    5. Mok R
    6. Totz JF
    7. Vidakovic L
    8. Eckhardt B
    9. Dunkel J
    10. Drescher K

    (2019) Learning the space-time phase diagram of bacterial swarm expansion

    PNAS 116:1489–1494.

    https://doi.org/10.1073/pnas.1811722116

    • PubMed
    • Google Scholar
24. 1. Lopez-Garrido J
    2. Ojkic N
    3. Khanna K
    4. Wagner FR
    5. Villa E
    6. Endres RG
    7. Pogliano K

    (2018) Chromosome translocation inflates Bacillus forespores and impacts cellular morphology

    Cell 172:758–770.

    https://doi.org/10.1016/j.cell.2018.01.027

    • PubMed
    • Google Scholar
25. 1. Marguerat S
    2. Bähler J

    (2012) Coordinating genome expression with cell size

    Trends in Genetics 28:560–565.

    https://doi.org/10.1016/j.tig.2012.07.003

    • PubMed
    • Google Scholar
26. 1. Nonejuie P
    2. Burkart M
    3. Pogliano K
    4. Pogliano J

    (2013) Bacterial cytological profiling rapidly identifies the cellular pathways targeted by antibacterial molecules

    PNAS 110:16169–16174.

    https://doi.org/10.1073/pnas.1311066110

    • PubMed
    • Google Scholar
27. 1. Ojkic N
    2. López-Garrido J
    3. Pogliano K
    4. Endres RG

    (2014) Bistable forespore engulfment in Bacillus subtilis by a zipper mechanism in absence of the cell wall

    PLOS Computational Biology 10:e1003912.

    https://doi.org/10.1371/journal.pcbi.1003912

    • PubMed
    • Google Scholar
28. 1. Ojkic N
    2. López-Garrido J
    3. Pogliano K
    4. Endres RG

    (2016) Cell-wall remodeling drives engulfment during Bacillus subtilis sporulation

    eLife 5:e18657.

    https://doi.org/10.7554/eLife.18657

    • PubMed
    • Google Scholar
29. 1. Pelling AE
    2. Li Y
    3. Shi W
    4. Gimzewski JK

    (2005) Nanoscale visualization and characterization of Myxococcus xanthus cells with atomic force microscopy

    PNAS 102:6484–6489.

    https://doi.org/10.1073/pnas.0501207102

    • PubMed
    • Google Scholar
30. 1. Potluri LP
    2. de Pedro MA
    3. Young KD

    (2012) Escherichia coli low-molecular-weight penicillin-binding proteins help orient septal FtsZ, and their absence leads to asymmetric cell division and branching

    Molecular Microbiology 84:203–224.

    https://doi.org/10.1111/j.1365-2958.2012.08023.x

    • PubMed
    • Google Scholar
31. 1. Quach DT
    2. Sakoulas G
    3. Nizet V
    4. Pogliano J
    5. Pogliano K

    (2016) Bacterial cytological profiling (BCP) as a rapid and accurate antimicrobial susceptibility testing method for Staphylococcus aureus

    EBioMedicine 4:95–103.

    https://doi.org/10.1016/j.ebiom.2016.01.020

    • PubMed
    • Google Scholar
32. 1. Sato F

    (2000)

    Helicobacter pylori in culture: an ultrastructural study

    The Hokkaido Journal of Medical Science 75:.

    • Google Scholar
33. 1. Schaechter M
    2. Maaløe O
    3. Kjeldgaard NO

    (1958) Dependency on medium and temperature of cell size and chemical composition during balanced growth of Salmonella typhimurium

    Microbiology 19:.

    https://doi.org/10.1099/00221287-19-3-592

    • Google Scholar
34. 1. Sharpe ME
    2. Hauser PM
    3. Sharpe RG
    4. Errington J

    (1998)

    Bacillus subtilis cell cycle as studied by fluorescence microscopy: constancy of cell length at initiation of DNA replication and evidence for active nucleoid partitioning

    Journal of Bacteriology 180:547–555.

    • PubMed
    • Google Scholar
35. 1. Shi H
    2. Colavin A
    3. Bigos M
    4. Tropini C
    5. Monds RD
    6. Huang KC

    (2017) Deep phenotypic mapping of bacterial cytoskeletal mutants reveals physiological robustness to cell size

    Current Biology 27:3419–3429.

    https://doi.org/10.1016/j.cub.2017.09.065

    • PubMed
    • Google Scholar
36. 1. Shi H
    2. Bratton BP
    3. Gitai Z
    4. Huang KC

    (2018) How to build a bacterial cell: mreb as the foreman of E. coli Construction

    Cell 172:1294–1305.

    https://doi.org/10.1016/j.cell.2018.02.050

    • PubMed
    • Google Scholar
37. 1. Si F
    2. Li D
    3. Cox SE
    4. Sauls JT
    5. Azizi O
    6. Sou C
    7. Schwartz AB
    8. Erickstad MJ
    9. Jun Y
    10. Li X
    11. Jun S

    (2017) Invariance of initiation mass and predictability of cell size in Escherichia coli

    Current Biology 27:1278–1287.

    https://doi.org/10.1016/j.cub.2017.03.022

    • PubMed
    • Google Scholar
38. 1. Si F
    2. Le Treut G
    3. Sauls JT
    4. Vadia S
    5. Levin PA
    6. Jun S

    (2019) Mechanistic origin of Cell-Size control and homeostasis in Bacteria

    Current Biology 29:1760–1770.

    https://doi.org/10.1016/j.cub.2019.04.062

    • PubMed
    • Google Scholar
39. 1. Söderström B
    2. Badrutdinov A
    3. Chan H
    4. Skoglund U

    (2018) Cell shape-independent FtsZ dynamics in synthetically remodeled bacterial cells

    Nature Communications 9:.

    https://doi.org/10.1038/s41467-018-06887-7

    • PubMed
    • Google Scholar
40. 1. Taheri-Araghi S
    2. Bradde S
    3. Sauls JT
    4. Hill NS
    5. Levin PA
    6. Paulsson J
    7. Vergassola M
    8. Jun S

    (2015) Cell-size control and homeostasis in Bacteria

    Current Biology 25:385–391.

    https://doi.org/10.1016/j.cub.2014.12.009

    • PubMed
    • Google Scholar
41. 1. Trachtenberg S

    (2004) Shaping and moving a Spiroplasma

    Journal of Molecular Microbiology and Biotechnology 7:78–87.

    https://doi.org/10.1159/000077872

    • PubMed
    • Google Scholar
42. 1. Vadia S
    2. Tse JL
    3. Lucena R
    4. Yang Z
    5. Kellogg DR
    6. Wang JD
    7. Levin PA

    (2017) Fatty acid availability sets cell envelope capacity and dictates microbial cell size

    Current Biology 27:1757–1767.

    https://doi.org/10.1016/j.cub.2017.05.076

    • PubMed
    • Google Scholar
43. 1. Volkmer B
    2. Heinemann M

    (2011) Condition-dependent cell volume and concentration of Escherichia coli to facilitate data conversion for systems biology modeling

    PLOS ONE 6:e23126.

    https://doi.org/10.1371/journal.pone.0023126

    • PubMed
    • Google Scholar
44. 1. Wallden M
    2. Fange D
    3. Lundius EG
    4. Baltekin Ö
    5. Elf J

    (2016) The synchronization of replication and division cycles in individual E. coli Cells

    Cell 166:729–739.

    https://doi.org/10.1016/j.cell.2016.06.052

    • PubMed
    • Google Scholar
45. 1. Weart RB
    2. Lee AH
    3. Chien AC
    4. Haeusser DP
    5. Hill NS
    6. Levin PA

    (2007) A metabolic sensor governing cell size in Bacteria

    Cell 130:335–347.

    https://doi.org/10.1016/j.cell.2007.05.043

    • PubMed
    • Google Scholar
46. 1. Wright CS
    2. Banerjee S
    3. Iyer-Biswas S
    4. Crosson S
    5. Dinner AR
    6. Scherer NF

    (2015) Intergenerational continuity of cell shape dynamics in Caulobacter crescentus

    Scientific Reports 5:.

    https://doi.org/10.1038/srep09155

    • Google Scholar
47. 1. Young KD

    (2006) The selective value of bacterial shape

    Microbiology and Molecular Biology Reviews 70:660–703.

    https://doi.org/10.1128/MMBR.00001-06

    • PubMed
    • Google Scholar
48. 1. Zaritsky A

    (1975) Dimensional rearrangement of Escherichia coli B/r cells during a nutritional shift-down

    Journal of Theoretical Biology 54:.

    https://doi.org/10.1099/00221287-139-11-2711

    • Google Scholar
49. 1. Zaritsky A

    (2015) Cell-shape homeostasis in Escherichia coli is driven by growth, division, and nucleoid complexity

    Biophysical Journal 109:178–181.

    https://doi.org/10.1016/j.bpj.2015.06.026

    • PubMed
    • Google Scholar
50. 1. Zheng H
    2. PY H
    3. Jiang M
    4. Tang B
    5. Liu W
    6. Li D
    7. Yu X
    8. Kleckner NE
    9. Amir A
    10. Liu C

    (2016) Interrogating the Escherichia coli cell cycle by cell dimension perturbations

    PNAS 113:15000–15005.

    https://doi.org/10.1073/pnas.1617932114

    • Google Scholar

[Editors’ note: this article was originally rejected after discussions between the reviewers, but the authors were invited to resubmit after an appeal against the decision.]

Thank you for submitting your work entitled "Universal surface-to-volume scaling and aspect ratio homeostasis in rod-shaped bacteria" for consideration by eLife. Your article has been reviewed by two peer reviewers, and the evaluation has been overseen by a Reviewing Editor and a Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Charles W Wolgemuth (Reviewer #1).

The reviewers had no doubt that the subject matter of your paper is both interesting and timely, but both have raised questions about its significance, the nature of the underlying assumptions, and the level of "universality" that can be claimed. Our decision on your manuscript has been reached after consultation between the reviewers. Based on these discussions and the individual reviews below, we regret to inform you that your work will not be considered further for publication in eLife.

Reviewer #1:

This manuscript compiles data on the length, width, and growth rate of E. coli under a number of experimental perturbations, such as changes in growth medium, incubation in antibiotics, inhibited protein synthesis, etc. and shows that the surface area to volume ratio is strongly conserved. This result is consistent with recent work from Julie Theriot's group (which is cited in Harris and Theriot, 2016; 2018). In this manuscript, the author's also add data from other rod-shaped bacteria that show similar behavior (Figure 1G). The authors use this result to develop a model for aspect ratio regulation that is based on exponential growth of the bacterial length at fixed width, FtsZ production at a rate proportional to volume growth rate (which by assuming constant width also assumes that the FtsZ production is proportional to the length growth rate), and division that occurs when FtsZ production reaches a critical value that is proportional to the width. This model predicts a constant aspect ratio and the authors then go on to predict the dynamics under impulse type perturbations.

I have two main concerns with the manuscript.

First, the only novelty of the model is the assumption that there is a critical amount of FtsZ required to divide the cell and that this depends on the width. I think that this is a reasonable assumption, but I also feel that the overall results are fairly obvious. That is, it is not clear that the model provides a significant advance in our understanding. That said, there is also a small problem with the model, in that we would expect the binding rate of FtsZ to depend on the surface to volume ratio (which turns out not to matter, because the authors end up making assumptions that the rate the ring is built is equal to the protein production rate). Note also that there is a typo in the equation fordPrdt, which has k_d multiplying both rates.

Second, and more important, is that while the results match well with the data, there are a number of aspects of the presentation that are misleading. The title claims that the results/model presented here are universal. In Figure 1G, the authors select 7 bacteria to claim that the scaling of SA = 2 π V^(2/3) is ubiquitous among bacteria. As noted, this also suggests that an aspect ratio of ~4 is the rule for rod-shaped bacteria. This is not true. As an example, myxococcus xanthus has an aspect ratio around 7-8 and spirochetes have aspect ratios of ~30! Even more importantly, single species don't always maintain the same aspect ratio. In B. subtilis, the aspect ratio can vary between at least 3.8 – 8 (see Ilkanaiv et al., 2017). Therefore, this model may be applicable to E. coli (and possibly some other bacteria), but it is not universal.

Reviewer #2:

In this study the authors set out to study the size and shape of a wide range of 'rod-shaped' cells by collecting image data from at least seven different species and thousands to total conditions (genotype x nutrients x antibiotics). Across all these conditions, the authors find a simple scaling law to the surface area/volume ratio, namely a scaling that preserves the aspect ratio of the cell at roughly 4:1. Given this observation, they build a simple, mechanistically inspired, quantitative model for the growth of the cell. Using this model, they are able to tune three parameters (k, k_p, β) to match a collection of genetic knockdown and antibiotic treatment experiments.

Understanding how cell size and shape homeostasis is maintained throughout the bacterial kingdom is a very interesting and important problem and these authors should be commended for pushing the community to consider that these mechanisms may be conserved across a wide phylogenetic range. However, given the extensive body of literature already available about cell size/shape homeostasis, and, in particular the review mentioned by the authors by Harris and Theriot, the scientific bar for productive engagement on this topic is already quite high. Much of the intellectual driving force for this work seems to follow directly from the hypothesis from Harris and Theriot that "While many studies have treated volume as the actively controlled parameter in this scenario, our recent work suggests that it is likely the other way around, and that SA/V is the actively regulated variable, with size following along as necessary [13]." The current work seeks to extend or provide alternatives for the mechanistic models presented in Harris and Theriot as well as integrate additional data in other species. However, given that the idea of SA/V scaling conservation is not new, appealing to a broad audience such as that of eLife would require experimental validation of their mechanistic model.

In addition to the core concerns around novelty of the central hypothesis and validity of the mechanistic model, there are a few issues the authors might choose to consider:

Major points:

1) The authors should clearly explain how their mechanistic model contrasts with the cell wall-focused model proposed by Harris and Theriot and should strive to propose experiments with predicted outcomes that would differentiate a peptidoglycan centric model from an FtsZ centric model. If the data already exist to rule out one of them, this should be clearly presented. As one such example, the authors show that tuning one parameter (k_p) is consistent with the experimental notion of knocking down the production of FtsZ. However, they fail to show if there is quantitative agreement between the production rate of FtsZ and the amount they expect to need to change k_p(40%).

2) The use of 'universal' in the paper's title significantly oversells of the breadth of species included in the observations and a power law describing data which span roughly one order of magnitude. While the authors do include a large collection of data, the collection is far from comprehensive for all size/shape data available and the authors do not clearly indicate why they limited themselves to the data they did. A quick literature search reveals anecdotal evidence of bacterial sizes that are much smaller than a micron such as Brevundimonas (PDA J Pharm Sci Technol. 2002 Mar-Apr;56(2):99-108.) to nearly a millimeter in length Epulopiscium (J. Protozoal., 35(4), 1988, pp. 565-569). Granted, these publications may not have the same type of data necessary to integrate it directly into their model, but for a discussion of the 'universal scaling', the authors should push themselves to cover as large of a length-scale as possible. When choosing a set of species for inclusion in this study, it seems like the microbiology community may have already picked an aspect ratio of about 4:1 in its definition of rod-shaped bacteria. For example, cells that have a much shorter aspect ratio are given the term ovoid or lancet (Streptococcus pneumoniae) or spherical (Staphylococcus aureus included here) and ones that are much longer are called filamentous (Sphaerotilus natans). Confusingly, these authors do not include species that have been traditionally classified as rod-shaped cells with a longer aspect ratio such as (Helicobacter, Spiroplasma, Spirochetes, Myxobacter).

3) I'm not entirely convinced that the universal scaling applies within the single cell data (Figure 1D). By plotting the single cell data from a variety of experiments, the range of the data seems to put a larger priority on the averages. However, within each condition there seems to be clear deviations from the 'single aspect ratio', consistent with the author's single cell growth model that cells grow without changing their diameter before dividing. This should result in a roughly factor of two change in aspect ratio from birth to division. I think this is what the authors refer to in the fourth paragraph of the Introduction but should discuss more fully.

4) I do not understand Figure 2B at all. In particular, the binning of the data that I have been able to find in Taheri-Araghi et al., 2015, is binned by the size of cells at birth, not the individual cell growth rate. Further, the authors do not describe how they go from the data in Taheri-Araghi et al., 2015, to the data in Figure 2B, but it could be that they obtained the raw data from the authors and performed a new type of analysis. If so, a description of this process should be included.

5) I am unclear about why the MreB and FtsZ knockdown data from Si et al. is included in the bulk Figure 1A data but the MreB and FtsZ knockdown data from Zheng et al. is treated as a completely separate experiment. If the approach that these two studies used was different, it may be helpful to explain why some data are included in one place and others are not.

[Editors’ note: the author responses to the first round of peer review follow.]

The mechanistic or molecular origin of bacterial aspect ratio control has remained an unsolved problem for more than four decades (see e.g. Zaritsky,2015; Zaritsky, 1975). In our manuscript we provide the first biophysical model for aspect-ratio homeostasis in rod-like bacteria and elucidate the underlying molecular mechanism that will inform future experimental studies. Our findings push the field of ‘bacterial cell size control’ to a new direction, which has so far focused on the individual control of cell volume, length or width, neglecting how bacterial length and width are coupled to give rise to rod-like cell shapes.

To support our model, we collected a large number of cell shape data (~5000 conditions) from many different laboratories, which indeed confirm that aspect ratio is conserved in E. coli (and 7 other organisms) under many different perturbations to the growth conditions (Figure 1). Importantly, our model also predicts under what conditions E. coli cells may deviate from their homeostatic aspect ratio of 4:1, and we tested our quantitative predictions for filamentous and spherical cell shapes against experimental data (Figure 2). Therefore, our thesis is not solely about the maintenance of 4:1 aspect ratio in E. coli, but more broadly about the control of bacterial cell shapes under many different perturbations.

It is evident from the reviewer comments that the principal issues with our manuscript lie in the presentation of results (e.g. claim of ‘universality’), and concerns about the novelty of our model in the context of previous studies. This is in part due to inadequate communication on our part. Having carefully read and deliberated on the reviewer comments, we believe that an improvement in the presentation of our results, an increased clarity of writing, and expanded description of the model will address all reviewer comments thoroughly and completely.

Reviewer #1:

This manuscript compiles data on the length, width, and growth rate of E. coli under a number of experimental perturbations, such as changes in growth medium, incubation in antibiotics, inhibited protein synthesis, etc. and shows that the surface area to volume ratio is strongly conserved. This result is consistent with recent work from Julie Theriot's group (which is cited in Harris and Theriot, 2016; 2018). In this manuscript, the author's also add data from other rod-shaped bacteria that show similar behavior (Figure 1G). The authors use this result to develop a model for aspect ratio regulation that is based on exponential growth of the bacterial length at fixed width, FtsZ production at a rate proportional to volume growth rate (which by assuming constant width also assumes that the FtsZ production is proportional to the length growth rate), and division that occurs when FtsZ production reaches a critical value that is proportional to the width. This model predicts a constant aspect ratio and the authors then go on to predict the dynamics under impulse type perturbations.

We appreciate the succinct summary of our work. While our findings are consistent with a recent model proposed by Julie Theriot’s group, it is important to note the key difference: Harris and Theriot showed that rod-shaped bacteria (E. coli, C. crescentus) maintain a homeostatic surface-to-volume ratio (S/V) in a growth-rate dependent manner. Here we uncover a much stronger geometric constraint that bacteria (especially E. coli) maintain the relationship S=γV^2/3 (with a constant γ), independent of growth rate. Furthermore Harris and Theriot model does not lead to aspect ratio control, as addressed in response to a comment from Reviewer #2.

I have two main concerns with the manuscript.

First, the only novelty of the model is the assumption that there is a critical amount of FtsZ required to divide the cell and that this depends on the width. I think that this is a reasonable assumption, but I also feel that the overall results are fairly obvious. That is, it is not clear that the model provides a significant advance in our understanding.

Our model and analysis expand the current state of understanding in many ways:

– We provide the first biophysical model for aspect ratio control in bacteria and identify the molecular origins. We support our model by large population measurements (~5011 growth conditions in E. coli, 50 different bacterial species) and single cell measurements in mother machine (n~80,000). Why E. coli cells maintain a constant aspect ratio has been puzzling for more than half a century (Zaritsky, 1975; Zaritsky, 2015), with no prior existing model.

– The model for aspect ratio homeostasis provides a conceptual leap in the field of bacterial cell size control, by showing that added cell length for rod-like cells is coupled to their diameter. Previous phenomenological models treated cell length and diameter as independent control variables (Taheri-Araghi et al., 2015, Harris and Theriot, 2016).

– We think it’s not obvious that E. coli cells preserve their aspect ratios under multitude of perturbations to the nutrient conditions, ribosomes, protein overexpression or deletion (Figure 1A). Our model not only identifies that the maintenance of aspect ratio emerges from balanced biosynthesis of growth and division proteins (k/k_P constant), but also predicts under what perturbations cells may deviate from their homeostatic aspect ratio of 4:1 (Figure 4). We provide a quantitative, experimentally testable model for cell shape control that goes beyond just the regulation of FtsZ kinetics.

That said, there is also a small problem with the model, in that we would expect the binding rate of FtsZ to depend on the surface to volume ratio (which turns out not to matter, because the authors end up making assumptions that the rate the ring is built is equal to the protein production rate).

It is unclear why the FtsZ binding rate should depend on S/V. The rate equations are formulated in terms of the amount of the surface-bound and cytoplasmic FtsZ, and not their concentrations. If only the rate equations were formulated in terms of the concentration of cytoplasmic (c) and surface-bound proteins (c_r), then the rate of increase in surface bound concentration of FtsZ would naturally depend on S/V.

dcrdt=kcr+kbVSc-kdcr

Second, if we were to consider changes in FtsZ binding rate with S/V, our simulations show that this has negligible effect on aspect ratio. For rod-shaped cells, S/V ~1/w in the first approximation, where w is cell width. Since during one cell cycle width doesn’t change, S/V stays approximately constant (new Figure 4—figure supplement 1A). If the width of the bacteria is changing due to changes in growth conditions, overall binding rate may be affected by (S/V) since the area of the Z-ring = δw ~δ/(S/V), where δ is the lateral width of the FtsZ ring. E. coli width changes in different growth conditions from approximately 0.5 till 1 μm (Taheri-Araghi et al., 2015), so S/V can change by a maximum factor of 2. To address the effect of changes in binding rate, we simulated our dynamic model by changing the ratio k_b/k_d across 4 order of magnitude. The figure in Figure 4-figure supplement 1B shows the dependence of cellular newborn aspect ratio (n = 10000, during steady-state growth) on k_b/k_d. In the limit k_b>>k_d, aspect ratio~4 as expected. However, the factor of 2 change even for the border line case of k_b/k_d = 10, has negligible impact on cell aspect ratio.

However, as noted by the reviewer, the rate of FtsZ recruitment to the Z-ring (~10s, Soderstrom et al., Nat Commun 2018) is much faster than the growth rate. As a result, the rate at which the ring is built is determined by the production rate of FtsZ in the cytoplasm.

Note also that there is a typo in the equation fordPrdt, which has kd multiplying both rates.

We corrected the typo in the manuscript and thank the reviewer for pointing this out.

Second, and more important, is that while the results match well with the data, there are a number of aspects of the presentation that are misleading. The title claims that the results/model presented here are universal. In Figure 1G, the authors select 7 bacteria to claim that the scaling of SA = 2 π V^(2/3) is ubiquitous among bacteria. As noted, this also suggests that an aspect ratio of ~4 is the rule for rod-shaped bacteria. This is not true. As an example, myxococcus xanthus has an aspect ratio around 7-8 and spirochetes have aspect ratios of ~30! Even more importantly, single species don't always maintain the same aspect ratio. In B. subtilis, the aspect ratio can vary between at least 3.8 – 8 (see Ilkanaiv et al., 2017). Therefore, this model may be applicable to E. coli (and possibly some other bacteria), but it is not universal.

We apologize for the misunderstanding, which may have been triggered by a lack of clarity in our presentation. In our original submission, we did not claim that the 4:1 aspect ratio, or equivalently S = 2πV ^2/3, is universal. Instead, we found that a ‘universal’ scaling law S=γV^2/3 is conserved amongst rod-shaped or coccoid bacterial species, implying the maintenance of a fixed aspect ratio (Figure 1A and E, dataset expanded). It is indeed possible for different bacteria to have different values for γ. For instance, in Figure 1E (previously 1G) we show that the coccoid S. aureus under different perturbations maintain the relation S = 4.92 V ^2/3, implying preservation of the same scaling factor (2/3) while maintaining a different aspect ratio (1.38 +/- 0.18). In the same figure (1E), we now show data for a total of 48 different rod-shaped bacteria, and 1 rod-shaped Archaea (H. vulcanii), all of which remarkably follow the curve S = 2πV^2/3.

Furthermore, our model also predicts how the aspect ratio and cell width can be altered by changing (k/k_p) and (k/β), leading to filamentous or spherical cells, in agreement with available experimental data. In Figure 4 (previously Figure 2) we show that our model indeed predicts the breakdown of the 4:1 aspect ratio in E. coli under FtsZ or MreB perturbations.

However, the reviewer has made an excellent point that long filamentous cells, such as Spirochetes, do not necessarily conserve their aspect ratios. In Figure 1E, we now also include the data for Spirochetes, as one of the exceptions to the rule S=γV^2/3. We have therefore removed the term ‘universal’ from the title and Abstract of our paper. The fact, however, remains that E. coli remarkably conserve their aspect ratios under diverse size perturbations spanning two orders of magnitude (Figure 1A), and so do 50 other cell types (Figure 1E).

Motivated by the comments of reviewers 1 and 2, we now include a schematic in Figure 1D to illustrate the expected scaling relations for different bacterial shapes. Filamentous cells (Helicobacter, Spiroplasma, Spirochetes, Myxobacter) would likely follow the relation SμV, whereas coccoid or the rod-shaped cells follow the scaling law: SμV^2/3.

Reviewer #2:

In this study the authors set out to study the size and shape of a wide range of 'rod-shaped' cells by collecting image data from at least seven different species and thousands to total conditions (genotype x nutrients x antibiotics). Across all these conditions, the authors find a simple scaling law to the surface area/volume ratio, namely a scaling that preserves the aspect ratio of the cell at roughly 4:1. Given this observation, they build a simple, mechanistically inspired, quantitative model for the growth of the cell. Using this model, they are able to tune three parameters (k, k_p, β) to match a collection of genetic knockdown and antibiotic treatment experiments.

Understanding how cell size and shape homeostasis is maintained throughout the bacterial kingdom is a very interesting and important problem and these authors should be commended for pushing the community to consider that these mechanisms may be conserved across a wide phylogenetic range. However, given the extensive body of literature already available about cell size/shape homeostasis, and, in particular the review mentioned by the authors by Harris and Theriot, the scientific bar for productive engagement on this topic is already quite high. Much of the intellectual driving force for this work seems to follow directly from the hypothesis from Harris and Theriot that "While many studies have treated volume as the actively controlled parameter in this scenario, our recent work suggests that it is likely the other way around, and that SA/V is the actively regulated variable, with size following along as necessary [13]." The current work seeks to extend or provide alternatives for the mechanistic models presented in Harris and Theriot as well as integrate additional data in other species. However, given that the idea of SA/V scaling conservation is not new, appealing to a broad audience such as that of eLife would require experimental validation of their mechanistic model.

In addition to the core concerns around novelty of the central hypothesis and validity of the mechanistic model, there are a few issues the authors might choose to consider:

We thank the reviewer for summarising the key aspects of our work and recognizing the importance of the field of study. Below we address some of the key comments raised above. “extensive body of literature already available about cell size/shape homeostasis” –A lot of work has been done over the past five years on developing phenomenological models for cell size control. Phenomenological models for the homeostasis of bacterial cell shape are treated the control of cell length separately from the control of cell width in rod-shaped bacteria. We provide a molecular model to show for the first time that bacterial cell dimensions are coupled to preserve aspect ratio, thereby linking the field of cell size and shape homeostasis.

“Much of the intellectual driving force for this work seems to follow directly from the hypothesis from Harris and Theriot” – Our model draws evidence from multiple recent experimental studies, while questioning the Harris and Theriot (HT) hypothesis. It is important to recognize the key distinctions between the two models. HT model does not lead to conservation of S-to-V scaling or aspect ratio, instead it leads to a model for the control of cell width (Eq. 3). HT model infers that S/V ratio is a function of growth media, such that cells reach a new homeostatic value of S/V upon perturbations in the growth rate. Here instead we propose a much stronger constraint that cells preserve the scaling relation, S = μV ^2/3 (μ a constant) under diverse growth perturbations (~5000 conditions) across ~50 different bacterial species. Furthermore, HT model is agnostic about molecular mechanisms. Here we provide an explicit molecular candidate (FtsZ) for bacterial shape control, in agreement with exciting new evidence from Si et al., 2019. Taken together, our model integrates the adder model for cell size homeostasis with the regulation of S/V ratio and FtsZ, providing an integrative framework that successfully predicts bacterial shape control with only three physiological parameters.

“the idea of SA/V scaling conservation is not new” – We are not aware of any other studies that propose conservation of the scaling relation S = μV ^2/3 across growth conditions, nor provide a model for it. Others have only shown evidence for the regulation of surface-to volume ratio by growth rate, which is a natural consequence of our model (Figure 1C).

“appealing to a broad audience such as that of eLife would require experimental validation of their mechanistic model” – Our model is tightly grounded in experimental data (see Figures 1-4), and we compare our model predictions extensively to experimental data, throughout the manuscript. As we are not an experimental lab, we have compiled data from a number of different laboratories to show that our model is consistent with all the available cell shape data across ~50 bacterial species and ~5000 growth conditions for E. coli. We definitely welcome suggestions to test our model further.

Major points:

1) The authors should clearly explain how their mechanistic model contrasts with the cell wall-focused model proposed by Harris and Theriot and should strive to propose experiments with predicted outcomes that would differentiate a peptidoglycan centric model from an FtsZ centric model. If the data already exist to rule out one of them, this should be clearly presented.

We agree with the reviewer that clearer discussion of the contrast between our model and that of Harris/Theriot should be articulated in the manuscript. In the revised manuscript, we have expanded the discussion to highlight the key differences between these two models.

Foremost amongst the comparison is that Harris and Theriot proposes a homeostatic regulation of S/V in a growth-rate dependent manner. Whereas we propose a much stronger geometric constraint that the scaling relation S = μV ^2/3 is preserved independent of growth rate. This result, however, does not contradict the model of Harris/Theriot.

Second, Harris and Theriot proposed a model where cells divide once a threshold amount of excess surface area material, ΔA, is accumulated in the cell. From this model it follows that, ΔA = ΔV (β/k – 2/r) = constant, where r is the cell radius of cross-section. This in turn 1, which is in contradiction to experimental data (Figure 1).

Third, we can indeed propose several experimental tests for our model, as highlighted in the revised manuscript:

– Our model would predict that FtsZ overexpression leads to minicells while FtsZ deletion would induce elongated phenotypes (Figure 4A). These predictions are consistent with data from Potluri et al., 1999, and Zheng et al., 2016.

– Oscillations in FtsZ amount would lead to cell size oscillations, in agreement with new data from Si et al., 2019.

– Total abundance of FtsZ scales with cell diameter, in agreement with data from Shi et al., 2017.

– We further predict that FtsZ knockdown would break aspect ratio preservation, whereas targeting cell wall precursors would change growth rate, but not alter aspect ratio or the scaling relation S = μV ^2/3. Figure 4—figure supplement 1C shows surface-to-volume scaling for E. coli cells treated with Fosfomycin that target MurA (affecting cell wall biogenesis) and FtsZ depletion. We find that Fosfomycin treated cells preserve the S~V^2/3 scaling, whereas FtsZ depletion breaks the S~V^2/3 scaling. This is a clear contrast between the role of cell wall precursors and FtsZ on cell shape control, implying that a cell wall precursor-based model alone is not sufficient to account for shape changes.

As one such example, the authors show that tuning one parameter (k_p) is consistent with the experimental notion of knocking down the production of FtsZ. However, they fail to show if there is quantitative agreement between the production rate of FtsZ and the amount they expect to need to change k_p (40%).

Our model predicts that a reduction in FtsZ production rate to 40% of WT leads to observed phenotype in Zheng et al., 2016. This is consistent with reduction of relative mRNA to ~ 40% corresponding to addition of 3 ng/ml of aTc (Figure 2B of the Zheng et al.). We comment on this in our manuscript and thank the reviewer for pointing this out.

2) The use of 'universal' in the paper's title significantly oversells of the breadth of species included in the observations and a power law describing data which span roughly one order of magnitude. While the authors do include a large collection of data, the collection is far from comprehensive for all size/shape data available and the authors do not clearly indicate why they limited themselves to the data they did. A quick literature search reveals anecdotal evidence of bacterial sizes that are much smaller than a micron such as Brevundimonas (PDA J Pharm Sci Technol. 2002 Mar-Apr;56(2):99-108.) to nearly a millimeter in length Epulopiscium (J. Protozoal., 35(4), 1988, pp. 565-569). Granted, these publications may not have the same type of data necessary to integrate it directly into their model, but for a discussion of the 'universal scaling', the authors should push themselves to cover as large of a length-scale as possible. When choosing a set of species for inclusion in this study, it seems like the microbiology community may have already picked an aspect ratio of about 4:1 in its definition of rod-shaped bacteria. For example, cells that have a much shorter aspect ratio are given the term ovoid or lancet (Streptococcus pneumoniae) or spherical (Staphylococcus aureus included here) and ones that are much longer are called filamentous (Sphaerotilus natans). Confusingly, these authors do not include species that have been traditionally classified as rod-shaped cells with a longer aspect ratio such as (Helicobacter, Spiroplasma, Spirochetes, Myxobacter).

We address this point in response to the first reviewer. Both the reviewers have raised a pertinent point that long filamentous cells, such as Spirochetes, do not necessarily conserve their aspect ratios. In Figure 1E, we now include available shape data for Spirochetes, as one of the exceptions to the rule S=γV ^2/3. We have therefore removed the term ‘universal’ from the title and Abstract of our paper. The fact, however, remains that rod-shaped bacteria (E. coli) remarkably conserve their aspect ratios under diverse size perturbations spanning two orders of magnitude (Figure 1A).

In Figure 1E we have now expanded the dataset to cover two orders of magnitude by including 49 different rod-shaped bacterial species and 1 rod-shaped Archea. All of them lie on the curve S=γV ^2/3, confirming our predictions. In addition, we have also expanded the E. coli dataset by 30 more nutrient growth conditions (Gray et al., 2019), confirming our initial statement of aspect-ratio homeostasis.

We are grateful to the reviewer for providing the papers reporting drastic volume range in bacteria spanning 2 orders of magnitude. Bacteria that we include in Figure 1E are those that are known divide using FtsZ machinery during binary fission. This is to maintain consistency with our model that is based on FtsZ regulation. For this reason, we did not include Epulopiscium in our analysis. We also did not include Sphaerotilus natans in our graph as we could not find good shape measurements for it. In line with reviewer’s comments we have now included longer filamentous cells in Figure 1E. We have also introduced a new cartoon in Figure 1D showing how long filamentous cells that keep their width constant, would have a different scaling law S~V.

3) I'm not entirely convinced that the universal scaling applies within the single cell data (Figure 1D). By plotting the single cell data from a variety of experiments, the range of the data seems to put a larger priority on the averages. However, within each condition there seems to be clear deviations from the 'single aspect ratio', consistent with the author's single cell growth model that cells grow without changing their diameter before dividing. This should result in a roughly factor of two change in aspect ratio from birth to division. I think this is what the authors refer to in the fourth paragraph of the Introduction but should discuss more fully.

In our original submission, we had already explored in detail the deviation from 2/3 scaling in the single-cell data (Figure 2—figure supplement 1A-B). The main reason for the deviation from 2/3 scaling comes from large fluctuations in newborn length for a given width of bacteria. Using our model, we can quantitatively explain the deviation from the universal scaling by incorporating experimentally measured fluctuations in cell width and length, in agreement with experimental data. We have now attempted to explain this better in the manuscript and in the supplementary figure caption.

4) I do not understand Figure 2B at all. In particular, the binning of the data that I have been able to find in Taheri-Araghi et al., 2015, is binned by the size of cells at birth, not the individual cell growth rate. Further, the authors do not describe how they go from the data in Taheri-Araghi et al., 2015, to the data in Figure 2B, but it could be that they obtained the raw data from the authors and performed a new type of analysis. If so, a description of this process should be included.

We were kindly provided with the raw data for single-cell width and length at various growth rates (conditions) by the Suckjoon Jun lab. We reanalysed the data, performed the necessary binning and analysis. We have clearly stated this in the Appendix and in each figure caption.

5) I am unclear about why the MreB and FtsZ knockdown data from Si et al. is included in the bulk Figure 1A data but the MreB and FtsZ knockdown data from Zheng et al. is treated as a completely separate experiment. If the approach that these two studies used was different, it may be helpful to explain why some data are included in one place and others are not.

For consistency, we now plot the MreB and FtsZ knockdown data from Si et al. in Figure 4B. The knockdown data from Si et al. cover a small dynamic range so it is hard to extract a clear trend from these data alone. This is presumably because cells in those knockdown experiments were grown in slow growth media (MOPS glucose + 6 a. a., with growth rate ~0.75 h-1) and small perturbations, whereas the data from Zheng et al. that show drastic cell shape changes (Figure 4B) are obtained from experiments on rich media (RDM + glucose, with growth rate 1.6 h-1) and large perturbations. The trend in Si et al. seems to be consistent with those in Zheng et al.

Author details

1. Nikola Ojkic

   Department of Physics and Astronomy, Institute for the Physics of Living Systems, University College London, London, United Kingdom

   Contribution

   Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—original draft

   Competing interests

   No competing interests declared

2. Diana Serbanescu

   Department of Physics and Astronomy, Institute for the Physics of Living Systems, University College London, London, United Kingdom

   Contribution

   Data curation, Investigation, Methodology

   Competing interests

   No competing interests declared

3. Shiladitya Banerjee

   Department of Physics and Astronomy, Institute for the Physics of Living Systems, University College London, London, United Kingdom

   Present address

   Department of Physics, Carnegie Mellon University, Pittsburgh, United States

   Contribution

   Conceptualization, Resources, Supervision, Funding acquisition, Validation, Investigation, Visualization, Writing—original draft, Project administration

   For correspondence

   shiladitya.banerjee@ucl.ac.uk

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-8000-2556

• 8,912

  Page views

• 614

  Downloads

• 37

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.