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.

The transcription factor network of E. coli steers global responses to shifts in RNAP concentration.

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

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.

A library of reporters of the global regulators of gene expression in Escherichia coli.

The topology of the transcription factor network (TFN) of Escherichia coli is far from uniform, with 22 global regulator (GR) proteins controlling one-third of all genes. So far, their production rates cannot be tracked by comparable fluorescent proteins. We developed a library of fluorescent reporters for 16 GRs for this purpose. Each consists of a single-copy plasmid coding for green fluorescent protein (GFP) fused to the full-length copy of the native promoter. We tracked their activity in exponential and stationary growth, as well as under weak and strong stresses. We show that the reporters have high sensitivity and specificity to all stresses tested and detect single-cell variability in transcription rates. Given the influence of GRs on the TFN, we expect that the new library will contribute to dissecting global transcriptional stress-response programs of E. coli. Moreover, the library can be invaluable in bioindustrial applications that tune those programs to, instead of cell growth, favor productivity while reducing energy consumption.IMPORTANCECells contain thousands of genes. Many genes are involved in the control of cellular activities. Some activities require a few hundred genes to run largely synchronous transcriptional programs. To achieve this, cells have evolved global regulator (GR) proteins that can influence hundreds of genes simultaneously. We have engineered a library of Escherichia coli strains to track the levels over time of these, phenotypically critical, GRs. Each strain has a single-copy plasmid coding for a fast-maturing green fluorescent protein whose transcription is controlled by a copy of the natural GR promoter. By allowing the tracking of GR levels, with sensitivity and specificity, this library should become of wide use in scientific research on bacterial gene expression (from molecular to synthetic biology) and, later, be used in applications in therapeutics and bioindustries.

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.

Dissecting the in vivo dynamics of transcription locking due to positive supercoiling buildup.

Positive supercoiling buildup (PSB) is a pervasive phenomenon in the transcriptional programs of Escherichia coli. After finding a range of Gyrase concentrations where the inverse of the transcription rate of a chromosome-integrated gene changes linearly with the inverse of Gyrase concentration, we apply a LineWeaver-Burk plot to dissect the expected in vivo transcription rate in absence of PSB. We validate the estimation by time-lapse microscopy of single-RNA production kinetics of the same gene when single-copy plasmid-borne, shown to be impervious to Gyrase inhibition. Next, we estimate the fraction of time in locked states and number of transcription events prior to locking, which we validate by measurements under Gyrase inhibition. Replacing the gene of interest by one with slower transcription rate decreases the fraction of time in locked states due to PSB. Finally, we combine data from both constructs to infer a range of possible transcription initiation locking kinetics in a chromosomal location, obtainable by tuning the transcription rate. We validate with measurements of transcription activity at different induction levels. This strategy for dissecting transcription initiation locking kinetics due to PSB can contribute to resolve the transcriptional programs of E. coli and in the engineering of synthetic genetic circuits.