Abstract
Escherichia coli cells respond to a period of famine by globally reorganizing their gene expression. The changes are known as the stringent response, which is orchestrated by the alarmone ppGpp that binds directly to RNA polymerase. The resulting changes in gene expression are particularly well studied in the case of amino acid starvation. We used deep RNA sequencing in combination with spike-in cells to measure global changes in the transcriptome after valine-induced isoleucine starvation of a standard E. coli K12 strain. Owing to the whole-cell spike-in method that eliminates variations in RNA extraction efficiency between samples, we show that ribosomal RNA levels are reduced during isoleucine starvation and we quantify how the change in cellular RNA content affects estimates of gene regulation. Specifically, we show that standard data normalization relying on sample sequencing depth underestimates the number of down-regulated genes in the stringent response and overestimates the number of up-regulated genes by approximately 40%. The whole-cell spike-in method also made it possible to quantify how rapidly the pool of total messenger RNA (mRNA) decreases upon amino acid starvation. A principal component analysis showed that the first two components together described 69% of the variability of the data, underlining that large and highly coordinated regulons are at play in the stringent response. The induction of starvation by sudden addition of high valine concentrations provoked prominent regulatory responses outside of the expected ppGpp, RpoS, and Lrp regulons. This underlines the notion that with the high resolution possible in deep RNA sequencing analysis, any different starvation method (e.g., nitrogen-deprivation, removal of an amino acid from an auxotroph strain, or valine addition to E. coli K12 strains) will produce measurable variations in the stress response produced by the cells to cope with the specific treatment.