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

  1. Jarno Mäkelä
  2. Alexandros Papagiannakis
  3. Wei-Hsiang Lin
  4. Michael Charles Lanz
  5. Skye Glenn
  6. Matthew Swaffer
  7. Georgi K Marinov
  8. Jan M Skotheim
  9. Christine Jacobs-Wagner  Is a corresponding author
  1. Howard Hughes Medical Institute, Stanford University, United States
  2. Sarafan Chemistry, Engineering, and Medicine for Human Health Institute, Stanford University, United States
  3. Institute of Biotechnology, Helsinki Institute of Life Science, University of Helsinki, Finland
  4. Department of Biology, Stanford University, United States
  5. Chan Zuckerberg Biohub, United Kingdom
  6. Department of Genetics, Stanford University, United States
  7. Department of Microbiology and Immunology, Stanford School of Medicine, United States
14 figures and 13 additional files

Figures

Figure 1 with 8 supplements
Growth rate and genome copy number in E. coli growing in M9glyCAAT.

(A) Illustration of 1N (CRISPR interference [CRISPRi] oriC, CJW7457) and multi-N (CRISPRi ftsZ, CJW7576) cells with different numbers of chromosomes along with representative microscopy images at …

Figure 1—figure supplement 1
The relative growth rate of wild-type (WT) (CJW7339) cells in M9glyCAAT grown after placing them on an agarose pad.

The plot includes 61,562 datapoints from 11,907 cells. Lines and shaded areas denote mean ± SD from three biological replicates.

Figure 1—figure supplement 2
Characterization of ploidy in CRISPR interference (CRISPRi) oriC cells.

Representative microscopy images of CRISPRi oriC cells expressing HU-CFP and ParB-mCherry with parS site at ori1 (CJW7517) in M9glyCAAT 210 min after addition of 0.2% L-arabinose. Scale bar: 1 µm.

Figure 1—figure supplement 3
Effects of cell area overestimation on the relative growth rate calculation.

The simulated relative growth rate for increasing cell area considering different levels of cell area overestimation (0–0.2 μm2) and a constant relative growth rate of 0.014 min–1. In the absence of …

Figure 1—figure supplement 4
Relationship between growth rate and cell area in cell types of different ploidy.

(A) Absolute and (B) relative growth rate in M9glyCAAT in 1N (32,735 datapoints from 1568 cells, CJW7457), multi-N (14,006 datapoints from 916 cells, CJW7576), and wild-type (WT) (19,495 datapoints …

Figure 1—figure supplement 5
Validation of stable growth under microscope observation and absolute growth rate determination of ppGpp0 and ∆recA cells.

(A) Plot showing the absolute growth rate of 1N (CJW7457) cells imaged for 2 hr (dotted line). Prior to spotting cells on the agarose pad containing medium, cells were grown for 90 min in a liquid …

Figure 1—figure supplement 6
Characterization of ploidy in dnaC2 cells.

(A) Representative microscopy images of wild-type (WT) (top) and dnaC2 (bottom) cells expressing HU-mCherry growing under microscope observation in M9glyCAAT. DNA replication in dnaC2 cells was …

Figure 1—figure supplement 7
Relationships between growth rate and cell volume across cell types of different ploidy in M9glyCAAT.

(A) Absolute and (B) relative growth rate based on cell volume in 1N (32,735 datapoints from 1568 cells, CJW7457), multi-N (14,006 datapoints from 916 cells, CJW7576), and dnaC2 1N (13,933 …

Figure 1—figure supplement 8
DNA-dependent growth in C. crescentus.

(A) Plot showing the absolute and (B) relative growth rate of 1N (DnaA depletion, strain CJW4823, 87 cells) and multi-N (FtsZ depletion, strain CJW3673, 181 cells) C. crescentus cells as a function …

Figure 2 with 2 supplements
Lower ribosome activity explains the reduced growth rate of 1N cells growing in M9glyCAAT.

(A) RpsB-msfGFP fluorescence concentration in 1N (6542 cells, CJW7478) and multi-N (10,537 cells, CJW7564) cells as a function of cell area. Lines and shaded areas denote mean ± SD from three …

Figure 2—figure supplement 1
Diffusive characteristics of labeled ribosomes in rifampicin-treated wild-type (WT) cells.

Plot showing the probability density of apparent diffusion coefficients (Da) of JF549-labeled RpsB-HaloTag in WT cells (CJW7528) treated with 200 µg/mL rifampicin for 30 min. Only tracks of length …

Figure 2—figure supplement 2
Diffusive characteristics of labeled ribosomes in 1N cells as a function of cell area.

Plots showing the probability density of apparent diffusion coefficients (Da) of JF549-labeled RpsB-HaloTag in 1N cells (CJW7529) in M9glyCAAT at different cell areas. Also shown is Da fitted by …

Figure 3 with 2 supplements
RNA polymerase (RNAP) activity is reduced in 1N cells growing in M9glyCAAT.

(A) RpoC-YFP fluorescence concentration in 1N (3580 cells, CJW7477) and multi-N (5554 cells, CJW7563) cells as a function of cell area. Lines and shaded areas denote mean ± SD from three …

Figure 3—figure supplement 1
Diffusive characteristics of labeled RNA polymerases (RNAPs) in rifampicin-treated cells.

Plot showing the probability density of apparent diffusion coefficients (Da) of JF549-labeled RpoC-HaloTag in CJW7519 cells treated with 200 µg/mL rifampicin for 30 min. Also shown is Da fitted by …

Figure 3—figure supplement 2
Determination of the relative Rsd concentration in 1N-rich and multi-N cells as a function of cell area.

Plots showing the relative protein concentration of Rsd in 1N-rich (SJ_XTL676) and multi-N (SJ_XTL229) cells, as determined by tandem-mass-tag (TMT)-mass spectrometry (MS). 1N-rich cells grown in …

Figure 4 with 3 supplements
RNASelect and EUB338 concentration measurements in 1N and multi-N cells.

(A) Images of representative cells from a mixed population of 1N (CRISPR interference [CRISPRi] oriC) and multi-N (CRISPRi ftsZ) cells. Strains CJW7457 and CJW7576 carrying HU-mCherry were used for …

Figure 4—figure supplement 1
Cell sampling to match cell size distribution in mixed populations of 1N and multi-N cells.

(A) Cell area distributions from mixed CJW7457 (CRISPR interference [CRISPRi] oriC) and CJW7576 (CRISPRi ftsZ) populations stained with SYTO RNASelect (aggregated data from five biological …

Figure 4—figure supplement 2
Comparison of RpoC-HaloTag-JF549 labeling between 1N and multi-N cells.

(A) Phase contrast (left) and RpoC-HaloTag-JF549 fluorescence (right) images of two representative cells from a mixed population of 1N (CRISPR interference [CRISPRi] oriC, CJW7520) and multi-N …

Figure 4—figure supplement 3
EUB338 staining comparison between fast (M9glyCAAT) and slow (M9gly) growing populations.

The time of dCas9 induction was adjusted to obtain comparable cell area distribution between 1N and multi-N cells (see Materials and methods). The cell areas were then sampled to achieve a perfect …

Figure 5 with 2 supplements
Mathematical modeling of DNA limitation.

(A–C) Plots comparing simulation results of model A (solid lines) with experimental data points (dots) and averages (open squares) in the M9glyCAAT condition. The multi-N and 1N cells are indicated …

Figure 5—figure supplement 1
Comparison between experimental results from the M9glyCAAT condition and simulation results using model B.

(A–C) Plots comparing simulation results of model B (solid lines) with experimental data (dots) and averages (open squares). The multi-N and 1N cells are indicated as blue and yellow, respectively: …

Figure 5—figure supplement 2
Model A-based simulations examining the effects of varying the rates in either mRNA synthesis or mRNA degradation on the relative growth rate of 1N cells as a function of cell area.

(A) Plot showing the decay of DNA concentration (black) and of the relative growth rate (blue) in 1N cells when the rate of bulk mRNA synthesis (r1) increases or decreases by 10-fold. Each quantity …

Figure 6 with 2 supplements
Scaling of the total active RNA polymerases (RNAPs), total active ribosomes, and growth rate with cell area during genome dilution in nutrient-poor media.

(A) Plot showing the total amount of active RNAPs (calculated by multiplying the total amount of RNAPs by the fraction of active RNAPs from Figure 6—figure supplement 1A and G) in wild-type (WT) …

Figure 6—figure supplement 1
Characterization of RNA polymerase (RNAP) diffusion and active fraction in poor media conditions.

(A) Plot showing the RpoC-YFP fluorescence concentration in 1N (CJW7477) cells grown in M9gly (three experiments) as a function of cell area. Lines and shaded areas denote mean ± SD between …

Figure 6—figure supplement 2
Characterization of ribosomal diffusion and active fraction in poor media conditions.

(A) Plot showing the RpsB-msfGFP fluorescence concentration in 1N (CJW7478) cells grown in M9gly (from three experiments) as a function of cell area. Lines and shaded areas denote mean ± SD between …

Figure 7 with 4 supplements
Proteome and transcriptome remodeling in 1N-rich cells.

(A) Schematic explaining the calculation of the protein slopes, which describes the scaling of the relative protein concentration (concentration of a given protein relative to the proteome) with …

Figure 7—figure supplement 1
Comparison of protein and mRNA scaling between biological replicates.

(A) Correlation of protein slopes across the proteome (2360 proteins) between two biological replicates for 1N cells. (B) Same as (A) but for multi-N cells. (C) Correlation of RNA slopes across the …

Figure 7—figure supplement 2
Comparison of our data with reference datasets.

(A) Comparison between the concentrations of 3446 RNAs present in both our dataset and that of Balakrishnan et al., 2022. The first time point (60 min) after CRISPR interference (CRISPRi) induction …

Figure 7—figure supplement 3
Protein slopes relative to the chromosome position of their gene in 1N-rich and multi-N cells.

(A) Relationship between the absolute distance of a gene from oriC and the slope of the protein it encodes. Data from 2268 proteins are shown (colormap: Gaussian kernel density estimation), as well …

Figure 7—figure supplement 4
Protein slopes relative to protein ion intensity for 1N-rich cells.

The summed ion intensity of each protein was divided by the protein sequence length and the quotient was log-transformed. The locations of selected proteins are annotated.

Appendix 1—figure 1
Estimation of DNA concentration.

(A) DNA replication pattern for different medium conditions. (B) Extraction of the genome content, cell volume, and genome concentration along the division cycle. (C) Genome copy, cell volume, and …

Appendix 2—figure 1
Interpolation and extrapolation of parameters.

The parameters of the ODE models were calculated using the fitted formulae. The obtained values are summarized in Supplementary file 9.

Appendix 3—scheme 1
Topology of the three RNA-polymerase states and their transition fluxes.
Appendix 3—figure 1
Comparing changes in active RNAP fraction between ODE models A and B.

The two-dimensional colormaps show the values of active RNAP fraction (αRNAP) under different promoter and RNAP concentrations. (A) Using the formula of model A with M9glyCAAT parameters. (B) Using …

Appendix 4—figure 1
Comparison between initial (Ini) and optimized (Opt) parameters.

(A) Parameters used in model A. (B) Parameters used in model B.

Author response image 1

The mRNA production rate equivalent (mRNA abundance at the first time point after CRISPRi oriC induction multiplied by the mRNA degradation rate measured by Balakrishnan et al., 2022, PMID: …

Author response image 2

The relative fitness of each gene (data by Hawkins et al., 2020, PMID: 33080209, median fitness from the highest sgRNA activity bin) plotted versus the gene-specific RNA and protein slopes that we …

Additional files

Supplementary file 1

Gene-specific protein slopes calculated from the tandem-mass-tag mass spectrometry measurements in growing 1N-rich or mutli-N cell populations.

https://cdn.elifesciences.org/articles/97465/elife-97465-supp1-v1.xlsx
Supplementary file 2

Description of the model parameters.

https://cdn.elifesciences.org/articles/97465/elife-97465-supp2-v1.docx
Supplementary file 3

Initial and optimized model parameters.

https://cdn.elifesciences.org/articles/97465/elife-97465-supp3-v1.docx
Supplementary file 4

Gene-specific RNA slopes calculated using RNA sequencing in growing 1N-rich cell populations.

The RNA slopes are also compared with the average protein slopes.

https://cdn.elifesciences.org/articles/97465/elife-97465-supp4-v1.xlsx
Supplementary file 5

Comparison between transcriptome and proteome remodeling statistics (RNA and protein slopes, respectively) with previously published gene expression statistics (Balakrishnan et al., 2022), or gene essentiality data (Gerdes et al., 2003; Goodall et al., 2018; Hashimoto et al., 2005).

https://cdn.elifesciences.org/articles/97465/elife-97465-supp5-v1.xlsx
Supplementary file 6

Strains used in this study.

The abbreviations kan, cat, and spec refer to gene cassette insertions conferring resistance to kanamycin, chloramphenicol, and spectinomycin, respectively. These insertions are flanked by Flp site-specific recombination sites (frt) that allow the removal of the insertion using Flp recombinase from plasmid pCP20 (Cherepanov and Wackernagel, 1995).

https://cdn.elifesciences.org/articles/97465/elife-97465-supp6-v1.docx
Supplementary file 7

Oligonucleotides used in this study.

https://cdn.elifesciences.org/articles/97465/elife-97465-supp7-v1.docx
Supplementary file 8

Sizes, DNA concentrations, and growth rates of cells in different growth media.

https://cdn.elifesciences.org/articles/97465/elife-97465-supp8-v1.docx
Supplementary file 9

mRNA and protein numbers per cells and bulk rates of transcription and translation.

This table describes how kinetic constants were estimated from the literature. Parameters Xini, Yini, r1,r2 were used in ordinary differential equation (ODE) simulations. We assumed exponential growth and used the relation Mini=M/2log2, where Mini is the biomass of a newborn cell and M is the average cellular biomass in the population (Koch and Schaechter, 1962). The bulk transcription rate r1 was defined as r1=mRNAsynthesisrateofproteinsthecell , and the bulk translation rate r2 was defined as r2=proteinsynthesisrateofproteinsthecell . Values were estimated from a previous study (Bremer and Dennis, 2008). For our estimations, we assumed that the average protein length is 310 amino acids and that the average mRNA length is about 1 kb (Ishihama et al., 2008).

https://cdn.elifesciences.org/articles/97465/elife-97465-supp9-v1.docx
Supplementary file 10

DNA and mRNA affinity constants (K1, K2).

The active fraction of RNA polymerases (RNAPs) and ribosomes, αRNAP and αribo, are given by the formulae αRNAP=ZK1+Z and αribo=XK2+X, where [Z] and [X] are the DNA concentration and the mRNA concentration in the cells, respectively.

To infer the parameters K1, K2, we used the values of αRNAP, αribo, [Z], and X of wild-type cells determined in our study and back-calculated the values of K1, K2.

https://cdn.elifesciences.org/articles/97465/elife-97465-supp10-v1.docx
Supplementary file 11

Estimation of the mRNA degradation rate.

The mRNA degradation data were obtained from an experimental study (Balakrishnan et al., 2022).

https://cdn.elifesciences.org/articles/97465/elife-97465-supp11-v1.docx
Supplementary file 12

Comparison between the model formulation with constant RNA polymerase (RNAP) concentration (model A) and the more complex model with increasing RNAP concentration (model B).

https://cdn.elifesciences.org/articles/97465/elife-97465-supp12-v1.docx
MDAR checklist
https://cdn.elifesciences.org/articles/97465/elife-97465-mdarchecklist1-v1.docx

Download links