The CODOs ("
coherent
domains of gene expression" denoted formerly as “functional domains”)
[7] emerge as growth condition-dependent, transient, spatially extended gene expression patterns. These spatially confined patterns are thought to manifest the structural-functional organization of the bacterial genome. In
D. dadantii, the CODOs comprising thermodynamically variable DNA sequences are associated with distinct genetic traits; in other words, the CODOs, which emerge in various constellations depending on the applied stress, integrate the chromosomal transcriptional response to the stress-induced peculiar changes of supercoiling with the expression of stress-specific virulence and adaptation genes
[15][16][17]. This coupling of the DNA physicochemical properties and the supercoiling response with particular genetic function within the CODOs, underscores once more the role of genomic sequence organization (i.e., the spatially ordered distribution of DNA thermodynamic stability) in coordinating the global transcriptional response. Also, in
Streptococcus pneumoniae the sequence composition of chromosomal domains harboring the distinct adaptation and virulence traits was found to be determinative for peculiar supercoiling response
[18][19], again underscoring the role of genomic sequence organization in coordinating the bacterial genetic response. In
E. coli, this strategic coupling of the supercoiling response and the genetic function in genomic sequence organization is made conspicuous in the selective enrichment of the relatively G/C-rich OriC and relatively A/T-rich Ter chromosomal poles for anabolic and catabolic genes, respectively
[20][21]. More compellingly, the bacterial genomic sequence organization reveals how during the growth cycle the environmentally determined availability of superhelical energy can be fittingly coupled to a genetic adaptive response. While in
E. coli this directional supercoiling response is encoded in and enabled by the peculiar genomic sequence organization along the OriC-Ter axis, in other examined bacteria (especially in pathogens having a relatively complex lifestyle) the genomic organization patterns may be more nuanced yet represent variations on a common theme
[16][17]. As mentioned above, in general, both the local and global transcriptional responses to supercoiling vary, with relatively G/C-rich and G/C-poor sequences responding to high and low levels of negative superhelicity, respectively
[22][23].
3. Coupling of DNA Topology, NAP Binding Effects, and Holoenzyme Sigma-Factor Composition: Major Regulatory Events during the E. coli Growth Cycle
On nutritional shift-up, the expression of the
gyrA and
gyrB genes encoding DNA gyrase—an enzyme introducing negative supercoils into the DNA in an ATP-dependent manner,—increases
[30][31] concomitantly with an increase in the ATP/ADP ratio
[32][33][34]. It is likely that at the chromosomal OriC pole negative superhelicity rapidly attains high densities of σ~−0.07 to −0.08, consistent with both the enrichment of the OriC pole for the gyrase binding sites and the preferential binding of gyrase downstream of the strongly transcribed genes organized predominantly around OriC including the exceptionally strong rRNA operons
[35][36]. Increased negative superhelicity both facilitates replication initiation and strongly activates the OriC-proximal
fis gene expression
[37][38][39], whereas the accumulation of FIS in turn maintains activated ribosomal RNA transcription
[40]. Structurally, the increase in negative superhelicity leads to a branching of plectonemically coiled DNA, thus multiplying the tightly bent apical loops and facilitating the wrapping of DNA by the RNAP σ
70 holoenzyme
[41][42][43][44], which prefers highly supercoiled templates for transcription
[45]. Notably, the RNAP σ
70 holoenzyme demonstrates an OriC-Ter gradient of binding site frequency distribution correlating with that of the gyrase binding sites
[1]. Concomitantly with activation of the
fis gene, the expression of the
rpoZ gene encoding the RpoZ (ω)subunit of RNA polymerase is also strongly increased
[20][31]. RpoZ stabilizes the polymerase σ
70 holoenzyme assembly
[46] and is also directly involved in mediating the response to the "alarmone" nucleotide ppGpp
[47][48]. FIS activates the
hupA gene
[49] such that at the early growth stage, the high levels of negative superhelicity and the RNAP σ
70 holoenzyme coexist with increased levels of the two major “early” NAPs—FIS and HUα. High FIS levels boost rRNA synthesis, whereas HU supports the organization of rRNA operons in “transcription factories” engaging hundreds of RNAP molecules
[35][50]. The latter effect requires high levels of negative superhelicity
[51]. Accordingly, both FIS and HU stabilize negative supercoils
[52][53][54].
On transition to the stationary phase associated with a shortage of resources, the rapid increase in ppGpp levels (ppGpp spike) decreases both rRNA transcription and the supercoil density around OriC, thus precluding replication initiation
[55]. The
gyrA/
B gene expression and the gyrase levels subside
[30] and the total superhelicity drops to σ~−0.03. The down-regulation of the chromosomal OriC pole at this stage is likely supported by the binding of the global repressor H-NS and the “late” NAPs such as Lrp and Dps
[1][56], whereas the concomitant activation of the Ter pole is attributable to the accumulation of the RNAPσ
S holoenzyme and IHF, the binding sites for both of which are strongly enriched around the chromosomal Ter pole
[1]. Notably, the “late” NAPs as well as the RNAPσ
S holoenzyme, preferentially bind relaxed DNA
[20]. Thus, in the stationary phase, the RNAPσ
S holoenzyme responsible for the maintenance function predominates, whereas the globally relaxed DNA is organized by the binding effects of three abundant NAPs including IHF, which stabilizes planar bends without the substantial constraints of negative superhelicity, Lrp, which constrains positive supercoils, and Dps, which packages the relaxed DNA in a protective crystalline lattice
[57][58].
What happens to the nucleoids of cells enriched for these abundant stationary-phase NAPs on nutritional shift-up? After nutritional shift-up, the compaction state of the nucleoid undergoes marked alterations adopting a more open conformation
[59], presumably due to competition between the early and late growth stage NAPs, respectively, FIS and Dps
[56][60]. The DNA binding of IHF and Dps was shown to depend on environmental factors
[61]. Under conditions of the increase in DNA negative superhelicity after nutritional shift-up, both IHF and Dps likely start to dissociate from genomic binding sites, perhaps accelerated in part by “facilitated dissociation”—a general mechanism thought to modulate gene expression by assisting in the local removal of DNA binding proteins from cognate sites
[62][63][64]. However, IHF remains stably bound at the chromosomal origin of replication, its dissociation being prevented by the replication initiator protein DnaA cooperating with IHF in initiating chromosomal replication
[65]. This latter process is facilitated by FIS protein produced at high levels on nutritional shift-up
[66].
4. Genomic Transcription during the Bacterial Growth Cycle Is Steered by Supercoil Energy
Changes in DNA supercoiling are assumed to mediate the transmission of environmental changes to the chromosome, responding not only to the altered availability of metabolic resources during the growth cycle
[67] but also to the various stress factors including suboptimal oxygen tension, temperature, and osmolarity
[32][68][69][70]. This mediation by supercoiling is due to the capacity of the DNA double helix to sense the environmental conditions and to respond by adjusting accordingly the superhelical density on the one hand, and by channeling the available superhelical energy into corresponding gene expression patterns on the other. During the bacterial growth cycle, changes in DNA topology modulate the binding of DNA architectural proteins and the activity of transcription machinery concertedly, resulting in orchestrated directional alterations of chromatin architecture and gene expression. Several lines of evidence are consistent with this notion.
First, the continuous alteration of DNA superhelical density during the growth cycle is dependent on the availability of nutritional resources
[67]. As in the progression of the population along the growth cycle, the superhelical density decreases, the tightness of apical loops and DNA interwindings also decreases, and so the accessibility to the RNA polymerase changes. As mentioned above, from the early exponential phase to the early stationary phase dominated by the RNAPσ
70 holoenzyme, the number of total superhelical turns per 1 kb would reduce from ~8 to ~3, with a reduction in unconstrained turns likely to be greater and hence, the available DNA superhelical torque lower. Under these conditions, the RNAPσ
S holoenzyme preferring relaxed DNA substrates for transcription becomes active
[45][71]. At the early growth stage the RpoZ-dependent stabilization of the RNAPσ
70 holoenzyme facilitates the utilization of negatively supercoiled templates, whereas the lack of RpoZ leads to a global DNA relaxation and an increased activity of the RNAPσ
S holoenzyme, switching the global transcription preferences to the utilization of the relaxed DNA templates
[72]. Importantly, the impact of accumulated σ
S could not be detected until the superhelical density of the reporter plasmids subsided to relatively low levels (∆Lk of ~3–4 between the late exponential phase and the early stationary phase, where the σ
S impact was observed). It was thus inferred that the variations in DNA supercoiling as a function of the growth phase act as a checkpoint, precluding the shift from the RNAPσ
70 to the RNAPσ
S transcriptional machinery until the growth conditions become unfavorable enough to cause entry into the stationary phase
[73].
Second, the superhelical density and the supercoiling response of the genomic transcription change coordinately both during the growth cycle of
E. coli [74] and during the circadian cycle of gene expression in cyanobacteria. A combination of topological analyses with transcriptomics data in
Synechococcus elongatus suggested that each topological state corresponded to a unique state of gene expression, indicating that supercoiling plays a primary role in regulating circadian gene expression
[75]. Notably, the DNA sequence characteristics of genes monotonically activated and repressed by chromosomal relaxation during the circadian cycle were similar to those of the supercoiling-responsive genes in
E. coli.
Third, the sequential activation of sets of primary and downstream regulatory genes was observed in response to the long-term supercoiling imbalance achieved by modulating the
topA gene expression in
Streptomyces coelicolor, whereby increased negative superhelicity modified the levels of topoisomerases and NAPs coordinately
[76]. Another relevant observation comes from experimental evolution studies, which identified mutations in genes encoding the NAPs and topoisomerases. These mutations induced inheritable adaptive changes of supercoiling and also provided fitness gains, thus revealing the pivotal role of NAPs and topoisomerases in organizing the global transcription program during adaptation
[77][78].
Fourth, supercoiling can impose directionality by rendering the structural transitions in DNA both deterministic
[79] and coordinated
[80][81][82][83], whereas the nucleoprotein structures stabilized by NAPs have been implicated in the directional channeling of torsional energy toward the transcription initiation sites
[84][85][86]. Available data indicate that the NAP-dependent alterations of gene expression during the growth cycle involve directional and coordinated transitions in the composition of regulatory nucleoprotein structures associated with gene promoter regions
[87][88][89][90]. These directional transitions of nucleoprotein complexes associated with individual gene promoters are paralleled by the coordinated redistribution of transcription machinery during the growth cycle at the global scale of the entire chromosome
[51][91].
Finally, it is noteworthy that the regulation of gene transcription by supercoil energy is also modulated by gene organization in the genome. The processes of DNA supercoiling and transcription are interdependent
[92] and in addition to the specific regulation of gene promoters, global transcriptional responses to changes in DNA supercoiling depend on constraints imposed by the local orientation of genes and the supercoil diffusion induced by the transcription process itself
[93][94][95][96]. Ongoing transcription generates negative and positive supercoils, respectively, upstream and downstream of the translocating RNAP and this transcription-coupled diffusion of supercoils (TCDS) modulates the activity of neighbouring gene promoters. The TCDS effect exerted on a particular gene depends on the mutual orientation (either convergent, divergent, or tandem) of its surrounding transcription units. For example, a gene embedded between two divergently oriented transcription units will experience high levels of negative superhelicity (negative TCDS), and vice versa in the case of flanking convergent units, high positive superhelicity (positive TCDS). Although the transcriptional regulation shaped by local genomic architecture can be modulated by changes in global supercoiling and by NAPs
[97], these local constraints provide additional means for fine-tuning the supercoiling-dependent impacts on gene expression during the growth cycle. In particular, the measurement of TCDS during the growth cycle of
E. coli shows that the genes responding to high negative superhelicity (
hyp genes) also experience negative TCDS from their neighbors, whereas the genes responding to DNA relaxation (
rel genes) experience positive TCDS. TCDS varies noticeably with growth time and chromosomal region. These observations are fully consistent with the notion of local genomic architecture providing additional means for modulating and fine-tuning the effects of changing DNA superhelicity on gene expression during the bacterial growth cycle.
5. Spatially Shifting Superhelicity Optimum Determines the Temporal Gene Expression
There is little doubt that supercoiling regulates global gene expression, but this by itself does not explain how the spatiotemporal order of gene expression is established. The latter has been correlated with gene replication
[98] and thus with the cell cycle, but what about the gene expression order during the growth cycle of the bacterial population?
It is conceivable that the genomic sequence organization, in conjunction with the temporal gradient of supercoiling, determines the spatiotemporal gene expression order during the bacterial growth cycle. As mentioned above, the genomic sequence of
E. coli, and Gamma-proteobacteria in general, demonstrates a gradient of DNA thermodynamic stability (approximated by G/C-richness) along the OriC-Ter axis of the chromosome, whereas in general, the relatively more A/T-rich and more G/C-rich sequences respond to low and high levels of superhelicity, respectively
[22][23][99].
In particular, the G/C composition of the promoter discriminator sequence (the sequence between the -10 promoter element and the transcription startpoint) was shown to be determinative for the supercoiling response
[99][100]. Individual gene promoters for which the supercoiling response optima were studied in detail, are shown in
Table 1. Although the sample of supercoiling-sensitive promoters provided in
Table 1 is far from comprehensive, it shows a trend suggesting that the closer the promoter is located to OriC, the higher the optimum superhelical density for its transcription, consistent with the gradual decline in the average DNA thermodynamic stability (~G/C content) as a function of distance from OriC. This general trend—a positive correlation between the transcriptional response to high negative superhelicity, high negative melting energy content (~G/C-richness), and proximity to OriC of the genes—is corroborated by time-resolved analyses of global genomic transcription during the growth cycle.
Researchers have argued that the process of replication as well as of ribosomal gene transcription proceeding directionally from OriC toward Ter imposes a supercoiling asymmetry with higher negative superhelicity accumulated at the OriC pole of the chromosome
[1][35]. At the early growth stage under conditions of high negative superhelical density, the OriC proximal G/C-rich sequences are likely to be transcribed optimally. However, it is conceivable that the gradual decrease in the global superhelical density during the growth cycle
[67] would shift the optimum of the promoter supercoiling response toward genomic regions with lower average G/C content and hence, away from OriC toward Ter. Thus, with the passage of time and the associated decrease in global negative superhelicity, progressively more OriC-distal genes with lower G/C content and lower optima of superhelical density for transcription would become maximally active. Researchers propose that the coupling of the temporal growth-phase-dependent gradient of chromosomal DNA supercoiling with the spatial gradient of genomic DNA thermodynamic stability along the chromosomal OriC-Ter axis acts as a timing chain determining the spatiotemporal order of gene transcription during the bacterial growth cycle.
6. The Temporal Gradient of Superhelicity Reflects the Gradient of Ion Composition and Intracellular pH
The temporal gradient of superhelicity apparent during the growth cycle and its correlation with the energy requirement for promoter opening, suggest that both these effects are a response to changes in energy availability. A primary response to shift-up is an almost immediate doubling of the intracellular energy charge from the stationary phase level
[103]. On a similar time scale, there is a rapid influx of K
+ mediated by the proton-dependent K
+ transporting P-type ATPase accompanied by an efflux of H
+ and Na
+ resulting in an excess of K
+ [104]. The net effect of these changes in the ion composition is an increase in intracellular pH
[105]. The change in energy charge could also promote DNA gyrase activity
[33]. During the subsequent growth cycle, the relative concentrations of K
+ and Na
+ are rebalanced
[104] whereas the energy charge falls from the onset of the stationary phase
[103], and the intracellular H
+ concentration increases
[106]. The observed changes in the intracellular ion composition during the growth cycle are likely relevant to transcriptional regulation and DNA compaction. The efficient in vitro transcription of supercoiled DNA is K
+ (and not Na
+) dependent
[107], whereas during the stationary phase the relative DNA binding of the major NAPs, IHF and Dps, depends on the pH, K
+, and Mg
++ concentrations, with that of Dps being favored by low pH and low [K
+]
[61]. Researchers suggest that the temporal gradient of superhelicity reflects a corresponding temporal gradient in the ion composition and intracellular pH, which in turn would influence energy availability via the F
0F
1 ATPase.