DNA Supercoiling and Genomic Sequence Organization: History
Please note this is an old version of this entry, which may differ significantly from the current revision.
Contributor: , ,

The bacterial growth cycle was described  as a closed, self-reproducing, or autopoietic circuit, reestablishing the physiological state of stationary cells initially inoculated in the growth medium. In batch culture, this process of self-reproduction is associated with the gradual decline in available metabolic energy and corresponding change in the physiological state of the population as a function of “travelled distance” along the autopoietic path. This directional alteration of cell physiology is both reflected in and supported by sequential gene expression along the chromosomal OriC-Ter axis. During the E. coli growth cycle, the spatiotemporal order of gene expression is established by coupling the temporal gradient of supercoiling energy to the spatial gradient of DNA thermodynamic stability along the chromosomal OriC-Ter axis.

  • bacterial growth cycle
  • growth phase-dependent gene expression
  • DNA supercoiling

1. Coupling of Chromosomal Gene Order and Transcription

The continuous information underpinning the operation of the bacterial growth cycle is apparently of a hereditary nature and, therefore, it must be encoded (either itself or its mechanism of production) in the cellular genome. In order to vary concertedly with bacterial growth, this information also has to be dynamic. The genetic information encoded in the form of unique genes is of a discontinuous or digital nature and is static. Therefore, it cannot readily satisfy the demands of continuity and dynamics. However, during the E. coli growth cycle the genes, especially those encoding the major regulators of growth and adaptation, are expressed in temporal succession corresponding to their spatial order along the chromosomal replication origin-to-terminus (OriC-Ter) axis. The latter was shown to serve as a coordinate system for genetic regulation [6,7]. Thus, although the underlying gene order is static, the gene expression is subject to control by a dynamic, continuous variable determining the sequential, chromosomal position-dependent gene expression pattern. In this article, we argue that the continuous variable governing the patterns of gene expression at the most fundamental level is the changing genomic distribution of the effective DNA superhelicity [8], defined as the torsional energy available for the untwisting of gene promoters and transcription initiation.
The bacterial chromosome (aka nucleoid) has been proposed to behave as a “smart polymer”, capable of undergoing large conformational transitions (e.g., reversible collapse) in response to small changes in environmental factors such as pH, temperature, or ionic strength [9]. Furthermore, the nucleoid was shown to undergo cycles of compaction-decompaction under the manipulation of crowding conditions [10], in line with the proposed switch-like conformational transition model for chromatin folding [11]. The regulatory role for “pre-programmed” phase transitions [9] was suggested by observations indicating that the nucleoid structure and gene expression are interdependent [6,12,13,14,15,16]. However, the tight coupling of chromosome structural transitions with gene expression dynamics appears at variance with the notion of gene expression as a fundamentally “noisy” stochastic process [17]. The deterministic average gene expression pattern observed at the cell population level emerges gradually from noisy single-cell expressions [18]. Yet, even though the gene expression patterns in individual bacterial cells might appear stochastic, the chromosomal gene order is not. Furthermore, temporally the gene transcription is correlated with the gene order along the OriC-Ter axis at the population level [12,19] and with the timing of gene replication at a single-cell level (see below).

2. Coupling of DNA Structure to Function Using Two Types of DNA Information

The CODOs ("coherent domains of gene expression" denoted formerly as “functional domains”) [12] 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 [96,103,104]. 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 [105,106], 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 [3,107]. 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 [103,104]. 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 [42,84].
The patterns of the sequence-dependent supercoiling responses of genomic transcription demonstrate that in addition to static genetic information (digital code), the chromosomal DNA polymer also provides dynamic information manifested in the variable distribution of the available DNA torque. The latter determines the organization of functionally meaningful spatial transcript patterns in the genome. Essentially, this dynamic information depends on the distinct thermodynamic stabilities of the consecutive DNA base steps, which overlap and thus provide information of a continuous or analogue type [52,107]. The DNA analogue information encoded in various arrangements of base steps stabilizing distinct DNA conformations facilitates the binding of proteins involved in DNA transactions including transcriptional control [108,109,110,111]. More specifically, the genome encodes both the amino acid sequences of DNA-binding proteins and enzymes governing DNA transactions, as well as the dynamic, supercoiling-dependent structures, serving as recognition signatures recruiting these enzymes and DNA-binding proteins to particular genomic loci [3]. This notion is consistent with the proposed evolutionary “multiplexing” of DNA genetic and structural information into the same molecular context [112]. Thus, by operating with these two—discontinuous (digital) and continuous (analogue)—types of information interwoven in the very same genomic sequence, the DNA appears to communicate with itself. We propose that this “self-communication” represents the global feedback mechanism endowing the growing cell population with the capacity to monitor its status quo, a capacity that is reflected, in part, in the interdependence of chromosome structural dynamics and spatial patterns of gene expression [6,15,16].

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 [123,124] concomitantly with an increase in the ATP/ADP ratio [125,126,127]. 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 [66,128]. Increased negative superhelicity both facilitates replication initiation and strongly activates the OriC-proximal fis gene expression [37,39,129], whereas the accumulation of FIS in turn maintains activated ribosomal RNA transcription [130]. 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 [88,131,132,133], which prefers highly supercoiled templates for transcription [134]. Notably, the RNAP σ70 holoenzyme demonstrates an OriC-Ter gradient of binding site frequency distribution correlating with that of the gyrase binding sites [6]. Concomitantly with activation of the fis gene, the expression of the rpoZ gene encoding the RpoZ (ω)subunit of RNA polymerase is also strongly increased [3,124]. RpoZ stabilizes the polymerase σ70 holoenzyme assembly [135] and is also directly involved in mediating the response to the "alarmone" nucleotide ppGpp [23,136]. FIS activates the hupA gene [137] 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 [66,138]. The latter effect requires high levels of negative superhelicity [121]. Accordingly, both FIS and HU stabilize negative supercoils [43,139,140].
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 [141]. The gyrA/B gene expression and the gyrase levels subside [123] 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 [6,142], 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 [6]. Notably, the “late” NAPs as well as the RNAPσS holoenzyme, preferentially bind relaxed DNA [3]. 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 [115,143].
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 [144], presumably due to competition between the early and late growth stage NAPs, respectively, FIS and Dps [142,145]. The DNA binding of IHF and Dps was shown to depend on environmental factors [146]. 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 [147,148,149]. 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 [150]. This latter process is facilitated by FIS protein produced at high levels on nutritional shift-up [151].

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 [68] but also to the various stress factors including suboptimal oxygen tension, temperature, and osmolarity [125,152,153,154]. 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 [68]. 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 [25,134]. 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 [155]. 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 [156].
Second, the superhelical density and the supercoiling response of the genomic transcription change coordinately both during the growth cycle of E. coli [31] 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 [157]. 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 [40]. 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 [158,159].
Fourth, supercoiling can impose directionality by rendering the structural transitions in DNA both deterministic [160] and coordinated [72,73,74,75], whereas the nucleoprotein structures stabilized by NAPs have been implicated in the directional channeling of torsional energy toward the transcription initiation sites [55,161,162]. 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 [163,164,165,166]. 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 [118,121].
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 [167] 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 [168,169,170,171]. 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 [57], 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 [59] 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 [42,84,85].
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 [83,85]. 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.
We 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 [6,66]. 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 [68] 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. We 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 [174]. 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+ [175]. The net effect of these changes in the ion composition is an increase in intracellular pH [176]. The change in energy charge could also promote DNA gyrase activity [126]. During the subsequent growth cycle, the relative concentrations of K+ and Na+ are rebalanced [175] whereas the energy charge falls from the onset of the stationary phase [174], and the intracellular H+ concentration increases [177]. 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 [178], 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+] [146]. We 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 F0F1 ATPase.

This entry is adapted from the peer-reviewed paper 10.3390/biom12060831

This entry is offline, you can click here to edit this entry!
Video Production Service