Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Editor's Pick Research Article | Molecular Biology and Physiology

Escherichia coli Has a Unique Transcriptional Program in Long-Term Stationary Phase Allowing Identification of Genes Important for Survival

Karin E. Kram, Autumn L. Henderson, Steven E. Finkel
Xiaoxia Lin, Editor
Karin E. Kram
aDepartment of Biology, California State University, Dominguez Hills, Carson, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Karin E. Kram
Autumn L. Henderson
aDepartment of Biology, California State University, Dominguez Hills, Carson, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Steven E. Finkel
bMolecular and Computational Biology Section, Department of Biological Sciences, University of Southern California, Los Angeles, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiaoxia Lin
University of Michigan—Ann Arbor
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00364-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Microbes live in complex and constantly changing environments, but it is difficult to replicate this in the laboratory. Escherichia coli has been used as a model organism in experimental evolution studies for years; specifically, we and others have used it to study evolution in complex environments by incubating the cells into long-term stationary phase (LTSP) in rich media. In LTSP, cells experience a variety of stresses and changing conditions. While we have hypothesized that this experimental system is more similar to natural environments than some other lab conditions, we do not yet know how cells respond to this environment biochemically or physiologically. In this study, we began to unravel the cells’ responses to this environment by characterizing the transcriptome of cells during LTSP. We found that cells in LTSP have a unique transcriptional program and that several genes are uniquely upregulated or downregulated in this phase. Further, we identified two genes, cspB and cspI, which are most highly expressed in LTSP, even though these genes are primarily known to respond to cold shock. By competing cells lacking these genes with wild-type cells, we show that these genes are also important for survival during LTSP. These data can help identify gene products that may play a role in survival in this complex environment and lead to identification of novel functions of proteins.

IMPORTANCE Experimental evolution studies have elucidated evolutionary processes, but usually in chemically well-defined and/or constant environments. Using complex environments is important to begin to understand how evolution may occur in natural environments, such as soils or within a host. However, characterizing the stresses that cells experience in these complex environments can be challenging. One way to approach this is by determining how cells biochemically acclimate to heterogenous environments. In this study, we began to characterize physiological changes by analyzing the transcriptome of cells in a dynamic complex environment. By characterizing the transcriptional profile of cells in long-term stationary phase, a heterogenous and stressful environment, we can begin to understand how cells physiologically and biochemically react to the laboratory environment, and how this compares to more-natural conditions.

INTRODUCTION

Experimental evolution studies of bacteria have revealed many insights into evolutionary processes (1–3). Many of these experiments have been performed in environments where only one factor is being experimentally manipulated, for instance, starvation for one nutrient (4), heat shock (5), antibiotic stress (6), etc. However, in natural environments, cells are likely experiencing multiple, as well as differing, stresses. We and others have previously used long-term batch culture experimental evolution to explore how evolutionary processes work in a heterogenous environment, where cells might experience some of the same environmental factors such as oxidative stress, anaerobic conditions, or low availability of certain nutrients (7–9). In these experiments, we allow cells to move through the entire life cycle of Escherichia coli in complex media, in order for cells to experience multiple types of stresses, which change throughout the incubation period.

During the E. coli life cycle in the laboratory, cells transition through the lag, log, and stationary phases and then into the death phase, where ∼99% of cells die, lyse, and release their cellular contents into the medium (10). This allows the remaining ∼1% of cells to enter long-term stationary phase (LTSP), where they can use the detritus of lysed cells as carbon and energy sources (10–12). The stresses of long-term stationary phase likely include high pH, high oxidative stress but low oxygen, and a lack of readily metabolized nutrients; however, none of these have been characterized completely (10, 13, 14). We previously showed that cells with mutations that may help cope with these stresses are selected for during this phase and that as populations continue further into LTSP, there is turnover of different mutant genotypes depending on the medium conditions at any given time (7, 9).

While we have hypothesized that this dynamism of genotypes indicates that different subpopulations of cells are growing and dying within LTSP cultures, we do not know if this phase consists of a collection of cells reflecting the four previous phases (lag, log, stationary, and/or death) or if LTSP consists of cells acting uniquely. We know that cells in lag, log, and stationary phases have unique transcriptional programs (15), but while expression changes due to genetic effects have been analyzed (16), to our knowledge, the transcriptional profile of LTSP in E. coli has not been studied. Elucidating the transcriptional profile of LTSP could help determine whether it is a phase in and of itself or whether it is more representative of an amalgam of cells experiencing the other phases.

In this study, we analyzed the transcriptome of E. coli throughout its laboratory life cycle, including death phase and LTSP. We have determined that the expression of a small set of genes is uniquely regulated during LTSP, indicating that this phase does have its own transcriptional program. Further, we hypothesize that the genes which are upregulated only in LTSP compared to other phases may be important for survival during this phase. We identify three genes in the cold shock protein family, cspB, cspF, and cspI, which are upregulated in LTSP compared to the other four phases. Loss of two of these genes (cspB and cspI) affects survival in LTSP during competition with wild-type cells. These results further indicate that (i) we may be able to use these types of data to characterize the biochemical and physiological responses to long-term stationary-phase stresses and (ii) we may be able to identify novel functions for proteins which have been characterized only in earlier phases of the life cycle. Understanding how cells respond biochemically to LTSP may allow us to determine how applicable standard laboratory conditions are as models to better understand natural environments.

RESULTS

The transcriptome of cells in LTSP can be distinguished from that of cells in other phases.In order to determine if there is a unique transcriptional program in LTSP, we incubated cells in triplicate into long-term stationary phase, extracted whole-cell RNA at six time points throughout the E. coli life cycle, and performed RNA sequencing (Fig. 1). Time points included log phase (4 h), late-log phase (8 h), stationary phase (24 h), death phase (72 h), and two time points in LTSP (144 and 192 h). We sequenced rRNA-depleted RNA from each cell population and analyzed at least 6 million reads per sample using HTSeq and DESeq2 (see Materials and Methods for a more detailed explanation) (17, 18).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Average cell density during RNA-extraction time points. Three cultures were used for RNA extraction and for counting viable cells at 4, 8, 24, 72, 144, and 192 h postinoculation. Points represent averages of results from the three replicate cultures. Error bars represent standard deviations. In some cases, the error bar is so small that the point covers it. Shading represents a new phase in the E. coli life cycle (lag, log, stationary, death, long-term stationary).

Using the gene expression data, we performed a principal-component analysis (PCA) to determine if there were features that distinguished the LTSP transcriptome from those of the other phases (Fig. 2). PCA ordinates the gene expression data across all 4,437 genes included in the analysis to identify principal-component axes (PCs) that capture the majority of the variance across all genes. We focused on the first 3 PCs as they capture 73.4% of the variation between samples (PC1 = 31.9%, PC2 = 22.8%, and PC3 = 18.7%). We can plot each sample across any combination of these axes to observe how similar each sample is to another sample with regard to the expression data and those principal components. Comparing PC1 to PC2 (Fig. 2A), we observed that cells in death phase and LTSP have similar positions along the two axes. However, cells in death phase were separated from cells in LTSP along PC3 (Fig. 2B), indicating that the expression patterns exhibited by cells in these two phases were similar with respect to PC1 and PC2 but different with respect to PC3. Cells in stationary phase were found to have a position similar to that seen with those in LTSP along PC1 and PC3 (Fig. 2B), indicating that there were similarities among their expression profiles for components 1 and 3 but that they were positioned uniquely along PC2 (Fig. 2A). Cells in log and late-log phase were positioned similarly along PC2 (as were cells in death phase; Fig. 2A) but were positioned uniquely along PC1 and PC3 (Fig. 2B). These data indicate that cells in the log and late-log phases had expression profiles that are similar to each other with respect to PC2 but different from each other with respect to PC1 and PC3. Interestingly, the expression profiles of cells in log and death phase were also similar along PC3 (Fig. 2B), indicating that there were similarities in expression patterns between these populations of cells with respect to this component. The fact that cells in LTSP had some gene expression patterns in common with those in stationary phase (PC1 and PC3) and in death phase (PC1 and PC2) supports the hypothesis that there are subpopulations of cells experiencing other phases during LTSP. However, while showing similarities in gene expression along some component axes, comparing each phase across all three axes simultaneously shows that each phase had a distinct expression profile (Fig. 2C).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

PCA of transcriptomes. Shaded circles represent 95% confidence intervals based on correlation matrices of the three replicates for each time point as follows: red, 4 h; orange, 8 h; yellow, 24 h; green, 72 h; blue, 144 h; purple, 192 h. (A) PC1 versus PC2. (B) PC1 versus PC3. (C) PC1, PC2, and PC3 plotted on the x, z, and y axes, respectively. The PCA results suggest that the cells had distinct gene expression patterns in each phase.

We also plotted a 95% confidence interval (CI) ellipse for each time point using a correlation matrix representing our three replicates (Fig. 2A and B). By comparing overlaps among ellipses for different time points, we were then able to determine how similar the transcriptomes were across cells at different time points. Where two CIs overlap, this means that, with respect to the two plotted PCs, the cells in those populations have similar expression patterns. These data indicate that cells mostly have unique transcriptional profiles in each stage of their life cycle (they stay within their own CI, which does not overlap another population’s CI) but that there may be time points throughout the life cycle when cells in one phase may also have subsets of gene expression in common with cells in other phases.

We noted that the 3 replicates for each time point before the cells reached LTSP clustered within their 95% CI, indicating that biological replicates respond very similarly to their environment even as it changes. Further, we observed that populations in LTSP are not as similar to each other as those in other phases, as these replicates did not all appear within the 95% confidence interval. We also observed that the 95% CIs for each time point overlapped only in the following cases: each LTSP time point overlapped the other LTSP time point’s CI in each PC comparison (Fig. 2A and B), and the LTSP CIs overlapped stationary-phase time points in the comparison of PC1 to PC3 (Fig. 2B).

The PC analysis indicates that cells harvested at the two LTSP time points were very similar to each other. In order to determine how similar the expression profiles were at these time points, we identified the genes that were differentially expressed (DE) between cultures incubated for 144 and 192 h. Only 13 genes in 4 operons were determined to be significantly different (false-discovery-rate [q] value < 0.05; fold change > 2) between cells at these two time points (see Table S1 in the supplemental material). Since this is a relatively small number compared to differences during other phases (for instance, 567 genes were significantly DE between log-phase and h 144 cells), we treated the cells from the two LTSP time points as one group (LTSP) for the remainder of the analyses.

TABLE S1

Genes differentially expressed between LTSP time points. Download Table S1, DOCX file, 0.01 MB.
Copyright © 2020 Kram et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Several genes are uniquely expressed in cells experiencing LTSP.While the PC analysis indicates that the cell population in LTSP was transcriptionally distinct from those in the other phases, it was not clear whether there were genes uniquely expressed in this phase or whether this distinct profile was due to a mixture of cells in one of the other four stages.

In order to address this issue, we compared the expression levels for each gene between LTSP and each of the four other phases. All of the statistically significant DE genes are represented in Table S2. A heat map of each comparison representing all genes in the E. coli genome is shown in Fig. 3. Of genes DE in LTSP in comparison to log phase, 755 were found to be expressed at higher levels and 823 at lower levels during LTSP; of genes DE in LTSP in comparison to late-log phase, 669 were expressed at higher levels and 463 at lower levels in LTSP; of genes DE in LTSP in comparison to stationary phase, 632 were expressed at higher levels and 718 at lower levels in LTSP; and of genes DE in LTSP in comparison to death phase, 311 were expressed at higher levels and 499 at lower levels in LTSP. In each comparison, with the exception of late-log phase, more genes were expressed at lower levels in LTSP than were expressed at higher levels. Further, the number of genes expressed at higher levels than populations in LTSP decreased as populations approach LTSP.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Heat map of the differences in expression between LTSP and other time points. Genes are indicated along the x axis in genome order. FC, log2 fold change.

TABLE S2

Genes differentially expressed between LTSP and at least one other phase. Download Table S2, XLSX file, 0.1 MB.
Copyright © 2020 Kram et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

While some genes were consistently DE from LTSP, there are many sets of genes that were DE in only one or two other phases. These data support the hypothesis that the population of cells in LTSP likely consists of cells experiencing one of the other phases but that LTSP may also have a unique transcriptional program that cells in some subset are expressing.

In order to identify gene expression patterns unique to LTSP, we identified 38 genes that were significantly DE in LTSP compared to all other phases (Table 1); 32 of these genes were expressed at higher levels, and 6 of these genes were expressed at lower levels.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Log2 fold change of genes differentially expressed in LTSP compared to all other phases

Some genes uniquely expressed during LTSP affect survival during this phase.We hypothesize that genes which are expressed most highly in LTSP compared to other time points may play a role in cell survival. Interestingly, 3 genes identified as expressed at higher levels in LTSP than in all other phases were csp genes: cspB, cspF, and cspI (Table 1). csp genes were initially characterized as important for physiological adaptation of cells to cold temperatures (19), which cells in our cultures did not experience. Therefore, we were especially interested in this set of genes as they may be playing a role other than in response to cold shock in long-term batch cultures. The expression patterns of all 9 known csp genes during the incubation period are shown in Fig. 4. We grouped these genes into categories I to IV. Category I included the genes mentioned above which were expressed at higher levels during LTSP than during all other phases (cspB, cspF, and cspI). Category II included the genes which were highly expressed (normalized count > 1,000) at all time points (cspC, cspD, and cspE). Category III included the genes whose expression changed throughout the time course but whose expression in LTSP was not different from that in all other phases (cspA and cspG). Category IV included the one gene which showed low expression (normalized count < 15) at all time points (cspH). While the levels of cspF and cspG expression were similar in LTSP, cspF expression was lower earlier in the life cycle, which is why it is statistically placed in category I, whereas cspG expression started higher earlier in the life cycle, which is why it is statistically grouped in category III (Fig. 4). Each csp gene was expressed from its own operon—none shared promoters, although some might share regulatory regions—and are therefore likely regulated independently of one another.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Expression pattern of csp genes. Normalized counts were outputted from DESeq2 analysis. Each point represents averages of results from the three cultures. We hypothesized that genes in category 1 would be important for survival in LTSP, genes in categories 2 and 3 might play a role in other phases of the life cycle, and the gene in category 4 would have no effect on cell growth or survival.

We predicted that all genes in category I would be important for survival in LTSP specifically, that genes in categories II and III may be important for growth or survival during other phases, and that genes in category IV would not affect cells under these culture conditions. To test these hypotheses, we deleted each gene from the PFM2 parental strain using P1 transduction from the Keio collection of gene knockouts (20) and incubated these cells into LTSP both alone in monoculture and in competition with wild-type cells. The data show that loss of genes in categories II, III, and IV had a negligible effect, if any, on growth of E. coli in competition with wild-type cells in long-term cultures (Fig. 5). Loss of cspC or cspE (both category II) might have given cells a slight (∼5-fold to 10-fold) advantage in early LTSP, but the advantage was gone 48 h later. Loss of the cspF category 1 gene also had no effect (Fig. 5F).

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Competitions (A to I) and monocultures (J) of each csp mutant versus wild-type cells. Wild-type cells are shown in black in each graph. Letters correspond to the csp gene that is missing in that mutant cell (i.e., panel A shows the competition between the wild-type and cspA::Kanr strains), with the exception of panel J, which shows monocultures of the cspB::Kanr (red), cspF::Kanr (green), and cspI::Kanr (blue) mutants and the WT strain (black). Note that panel J has a different y axis scale. For panels B, F, I, and J, the points represent averages of results from three cultures. For A, C, D, E, G, and H, the points represent averages of results from two cultures. Error bars represent standard deviations. In some cases, the error bar is so small that the point covers it.

However, loss of two of the other category I genes, cspB and cspI, did affect the cell’s ability to compete with the wild-type strain, to various degrees (Fig. 5B and I). Cells mutant for cspB grew similarly to wild-type cells and experienced a similar death phase but then did not recover to the same level as the wild-type cells at 144 h (∼80-fold lower CFU/ml). In this case, wild-type cell density decreased over the next 48 h, which brought them to the same level as the mutant cells (Fig. 5B). Cells lacking cspI again grew similarly to wild-type cells but then had a slightly more severe death phase (∼10-fold-lower CFU/ml at 72 h) and never recovered fully from death phase, leading to complete loss of viable mutant cells by 192 h (Fig. 5I). Interestingly, when cells missing each of these genes were incubated in monoculture, without having to compete with wild-type cells, they did not show any defects in survival during death phase or LTSP (Fig. 5J).

DISCUSSION

We have determined that each phase in the E. coli life cycle has a unique transcriptional profile (Fig. 2). This has previously been observed for log phase and stationary phase, but to our knowledge this is the first indication that gene expression in death phase and LTSP is significantly different from gene expression in other phases. Further, we have also shown that there are genes which are differentially expressed in LTSP compared to all other phases, indicating that this phase has a unique transcriptional program and does not represent only a combination of cells in other phases (although the populations are almost certainly also made up of cells in the other four phases, as well as of those experiencing this unique transcriptional program).

At least two of the genes that are expressed at higher levels in LTSP than in other phases, cspB and cspI, are important for survival in LTSP in competition against wild-type cells (Fig. 4 and 5). CspB and CspI are both cold shock proteins, which were initially discovered as induced at low temperatures (19, 21, 22). Generally, Csp proteins help to recover translation levels of proteins by binding to and relaxing mRNA that becomes too structured to be translated appropriately in colder temperatures, although some Csp proteins, including CspB and CspI, may be chaperones for single-stranded DNA (ssDNA) as well (19, 23, 24). Because cultures in these experiments are consistently incubated at 37°C, it is likely that these proteins are actually performing previously unidentified functions. Interestingly, both cspB and cspI are expressed in stationary phase even without cold shock (21), further supporting the idea that these genes are likely being expressed in response to stresses other than cold shock.

Other data indicate that these particular genes, as well as csp genes in general, may respond to other types of stresses. For instance, neither of these genes responded to cold stress in a pathogenic strain of E. coli (25). Further, Brandi et al. have shown that CspA, the originally identified cold shock protein, can also be regulated by the global regulators Fis and H-NS, indicating that it likely responds to stresses other than cold shock (26). cspC, cspD, and cspE are regulated by growth arrest during diauxic shift, H2O2 stress, and the transition to stationary phase (15). CspC and CspE are important for virulence in Salmonella, again indicating that they likely respond to a signal other than cold shock (27). Future studies will elucidate the roles of these and other genes during LTSP.

It is unsurprising that not all of the LTSP-specific genes identified play a specific role in LTSP—a similar finding was noted with transcriptome data in long-term cultures of Rhodopseudomonas palustris (28). However, data presented here indicate that, in fact, we can use expression patterns in LTSP to identify at least some genes that are important for survival in this environment and, further, possibly determine the growth state of cells within particular environments.

Determining which genes are expressed or repressed in lower-nutrient environments, as well as which proteins are essential for survival, will allow us to characterize the biochemical and physiological responses to stress in LTSP. We can then assess whether LTSP (or other laboratory conditions) may be a good model for natural environments, such as in a host or soil, by comparing those biochemical responses due to similar environmental factors, such as oxygen level and access to certain nutrients, to those in cells in natural environments.

MATERIALS AND METHODS

Bacterial strains, growth conditions, and viable cell counts.All strains in this study originated from strain PFM2, derived from the E. coli K-12 lineage strain, MG1655 (29) kanamycin-resistant (Kanr) and chloramphenicol-resistant (Camr) derivatives, or either mutant or “wild-type” cells, as appropriate. Before any experiment was initiated, overnight cultures were inoculated from frozen stocks into 5 ml of Luria-Bertani (Lennox) medium (LB) (Difco) in 18-mm-by-150-mm borosilicate test tubes, which were incubated with aeration in a TC-7 rolling drum (New Brunswick Scientific, Edison, NJ) at 37°C. We monitored cell growth and survival as described previously (13). Briefly, after inoculation of cells at a 1:1,000 (vol/vol) dilution from overnight growth frozen stocks into 5 ml of LB, we determined viable cell counts by serially plating dilutions of cultures on LB agar plates supplemented with kanamycin (50 μg/ml) or chloramphenicol (50 μg/ml), as appropriate.

RNA preparation and sequencing.We inoculated all cultures for transcriptome analysis from a single PFM2 overnight culture. At each time point, we sacrificed 3 cultures to remove 1 ml (4 h), 0.5 ml (8 h), or 2 ml (24, 72, 144, and 196 h) of cells, and viable cell counts were determined by serial dilution as described above. RNAprotect reagent (Qiagen) was added per the manufacturer’s instructions. Cells were pelleted and frozen at −80°C until total cellular RNA extraction was performed using an RNeasy Protect Bacteria minikit (Qiagen). Total RNA was processed, assessed for quality, and sequenced by the University of Southern California (USC) Genomics Core, where rRNA was removed using a Ribo-Zero kit (Illumina, Inc.). Single-end reads (75 bp) were generated from each sample on an Illumina HiSeq 2500 platform. We received an average of 7,882,901 reads per sample, with a range of ∼6.5 to ∼9.5 million reads per sample.

Transcriptome analysis.Using a custom Python pipeline, we used TopHat2 to align the raw reads to the E. coli MG1655 genome (30), SAMtools to generate binary alignment files (BAM) (31), and HTSeq to calculate read counts per gene and per sample (18). An average of 89.4% of reads across all samples mapped to our reference genome (Escherichia coli K-12 MG1655) (32), with a range of ∼66% to ∼96%. All rRNA gene reads, as well as reads of any genes that showed no expression (sum of gene counts in all samples < 1), were removed from the data set. We used the HTSeq output files as the input to analyze gene expression levels using DESeq2 (17), with comparisons between samples made based on the time of incubation. Once we had determined that expression differences between h 144 and h 192 were minimal, we marked these data points as the same time (“LTSP”) in our data frame. We identified differentially expressed genes using pairwise comparisons between time points. We considered genes to be differentially expressed between cells in LTSP and cells at other time points if the log2 fold change value was greater than or equal to 1 and if the q value was less than or equal to 0.05. We summarized the data using principal-component analysis (PCA) and regularized log scale (rlog) values of the normalized counts of each gene (derived from DESeq2) using the prcomp function with centering and scaling in base R. We added 95% confidence intervals by calculating correlation matrices for each the samples in each time point and then adding these intervals to our plot using the polygon function in the ellipse package in R (33). The three-dimensional (3D) PCA plot was created using the same values from prcomp, and the values were plotted using the Scatterplot3D package in R (34). We also created a heat map comparing the log2 fold change between the average normalized counts for each gene in LTSP versus each other time point using ggplot2 (35).

Mutagenesis of csp genes.We constructed in-frame knockout mutations of each csp gene individually (cspA to cspI) to create nine mutant strains by the use of P1 transduction from strains in the Keio collection (20). Briefly, we inoculated overnight cultures from a frozen stock in the Keio collection. The following day, we inoculated cultures at a 1:100 dilution into fresh medium and allowed them to grow to mid-log phase. Lysates were made from the donor cells using P1 stock. We then used the lysates to transduce wild-type PFM2 and selected for mutants on LB agar plates supplemented with kanamycin. We confirmed gene replacement using PCR.

Competitions and monocultures with wild-type and mutant strains.For monocultures, we inoculated cultures as described above. We also performed competitions between mutant strains and parental strains. Parental strains were marked with a chloramphenicol cassette replacing lacZ (lacZ::Camr) (7), and mutant strains were marked with a kanamycin cassette replacing the gene of interest (cspA::Kanr to cspI::Kanr). We inoculated wild-type cells (lacZ::Camr) and cells of each of the mutants (csp*::Kanr) from frozen stocks as described above, and overnight cultures were each coinoculated at 1:1,000 (vol/vol) into LB broth in test tubes. We assayed both monoculture and competition populations for viable counts throughout the experiment (0, 4, 8, 24, 72, 144, and 192 h after inoculation) as described above.

Data availability.The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus (36) and are accessible through GEO series accession number GSE152619.

ACKNOWLEDGMENTS

We thank Lacey L. Westphal for help with data analysis and Daniel Stoebel, Sonal Singhal, and members of the Kram lab for helpful discussion of the project and manuscript.

This work was supported in part by U.S. Army Research Office grants W911NF1010444 and W911NF1210321 (to S.E.F.) and by National Science Foundation RUI 1715066 (to K.E.K.). K.E.K. was also partially supported by the CSUDH Research, Scholarship, and Creative Activity award, the CSUDH Norris Foundation Award, and the CSUDH Emeritus Faculty Fund Award. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

FOOTNOTES

    • Received April 27, 2020.
    • Accepted July 23, 2020.
  • Copyright © 2020 Kram et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Lenski RE,
    2. Wiser MJ,
    3. Ribeck N,
    4. Blount ZD,
    5. Nahum JR,
    6. Morris JJ,
    7. Morris JJ,
    8. Zaman L,
    9. Turner CB,
    10. Wade BD,
    11. Maddamsetti R,
    12. Burmeister AR,
    13. Baird EJ,
    14. Bundy J,
    15. Grant NA,
    16. Card KJ,
    17. Rowles M,
    18. Weatherspoon K,
    19. Papoulis SE,
    20. Sullivan R,
    21. Clark C,
    22. Mulka JS,
    23. Hajela N
    . 2015. Sustained fitness gains and variability in fitness trajectories in the long-term evolution experiment with Escherichia coli. Proc R Soc Lond B Biol Sci 282:20152292. doi:10.1098/rspb.2015.2292.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Maharjan RP,
    2. Liu B,
    3. Feng L,
    4. Ferenci T,
    5. Wang L
    . 2015. Simple phenotypic sweeps hide complex genetic changes in populations. Genome Biol Evol 7:531–544. doi:10.1093/gbe/evv004.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Hindré T,
    2. Knibbe C,
    3. Beslon G,
    4. Schneider D
    . 2012. New insights into bacterial adaptation through in vivo and in silico experimental evolution. Nat Rev Microbiol 10:352–365. doi:10.1038/nrmicro2750.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Blount ZD,
    2. Borland CZ,
    3. Lenski RE
    . 2008. Historical contingency and the evolution of a key innovation in an experimental population of Escherichia coli. Proc Natl Acad Sci U S A 105:7899–7906. doi:10.1073/pnas.0803151105.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Rodríguez-Verdugo A,
    2. Tenaillon O,
    3. Gaut BS
    . 2016. First-step mutations during adaptation restore the expression of hundreds of genes. Mol Biol Evol 33:25–39. doi:10.1093/molbev/msv228.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Toprak E,
    2. Veres A,
    3. Michel J-B,
    4. Chait R,
    5. Hartl DL,
    6. Kishony R
    . 2011. Evolutionary paths to antibiotic resistance under dynamically sustained drug selection. Nat Genet 44:101–105. doi:10.1038/ng.1034.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Westphal LL,
    2. Lau J,
    3. Negro Z,
    4. Moreno IJ,
    5. Ismail Mohammed W,
    6. Lee H,
    7. Tang H,
    8. Finkel SE,
    9. Kram KE
    . 2018. Adaptation of Escherichia coli to long-term batch culture in various rich media. Res Microbiol 169:145–156. doi:10.1016/j.resmic.2018.01.003.
    OpenUrlCrossRef
  8. 8.↵
    1. Chib S,
    2. Ali F,
    3. Seshasayee ASN
    . 2017. Genomewide mutational diversity in Escherichia coli population evolving in prolonged stationary phase. mSphere 2:e00059-17. doi:10.1128/mSphere.00059-17.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Kram KE,
    2. Geiger C,
    3. Ismail WM,
    4. Lee H,
    5. Tang H,
    6. Foster PL,
    7. Finkel SE
    . 2017. Adaptation of Escherichia coli to long-term serial passage in complex medium: evidence of parallel evolution. mSystems 2:e00192-16. doi:10.1128/mSystems.00192-16.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Finkel SE
    . 2006. Long-term survival during stationary phase: evolution and the GASP phenotype. Nat Rev Microbiol 4:113–120. doi:10.1038/nrmicro1340.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Takano S,
    2. Pawlowska BJ,
    3. Gudel I,
    4. Yomo T,
    5. Tsuru S
    . 2017. Density-dependent recycling promotes the long-term survival of bacterial populations during periods of starvation. mBio 8:e02336-16. doi:10.1128/mBio.02336-16.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Schink SJ,
    2. Biselli E,
    3. Ammar C,
    4. Gerland U
    . 2019. Death rate of E. coli during starvation is set by maintenance cost and biomass recycling. Cell Syst 9:64–73. doi:10.1016/j.cels.2019.06.003.
    OpenUrlCrossRef
  13. 13.↵
    1. Kram KE,
    2. Finkel SE
    . 2015. Rich medium composition affects Escherichia coli survival, glycation, and mutation frequency during long-term batch culture. Appl Environ Microbiol 81:4442–4450. doi:10.1128/AEM.00722-15.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Kram KE,
    2. Finkel SE
    . 2014. Culture volume and vessel affect long-term survival, mutation frequency, and oxidative stress in. Appl Environ Microbiol 80:1732–1738. doi:10.1128/AEM.03150-13.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Chang DE,
    2. Smalley DJ,
    3. Conway T
    . 2002. Gene expression profiling of Escherichia coli growth transitions: an expanded stringent response model. Mol Microbiol 45:289–306. doi:10.1046/j.1365-2958.2002.03001.x.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Arunasri K,
    2. Adil M,
    3. Khan PA,
    4. Shivaji S
    . 2014. Global gene expression analysis of long-term stationary phase effects in E. coli K12 MG1655. PLoS One 9:e96701. doi:10.1371/journal.pone.0096701.
    OpenUrlCrossRef
  17. 17.↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    . 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15:550. doi:10.1186/s13059-014-0550-8.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Anders S,
    2. Pyl PT,
    3. Huber W
    . 2015. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics 31:166–169. doi:10.1093/bioinformatics/btu638.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    1. Zhang Y,
    2. Burkhardt DH,
    3. Rouskin S,
    4. Li GW,
    5. Weissman JS,
    6. Gross CA
    . 2018. A stress response that monitors and regulates mRNA structure is central to cold shock adaptation. Mol Cell 70:274–286.e7. doi:10.1016/j.molcel.2018.02.035.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Baba T,
    2. Ara T,
    3. Hasegawa M,
    4. Takai Y,
    5. Okumura Y,
    6. Baba M,
    7. Datsenko KA,
    8. Tomita M,
    9. Wanner BL,
    10. Mori H
    . 2006. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol 2:2006.0008. doi:10.1038/msb4100050.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Czapski TR,
    2. Trun N
    . 2014. Expression of csp genes in E. coli K-12 in defined rich and defined minimal media during normal growth, and after cold-shock. Gene 547:91–97. doi:10.1016/j.gene.2014.06.033.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Wang N,
    2. Yamanaka K,
    3. Inouye M
    . 1999. CspI, the ninth member of the CspA family of Escherichia coli, is induced upon cold shock. J Bacteriol 181:1603–1609. doi:10.1128/JB.181.5.1603-1609.1999.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Yamanaka K,
    2. Fang L,
    3. Inouye M
    . 1998. The CspA family in Escherichia coli: multiple gene duplication for stress adaptation. Mol Microbiol 27:247–255. doi:10.1046/j.1365-2958.1998.00683.x.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. Yu T,
    2. Keto-Timonen R,
    3. Jiang X,
    4. Virtanen JP,
    5. Korkeala H
    . 2019. Insights into the phylogeny and evolution of cold shock proteins: from enteropathogenic Yersinia and Escherichia coli to eubacteria. Int J Mol Sci 20:4059. doi:10.3390/ijms20164059.
    OpenUrlCrossRef
  25. 25.↵
    1. Li Y,
    2. Zhou D,
    3. Hu S,
    4. Xiao X,
    5. Yu Y,
    6. Li X
    . 2018. Transcriptomic analysis by RNA-seq of Escherichia coli O157:H7 response to prolonged cold stress. LWT 97:17–24. doi:10.1016/j.lwt.2018.06.025.
    OpenUrlCrossRef
  26. 26.↵
    1. Brandi A,
    2. Giangrossi M,
    3. Giuliodori AM,
    4. Falconi M
    . 2016. An interplay among FIS, H-NS, and guanosine tetraphosphate modulates transcription of the Escherichia coli cspA gene under physiological growth conditions. Front Mol Biosci 3:19. doi:10.3389/fmolb.2016.00019.
    OpenUrlCrossRef
  27. 27.↵
    1. Michaux C,
    2. Holmqvist E,
    3. Vasicek E,
    4. Sharan M,
    5. Barquist L,
    6. Westermann AJ,
    7. Gunn JS,
    8. Vogel J
    . 2017. RNA target profiles direct the discovery of virulence functions for the cold-shock proteins CspC and CspE. Proc Natl Acad Sci U S A 114:6824–6829. doi:10.1073/pnas.1620772114.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Pechter KB,
    2. Yin L,
    3. Oda Y,
    4. Gallagher L,
    5. Yang J,
    6. Manoil C,
    7. Harwood CS
    . 2017. Molecular basis of bacterial longevity. mBio 8:e01726-17. doi:10.1128/mBio.01726-17.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    1. Foster PL,
    2. Lee H,
    3. Popodi E,
    4. Townes JP,
    5. Tang H
    . 2015. Determinants of spontaneous mutation in the bacterium Escherichia coli as revealed by whole-genome sequencing. Proc Natl Acad Sci U S A 112:E5990–E5999. doi:10.1073/pnas.1512136112.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Kim D,
    2. Pertea G,
    3. Trapnell C,
    4. Pimentel H,
    5. Kelley R,
    6. Salzberg SL
    . 2013. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14:R36. doi:10.1186/gb-2013-14-4-r36.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Li H,
    2. Handsaker B,
    3. Wysoker A,
    4. Fennell T,
    5. Ruan J,
    6. Homer N,
    7. Marth G,
    8. Abecasis G,
    9. Durbin R
    , 1000 Genome Project Data Processing Subgroup. 2009. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079. doi:10.1093/bioinformatics/btp352.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Riley M,
    2. Abe T,
    3. Arnaud MB,
    4. Berlyn MKB,
    5. Blattner FR,
    6. Chaudhuri RR,
    7. Glasner JD,
    8. Horiuchi T,
    9. Keseler IM,
    10. Kosuge T,
    11. Mori H,
    12. Perna NT,
    13. Iii GP,
    14. Rudd KE,
    15. Serres MH,
    16. Thomas GH,
    17. Thomson NR,
    18. Wishart D,
    19. Wanner BL
    . 2006. Escherichia coli K-12: a cooperatively developed annotation snapshot–2005. Nucleic Acids Res 34:1–9. doi:10.1093/nar/gkj405.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Murdoch DJ,
    2. Chow ED
    . 1996. A graphical display of large correlation matrices. Am Stat 50:178–180. doi:10.2307/2684435.
    OpenUrlCrossRefWeb of Science
  34. 34.↵
    1. Ligges U,
    2. Mächler M
    . 2003. Scatterplot3d - an R package for visualizing multivariate data. J Stat Soft 8:1–20. doi:10.18637/jss.v008.i11.
    OpenUrlCrossRef
  35. 35.↵
    1. Wickham H
    . 2016. ggplot2: elegant graphics for data analysis. Springer-Verlag, New York, NY.
  36. 36.↵
    1. Edgar R,
    2. Domrachev M,
    3. Lash AE
    . 2002. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 30:207–210. doi:10.1093/nar/30.1.207.
    OpenUrlCrossRefPubMedWeb of Science
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Escherichia coli Has a Unique Transcriptional Program in Long-Term Stationary Phase Allowing Identification of Genes Important for Survival
Karin E. Kram, Autumn L. Henderson, Steven E. Finkel
mSystems Aug 2020, 5 (4) e00364-20; DOI: 10.1128/mSystems.00364-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Escherichia coli Has a Unique Transcriptional Program in Long-Term Stationary Phase Allowing Identification of Genes Important for Survival
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Escherichia coli Has a Unique Transcriptional Program in Long-Term Stationary Phase Allowing Identification of Genes Important for Survival
Karin E. Kram, Autumn L. Henderson, Steven E. Finkel
mSystems Aug 2020, 5 (4) e00364-20; DOI: 10.1128/mSystems.00364-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

cold shock
experimental evolution
genetics
long-term stationary phase
transcriptomics

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077