Effects of σ factor competition are promoter initiation kinetics dependent

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.

Construction of Strong Promoters by Assembling Sigma Factor Binding Motifs

Development of strong promoters is of growing interest in the field of biotechnology and synthetic biology. Here we present a protocol for the construction of strong prokaryotic promoters that can be recognized by designated multiple sigma factors by interlocking their cognate binding motifs on DNA strands. Strong and stress responsive promoters for Escherichia coli and Bacillus subtilis have been created following the presented protocol. Customized promoters could be easily developed for fine-tuning gene expression or overproducing enzymes with prokaryotic cell factories.

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 ‘d_TSS’, 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 d_TSS are the key regulators of transcription interference due to TSS occlusion by RNAP at one of the promoters (when d_TSS ≤ 35 bp) and RNAP occupancy of the downstream promoter (when d_TSS > 35 bp). Occlusion and downstream promoter occupancy are modeled as linear functions of occupancy time, while the influence of d_TSS 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 d_TSS. 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.

Tandem promoters are common in nature, but investigations on their dynamics have so far largely relied on synthetic constructs. Thus, their regulation and potentially unique dynamics remain unexplored. We first performed a comprehensive exploration of the conservation of genes regulated by these promoters in E. coli and the properties of their input transcription factors. We then measured protein and RNA levels expressed by 30 Escherichia coli tandem promoters, to establish an analytical model of the expression dynamics of genes controlled by such promoters. We show that start site occlusion and downstream RNAP occupancy can be realistically captured by a model with RNAP binding affinity, the time length of open complex formation, and the nucleotide distance between transcription start sites. This study contributes to a better understanding of the unique dynamics tandem promoters can bring to the dynamics of gene networks and will assist in their use in synthetic genetic circuits.

Involvement of three FliA-family sigma factors in the sporangium formation, spore dormancy and sporangium dehiscence in Actinoplanes missouriensis

The rare actinomycete Actinoplanes missouriensis forms sporangia, which open up and release zoospores in response to water. Here, we report a genetic and functional analysis of four FliA-family sigma factors, FliA1, FliA2, FliA3 and FliA4. Transcription of fliA1, fliA2 and fliA3 was directly activated by the global transcriptional activator TcrA during sporangium formation and dehiscence, while fliA4 was almost always transcribed at low levels. Gene disruption analysis showed that (a) deletion of fliA2 reduced the zoospore swimming speed by half, (b) the fliA1-fliA2 double-deletion mutant formed abnormal sporangia in which mutant spores ectopically germinated and (c) deletion of fliA3 induced no phenotypic changes in the wild-type and mutant strains of fliA1 and/or fliA2. Comparative RNA-Seq analyses among the wild-type and gene deletion mutant strains showed probable targets of each FliA-family sigma factor, indicating that FliA1- and FliA2-dependent promoters are quite similar to each other, while the FliA3-dependent promoter is somewhat different. Gene complementation experiments also indicated that the FliA1 regulon overlaps with the FliA2 regulon. These results demonstrate that A. missouriensis has developed a complex transcriptional regulatory network involving multiple FliA-family sigma factors for the accomplishment of its characteristic reproduction process, including sporangium formation, spore dormancy and sporangium dehiscence.

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.

Transcription closed and open complex formation coordinate expression of genes with a shared promoter region

Many genes are spaced closely, allowing coordination without explicit control through shared regulatory elements and molecular interactions. We study the dynamics of a stochastic model of a gene-pair in a head-to-head configuration, sharing promoter elements, which accounts for the rate-limiting steps in transcription initiation. We find that only in specific regions of the parameter space of the rate-limiting steps is orderly coexpression exhibited, suggesting that successful cooperation between closely spaced genes requires the coevolution of compatible rate-limiting step configuration. The model predictions are validated using in vivo single-cell, single-RNA measurements of the dynamics of pairs of genes sharing promoter elements. Our results suggest that, in E. coli, the kinetics of the rate-limiting steps in active transcription can play a central role in shaping the dynamics of gene-pairs sharing promoter elements.

Estimating RNA numbers in single cells by RNA fluorescent tagging and flow cytometry

Estimating the statistics of single-cell RNA numbers has become a key source of information on gene expression dynamics. One of the most informative methods of in vivo single-RNA detection is MS2d-GFP tagging. So far, it requires microscopy and laborious semi-manual image analysis, which hampers the amount of collectable data. To overcome this limitation, we present a new methodology for quantifying the mean, standard deviation, and skewness of single-cell distributions of RNA numbers, from flow cytometry data on cells expressing RNA tagged with MS2d-GFP. The quantification method, based on scaling flow-cytometry data from microscopy single-cell data on integer-valued RNA numbers, is shown to readily produce precise, big data on in vivo single-cell distributions of RNA numbers and, thus, can assist in studies of transcription dynamics.

Diversity, versatility and complexity of bacterial gene regulation mechanisms: opportunities and drawbacks for applications in synthetic biology

Gene expression occurs in two essential steps: transcription and translation. In bacteria, the two processes are tightly coupled in time and space, and highly regulated. Tight regulation of gene expression is crucial. It limits wasteful consumption of resources and energy, prevents accumulation of potentially growth inhibiting reaction intermediates, and sustains the fitness and potential virulence of the organism in a fluctuating, competitive and frequently stressful environment. Since the onset of studies on regulation of enzyme synthesis, numerous distinct regulatory mechanisms modulating transcription and/or translation have been discovered. Mostly, various regulatory mechanisms operating at different levels in the flow of genetic information are used in combination to control and modulate the expression of a single gene or operon. Here, we provide an extensive overview of the very diverse and versatile bacterial gene regulatory mechanisms with major emphasis on their combined occurrence, intricate intertwinement and versatility. Furthermore, we discuss the potential of well-characterized basal expression and regulatory elements in synthetic biology applications, where they may ensure orthogonal, predictable and tunable expression of (heterologous) target genes and pathways, aiming at a minimal burden for the host.

Chromosome and plasmid-borne PLacO3O1 promoters differ in sensitivity to critically low temperatures

Temperature shifts trigger genome-wide changes in Escherichia coli's gene expression. We studied if chromosome integration impacts on a gene's sensitivity to these shifts, by comparing the single-RNA production kinetics of a PLacO3O1 promoter, when chromosomally-integrated and when single-copy plasmid-borne. At suboptimal temperatures their induction range, fold change, and response to decreasing temperatures are similar. At critically low temperatures, the chromosome-integrated promoter becomes weaker and noisier. Dissection of its initiation kinetics reveals longer lasting states preceding open complex formation, suggesting enhanced supercoiling buildup. Measurements with Gyrase and Topoisomerase I inhibitors suggest hindrance to escape supercoiling buildup at low temperatures. Consistently, similar phenomena occur in energy-depleted cells by DNP at 30 °C. Transient, critically-low temperatures have no long-term consequences, as raising temperature quickly restores transcription rates. We conclude that the chromosomally-integrated PLacO3O1 has higher sensitivity to low temperatures, due to longer-lasting super-coiled states. A lesser active, chromosome-integrated native lac is shown to be insensitive to Gyrase overexpression, even at critically low temperatures, indicating that the rate of escaping positive supercoiling buildup is temperature and transcription rate dependent. A genome-wide analysis supports this, since cold-shock genes exhibit atypical supercoiling-sensitivities. This phenomenon might partially explain the temperature-sensitivity of some transcriptional programs of E. coli.

Added value of autoregulation and multi-step kinetics of transcription initiation

Bacterial gene expression regulation occurs mostly during transcription, which has two main rate-limiting steps: the close complex formation, when the RNA polymerase binds to an active promoter, and the subsequent open complex formation, after which it follows elongation. Tuning these steps' kinetics by the action of e.g. transcription factors, allows for a wide diversity of dynamics. For example, adding autoregulation generates single-gene circuits able to perform more complex tasks. Using stochastic models of transcription kinetics with empirically validated parameter values, we investigate how autoregulation and the multi-step transcription initiation kinetics of single-gene autoregulated circuits can be combined to fine-tune steady state mean and cell-to-cell variability in protein expression levels, as well as response times. Next, we investigate how they can be jointly tuned to control complex behaviours, namely, time counting, switching dynamics and memory storage. Overall, our finding suggests that, in bacteria, jointly regulating a single-gene circuit's topology and the transcription initiation multi-step dynamics allows enhancing complex task performance.

Temperature-dependence of the single-cell variability in the kinetics of transcription activation in Escherichia coli

From in vivo single-cell, single-RNA measurements of the activation times and subsequent steady-state active transcription kinetics of a single-copy Lac-ara-1 promoter in Escherichia coli, we characterize the intake kinetics of the inducer (IPTG) from the media, following temperature shifts. For this, for temperature shifts of various degrees, we obtain the distributions of transcription activation times as well as the distributions of intervals between consecutive RNA productions following activation in individual cells. We then propose a novel methodology that makes use of deconvolution techniques to extract the mean and the variability of the distribution of intake times. We find that cells, following shifts to low temperatures, have higher intake times, although, counter-intuitively, the cell-to-cell variability of these times is lower. We validate the results using a new methodology for direct estimation of mean intake times from measurements of activation times at various inducer concentrations. The results confirm that E. coli's inducer intake times from the environment are significantly higher following a shift to a sub-optimal temperature. Finally, we provide evidence that this is likely due to the emergence of additional rate-limiting steps in the intake process at low temperatures, explaining the reduced cell-to-cell variability in intake times.

Rate-limiting steps in transcription dictate sensitivity to variability in cellular components

Cell-to-cell variability in cellular components generates cell-to-cell diversity in RNA and protein production dynamics. As these components are inherited, this should also cause lineage-to-lineage variability in these dynamics. We conjectured that these effects on transcription are promoter initiation kinetics dependent. To test this, first we used stochastic models to predict that variability in the numbers of molecules involved in upstream processes, such as the intake of inducers from the environment, acts only as a transient source of variability in RNA production numbers, while variability in the numbers of a molecular species controlling transcription of an active promoter acts as a constant source. Next, from single-cell, single-RNA level time-lapse microscopy of independent lineages of Escherichia coli cells, we demonstrate the existence of lineage-to-lineage variability in gene activation times and mean RNA production rates, and that these variabilities differ between promoters and inducers used. Finally, we provide evidence that this can be explained by differences in the kinetics of the rate-limiting steps in transcription between promoters and induction schemes. We conclude that cell-to-cell and consequent lineage-to-lineage variability in RNA and protein numbers are both promoter sequence-dependent and subject to regulation.

Two Distinct Mechanisms Govern RpoS-Mediated Repression of Tick-Phase Genes during Mammalian Host Adaptation by Borrelia burgdorferi, the Lyme Disease Spirochete

The alternative sigma factor RpoS plays a key role modulating gene expression in Borrelia burgdorferi, the Lyme disease spirochete, by transcribing mammalian host-phase genes and repressing σ⁷⁰-dependent genes required within the arthropod vector. To identify cis regulatory elements involved in RpoS-dependent repression, we analyzed green fluorescent protein (GFP) transcriptional reporters containing portions of the upstream regions of the prototypical tick-phase genes ospAB, the glp operon, and bba74. As RpoS-mediated repression occurs only following mammalian host adaptation, strains containing the reporters were grown in dialysis membrane chambers (DMCs) implanted into the peritoneal cavities of rats. Wild-type spirochetes harboring ospAB- and glp-gfp constructs containing only the minimal (−35/−10) σ⁷⁰ promoter elements had significantly lower expression in DMCs relative to growth in vitro at 37°C; no reduction in expression occurred in a DMC-cultivated RpoS mutant harboring these constructs. In contrast, RpoS-mediated repression of bba74 required a stretch of DNA located between −165 and −82 relative to its transcriptional start site. Electrophoretic mobility shift assays employing extracts of DMC-cultivated B. burgdorferi produced a gel shift, whereas extracts from RpoS mutant spirochetes did not. Collectively, these data demonstrate that RpoS-mediated repression of tick-phase borrelial genes occurs by at least two distinct mechanisms. One (e.g., ospAB and the glp operon) involves primarily sequence elements near the core promoter, while the other (e.g., bba74) involves an RpoS-induced transacting repressor. Our results provide a genetic framework for further dissection of the essential “gatekeeper” role of RpoS throughout the B. burgdorferi enzootic cycle.

Borrelia burgdorferi, the Lyme disease spirochete, modulates gene expression to adapt to the distinctive environments of its mammalian host and arthropod vector during its enzootic cycle. The alternative sigma factor RpoS has been referred to as a “gatekeeper” due to its central role in regulating the reciprocal expression of mammalian host- and tick-phase genes. While RpoS-dependent transcription has been studied extensively, little is known regarding the mechanism(s) of RpoS-mediated repression. We employed a combination of green fluorescent protein transcriptional reporters along with an in vivo model to define cis regulatory sequences responsible for RpoS-mediated repression of prototypical tick-phase genes. Repression of ospAB and the glp operon requires only sequences near their core promoters, whereas modulation of bba74 expression involves a putative RpoS-dependent repressor that binds upstream of the core promoter. Thus, Lyme disease spirochetes employ at least two different RpoS-dependent mechanisms to repress tick-phase genes within the mammal.