Alteration of DNA supercoiling serves as a trigger of short-term cold shock repressed genes of E. coli.

Cold shock adaptability is a key survival skill of gut bacteria of warm-blooded animals. Escherichia coli cold shock responses are controlled by a complex multi-gene, timely-ordered transcriptional program. We investigated its underlying mechanisms. Having identified short-term, cold shock repressed genes, we show that their responsiveness is unrelated to their transcription factors or global regulators, while their single-cell protein numbers' variability increases after cold shock. We hypothesized that some cold shock repressed genes could be triggered by high propensity for transcription locking due to changes in DNA supercoiling (likely due to DNA relaxation caused by an overall reduction in negative supercoiling). Concomitantly, we found that nearly half of cold shock repressed genes are also highly responsive to gyrase inhibition (albeit most genes responsive to gyrase inhibition are not cold shock responsive). Further, their response strengths to cold shock and gyrase inhibition correlate. Meanwhile, under cold shock, nucleoid density increases, and gyrases and nucleoid become more colocalized. Moreover, the cellular energy decreases, which may hinder positive supercoils resolution. Overall, we conclude that sensitivity to diminished negative supercoiling is a core feature of E. coli's short-term, cold shock transcriptional program, and could be used to regulate the temperature sensitivity of synthetic circuits.

Sequence-dependent model of genes with dual σ factor preference.

Escherichia coli uses σ factors to quickly control large gene cohorts during stress conditions. While most of its genes respond to a single σ factor, approximately 5% of them have dual σ factor preference. The most common are those responsive to both σ70, which controls housekeeping genes, and σ38, which activates genes during stationary growth and stresses. Using RNA-seq and flow-cytometry measurements, we show that 'σ70+38 genes' are nearly as upregulated in stationary growth as 'σ38 genes'. Moreover, we find a clear quantitative relationship between their promoter sequence and their response strength to changes in σ38 levels. We then propose and validate a sequence dependent model of σ70+38 genes, with dual sensitivity to σ38 and σ70, that is applicable in the exponential and stationary growth phases, as well in the transient period in between. We further propose a general model, applicable to other stresses and σ factor combinations. Given this, promoters controlling σ70+38 genes (and variants) could become important building blocks of synthetic circuits with predictable, sequence-dependent sensitivity to transitions between the exponential and stationary growth phases.

Analytical kinetic model of native tandem promoters in E. coli.

Closely spaced promoters in tandem formation are abundant in bacteria. We investigated the evolutionary conservation, biological functions, and the RNA and single-cell protein expression of genes regulated by tandem promoters in E. coli. We also studied the sequence (distance between transcription start sites 'dTSS', pause sequences, and distances from oriC) and potential influence of the input transcription factors of these promoters. From this, we propose an analytical model of gene expression based on measured expression dynamics, where RNAP-promoter occupancy times and dTSS are the key regulators of transcription interference due to TSS occlusion by RNAP at one of the promoters (when dTSS ≤ 35 bp) and RNAP occupancy of the downstream promoter (when dTSS > 35 bp). Occlusion and downstream promoter occupancy are modeled as linear functions of occupancy time, while the influence of dTSS is implemented by a continuous step function, fit to in vivo data on mean single-cell protein numbers of 30 natural genes controlled by tandem promoters. The best-fitting step is at 35 bp, matching the length of DNA occupied by RNAP in the open complex formation. This model accurately predicts the squared coefficient of variation and skewness of the natural single-cell protein numbers as a function of dTSS. Additional predictions suggest that promoters in tandem formation can cover a wide range of transcription dynamics within realistic intervals of parameter values. By accurately capturing the dynamics of these promoters, this model can be helpful to predict the dynamics of new promoters and contribute to the expansion of the repertoire of expression dynamics available to synthetic genetic constructs.