Transcriptional regulation shapes the organization of genes on bacterial chromosomes

An ancient bacterial zinc acquisition system identified from a cyanobacterial exoproteome

Bacteria have developed fine-tuned responses to cope with potential zinc limitation. The Zur protein is a key player in coordinating this response in most species. Comparative proteomics conducted on the cyanobacterium Anabaena highlighted the more abundant proteins in a zur mutant compared to the wild type. Experimental evidence showed that the exoprotein ZepA mediates zinc uptake. Genomic context of the zepA gene and protein structure prediction provided additional insights on the regulation and putative function of ZepA homologs. Phylogenetic analysis suggests that ZepA represents a primordial system for zinc acquisition that has been conserved for billions of years in a handful of species from distant bacterial lineages. Furthermore, these results show that Zur may have been one of the first regulators of the FUR family to evolve, consistent with the scarcity of zinc in the ecosystems of the Archean eon.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

Composition of Transcription Machinery and Its Crosstalk with Nucleoid-Associated Proteins and Global Transcription Factors

The coordination of bacterial genomic transcription involves an intricate network of interdependent genes encoding nucleoid-associated proteins (NAPs), DNA topoisomerases, RNA polymerase subunits and modulators of transcription machinery. The central element of this homeostatic regulatory system, integrating the information on cellular physiological state and producing a corresponding transcriptional response, is the multi-subunit RNA polymerase (RNAP) holoenzyme. In this review article, we argue that recent observations revealing DNA topoisomerases and metabolic enzymes associated with RNAP supramolecular complex support the notion of structural coupling between transcription machinery, DNA topology and cellular metabolism as a fundamental device coordinating the spatiotemporal genomic transcription. We analyse the impacts of various combinations of RNAP holoenzymes and global transcriptional regulators such as abundant NAPs, on genomic transcription from this viewpoint, monitoring the spatiotemporal patterns of couplons—overlapping subsets of the regulons of NAPs and RNAP sigma factors. We show that the temporal expression of regulons is by and large, correlated with that of cognate regulatory genes, whereas both the spatial organization and temporal expression of couplons is distinctly impacted by the regulons of NAPs and sigma factors. We propose that the coordination of the growth phase-dependent concentration gradients of global regulators with chromosome configurational dynamics determines the spatiotemporal patterns of genomic expression.

Chromosomal origin of replication coordinates logically distinct types of bacterial genetic regulation

For a long time it has been hypothesized that bacterial gene regulation involves an intricate interplay of the transcriptional regulatory network (TRN) and the spatial organization of genes in the chromosome. Here we explore this hypothesis both on a structural and on a functional level. On the structural level, we study the TRN as a spatially embedded network. On the functional level, we analyze gene expression patterns from a network perspective (“digital control”), as well as from the perspective of the spatial organization of the chromosome (“analog control”). Our structural analysis reveals the outstanding relevance of the symmetry axis defined by the origin (Ori) and terminus (Ter) of replication for the network embedding and, thus, suggests the co-evolution of two regulatory infrastructures, namely the transcriptional regulatory network and the spatial arrangement of genes on the chromosome, to optimize the cross-talk between two fundamental biological processes: genomic expression and replication. This observation is confirmed by the functional analysis based on the differential gene expression patterns of more than 4000 pairs of microarray and RNA-Seq datasets for E. coli from the Colombos Database using complex network and machine learning methods. This large-scale analysis supports the notion that two logically distinct types of genetic control are cooperating to regulate gene expression in a complementary manner. Moreover, we find that the position of the gene relative to the Ori is a feature of very high predictive value for gene expression, indicating that the Ori–Ter symmetry axis coordinates the action of distinct genetic control mechanisms.

Chromosomal Organization and Regulation of Genetic Function in Escherichia coli Integrates the DNA Analog and Digital Information

In this article, we summarize our current understanding of the bacterial genetic regulation brought about by decades of studies using the Escherichia coli model. It became increasingly evident that the cellular genetic regulation system is organizationally closed, and a major challenge is to describe its circular operation in quantitative terms. We argue that integration of the DNA analog information (i.e., the probability distribution of the thermodynamic stability of base steps) and digital information (i.e., the probability distribution of unique triplets) in the genome provides a key to understanding the organizational logic of genetic control. During bacterial growth and adaptation, this integration is mediated by changes of DNA supercoiling contingent on environmentally induced shifts in intracellular ionic strength and energy charge. More specifically, coupling of dynamic alterations of the local intrinsic helical repeat in the structurally heterogeneous DNA polymer with structural-compositional changes of RNA polymerase holoenzyme emerges as a fundamental organizational principle of the genetic regulation system. We present a model of genetic regulation integrating the genomic pattern of DNA thermodynamic stability with the gene order and function along the chromosomal OriC-Ter axis, which acts as a principal coordinate system organizing the regulatory interactions in the genome.

Systematic discovery of uncharacterized transcription factors in Escherichia coli K-12 MG1655

Transcriptional regulation enables cells to respond to environmental changes. Of the estimated 304 candidate transcription factors (TFs) in Escherichia coli K-12 MG1655, 185 have been experimentally identified, but ChIP methods have been used to fully characterize only a few dozen. Identifying these remaining TFs is key to improving our knowledge of the E. coli transcriptional regulatory network (TRN). Here, we developed an integrated workflow for the computational prediction and comprehensive experimental validation of TFs using a suite of genome-wide experiments. We applied this workflow to (i) identify 16 candidate TFs from over a hundred uncharacterized genes; (ii) capture a total of 255 DNA binding peaks for ten candidate TFs resulting in six high-confidence binding motifs; (iii) reconstruct the regulons of these ten TFs by determining gene expression changes upon deletion of each TF and (iv) identify the regulatory roles of three TFs (YiaJ, YdcI, and YeiE) as regulators of l-ascorbate utilization, proton transfer and acetate metabolism, and iron homeostasis under iron-limited conditions, respectively. Together, these results demonstrate how this workflow can be used to discover, characterize, and elucidate regulatory functions of uncharacterized TFs in parallel.

Large-Scale Functional Analysis of CRP-Mediated Feed-Forward Loops

Feed-forward loops (FFLs) represent an important and basic network motif to understand specific biological functions. Cyclic-AMP (cAMP) receptor protein (CRP), a transcription factor (TF), mediates catabolite repression and regulates more than 400 genes in response to changes in intracellular concentrations of cAMP in Escherichia coli. CRP participates in some FFLs, such as araBAD and araFGH operons and adapts to fluctuating environmental nutrients, thereby enhancing the survivability of E. coli. Although computational simulations have been conducted to explore the potential functionality of FFLs, a comprehensive study on the functions of all structural types on the basis of in vivo data is lacking. Moreover, the regulatory role of CRP-mediated FFLs (CRP-FFLs) remains obscure. We identified 393 CRP-FFLs in E. coli using EcoCyc and RegulonDB. Dose⁻response genomic microarray of E. coli revealed dynamic gene expression of each target gene of CRP-FFLs in response to a range of cAMP dosages. All eight types of FFLs were present in CRP regulon with various expression patterns of each CRP-FFL, which were further divided into five functional groups. The microarray and reported regulatory relationships identified 202 CRP-FFLs that were directly regulated by CRP in these eight types of FFLs. Interestingly, 34% (147/432) of genes were directly regulated by CRP and CRP-regulated TFs, which indicates that these CRP-regulated genes were also regulated by other CRP-regulated TFs responding to environmental signals through CRP-FFLs. Furthermore, we applied gene ontology annotation to reveal the biological functions of CRP-FFLs.

Broad-scale redistribution of mRNA abundance and transcriptional machinery in response to growth rate in Salmonella enterica serovar Typhimurium

We have investigated the connection between the four-dimensional architecture of the bacterial nucleoid and the organism's global gene expression programme. By localizing the transcription machinery and the transcriptional outputs across the genome of the model bacterium Salmonella enterica serovar Typhimurium at different stages of the growth cycle, a surprising disconnection between gene dosage and transcriptional output was revealed. During exponential growth, gene output occurred chiefly in the Ori (origin), Ter (terminus) and NSL (non-structured left) domains, whereas the Left macrodomain remained transcriptionally quiescent at all stages of growth. The apparently high transcriptional output in Ter was correlated with an enhanced stability of the RNA expressed there during exponential growth, suggesting that longer mRNA half-lives compensate for low gene dosage. During exponential growth, RNA polymerase (RNAP) was detected everywhere, whereas in stationary phase cells, RNAP was concentrated in the Ter macrodomain. The alternative sigma factors RpoE, RpoH and RpoN were not required to drive transcription in these growth conditions, consistent with their observed binding to regions away from RNAP and regions of active transcription. Specifically, these alternative sigma factors were found in the Ter macrodomain during exponential growth, whereas they were localized at the Ori macrodomain in stationary phase.

Chromosomal organization of transcription: in a nutshell

Early studies of transcriptional regulation focused on individual gene promoters defined specific transcription factors as central agents of genetic control. However, recent genome-wide data propelled a different view by linking spatially organized gene expression patterns to chromosomal dynamics. Therefore, the major problem in contemporary molecular genetics concerned with transcriptional gene regulation is to establish a unifying model that reconciles these two views. This problem, situated at the interface of polymer physics and network theory, requires development of an integrative methodology. In this review, we discuss recent achievements in classical model organism E. coli and provide some novel insights gained from studies of a bacterial plant pathogen, D. dadantii. We consider DNA topology and the basal transcription machinery as key actors of regulation, in which activation of functionally relevant genes is coupled to and coordinated with the establishment of extended chromosomal domains of coherent transcription. We argue that the spatial organization of genome plays a fundamental role in its own regulation.

Sickeningly Sweet: L-rhamnose stimulates Flavobacterium columnare biofilm formation and virulence

Flavobacterium columnare, the causative agent of columnaris disease, causes substantial mortality worldwide in numerous freshwater finfish species. Due to its global significance and impact on the aquaculture industry continual efforts to better understand basic mechanisms that contribute to disease are urgently needed. The current work sought to evaluate the effect of L-rhamnose on the growth characteristics of F. columnare. While we initially did not observe any key changes during the total growth of F. columnare isolates tested when treated with L-rhamnose, it soon became apparent that the difference lies in the ability of this carbohydrate to facilitate the formation of biofilms. The addition of different concentrations of L-rhamnose consistently promoted the development of biofilms among different F. columnare isolates; however, it does not appear to be sufficient as a sole carbon source for biofilm growth. Our data also suggest that iron acquisition machinery is required for biofilm development. Finally, the addition of different concentrations of L-rhamnose to F. columnare prior to a laboratory challenge increased mortality rates in channel catfish (Ictalurus punctatus) as compared to controls. These results provide further evidence that biofilm formation is an integral virulence factor in the initiation of disease in fish.

Effect of genomic distance on coexpression of coregulated genes in E. coli

In prokaryotes, genomic distance is a feature that in addition to coregulation affects coexpression. Several observations, such as genomic clustering of highly coexpressed small regulons, support the idea that coexpression behavior of coregulated genes is affected by the distance between the coregulated genes. However, the specific contribution of distance in addition to coregulation in determining the degree of coexpression has not yet been studied systematically. In this work, we exploit the rich information in RegulonDB to study how the genomic distance between coregulated genes affects their degree of coexpression, measured by pairwise similarity of expression profiles obtained under a large number of conditions. We observed that, in general, coregulated genes display higher degrees of coexpression as they are more closely located on the genome. This contribution of genomic distance in determining the degree of coexpression was relatively small compared to the degree of coexpression that was determined by the tightness of the coregulation (degree of overlap of regulatory programs) but was shown to be evolutionary constrained. In addition, the distance effect was sufficient to guarantee coexpression of coregulated genes that are located at very short distances, irrespective of their tightness of coregulation. This is partly but definitely not always because the close distance is also the cause of the coregulation. In cases where it is not, we hypothesize that the effect of the distance on coexpression could be caused by the fact that coregulated genes closely located to each other are also relatively more equidistantly located from their common TF and therefore subject to more similar levels of TF molecules. The absolute genomic distance of the coregulated genes to their common TF-coding gene tends to be less important in determining the degree of coexpression. Our results pinpoint the importance of taking into account the combined effect of distance and coregulation when studying prokaryotic coexpression and transcriptional regulation.

Functional architecture and global properties of the Corynebacterium glutamicum regulatory network: Novel insights from a dataset with a high genomic coverage

Corynebacterium glutamicum is a Gram-positive, anaerobic, rod-shaped soil bacterium able to grow on a diversity of carbon sources like sugars and organic acids. It is a biotechnological relevant organism because of its highly efficient ability to biosynthesize amino acids, such as l-glutamic acid and l-lysine. Here, we reconstructed the most complete C. glutamicum regulatory network to date and comprehensively analyzed its global organizational properties, systems-level features and functional architecture. Our analyses show the tremendous power of Abasy Atlas to study the functional organization of regulatory networks. We created two models of the C. glutamicum regulatory network: all-evidences (containing both weak and strong supported interactions, genomic coverage=73%) and strongly-supported (only accounting for strongly supported evidences, genomic coverage=71%). Using state-of-the-art methodologies, we prove that power-law behaviors truly govern the connectivity and clustering coefficient distributions. We found a non-previously reported circuit motif that we named complex feed-forward motif. We highlighted the importance of feedback loops for the functional architecture, beyond whether they are statistically over-represented or not in the network. We show that the previously reported top-down approach is inadequate to infer the hierarchy governing a regulatory network because feedback bridges different hierarchical layers, and the top-down approach disregards the presence of intermodular genes shaping the integration layer. Our findings all together further support a diamond-shaped, three-layered hierarchy exhibiting some feedback between processing and coordination layers, which is shaped by four classes of systems-level elements: global regulators, locally autonomous modules, basal machinery and intermodular genes.

Conserved Units of Co-Expression in Bacterial Genomes: An Evolutionary Insight into Transcriptional Regulation

Genome-wide measurements of transcriptional activity in bacteria indicate that the transcription of successive genes is strongly correlated beyond the scale of operons. Here, we analyze hundreds of bacterial genomes to identify supra-operonic segments of genes that are proximal in a large number of genomes. We show that these synteny segments correspond to genomic units of strong transcriptional co-expression. Structurally, the segments contain operons with specific relative orientations (co-directional or divergent) and nucleoid-associated proteins are found to bind at their boundaries. Functionally, operons inside a same segment are highly co-expressed even in the apparent absence of regulatory factors at their promoter regions. Remote operons along DNA can also be co-expressed if their corresponding segments share a transcriptional or sigma factor, without requiring these factors to bind directly to the promoters of the operons. As evidence that these results apply across the bacterial kingdom, we demonstrate them both in the Gram-negative bacterium Escherichia coli and in the Gram-positive bacterium Bacillus subtilis. The underlying process that we propose involves only RNA-polymerases and DNA: it implies that the transcription of an operon mechanically enhances the transcription of adjacent operons. In support of a primary role of this regulation by facilitated co-transcription, we show that the transcription en bloc of successive operons as a result of transcriptional read-through is strongly and specifically enhanced in synteny segments. Finally, our analysis indicates that facilitated co-transcription may be evolutionary primitive and may apply beyond bacteria.

Engineering Synthetic cis-Regulatory Elements for Simultaneous Recognition of Three Transcriptional Factors in Bacteria

Recognition of cis-regulatory elements by transcription factors (TF) at target promoters is crucial to gene regulation in bacteria. In this process, binding of TFs to their cognate sequences depends on a set of physical interactions between these proteins and specific nucleotides in the operator region. Previously, we showed that in silico optimization algorithms are able to generate short sequences that are recognized by two different TFs of Escherichia coli, namely, CRP and IHF, thus generating an AND logic gate. Here, we expanded this approach in order to engineer DNA sequences that can be simultaneously recognized by three unrelated TFs (CRP, IHF, and Fis). Using in silico optimization and experimental validation strategies, we were able to obtain a candidate promoter (Plac-CFI1) regulated by only two TFs with an AND logic, thus demonstrating a limitation in the design. Subsequently, we modified the algorithm to allow the optimization of extended sequences, and were able to design two synthetic promoters (PCFI20-1 and PCFI22-5) that were functional in vivo. Expression assays in E. coli mutant strains for each TF revealed that while CRP positively regulates the promoter activities, IHF and Fis are strong repressors of both the promoter variants. Taken together, our results demonstrate the potential of in silico strategies in bacterial synthetic promoter engineering. Furthermore, the study also shows how small modifications in cis-regulatory elements can drastically affect the final logic of the resulting promoter.

Pharmacological manipulation of transcription factor protein-protein interactions: opportunities and obstacles

Much research on transcription factor biology and their genetic pathways has been undertaken over the last 30 years, especially in the field of developmental biology and cancer. Yet, very little is known about the molecular modalities of highly dynamic interactions between transcription factors, genomic DNA, and protein partners. Methodological breakthroughs such as RNA-seq (RNA-sequencing), ChIP-seq (chromatin immunoprecipitation sequencing), RIME (rapid immunoprecipitation mass spectrometry of endogenous proteins), and single-molecule imaging will dramatically accelerate the discovery rate of their molecular mode of action in the next few years. From a pharmacological viewpoint, conventional methods used to target transcription factor activity with molecules mimicking endogenous ligands fail to achieve high specificity and are limited by a lack of identification of new molecular targets. Protein-protein interactions are likely to represent one of the next major classes of therapeutic targets. Transcription factors, known to act mostly via protein-protein interaction, may well be at the forefront of this type of drug development. One hurdle in this field remains the difficulty to collate structural data into meaningful information for rational drug design. Another hurdle is the lack of chemical libraries meeting the structural requirements of protein-protein interaction disruption. As more attempts at modulating transcription factor activity are undertaken, valuable knowledge will be accumulated on the modality of action required to modulate transcription and how these findings can be applied to developing transcription factor drugs. Key discoveries will spawn into new therapeutic approaches not only as anticancer targets but also for other indications, such as those with an inflammatory component including neurodegenerative disorders, diabetes, and chronic liver and kidney diseases.

Expanding the logic of bacterial promoters using engineered overlapping operators for global regulators

The understanding of how the architecture of cis-regulatory elements at bacterial promoters determines their final output is of central interest in modern biology. In this work, we attempt to gain insight into this process by analyzing complex promoter architectures in the model organism Escherichia coli. By focusing on the relationship between different TFs at the genomic scale in terms of their binding site arrangement and their effect on the target promoters, we found no strong constraint limiting the combinatorial assembly of TF pairs in E. coli. More strikingly, overlapping binding sites were found equally associated with both equivalent (both TFs have the same effect on the promoter) and opposite (one TF activates while the other repress the promoter) effects on gene expression. With this information on hand, we set an in silico approach to design overlapping sites for three global regulators (GRs) of E. coli, specifically CRP, Fis, and IHF. Using random sequence assembly and an evolutionary algorithm, we were able to identify potential overlapping operators for all TF pairs. In order to validate our prediction, we constructed two lac promoter variants containing overlapping sites for CRP and IHF designed in silico. By assaying the synthetic promoters using a GFP reporter system, we demonstrated that these variants were functional and activated by CRP and IHF in vivo. Taken together, presented results add new information on the mechanisms of signal integration in bacterial promoters and provide new strategies for the engineering of synthetic regulatory circuits in bacteria.

Spatial organization of transcription in bacterial cells

Prokaryotic transcription has been extensively studied over the past half a century. However, there often exists a gap between the structural, mechanistic description of transcription obtained from in vitro biochemical studies, and the cellular, phenomenological observations from in vivo genetic studies. It is now accepted that a living bacterial cell is a complex entity; the heterogeneous cellular environment is drastically different from the homogenous, well-mixed situation in vitro. Where molecules are inside a cell may be important for their function; hence, the spatial organization of different molecular components may provide a new means of transcription regulation in vivo, possibly bridging this gap. In this review, we survey current evidence for the spatial organization of four major components of transcription [genes, transcription factors, RNA polymerase (RNAP) and RNAs] and critically analyze their biological significance.

Integration of syntactic and semantic properties of the DNA code reveals chromosomes as thermodynamic machines converting energy into information

Understanding genetic regulation is a problem of fundamental importance. Recent studies have made it increasingly evident that, whereas the cellular genetic regulation system embodies multiple disparate elements engaged in numerous interactions, the central issue is the genuine function of the DNA molecule as information carrier. Compelling evidence suggests that the DNA, in addition to the digital information of the linear genetic code (the semantics), encodes equally important continuous, or analog, information that specifies the structural dynamics and configuration (the syntax) of the polymer. These two DNA information types are intrinsically coupled in the primary sequence organisation, and this coupling is directly relevant to regulation of the genetic function. In this review, we emphasise the critical need of holistic integration of the DNA information as a prerequisite for understanding the organisational complexity of the genetic regulation system.

Effect of DNA conformation on facilitated diffusion

Within a living cell, site-specific DNA-binding proteins need to search the whole genome to find a target of ~10-20 bp. That they find the target, and do so quickly, is vital for the correct functioning of the DNA, and of the cell as a whole. The current understanding is that this search is performed via facilitated diffusion, i.e. by combining three-dimensional bulk diffusion within the cytoplasm or nucleoplasm, with one-dimensional diffusion along the DNA backbone, to which the protein binds non-specifically. After reviewing the standard theory of facilitated diffusion, we discuss in the present article the still rather rare direct computer simulations of this process, focusing on the three-dimensional part of the search, and the effect of DNA looping and the general DNA conformation on its efficiency. We close by highlighting some open questions in this field.

DNA-binding-protein inhomogeneity in E. coli modeled as biphasic facilitated diffusion

We have recently shown that nonspecifically bound lac repressors are spatially inhomogeneous in E. coli cells and depends upon the location of its encoding gene and the DNA compaction state [Kuhlman and Cox, Mol. Syst. Biol. 8, 610 (2012)]. Here we model this inhomogeneity as a consequence of diffusion within and exchange between two distinct intracellular phases: the condensed chromosomal DNA and an extrachromosomal compartment, the cytoplasm. We discuss the consequences of this model for the target search process.

Distance matters: the impact of gene proximity in bacterial gene regulation

Following recent discoveries of colocalization of downstream-regulating genes in living cells, the impact of the spatial distance between such genes on the kinetics of gene product formation is increasingly recognized. We here show from analytical and numerical analysis that the distance between a transcription factor (TF) gene and its target gene drastically affects the speed and reliability of transcriptional regulation in bacterial cells. For an explicit model system, we develop a general theory for the interactions between a TF and a transcription unit. The observed variations in regulation efficiency are linked to the magnitude of the variation of the TF concentration peaks as a function of the binding site distance from the signal source. Our results support the role of rapid binding site search for gene colocalization and emphasize the role of local concentration differences.

Transcriptional regulatory network shapes the genome structure of Saccharomyces cerevisiae

Among cellular processes gene transcription is central. More and more evidence is mounting that transcription is tightly connected with the spatial organization of the chromosomes. Spatial proximity of genes sharing transcriptional machinery is one of the consequences of this organization. Motivated by information on the physical relationship among genes identified via chromosomal conformation capture methods, we complement the spatial organization with the idea that genes under similar transcription factor control, but possible scattered throughout the genome, might be in physically proximity to facilitate the access of their commonly used transcription factors. Unlike the transcription factory model, "interacting" genes in our "Gene Proximity Model" are not necessarily immediate physical neighbors but are in spatial proximity. Considering the stochastic nature of TF-promoter binding, this local condensation mechanism could serve as a tie to recruit co-regulated genes to guarantee the swiftness of biological reactions. We tested this idea with a simple eukaryotic organism, Saccharomyces cerevisiae. Chromosomal interaction patterns and folding behavior generated by our model re-construct those obtained from experiments. We show that the transcriptional regulatory network has a close linkage with the genome organization in budding yeast, which is fundamental and instrumental to later studies on other more complex eukaryotes.

DNA thermodynamic stability and supercoil dynamics determine the gene expression program during the bacterial growth cycle

The chromosomal DNA polymer constituting the cellular genetic material is primarily a device for coding information. Whilst the gene sequences comprise the digital (discontinuous) linear code, physiological alterations of the DNA superhelical density generate in addition analog (continuous) three-dimensional information essential for regulation of both chromosome compaction and gene expression. Insight into the relationship between the DNA analog information and the digital linear code is of fundamental importance for understanding genetic regulation. Our previous study in the model organism Escherichia coli suggested that the chromosomal gene order and a spatiotemporal gradient of DNA superhelicity associated with DNA replication determine the growth phase-dependent gene transcription. In this study we reveal a general gradient of DNA thermodynamic stability correlated with the polarity of chromosomal replication and manifest in the spatiotemporal pattern of gene transcription during the bacterial growth cycle. Furthermore, by integrating the physical and dynamic features of the transcribed sequences with their functional content we identify spatiotemporal domains of gene expression encompassing different functions. We thus provide both an insight into the organisational principle of the bacterial growth program and a novel holistic methodology for exploring chromosomal dynamics.

Hierarchical dynamics of a transcription factors network in E. coli

In the present work we study the transcription-factor regulatory network that controls the synthesis of flagella in E. coli. Our objective is to address how the transcription-factor dynamics (in terms of their promoter activities and associated rates) correlate with their positions in the hierarchical organization of this regulatory network. Our results suggest that global-regulator promoters express at higher rates than those of local regulators, particularly when the bacterial populations are actively growing. Furthermore, promoter activity decreases together with the rate of cellular division. And finally, local-regulator promoters reach their maximal activity later than global-regulator promoters do. In summary, our results suggest a strong correlation between promoter activities and their hierarchical organization in this particular regulatory network.

Global genomic arrangement of bacterial genes is closely tied with the total transcriptional efficiency

The availability of a large number of sequenced bacterial genomes allows researchers not only to derive functional and regulation information about specific organisms but also to study the fundamental properties of the organization of a genome. Here we address an important and challenging question regarding the global arrangement of operons in a bacterial genome: why operons in a bacterial genome are arranged in the way they are. We have previously studied this question and found that operons of more frequently activated pathways tend to be more clustered together in a genome. Specifically, we have developed a simple sequential distance-based pseudo energy function and found that the arrangement of operons in a bacterial genome tend to minimize the clusteredness function (C value) in comparison with artificially-generated alternatives, for a variety of bacterial genomes. Here we extend our previous work, and report a number of new observations: (a) operons of the same pathways tend to group into a few clusters rather than one; and (b) the global arrangement of these operon clusters tend to minimize a new “energy” function (C⁺ value) that reflects the efficiency of the transcriptional activation of the encoded pathways. These observations provide insights into further study of the genomic organization of genes in bacteria.

Extensive cross-talk and global regulators identified from an analysis of the integrated transcriptional and signaling network in Escherichia coli

To understand the regulatory dynamics of transcription factors (TFs) and their interplay with other cellular components we have integrated transcriptional, protein-protein and the allosteric or equivalent interactions which mediate the physiological activity of TFs in Escherichia coli. To study this integrated network we computed a set of network measurements followed by principal component analysis (PCA), investigated the correlations between network structure and dynamics, and carried out a procedure for motif detection. In particular, we show that outliers identified in the integrated network based on their network properties correspond to previously characterized global transcriptional regulators. Furthermore, outliers are highly and widely expressed across conditions, thus supporting their global nature in controlling many genes in the cell. Motifs revealed that TFs not only interact physically with each other but also obtain feedback from signals delivered by signaling proteins supporting the extensive cross-talk between different types of networks. Our analysis can lead to the development of a general framework for detecting and understanding global regulatory factors in regulatory networks and reinforces the importance of integrating multiple types of interactions in underpinning the interrelationships between them.

Regulatory consequences of gene translocation in bacteria

Gene translocations play an important role in the plasticity and evolution of bacterial genomes. In this study, we investigated the impact on gene regulation of three genome organizational features that can be altered by translocations: (i) chromosome position; (ii) gene orientation; and (iii) the distance between a target gene and its transcription factor gene (‘target-TF distance’). Specifically, we quantified the effect of these features on constitutive expression, transcription factor binding and/or gene expression noise using a synthetic network in Escherichia coli composed of a transcription factor (LacI repressor) and its target gene (yfp). Here we show that gene regulation is generally robust to changes in chromosome position, gene orientation and target-TF distance. The only demonstrable effect was that chromosome position alters constitutive expression, due to changes in gene copy number and local sequence effects, and that this determines maximum and minimum expression levels. The results were incorporated into a mathematical model which was used to quantitatively predict the responses of a simple gene network to gene translocations; the predictions were confirmed experimentally. In summary, gene translocation can modulate constitutive gene expression levels due to changes in chromosome position but it has minimal impact on other facets of gene regulation.

The percentage of bacterial genes on leading versus lagging strands is influenced by multiple balancing forces

The majority of bacterial genes are located on the leading strand, and the percentage of such genes has a large variation across different bacteria. Although some explanations have been proposed, these are at most partial explanations as they cover only small percentages of the genes and do not even consider the ones biased toward the lagging strand. We have carried out a computational study on 725 bacterial genomes, aiming to elucidate other factors that may have influenced the strand location of genes in a bacterium. Our analyses suggest that (i) genes of some functional categories such as ribosome have higher preferences to be on the leading strands; (ii) genes of some functional categories such as transcription factor have higher preferences on the lagging strands; (iii) there is a balancing force that tends to keep genes from all moving to the leading and more efficient strand and (iv) the percentage of leading-strand genes in an bacterium can be accurately explained based on the numbers of genes in the functional categories outlined in (i) and (ii), genome size and gene density, indicating that these numbers implicitly contain the information about the percentage of genes on the leading versus lagging strand in a genome.

Prokaryotic regulatory systems biology: Common principles governing the functional architectures of Bacillus subtilis and Escherichia coli unveiled by the natural decomposition approach

Escherichia coli and Bacillus subtilis are two of the best-studied prokaryotic model organisms. Previous analyses of their transcriptional regulatory networks have shown that they exhibit high plasticity during evolution and suggested that both converge to scale-free-like structures. Nevertheless, beyond this suggestion, no analyses have been carried out to identify the common systems-level components and principles governing these organisms. Here we show that these two phylogenetically distant organisms follow a set of common novel biologically consistent systems principles revealed by the mathematically and biologically founded natural decomposition approach. The discovered common functional architecture is a diamond-shaped, matryoshka-like, three-layer (coordination, processing, and integration) hierarchy exhibiting feedback, which is shaped by four systems-level components: global transcription factors (global TFs), locally autonomous modules, basal machinery and intermodular genes. The first mathematical criterion to identify global TFs, the κ-value, was reassessed on B. subtilis and confirmed its high predictive power by identifying all the previously reported, plus three potential, master regulators and eight sigma factors. The functionally conserved cores of modules, basal cell machinery, and a set of non-orthologous common physiological global responses were identified via both orthologous genes and non-orthologous conserved functions. This study reveals novel common systems principles maintained between two phylogenetically distant organisms and provides a comparison of their lifestyle adaptations. Our results shed new light on the systems-level principles and the fundamental functions required by bacteria to sustain life.

Regulatory design governing progression of population growth phases in bacteria

It has long been noted that batch cultures inoculated with resting bacteria exhibit a progression of growth phases traditionally labeled lag, exponential, pre-stationary and stationary. However, a detailed molecular description of the mechanisms controlling the transitions between these phases is lacking. A core circuit, formed by a subset of regulatory interactions involving five global transcription factors (FIS, HNS, IHF, RpoS and GadX), has been identified by correlating information from the well- established transcriptional regulatory network of Escherichia coli and genome-wide expression data from cultures in these different growth phases. We propose a functional role for this circuit in controlling progression through these phases. Two alternative hypotheses for controlling the transition between the growth phases are first, a continuous graded adjustment to changing environmental conditions, and second, a discontinuous hysteretic switch at critical thresholds between growth phases. We formulate a simple mathematical model of the core circuit, consisting of differential equations based on the power-law formalism, and show by mathematical and computer-assisted analysis that there are critical conditions among the parameters of the model that can lead to hysteretic switch behavior, which--if validated experimentally--would suggest that the transitions between different growth phases might be analogous to cellular differentiation. Based on these provocative results, we propose experiments to test the alternative hypotheses.

A model for Escherichia coli chromosome packaging supports transcription factor-induced DNA domain formation

What physical mechanism leads to organization of a highly condensed and confined circular chromosome? Computational modeling shows that confinement-induced organization is able to overcome the chromosome's propensity to mix by the formation of topological domains. The experimentally observed high precision of separate subcellular positioning of loci (located on different chromosomal domains) in Escherichia coli naturally emerges as a result of entropic demixing of such chromosomal loops. We propose one possible mechanism for organizing these domains: regulatory control defined by the underlying E. coli gene regulatory network requires the colocalization of transcription factor genes and target genes. Investigating this assumption, we find the DNA chain to self-organize into several topologically distinguishable domains where the interplay between the entropic repulsion of chromosomal loops and their compression due to the confining geometry induces an effective nucleoid filament-type of structure. Thus, we propose that the physical structure of the chromosome is a direct result of regulatory interactions. To reproduce the observed precise ordering of the chromosome, we estimate that the domain sizes are distributed between 10 and 700 kb, in agreement with the size of topological domains identified in the context of DNA supercoiling.

Towards a molecular view of transcriptional control

The accumulation of experimental data over recent years has fueled theoretical work on how transcription factors (TFs) search for and recognise their DNA target sites, how they interact with one another, or with other DNA-binding proteins, and how they cope with the compaction of DNA within bacterial nucleoids or within eukaryotic chromatin. Many models have been built to study the kinetic, thermodynamic and mechanistic aspects of these questions. In some cases they have resulted in a relatively clear consensus view, but a number of questions remain controversial. We present an overview of recent work, with an emphasis on models that provide, or can inspire, a better understanding of transcriptional control at a detailed molecular level.

Genomic arrangement of regulons in bacterial genomes

Regulons, as groups of transcriptionally co-regulated operons, are the basic units of cellular response systems in bacterial cells. While the concept has been long and widely used in bacterial studies since it was first proposed in 1964, very little is known about how its component operons are arranged in a bacterial genome. We present a computational study to elucidate of the organizational principles of regulons in a bacterial genome, based on the experimentally validated regulons of E. coli and B. subtilis. Our results indicate that (1) genomic locations of transcriptional factors (TFs) are under stronger evolutionary constraints than those of the operons they regulate so changing a TF's genomic location will have larger impact to the bacterium than changing the genomic position of any of its target operons; (2) operons of regulons are generally not uniformly distributed in the genome but tend to form a few closely located clusters, which generally consist of genes working in the same metabolic pathways; and (3) the global arrangement of the component operons of all the regulons in a genome tends to minimize a simple scoring function, indicating that the global arrangement of regulons follows simple organizational principles.

Gene location and DNA density determine transcription factor distributions in Escherichia coli

The diffusion coefficient of the transcription factor LacI within living Escherichia coli has been measured directly by in vivo tracking to be D = 0.4 μm(2)/s. At this rate, simple models of diffusion lead to the expectation that LacI and other proteins will rapidly homogenize throughout the cell. Here, we test this expectation of spatial homogeneity by single-molecule visualization of LacI molecules non-specifically bound to DNA in fixed cells. Contrary to expectation, we find that the distribution depends on the spatial location of its encoding gene. We demonstrate that the spatial distribution of LacI is also determined by the local state of DNA compaction, and that E. coli can dynamically redistribute proteins by modifying the state of its nucleoid. Finally, we show that LacI inhomogeneity increases the strength with which targets located proximally to the LacI gene are regulated. We propose a model for intranucleoid diffusion that can reconcile these results with previous measurements of LacI diffusion, and we discuss the implications of these findings for gene regulation in bacteria and eukaryotes.

More than just a metabolic regulator--elucidation and validation of new targets of PdhR in Escherichia coli

BACKGROUND

The pyruvate dehydrogenase regulator protein (PdhR) of Escherichia coli acts as a transcriptional regulator in a pyruvate dependent manner to control central metabolic fluxes. However, the complete PdhR regulon has not yet been uncovered. To achieve an extended understanding of its gene regulatory network, we combined large-scale network inference and experimental verification of results obtained by a systems biology approach.

RESULTS

22 new genes contained in two operons controlled by PdhR (previously only 20 regulatory targets in eight operons were known) were identified by analysing a large-scale dataset of E. coli from the Many Microbes Microarray Database and novel expression data from a pdhR knockout strain, as well as a PdhR overproducing strain. We identified a regulation of the glycolate utilization operon glcDEFGBA using chromatin immunoprecipitation and gel shift assays. We show that this regulation could be part of a cross-induction between genes necessary for acetate and pyruvate utilisation controlled through PdhR. Moreover, a link of PdhR regulation to the replication machinery of the cell via control of the transcription of the dcw-cluster was verified in experiments. This augments our knowledge of the functions of the PdhR-regulon and demonstrates its central importance for further cellular processes in E. coli.

CONCLUSIONS

We extended the PdhR regulon by 22 new genes contained in two operons and validated the regulation of the glcDEFGBA operon for glycolate utilisation and the dcw-cluster for cell division proteins experimentally. Our results provide, for the first time, a plausible regulatory link between the nutritional status of the cell and cell replication mediated by PdhR.

The low cost of recombination in creating novel phenotypes: Recombination can create new phenotypes while disrupting well-adapted phenotypes much less than mutation

Recombination is often considered a disruptive force for well-adapted phenotypes, but recent evidence suggests that this cost of recombination can be small. A key benefit of recombination is that it can help create proteins and regulatory circuits with novel and useful phenotypes more efficiently than point mutation. Its effectiveness stems from the large-scale reorganization of genotypes that it causes, which can help explore far-flung regions in genotype space. Recent work on complex phenotypes in model gene regulatory circuits and proteins shows that the disruptive effects of recombination can be very mild compared to the effects of mutation. Recombination thus can have great benefits at a modest cost, but we do not understand the reasons well. A better understanding might shed light on the evolution of recombination and help improve evolutionary strategies in biochemical engineering.

Regulation and control of metabolic fluxes in microbes

After about ten years of research renaissance in metabolism, the present challenge is to understand how metabolic fluxes are controlled by a complex interplay of overlapping regulatory mechanisms. Reconstruction of various regulatory network topologies is steaming, illustrating that we underestimated the broad importance of post-translational modifications such as enzyme phosphorylation or acetylation for microbial metabolism. With the growing topological knowledge, the functional relevance of these regulatory events becomes an even more pressing need. A major knowledge gap resides in the regulatory network of protein-metabolite interactions, simply because we lacked pertinent methods for systematic analyses - but a start has now been made. Perhaps most dramatic was the conceptual shift in our perception of metabolism from an engine of cellular operation to a generator of input and feedback signals for regulatory circuits that govern many important decisions on cell proliferation, differentiation, death, and naturally metabolism.

Pervasive regulation of nucleoid structure and function by nucleoid-associated proteins

Bacterial DNA is organised in a compact nucleoid body that is tightly associated with the coupled transcription and translation of mRNAs. This structure contains abundant DNA-binding proteins which perform both structural and regulatory roles, and, in Escherichia coli, serve to buffer and organise pervasive DNA superhelicity. We argue that NAPs coordinate regulation of gene expression and superhelicity at the global (or chromosomal) and at local (corresponding to promoter activity and genetic recombination) levels.

Insights into plant biomass conversion from the genome of the anaerobic thermophilic bacterium Caldicellulosiruptor bescii DSM 6725

Caldicellulosiruptor bescii DSM 6725 utilizes various polysaccharides and grows efficiently on untreated high-lignin grasses and hardwood at an optimum temperature of ∼80°C. It is a promising anaerobic bacterium for studying high-temperature biomass conversion. Its genome contains 2666 protein-coding sequences organized into 1209 operons. Expression of 2196 genes (83%) was confirmed experimentally. At least 322 genes appear to have been obtained by lateral gene transfer (LGT). Putative functions were assigned to 364 conserved/hypothetical protein (C/HP) genes. The genome contains 171 and 88 genes related to carbohydrate transport and utilization, respectively. Growth on cellulose led to the up-regulation of 32 carbohydrate-active (CAZy), 61 sugar transport, 25 transcription factor and 234 C/HP genes. Some C/HPs were overproduced on cellulose or xylan, suggesting their involvement in polysaccharide conversion. A unique feature of the genome is enrichment with genes encoding multi-modular, multi-functional CAZy proteins organized into one large cluster, the products of which are proposed to act synergistically on different components of plant cell walls and to aid the ability of C. bescii to convert plant biomass. The high duplication of CAZy domains coupled with the ability to acquire foreign genes by LGT may have allowed the bacterium to rapidly adapt to changing plant biomass-rich environments.

Early Career Research Award Lecture. Structure, evolution and dynamics of transcriptional regulatory networks

The availability of entire genome sequences and the wealth of literature on gene regulation have enabled researchers to model an organism's transcriptional regulation system in the form of a network. In such a network, TFs (transcription factors) and TGs (target genes) are represented as nodes and regulatory interactions between TFs and TGs are represented as directed links. In the present review, I address the following topics pertaining to transcriptional regulatory networks. (i) Structure and organization: first, I introduce the concept of networks and discuss our understanding of the structure and organization of transcriptional networks. (ii) Evolution: I then describe the different mechanisms and forces that influence network evolution and shape network structure. (iii) Dynamics: I discuss studies that have integrated information on dynamics such as mRNA abundance or half-life, with data on transcriptional network in order to elucidate general principles of regulatory network dynamics. In particular, I discuss how cell-to-cell variability in the expression level of TFs could permit differential utilization of the same underlying network by distinct members of a genetically identical cell population. Finally, I conclude by discussing open questions for future research and highlighting the implications for evolution, development, disease and applications such as genetic engineering.

Direct and indirect effects of H-NS and Fis on global gene expression control in Escherichia coli

Nucleoid-associated proteins (NAPs) are global regulators of gene expression in Escherichia coli, which affect DNA conformation by bending, wrapping and bridging the DNA. Two of these--H-NS and Fis--bind to specific DNA sequences and structures. Because of their importance to global gene expression, the binding of these NAPs to the DNA was previously investigated on a genome-wide scale using ChIP-chip. However, variation in their binding profiles across the growth phase and the genome-scale nature of their impact on gene expression remain poorly understood. Here, we present a genome-scale investigation of H-NS and Fis binding to the E. coli chromosome using chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-seq). By performing our experiments under multiple time-points during growth in rich media, we show that the binding regions of the two proteins are mutually exclusive under our experimental conditions. H-NS binds to significantly longer tracts of DNA than Fis, consistent with the linear spread of H-NS binding from high- to surrounding lower-affinity sites; the length of binding regions is associated with the degree of transcriptional repression imposed by H-NS. For Fis, a majority of binding events do not lead to differential expression of the proximal gene; however, it has a significant indirect effect on gene expression partly through its effects on the expression of other transcription factors. We propose that direct transcriptional regulation by Fis is associated with the interaction of tandem arrays of Fis molecules to the DNA and possible DNA bending, particularly at operon-upstream regions. Our study serves as a proof-of-principle for the use of ChIP-seq for global DNA-binding proteins in bacteria, which should become significantly more economical and feasible with the development of multiplexing techniques.

Evolution of gene regulatory networks by fluctuating selection and intrinsic constraints

Various characteristics of complex gene regulatory networks (GRNs) have been discovered during the last decade, e.g., redundancy, exponential indegree distributions, scale-free outdegree distributions, mutational robustness, and evolvability. Although progress has been made in this field, it is not well understood whether these characteristics are the direct products of selection or those of other evolutionary forces such as mutational biases and biophysical constraints. To elucidate the causal factors that promoted the evolution of complex GRNs, we examined the effect of fluctuating environmental selection and some intrinsic constraining factors on GRN evolution by using an individual-based model. We found that the evolution of complex GRNs is remarkably promoted by fixation of beneficial gene duplications under unpredictably fluctuating environmental conditions and that some internal factors inherent in organisms, such as mutational bias, gene expression costs, and constraints on expression dynamics, are also important for the evolution of GRNs. The results indicate that various biological properties observed in GRNs could evolve as a result of not only adaptation to unpredictable environmental changes but also non-adaptive processes owing to the properties of the organisms themselves. Our study emphasizes that evolutionary models considering such intrinsic constraining factors should be used as null models to analyze the effect of selection on GRN evolution.

Various organismal traits, including the morphology of multicellular species and metabolism in unicellular species, are determined by the amount and combinations of proteins in the cell. The complex regulatory network plays an important role in controlling the protein profiles in a cell. Recent studies have revealed that gene regulatory networks have many interesting structural and mutational features such as their scale-free structure, mutational robustness, and evolvability. However, why and how these features have emerged from evolution is unknown. In this paper, we constructed an evolutionary model of gene regulatory networks and simulated its evolution under various environmental conditions. The results show that most features of known gene regulatory networks evolve as a result of adaptation to unpredictable environmental fluctuations. In addition, some internal organismal factors, such as mutational bias, gene expression costs, and constraints on expression dynamics, are also important for GRN evolution observed in real organisms. Thus, these GRN features appear to evolve as a result of not only adaptation to unpredictable environmental changes but also non-adaptive processes owing to the properties of the organisms themselves.

Dissecting the expression patterns of transcription factors across conditions using an integrated network-based approach

In prokaryotes, regulation of gene expression is predominantly controlled at the level of transcription. Transcription in turn is mediated by a set of DNA-binding factors called transcription factors (TFs). In this study, we map the complete repertoire of ∼300 TFs of the bacterial model, Escherichia coli, onto gene expression data for a number of nonredundant experimental conditions and show that TFs are generally expressed at a lower level than other gene classes. We also demonstrate that different conditions harbor varying number of active TFs, with an average of about 15% of the total repertoire, with certain stress and drug-induced conditions exhibiting as high as one-third of the collection of TFs. Our results also show that activators are more frequently expressed than repressors, indicating that activation of promoters might be a more common phenomenon than repression in bacteria. Finally, to understand the association of TFs with different conditions and to elucidate their dynamic interplay with other TFs, we develop a network-based framework to identify TFs which act as markers, defined as those which are responsible for condition-specific transcriptional rewiring. This approach allowed us to pinpoint several marker TFs as being central in various specialized conditions such as drug induction or growth condition variations, which we discuss in light of previously reported experimental findings. Further analysis showed that a majority of identified markers effectively control the expression of their regulons and, in general, transcriptional programs of most conditions can be effectively rewired by a very small number of TFs. It was also found that closeness is a key centrality measure which can aid in the successful identification of marker TFs in regulatory networks. Our results suggest the utility of the network-based approaches developed in this study to be applicable for understanding other interactomic data sets.

General organisational principles of the transcriptional regulation system: a tree or a circle?

Recent advances of systemic approaches to gene expression and cellular metabolism provide unforeseen opportunities for relating and integrating extensive datasets describing the transcriptional regulation system as a whole. However, due to the multifaceted nature of the phenomenon, these datasets often contain logically distinct types of information determined by underlying approach and adopted methodology of data analysis. Consequently, to integrate the datasets comprising information on the states of chromatin structure, transcriptional regulatory network and cellular metabolism, a novel methodology enabling interconversion of logically distinct types of information is required. Here we provide a holistic conceptual framework for analysis of global transcriptional regulation as a system coordinated by structural coupling between the transcription machinery and DNA topology, acting as interdependent sensors and determinants of metabolic functions. In this operationally closed system any transition in physiological state represents an emergent property determined by shifts in structural coupling, whereas genetic regulation acts as a genuine device converting one logical type of information into the other.

Identification and genomic analysis of transcription factors in archaeal genomes exemplifies their functional architecture and evolutionary origin

Archaea, which represent a large fraction of the phylogenetic diversity of organisms, are prokaryotes with eukaryote-like basal transcriptional machinery. This organization makes the study of their DNA-binding transcription factors (TFs) and their transcriptional regulatory networks particularly interesting. In addition, there are limited experimental data regarding their TFs. In this work, 3,918 TFs were identified and exhaustively analyzed in 52 archaeal genomes. TFs represented less than 5% of the gene products in all the studied species comparable with the number of TFs identified in parasites or intracellular pathogenic bacteria, suggesting a deficit in this class of proteins. A total of 75 families were identified, of which HTH_3, AsnC, TrmB, and ArsR families were universally and abundantly identified in all the archaeal genomes. We found that archaeal TFs are significantly small compared with other protein-coding genes in archaea as well as bacterial TFs, suggesting that a large fraction of these small-sized TFs could supply the probable deficit of TFs in archaea, by possibly forming different combinations of monomers similar to that observed in eukaryotic transcriptional machinery. Our results show that although the DNA-binding domains of archaeal TFs are similar to bacteria, there is an underrepresentation of ligand-binding domains in smaller TFs, which suggests that protein–protein interactions may act as mediators of regulatory feedback, indicating a chimera of bacterial and eukaryotic TFs’ functionality. The analysis presented here contributes to the understanding of the details of transcriptional apparatus in archaea and provides a framework for the analysis of regulatory networks in these organisms.

Structural and functional map of a bacterial nucleoid

Mapping global protein binding in the E. coli genome reveals extended domains of high protein occupancy.

Genome-wide mapping of transcription factor-DNA interactions in bacterial chromosomes in vivo has begun to reveal global zones occupied by these factors that serve two purposes: compacting the bacterial DNA and influencing global programs of gene transcription.

Mechanisms and evolution of control logic in prokaryotic transcriptional regulation

A major part of organismal complexity and versatility of prokaryotes resides in their ability to fine-tune gene expression to adequately respond to internal and external stimuli. Evolution has been very innovative in creating intricate mechanisms by which different regulatory signals operate and interact at promoters to drive gene expression. The regulation of target gene expression by transcription factors (TFs) is governed by control logic brought about by the interaction of regulators with TF binding sites (TFBSs) in cis-regulatory regions. A factor that in large part determines the strength of the response of a target to a given TF is motif stringency, the extent to which the TFBS fits the optimal TFBS sequence for a given TF. Advances in high-throughput technologies and computational genomics allow reconstruction of transcriptional regulatory networks in silico. To optimize the prediction of transcriptional regulatory networks, i.e., to separate direct regulation from indirect regulation, a thorough understanding of the control logic underlying the regulation of gene expression is required. This review summarizes the state of the art of the elements that determine the functionality of TFBSs by focusing on the molecular biological mechanisms and evolutionary origins of cis-regulatory regions.