Genome-wide analysis of Fis binding in Escherichia coli indicates a causative role for A-/AT-tracts

The interplay between DNA topology and accessory factors in site-specific recombination in bacteria and their bacteriophages

Site-specific recombination is employed widely in bacteria and bacteriophage as a basis for genetic switching events that control phenotypic variation. It plays a vital role in the life cycles of phages and in the replication cycles of chromosomes and plasmids in bacteria. Site-specific recombinases drive these processes using very short segments of identical (or nearly identical) DNA sequences. In some cases, the efficiencies of the recombination reactions are modulated by the topological state of the participating DNA sequences and by the availability of accessory proteins that shape the DNA. These dependencies link the molecular machines that conduct the recombination reactions to the physiological state of the cell. This is because the topological state of bacterial DNA varies constantly during the growth cycle and so does the availability of the accessory factors. In addition, some accessory factors are under allosteric control by metabolic products or second messengers that report the physiological status of the cell. The interplay between DNA topology, accessory factors and site-specific recombination provides a powerful illustration of the connectedness and integration of molecular events in bacterial cells and in viruses that parasitise bacterial cells.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Lack of the H-NS Protein Results in Extended and Aberrantly Positioned DNA during Chromosome Replication and Segregation in Escherichia coli

The architectural protein H-NS binds nonspecifically to hundreds of sites throughout the chromosome and can multimerize to stiffen segments of DNA as well as to form DNA-protein-DNA bridges. H-NS has been suggested to contribute to the orderly folding of the Escherichia coli chromosome in the highly compacted nucleoid. In this study, we investigated the positioning and dynamics of the origins, the replisomes, and the SeqA structures trailing the replication forks in cells lacking the H-NS protein. In H-NS mutant cells, foci of SeqA, replisomes, and origins were irregularly positioned in the cell. Further analysis showed that the average distance between the SeqA structures and the replisome was increased by ∼100 nm compared to that in wild-type cells, whereas the colocalization of SeqA-bound sister DNA behind replication forks was not affected. This result may suggest that H-NS contributes to the folding of DNA along adjacent segments. H-NS mutant cells were found to be incapable of adopting the distinct and condensed nucleoid structures characteristic of E. coli cells growing rapidly in rich medium. It appears as if H-NS mutant cells adopt a “slow-growth” type of chromosome organization under nutrient-rich conditions, which leads to a decreased cellular DNA content.

IMPORTANCE

It is not fully understood how and to what extent nucleoid-associated proteins contribute to chromosome folding and organization during replication and segregation in Escherichia coli. In this work, we find in vivo indications that cells lacking the nucleoid-associated protein H-NS have a lower degree of DNA condensation than wild-type cells. Our work suggests that H-NS is involved in condensing the DNA along adjacent segments on the chromosome and is not likely to tether newly replicated strands of sister DNA. We also find indications that H-NS is required for rapid growth with high DNA content and for the formation of a highly condensed nucleoid structure under such conditions.

FabR regulates Salmonella biofilm formation via its direct target FabB

Background

Biofilm formation is an important survival strategy of Salmonella in all environments. By mutant screening, we showed a knock-out mutant of fabR, encoding a repressor of unsaturated fatty acid biosynthesis (UFA), to have impaired biofilm formation. In order to unravel how this regulator impinges on Salmonella biofilm formation, we aimed at elucidating the S. Typhimurium FabR regulon. Hereto, we applied a combinatorial high-throughput approach, combining ChIP-chip with transcriptomics.

Results

All the previously identified E. coli FabR transcriptional target genes (fabA, fabB and yqfA) were shown to be direct S. Typhimurium FabR targets as well. As we found a fabB overexpressing strain to partly mimic the biofilm defect of the fabR mutant, the effect of FabR on biofilms can be attributed at least partly to FabB, which plays a key role in UFA biosynthesis. Additionally, ChIP-chip identified a number of novel direct FabR targets (the intergenic regions between hpaR/hpaG and ddg/ydfZ) and yet putative direct targets (i.a. genes involved in tRNA metabolism, ribosome synthesis and translation). Next to UFA biosynthesis, a number of these direct targets and other indirect targets identified by transcriptomics (e.g. ribosomal genes, ompA, ompC, ompX, osmB, osmC, sseI), could possibly contribute to the effect of FabR on biofilm formation.

Conclusion

Overall, our results point at the importance of FabR and UFA biosynthesis in Salmonella biofilm formation and their role as potential targets for biofilm inhibitory strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2387-x) contains supplementary material, which is available to authorized users.

DNA Sequence Determinants Controlling Affinity, Stability and Shape of DNA Complexes Bound by the Nucleoid Protein Fis

The abundant Fis nucleoid protein selectively binds poorly related DNA sequences with high affinities to regulate diverse DNA reactions. Fis binds DNA primarily through DNA backbone contacts and selects target sites by reading conformational properties of DNA sequences, most prominently intrinsic minor groove widths. High-affinity binding requires Fis-stabilized DNA conformational changes that vary depending on DNA sequence. In order to better understand the molecular basis for high affinity site recognition, we analyzed the effects of DNA sequence within and flanking the core Fis binding site on binding affinity and DNA structure. X-ray crystal structures of Fis-DNA complexes containing variable sequences in the noncontacted center of the binding site or variations within the major groove interfaces show that the DNA can adapt to the Fis dimer surface asymmetrically. We show that the presence and position of pyrimidine-purine base steps within the major groove interfaces affect both local DNA bending and minor groove compression to modulate affinities and lifetimes of Fis-DNA complexes. Sequences flanking the core binding site also modulate complex affinities, lifetimes, and the degree of local and global Fis-induced DNA bending. In particular, a G immediately upstream of the 15 bp core sequence inhibits binding and bending, and A-tracts within the flanking base pairs increase both complex lifetimes and global DNA curvatures. Taken together, our observations support a revised DNA motif specifying high-affinity Fis binding and highlight the range of conformations that Fis-bound DNA can adopt. The affinities and DNA conformations of individual Fis-DNA complexes are likely to be tailored to their context-specific biological functions.

Understanding the host-adapted state of Citrobacter rodentium by transcriptomic analysis

Citrobacter rodentium (Cr) is a mouse pathogen that mimics many aspects of enteropathogenic Escherichia coli infections including producing attaching and effacing (A/E) lesions. Host-adapted (HA) Cr cells that are shed at the peak of infection have been reported to be hyper-infective. The exact mechanism underlying this phenomenon has remained elusive since the pathogen loses its HA 'status' immediately upon subculturing in laboratory media. We sequenced the entire transcriptome of Cr directly from the feces of infected mice and analyzed the gene expression pattern. We observed that the entire transcriptional machinery as well as several transcriptional regulators to be differentially expressed when compared with the transcriptome of cells grown on laboratory media. Major adhesion and effector genes, tir and eae, were highly expressed in HA along with many genes located on all five loci of enterocyte effacement regions (LEE 1-5). Notable absent among the HA expressed genes were 19 fimbrial operons and non-fimbrial adhesions and several non-LEE encoded effectors. These results demonstrate that host-adapted Cr has a unique transcriptome that is associated with increased host transmission.

Development of a 3-hydroxypropionate resistant Escherichia coli strain

3-hydroxypropionate (3HP) is an important platform chemical, and its biosynthesis is severely restricted by the toxicity of 3HP on cell growth and survival. To improve Escherichia coli resistance to 3HP and reduce the total production cost in industrial applications, we have identified variations in protein expression level exposed to sub-lethal concentration of this chemical using 2-dimensional gel electrophoresis. Under 3HP stress, the amount of 46 proteins was increased while the amount of 23 proteins was reduced. According to the proteomic results, overexpression of some identified proteins significantly increased the E. coli survival rate under 3HP stress. This study shed light on clues for developing E. coli strains with higher resistance to 3HP toxicity and lower production cost for industrial applications.

Reconstruction and Use of Microbial Metabolic Networks: the Core Escherichia coli Metabolic Model as an Educational Guide

Biochemical network reconstructions have become popular tools in systems biology. Metabolicnetwork reconstructions are biochemically, genetically, and genomically (BiGG) structured databases of biochemical reactions and metabolites. They contain information such as exact reaction stoichiometry, reaction reversibility, and the relationships between genes, proteins, and reactions. Network reconstructions have been used extensively to study the phenotypic behavior of wild-type and mutant stains under a variety of conditions, linking genotypes with phenotypes. Such phenotypic simulations have allowed for the prediction of growth after genetic manipulations, prediction of growth phenotypes after adaptive evolution, and prediction of essential genes. Additionally, because network reconstructions are organism specific, they can be used to understand differences between organisms of species in a functional context.There are different types of reconstructions representing various types of biological networks (metabolic, regulatory, transcription/translation). This chapter serves as an introduction to metabolic and regulatory network reconstructions and models and gives a complete description of the core Escherichia coli metabolic model. This model can be analyzed in any computational format (such as MATLAB or Mathematica) based on the information given in this chapter. The core E. coli model is a small-scale model that can be used for educational purposes. It is meant to be used by senior undergraduate and first-year graduate students learning about constraint-based modeling and systems biology. This model has enough reactions and pathways to enable interesting and insightful calculations, but it is also simple enough that the results of such calculations can be understoodeasily.

The Nucleoid: an Overview

This review provides a brief review of the current understanding of the structure-function relationship of the Escherichia coli nucleoid developed after the overview by Pettijohn focusing on the physical properties of nucleoids. Isolation of nucleoids requires suppression of DNA expansion by various procedures. The ability to control the expansion of nucleoids in vitro has led to purification of nucleoids for chemical and physical analyses and for high-resolution imaging. Isolated E. coli genomes display a number of individually intertwined supercoiled loops emanating from a central core. Metabolic processes of the DNA double helix lead to three types of topological constraints that all cells must resolve to survive: linking number, catenates, and knots. The major species of nucleoid core protein share functional properties with eukaryotic histones forming chromatin; even the structures are different from histones. Eukaryotic histones play dynamic roles in the remodeling of eukaryotic chromatin, thereby controlling the access of RNA polymerase and transcription factors to promoters. The E. coli genome is tightly packed into the nucleoid, but, at each cell division, the genome must be faithfully replicated, divided, and segregated. Nucleoid activities such as transcription, replication, recombination, and repair are all affected by the structural properties and the special conformations of nucleoid. While it is apparent that much has been learned about the nucleoid, it is also evident that the fundamental interactions organizing the structure of DNA in the nucleoid still need to be clearly defined.

Reconstruction of novel transcription factor regulons through inference of their binding sites

Background

In most sequenced organisms the number of known regulatory genes (e.g., transcription factors (TFs)) vastly exceeds the number of experimentally-verified regulons that could be associated with them. At present, identification of TF regulons is mostly done through comparative genomics approaches. Such methods could miss organism-specific regulatory interactions and often require expensive and time-consuming experimental techniques to generate the underlying data.

Results

In this work, we present an efficient algorithm that aims to identify a given transcription factor’s regulon through inference of its unknown binding sites, based on the discovery of its binding motif. The proposed approach relies on computational methods that utilize gene expression data sets and knockout fitness data sets which are available or may be straightforwardly obtained for many organisms. We computationally constructed the profiles of putative regulons for the TFs LexA, PurR and Fur in E. coli K12 and identified their binding motifs. Comparisons with an experimentally-verified database showed high recovery rates of the known regulon members, and indicated good predictions for the newly found genes with high biological significance. The proposed approach is also applicable to novel organisms for predicting unknown regulons of the transcriptional regulators. Results for the hypothetical protein Dde0289 in D. alaskensis include the discovery of a Fis-type TF binding motif.

Conclusions

The proposed motif-based regulon inference approach can discover the organism-specific regulatory interactions on a single genome, which may be missed by current comparative genomics techniques due to their limitations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0685-y) contains supplementary material, which is available to authorized users.

Dynamic Transcriptional Regulation of Fis in Salmonella During the Exponential Phase

Fis is one of the most important global regulators and has attracted extensive research attention. Many studies have focused on comparing the Fis global regulatory networks for exploring Fis function during different growth stages, such as the exponential and stationary stages. Although the Fis protein in bacteria is mainly expressed in the exponential phase, the dynamic transcriptional regulation of Fis during the exponential phase remains poorly understood. To address this question, we used RNA-seq technology to identify the Fis-regulated genes in the S. enterica serovar Typhimurium during the early exponential phase, and qRT-PCR was performed to validate the transcriptional data. A total of 1495 Fis-regulated genes were successfully identified, including 987 Fis-repressed genes and 508 Fis-activated genes. Comparing the results of this study with those of our previous study, we found that the transcriptional regulation of Fis was diverse during the early- and mid-exponential phases. The results also showed that the strong positive regulation of Fis on Salmonella pathogenicity island genes in the mid-exponential phase transitioned into insignificant effect in the early exponential phase. To validate these results, we performed a cell infection assay and found that Δfis only exhibited a 1.49-fold decreased capacity compared with the LT2 wild-type strain, indicating a large difference from the 6.31-fold decrease observed in the mid-exponential phase. Our results provide strong evidence for a need to thoroughly understand the dynamic transcriptional regulation of Fis in Salmonella during the exponential phase.

Genome-wide Reconstruction of OxyR and SoxRS Transcriptional Regulatory Networks under Oxidative Stress in Escherichia coli K-12 MG1655

Three transcription factors (TFs), OxyR, SoxR, and SoxS, play a critical role in transcriptional regulation of the defense system for oxidative stress in bacteria. However, their full genome-wide regulatory potential is unknown. Here, we perform a genome-scale reconstruction of the OxyR, SoxR, and SoxS regulons in Escherichia coli K-12 MG1655. Integrative data analysis reveals that a total of 68 genes in 51 transcription units (TUs) belong to these regulons. Among them, 48 genes showed more than 2-fold changes in expression level under single-TF-knockout conditions. This reconstruction expands the genome-wide roles of these factors to include direct activation of genes related to amino acid biosynthesis (methionine and aromatic amino acids), cell wall synthesis (lipid A biosynthesis and peptidoglycan growth), and divalent metal ion transport (Mn(2+), Zn(2+), and Mg(2+)). Investigating the co-regulation of these genes with other stress-response TFs reveals that they are independently regulated by stress-specific TFs.

Modes of Escherichia coli Dps Interaction with DNA as Revealed by Atomic Force Microscopy

Multifunctional protein Dps plays an important role in iron assimilation and a crucial role in bacterial genome packaging. Its monomers form dodecameric spherical particles accumulating ~400 molecules of oxidized iron ions within the protein cavity and applying a flexible N-terminal ends of each subunit for interaction with DNA. Deposition of iron is a well-studied process by which cells remove toxic Fe²⁺ ions from the genetic material and store them in an easily accessible form. However, the mode of interaction with linear DNA remained mysterious and binary complexes with Dps have not been characterized so far. It is widely believed that Dps binds DNA without any sequence or structural preferences but several lines of evidence have demonstrated its ability to differentiate gene expression, which assumes certain specificity. Here we show that Dps has a different affinity for the two DNA fragments taken from the dps gene regulatory region. We found by atomic force microscopy that Dps predominantly occupies thermodynamically unstable ends of linear double-stranded DNA fragments and has high affinity to the central part of the branched DNA molecule self-assembled from three single-stranded oligonucleotides. It was proposed that Dps prefers binding to those regions in DNA that provide more contact pads for the triad of its DNA-binding bundle associated with one vertex of the protein globule. To our knowledge, this is the first study revealed the nucleoid protein with an affinity to branched DNA typical for genomic regions with direct and inverted repeats. As a ubiquitous feature of bacterial and eukaryotic genomes, such structural elements should be of particular care, but the protein system evolutionarily adapted for this function is not yet known, and we suggest Dps as a putative component of this system.

Nucleoid-associated proteins encoded on plasmids: Occurrence and mode of function

Nucleoid-associated proteins (NAPs) play a role in changing the shape of microbial DNA, making it more compact and affecting the regulation of transcriptional networks in host cells. Genes that encode NAPs include H-NS family proteins (H-NS, Ler, MvaT, BpH3, Bv3F, HvrA, and Lsr2), FIS, HU, IHF, Lrp, and NdpA, and are found in both microbial chromosomes and plasmid DNA. In the present study, NAP genes were distributed among 442 plasmids out of 4602 plasmid sequences, and many H-NS family proteins, and HU, IHF, Lrp, and NdpA were found in plasmids of Alpha-, Beta-, and Gammaproteobacteria, while HvrA, Lsr2, HU, and Lrp were found in other classes including Actinobacteria and Bacilli. Larger plasmids frequently carried multiple NAP genes. In addition, NAP genes were more frequently found in conjugative plasmids than non-transmissible plasmids. Several host cells carried the same types of H-NS family proteins on both their plasmids and chromosome(s), while this was not observed for other NAPs. Recent studies have shown that NAP genes on plasmids and chromosomes play important roles in the physical and regulatory integration of plasmids into the host cell.

The architecture of ArgR-DNA complexes at the genome-scale in Escherichia coli

DNA-binding motifs that are recognized by transcription factors (TFs) have been well studied; however, challenges remain in determining the in vivo architecture of TF-DNA complexes on a genome-scale. Here, we determined the in vivo architecture of Escherichia coli arginine repressor (ArgR)-DNA complexes using high-throughput sequencing of exonuclease-treated chromatin-immunoprecipitated DNA (ChIP-exo). The ChIP-exo has a unique peak-pair pattern indicating 5′ and 3′ ends of ArgR-binding region. We identified 62 ArgR-binding loci, which were classified into three groups, comprising single, double and triple peak-pairs. Each peak-pair has a unique 93 base pair (bp)-long (±2 bp) ArgR-binding sequence containing two ARG boxes (39 bp) and residual sequences. Moreover, the three ArgR-binding modes defined by the position of the two ARG boxes indicate that DNA bends centered between the pair of ARG boxes facilitate the non-specific contacts between ArgR subunits and the residual sequences. Additionally, our approach may also reveal other fundamental structural features of TF-DNA interactions that have implications for studying genome-scale transcriptional regulatory networks.

Gene regulation by H-NS as a function of growth conditions depends on chromosomal position in Escherichia coli

Cellular adaptation to changing environmental conditions requires the coordinated regulation of expression of large sets of genes by global regulatory factors such as nucleoid associated proteins. Although in eukaryotic cells genomic position is known to play an important role in regulation of gene expression, it remains to be established whether in bacterial cells there is an influence of chromosomal position on the efficiency of these global regulators. Here we show for the first time that genome position can affect transcription activity of a promoter regulated by the histone-like nucleoid-structuring protein (H-NS), a global regulator of bacterial transcription and genome organization. We have used as a local reporter of H-NS activity the level of expression of a fluorescent reporter protein under control of an H-NS−regulated promoter (Phns) at different sites along the genome. Our results show that the activity of the Phns promoter depends on whether it is placed within the AT-rich regions of the genome that are known to be bound preferentially by H-NS. This modulation of gene expression moreover depends on the growth phase and the growth rate of the cells, reflecting the changes taking place in the relative abundance of different nucleoid proteins and the inherent heterogeneous organization of the nucleoid. Genomic position can thus play a significant role in the adaptation of the cells to environmental changes, providing a fitness advantage that can explain the selection of a gene’s position during evolution.

Site-specific DNA Inversion by Serine Recombinases

Reversible site-specific DNA inversion reactions are widely distributed in bacteria and their viruses. They control a range of biological reactions that most often involve alterations of molecules on the surface of cells or phage. These programmed DNA rearrangements usually occur at a low frequency, thereby preadapting a small subset of the population to a change in environmental conditions, or in the case of phages, an expanded host range. A dedicated recombinase, sometimes with the aid of additional regulatory or DNA architectural proteins, catalyzes the inversion of DNA. RecA or other components of the general recombination-repair machinery are not involved. This chapter discusses site-specific DNA inversion reactions mediated by the serine recombinase family of enzymes and focuses on the extensively studied serine DNA invertases that are stringently controlled by the Fis-bound enhancer regulatory system. The first section summarizes biological features and general properties of inversion reactions by the Fis/enhancer-dependent serine invertases and the recently described serine DNA invertases in Bacteroides. Mechanistic studies of reactions catalyzed by the Hin and Gin invertases are then discussed in more depth, particularly with regards to recent advances in our understanding of the function of the Fis/enhancer regulatory system, the assembly of the active recombination complex (invertasome) containing the Fis/enhancer, and the process of DNA strand exchange by rotation of synapsed subunit pairs within the invertasome. The role of DNA topological forces that function in concert with the Fis/enhancer controlling element in specifying the overwhelming bias for DNA inversion over deletion and intermolecular recombination is emphasized.

Analysis of the Salmonella regulatory network suggests involvement of SsrB and H-NS in σ(E)-regulated SPI-2 gene expression

The extracytoplasmic functioning sigma factor σ^E is known to play an essential role for Salmonella enterica serovar Typhimurium to survive and proliferate in macrophages and mice. However, its regulatory network is not well-characterized, especially during infection. Here we used microarray to identify genes regulated by σ^E in Salmonella grown in three conditions: a nutrient-rich condition and two others that mimic early and late intracellular infection. We found that in each condition σ^E regulated different sets of genes, and notably, several global regulators. When comparing nutrient-rich and infection-like conditions, large changes were observed in the expression of genes involved in Salmonella pathogenesis island (SPI)-1 type-three secretion system (TTSS), SPI-2 TTSS, protein synthesis, and stress responses. In total, the expression of 58% of Salmonella genes was affected by σ^E in at least one of the three conditions. An important finding is that σ^E up-regulates SPI-2 genes, which are essential for Salmonella intracellular survival, by up-regulating SPI-2 activator ssrB expression at the early stage of infection and down-regulating SPI-2 repressor hns expression at a later stage. Moreover, σ^E is capable of countering the silencing of H-NS, releasing the expression of SPI-2 genes. This connection between σ^E and SPI-2 genes, combined with the global regulatory effect of σ^E, may account for the lethality of rpoE-deficient Salmonella in murine infection.

Global analysis of photosynthesis transcriptional regulatory networks

Photosynthesis is a crucial biological process that depends on the interplay of many components. This work analyzed the gene targets for 4 transcription factors: FnrL, PrrA, CrpK and MppG (RSP_2888), which are known or predicted to control photosynthesis in Rhodobacter sphaeroides. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) identified 52 operons under direct control of FnrL, illustrating its regulatory role in photosynthesis, iron homeostasis, nitrogen metabolism and regulation of sRNA synthesis. Using global gene expression analysis combined with ChIP-seq, we mapped the regulons of PrrA, CrpK and MppG. PrrA regulates ∼34 operons encoding mainly photosynthesis and electron transport functions, while CrpK, a previously uncharacterized Crp-family protein, regulates genes involved in photosynthesis and maintenance of iron homeostasis. Furthermore, CrpK and FnrL share similar DNA binding determinants, possibly explaining our observation of the ability of CrpK to partially compensate for the growth defects of a ΔFnrL mutant. We show that the Rrf2 family protein, MppG, plays an important role in photopigment biosynthesis, as part of an incoherent feed-forward loop with PrrA. Our results reveal a previously unrealized, high degree of combinatorial regulation of photosynthetic genes and significant cross-talk between their transcriptional regulators, while illustrating previously unidentified links between photosynthesis and the maintenance of iron homeostasis.

Photosynthetic organisms are among the most abundant life forms on earth. Their unique ability to harvest solar energy and use it to fix atmospheric carbon dioxide is at the foundation of the global food chain. This paper reports the first comprehensive analysis of networks that control expression of photosynthesis genes using Rhodobacter sphaeroides, a microbe that has been studied for decades as a model of solar energy capture and other aspects of the photosynthetic lifestyle. We find a previously unappreciated complexity in the level of control of photosynthetic genes, while identifying new links between photosynthesis and central processes like iron availability. This organism is an ancestor of modern day plants, so our data can inform studies in other photosynthetic organisms and improve our ability to harness solar energy for food and industrial processes.

Timely binding of IHF and Fis to DARS2 regulates ATP-DnaA production and replication initiation

In Escherichia coli, the ATP-bound form of DnaA (ATP-DnaA) promotes replication initiation. During replication, the bound ATP is hydrolyzed to ADP to yield the ADP-bound form (ADP-DnaA), which is inactive for initiation. The chromosomal site DARS2 facilitates the regeneration of ATP-DnaA by catalyzing nucleotide exchange between free ATP and ADP bound to DnaA. However, the regulatory mechanisms governing this exchange reaction are unclear. Here, using in vitro reconstituted experiments, we show that two nucleoid-associated proteins, IHF and Fis, bind site-specifically to DARS2 to activate coordinately the exchange reaction. The regenerated ATP-DnaA was fully active in replication initiation and underwent DnaA-ATP hydrolysis. ADP-DnaA formed heteromultimeric complexes with IHF and Fis on DARS2, and underwent nucleotide dissociation more efficiently than ATP-DnaA. Consistently, mutant analyses demonstrated that specific binding of IHF and Fis to DARS2 stimulates the formation of ATP-DnaA production, thereby promoting timely initiation. Moreover, we show that IHF-DARS2 binding is temporally regulated during the cell cycle, whereas Fis only binds to DARS2 in exponentially growing cells. These results elucidate the regulation of ATP-DnaA and replication initiation in coordination with the cell cycle and growth phase.

Deciphering Fur transcriptional regulatory network highlights its complex role beyond iron metabolism in Escherichia coli

The ferric uptake regulator (Fur) plays a critical role in the transcriptional regulation of iron metabolism. However, the full regulatory potential of Fur remains undefined. Here we comprehensively reconstruct the Fur transcriptional regulatory network in Escherichia coli K-12 MG1655 in response to iron availability using genome-wide measurements (ChIP-exo and RNA-seq). Integrative data analysis reveals that a total of 81 genes in 42 transcription units are directly regulated by three different modes of Fur regulation, including apo- and holo-Fur activation and holo-Fur repression. We show that Fur connects iron transport and utilization enzymes with negative-feedback loop pairs for iron homeostasis. In addition, direct involvement of Fur in the regulation of DNA synthesis, energy metabolism, and biofilm development is found. These results show how Fur exhibits a comprehensive regulatory role affecting many fundamental cellular processes linked to iron metabolism in order to coordinate the overall response of E. coli to iron availability.

Two Fis regulators directly repress the expression of numerous effector-encoding genes in Legionella pneumophila

Legionella pneumophila is an intracellular human pathogen that utilizes the Icm/Dot type IVB secretion system to translocate a large repertoire of effectors into host cells. For most of these effectors, there is no information regarding their regulation. Therefore, the aim of this study was to examine the involvement of the three L. pneumophila Fis homologs in the regulation of effector-encoding genes. Deletion mutants constructed in the genes encoding the three Fis regulators revealed that Fis1 (lpg0542 gene) and Fis3 (lpg1743) but not Fis2 (lpg1370) are partially required for intracellular growth of L. pneumophila in Acanthamoeba castellanii. To identify pathogenesis-related genes directly regulated by Fis, we established a novel in vivo system which resulted in the discovery of numerous effector-encoding genes directly regulated by Fis. Further examination of these genes revealed that Fis1 and Fis3 repress the level of expression of effector-encoding genes during exponential phase. Three groups of effector-encoding genes were identified: (i) effectors regulated mainly by Fis1, (ii) effectors regulated mainly by Fis3, and (iii) effectors regulated by both Fis1 and Fis3. Examination of the upstream regulatory region of all of these effector-encoding genes revealed multiple putative Fis regulatory elements, and site-directed mutagenesis confirmed that a few of these sites constitute part of a repressor binding element. Furthermore, gel mobility shift assays demonstrated the direct relation between the Fis1 and Fis3 regulators and these regulatory elements. Collectively, our results demonstrate for the first time that two of the three L. pneumophila Fis regulators directly repress the expression of Icm/Dot effector-encoding genes.

Chromosome position effects on gene expression in Escherichia coli K-12

In eukaryotes, the location of a gene on the chromosome is known to affect its expression, but such position effects are poorly understood in bacteria. Here, using Escherichia coli K-12, we demonstrate that expression of a reporter gene cassette, comprised of the model E. coli lac promoter driving expression of gfp, varies by ∼300-fold depending on its precise position on the chromosome. At some positions, expression was more than 3-fold higher than at the natural lac promoter locus, whereas at several other locations, the reporter cassette was completely silenced: effectively overriding local lac promoter control. These effects were not due to differences in gene copy number, caused by partially replicated genomes. Rather, the differences in gene expression occur predominantly at the level of transcription and are mediated by several different features that are involved in chromosome organization. Taken together, our findings identify a tier of gene regulation above local promoter control and highlight the importance of chromosome position effects on gene expression profiles in bacteria.

Intracellular concentrations of 65 species of transcription factors with known regulatory functions in Escherichia coli

The expression pattern of the Escherichia coli genome is controlled in part by regulating the utilization of a limited number of RNA polymerases among a total of its approximately 4,600 genes. The distribution pattern of RNA polymerase changes from modulation of two types of protein-protein interactions: the interaction of core RNA polymerase with seven species of the sigma subunit for differential promoter recognition and the interaction of RNA polymerase holoenzyme with about 300 different species of transcription factors (TFs) with regulatory functions. We have been involved in the systematic search for the target promoters recognized by each sigma factor and each TF using the newly developed Genomic SELEX system. In parallel, we developed the promoter-specific (PS)-TF screening system for identification of the whole set of TFs involved in regulation of each promoter. Understanding the regulation of genome transcription also requires knowing the intracellular concentrations of the sigma subunits and TFs under various growth conditions. This report describes the intracellular levels of 65 species of TF with known function in E. coli K-12 W3110 at various phases of cell growth and at various temperatures. The list of intracellular concentrations of the sigma factors and TFs provides a community resource for understanding the transcription regulation of E. coli under various stressful conditions in nature.

Beyond DnaA: the role of DNA topology and DNA methylation in bacterial replication initiation

The replication of chromosomal DNA is a fundamental event in the life cycle of every cell. The first step of replication, initiation, is controlled by multiple factors to ensure only one round of replication per cell cycle. The process of initiation has been described most thoroughly for bacteria, especially Escherichia coli, and involves many regulatory proteins that vary considerably between different species. These proteins control the activity of the two key players of initiation in bacteria: the initiator protein DnaA and the origin of chromosome replication (oriC). Factors involved in the control of the availability, activity, or oligomerization of DnaA during initiation are generally regarded as the most important and thus have been thoroughly characterized. Other aspects of the initiation process, such as origin accessibility and susceptibility to unwinding, have been less explored. However, recent findings indicate that these factors have a significant role. This review focuses on DNA topology, conformation, and methylation as important factors that regulate the initiation process in bacteria. We present a comprehensive summary of the factors involved in the modulation of DNA topology, both locally at oriC and more globally at the level of the entire chromosome. We show clearly that the conformation of oriC dynamically changes, and control of this conformation constitutes another, important factor in the regulation of bacterial replication initiation. Furthermore, the process of initiation appears to be associated with the dynamics of the entire chromosome and this association is an important but largely unexplored phenomenon.

Promoter islands as a platform for interaction with nucleoid proteins and transcription factors

Seventy-eight promoter islands with an extraordinarily high density of potential promoters have been recently found in the genome of Escherichia coli. It has been shown that RNA polymerase binds internal promoters of these islands and produces short oligonucleotides, while the synthesis of normal mRNAs is suppressed. This quenching may be biologically relevant, as most islands are associated with foreign genes, which expression may deplete cellular resources. However, a molecular mechanism of silencing with the participation of these promoter-rich regions remains obscure. It has been demonstrated that all islands interact with histone-like protein H-NS--a specific sentinel of foreign genes. In this study, we demonstrated the inhibitory effect of H-NS using Δhns mutant of Escherichia coli and showed that deletion of dps, encoding another protein of bacterial nucleoid, tended to decrease rather than increase the amount of island-specific transcripts. This observation precluded consideration of promoter islands as sites for targeted heterochromatization only and a computer search for the binding sites of 53 transcription factors (TFs) revealed six proteins, which may specifically regulate their transcriptional output.

Investigating the interplay between nucleoid-associated proteins, DNA curvature, and CRISPR elements using comparative genomics

Many prokaryotic and eukaryotic genomes feature a characteristic periodic signal in distribution of short runs of A or T (A-tracts) phased with the DNA helical period of ∼10-11 bp. Such periodic spacing of A-tracts has been associated with intrinsic DNA curvature. In eukaryotes, this periodicity is a major component of the nucleosome positioning signal but its physiological role in prokaryotes is not clear. One hypothesis centers on possible role of intrinsic DNA bends in nucleoid compaction. We use comparative genomics to investigate possible relationship between the A-tract periodicity and nucleoid-associated proteins in prokaryotes. We found that genomes with DNA-bridging proteins tend to exhibit stronger A-tract periodicity, presumably indicative of more prevalent intrinsic DNA curvature. A weaker relationship was detected for nucleoid-associated proteins that do not form DNA bridges. We consider these results an indication that intrinsic DNA curvature acts collaboratively with DNA-bridging proteins in maintaining the compact structure of the nucleoid, and that previously observed differences among prokaryotic genomes in terms DNA curvature-related sequence periodicity may reflect differences in nucleoid organization. We subsequently investigated the relationship between A-tract periodicity and presence of CRISPR elements and we found that genomes with CRISPR tend to have stronger A-tract periodicity. This result is consistent with our earlier hypothesis that extensive A-tract periodicity could help protect the chromosome against integration of prophages, possibly due to its role in compaction of the nucleoid.

The iron stimulon and fur regulon of Geobacter sulfurreducens and their role in energy metabolism

Iron plays a critical role in the physiology of Geobacter species. It serves as both an essential component for proteins and cofactors and an electron acceptor during anaerobic respiration. Here, we investigated the iron stimulon and ferric uptake regulator (Fur) regulon of Geobacter sulfurreducens to examine the coordination between uptake of Fe(II) and the reduction of Fe(III) at the transcriptional level. Gene expression studies across a variety of different iron concentrations in both the wild type and a Δfur mutant strain were used to determine the iron stimulon. The stimulon consists of a broad range of gene products, ranging from iron-utilizing to central metabolism and iron reduction proteins. Integration of gene expression and chromatin immunoprecipitation (ChIP) data sets assisted in the identification of the Fur transcriptional regulatory network and Fur's role as a regulator of the iron stimulon. Additional physiological and transcriptional analyses of G. sulfurreducens grown with various Fe(II) concentrations revealed the depth of Fur's involvement in energy metabolism and the existence of redundancy within the iron-regulatory network represented by IdeR, an alternative iron transcriptional regulator. These characteristics enable G. sulfurreducens to thrive in environments with fluctuating iron concentrations by providing it with a robust mechanism to maintain tight and deliberate control over intracellular iron homeostasis.

Genome-scale reconstruction of the sigma factor network in Escherichia coli: topology and functional states

Background

At the beginning of the transcription process, the RNA polymerase (RNAP) core enzyme requires a σ-factor to recognize the genomic location at which the process initiates. Although the crucial role of σ-factors has long been appreciated and characterized for many individual promoters, we do not yet have a genome-scale assessment of their function.

Results

Using multiple genome-scale measurements, we elucidated the network of σ-factor and promoter interactions in Escherichia coli. The reconstructed network includes 4,724 σ-factor-specific promoters corresponding to transcription units (TUs), representing an increase of more than 300% over what has been previously reported. The reconstructed network was used to investigate competition between alternative σ-factors (the σ⁷⁰ and σ³⁸ regulons), confirming the competition model of σ substitution and negative regulation by alternative σ-factors. Comparison with σ-factor binding in Klebsiella pneumoniae showed that transcriptional regulation of conserved genes in closely related species is unexpectedly divergent.

Conclusions

The reconstructed network reveals the regulatory complexity of the promoter architecture in prokaryotic genomes, and opens a path to the direct determination of the systems biology of their transcriptional regulatory networks.

Shared control of gene expression in bacteria by transcription factors and global physiology of the cell

A simple, parameterless mathematical model, in combination with real-time monitoring of promoter activities, shows how control of gene expression in bacteria is shared between transcription factors and global physiological effects.

We present an approach based on a simple, paramaterless mathematical model to analyze the control of gene expression by transcription factors and the global physiological state of the cell.

We illustrate the strength of this approach by means of time-resolved measurements of the transcriptional activities of genes in a central regulatory circuit in Escherichia coli.

We conclude that global physiological effects rather than transcription factors dominate the control of gene expression during a growth transition.

Our results call for a reappraisal of the role of transcription factors, which may be most appropriately viewed as complementing and finetuning global control exerted by the physiological state of the cell.

Gene expression is controlled by the joint effect of (i) the global physiological state of the cell, in particular the activity of the gene expression machinery, and (ii) DNA-binding transcription factors and other specific regulators. We present a model-based approach to distinguish between these two effects using time-resolved measurements of promoter activities. We demonstrate the strength of the approach by analyzing a circuit involved in the regulation of carbon metabolism in E. coli. Our results show that the transcriptional response of the network is controlled by the physiological state of the cell and the signaling metabolite cyclic AMP (cAMP). The absence of a strong regulatory effect of transcription factors suggests that they are not the main coordinators of gene expression changes during growth transitions, but rather that they complement the effect of global physiological control mechanisms. This change of perspective has important consequences for the interpretation of transcriptome data and the design of biological networks in biotechnology and synthetic biology.

Comparative genomics and transcriptomics analysis of experimentally evolved Escherichia coli MC1000 in complex environments

It has recently become feasible to study the basis and nature of evolutionary changes in bacteria in an experimental setting using defined media. However, assessment of adaptive changes in complex environments has been scarce. In an effort to describe the responses in such environments, we unravel, in a comparative approach, the transcriptional and genetic profiles of 19 Escherichia coli strains that evolved in Luria Bertani medium under three different oxygen regimes over 1000 generations. A positive relationship between upregulation of gene expression and the number of mutations was observed, suggesting that a number of metabolic pathways were activated. Phenotypic polymorphisms were observed in parallel cultures, of which some were related with mutations at the regulatory level. Non-parallel responses were observed at the intrapopulational level, which is indicative of diversifying selection. Parallel responses encompassed transcriptome diversity, and their effects were directly affected by differing genomic backgrounds. A fluctuating selective force produced higher phenotypic diversity compared with constant forces. This study demonstrates how phenotypic innovations may depend on the relationship between genomic changes and local ecological conditions. Using both comparative genomics and transcriptomics approaches, the results help elucidating various adaptive responses in cultures in unexplored complex environments.

Genome-scale analysis of escherichia coli FNR reveals complex features of transcription factor binding

FNR is a well-studied global regulator of anaerobiosis, which is widely conserved across bacteria. Despite the importance of FNR and anaerobiosis in microbial lifestyles, the factors that influence its function on a genome-wide scale are poorly understood. Here, we report a functional genomic analysis of FNR action. We find that FNR occupancy at many target sites is strongly influenced by nucleoid-associated proteins (NAPs) that restrict access to many FNR binding sites. At a genome-wide level, only a subset of predicted FNR binding sites were bound under anaerobic fermentative conditions and many appeared to be masked by the NAPs H-NS, IHF and Fis. Similar assays in cells lacking H-NS and its paralog StpA showed increased FNR occupancy at sites bound by H-NS in WT strains, indicating that large regions of the genome are not readily accessible for FNR binding. Genome accessibility may also explain our finding that genome-wide FNR occupancy did not correlate with the match to consensus at binding sites, suggesting that significant variation in ChIP signal was attributable to cross-linking or immunoprecipitation efficiency rather than differences in binding affinities for FNR sites. Correlation of FNR ChIP-seq peaks with transcriptomic data showed that less than half of the FNR-regulated operons could be attributed to direct FNR binding. Conversely, FNR bound some promoters without regulating expression presumably requiring changes in activity of condition-specific transcription factors. Such combinatorial regulation may allow Escherichia coli to respond rapidly to environmental changes and confer an ecological advantage in the anaerobic but nutrient-fluctuating environment of the mammalian gut.

Regulation of gene expression by transcription factors (TFs) is key to adaptation to environmental changes. Our comprehensive, genome-scale analysis of a prototypical global TF, the anaerobic regulator FNR from Escherichia coli, leads to several novel and unanticipated insights into the influences on FNR binding genome-wide and the complex structure of bacterial regulons. We found that binding of NAPs restricts FNR binding at a subset of sites, suggesting that the bacterial genome is not freely accessible for FNR binding. Our finding that less than half of the predicted FNR binding sites were occupied in vivo further challenges the utility of using bioinformatic searches alone to predict regulon structure, reinforcing the need for experimental determination of TF binding. By correlating the occupancy data with transcriptomic data, we confirm that FNR serves as a global signal of anaerobiosis but expression of some operons in the FNR regulon require other regulators sensitive to alternative environmental stimuli. Thus, FNR binding and regulation appear to depend on both the nucleoprotein structure of the chromosome and on combinatorial binding of FNR with other regulators. Both of these phenomena are typical of TF binding in eukaryotes; our results establish that they are also features of bacterial TF binding.

Variation of the folding and dynamics of the Escherichia coli chromosome with growth conditions

We examine whether the Escherichia coli chromosome is folded into a self-adherent nucleoprotein complex, or alternately is a confined but otherwise unconstrained self-avoiding polymer. We address this through in vivo visualization, using an inducible GFP fusion to the nucleoid-associated protein Fis to non-specifically decorate the entire chromosome. For a range of different growth conditions, the chromosome is a compact structure that does not fill the volume of the cell, and which moves from the new pole to the cell centre. During rapid growth, chromosome segregation occurs well before cell division, with daughter chromosomes coupled by a thin inter-daughter filament before complete segregation, whereas during slow growth chromosomes stay adjacent until cell division occurs. Image correlation analysis indicates that sub-nucleoid structure is stable on a 1 min timescale, comparable to the timescale for redistribution time measured for GFP-Fis after photobleaching. Optical deconvolution and writhe calculation analysis indicate that the nucleoid has a large-scale coiled organization rather than being an amorphous mass. Our observations are consistent with the chromosome having a self-adherent filament organization.

Genome-wide analysis of the salmonella Fis regulon and its regulatory mechanism on pathogenicity islands

Fis, one of the most important nucleoid-associated proteins, functions as a global regulator of transcription in bacteria that has been comprehensively studied in Escherichia coli K12. Fis also influences the virulence of Salmonella enterica and pathogenic E. coli by regulating their virulence genes, however, the relevant mechanism is unclear. In this report, using combined RNA-seq and chromatin immunoprecipitation (ChIP)-seq technologies, we first identified 1646 Fis-regulated genes and 885 Fis-binding targets in the S. enterica serovar Typhimurium, and found a Fis regulon different from that in E. coli. Fis has been reported to contribute to the invasion ability of S. enterica. By using cell infection assays, we found it also enhances the intracellular replication ability of S. enterica within macrophage cell, which is of central importance for the pathogenesis of infections. Salmonella pathogenicity islands (SPI)-1 and SPI-2 are crucial for the invasion and survival of S. enterica in host cells. Using mutation and overexpression experiments, real-time PCR analysis, and electrophoretic mobility shift assays, we demonstrated that Fis regulates 63 of the 94 Salmonella pathogenicity island (SPI)-1 and SPI-2 genes, by three regulatory modes: i) binds to SPI regulators in the gene body or in upstream regions; ii) binds to SPI genes directly to mediate transcriptional activation of themselves and downstream genes; iii) binds to gene encoding OmpR which affects SPI gene expression by controlling SPI regulators SsrA and HilD. Our results provide new insights into the impact of Fis on SPI genes and the pathogenicity of S. enterica.

Promoters of Escherichia coli versus promoter islands: function and structure comparison

Expression of bacterial genes takes place under the control of RNA polymerase with exchangeable σ-subunits and multiple transcription factors. A typical promoter region contains one or several overlapping promoters. In the latter case promoters have the same or different σ-specificity and are often subjected to different regulatory stimuli. Genes, transcribed from multiple promoters, have on average higher expression levels. However, recently in the genome of Escherichia coli we found 78 regions with an extremely large number of potential transcription start points (promoter islands, PIs). It was shown that all PIs interact with RNA polymerase in vivo and are able to form transcriptionally competent open complexes both in vitro and in vivo but their transcriptional activity measured by oligonucleotide microarrays was very low, if any. Here we confirmed transcriptional defectiveness of PIs by analyzing the 5'-end specific RNA-seq data, but showed their ability to produce short oligos (9-14 bases). This combination of functional properties indicated a deliberate suppression of transcriptional activity within PIs. According to our data this suppression may be due to a specific conformation of the DNA double helix, which provides an ideal platform for interaction with both RNA polymerase and the histone-like nucleoid protein H-NS. The genomic DNA of E.coli contains therefore several dozen sites optimized by evolution for staying in a heterochromatin-like state. Since almost all promoter islands are associated with horizontally acquired genes, we offer them as specific components of bacterial evolution involved in acquisition of foreign genetic material by turning off the expression of toxic or useless aliens or by providing optimal promoter for beneficial genes. The putative molecular mechanism underlying the appearance of promoter islands within recipient genomes is discussed.

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.

E. coli Fis protein insulates the cbpA gene from uncontrolled transcription

The Escherichia coli curved DNA binding protein A (CbpA) is a poorly characterised nucleoid associated factor and co-chaperone. It is expressed at high levels as cells enter stationary phase. Using genetics, biochemistry, and genomics, we have examined regulation of, and DNA binding by, CbpA. We show that Fis, the dominant growth-phase nucleoid protein, prevents CbpA expression in growing cells. Regulation by Fis involves an unusual “insulation” mechanism. Thus, Fis protects cbpA from the effects of a distal promoter, located in an adjacent gene. In stationary phase, when Fis levels are low, CbpA binds the E. coli chromosome with a preference for the intrinsically curved Ter macrodomain. Disruption of the cbpA gene prompts dramatic changes in DNA topology. Thus, our work identifies a novel role for Fis and incorporates CbpA into the growing network of factors that mediate bacterial chromosome structure.

Compaction of chromosomal DNA is a fundamental process that impacts on all aspects of cellular biology. However, our understanding of chromosome organisation in bacteria is poorly developed. Since bacteria are amongst the most abundant living organisms on the planet, this represents a startling gap in our knowledge. Despite our lack of understanding, it has long been known that Escherichia coli, and other bacteria, radically re-model their chromosomes in response to environmental stress. This is most notable during periods of starvation, when the E. coli chromosome is super compacted. In dissecting the molecular mechanisms that control this phenomenon, we have found that regulatory cross-talk between DNA–organising proteins plays an essential role. Thus, the major DNA folding protein from growing E. coli inhibits production of the major chromosome organisers in starved cells. Our findings illustrate the highly dynamic nature of bacterial chromosomes. Thus, DNA topology, gene transcription, and chromosome folding proteins entwine to create a web of interactions that define the properties of the chromosome.

Binding of nucleoid-associated protein fis to DNA is regulated by DNA breathing dynamics

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.

Cellular transcription factors (TFs) are proteins that regulate gene expression, and thereby cellular activity and fate, by binding to specific DNA segments. The physicochemical determinants of protein-DNA binding specificity are not completely understood. Here, we report that the propensity of transient opening and re-closing of the double helix, resulting from thermal fluctuations, aka “DNA breathing” or “DNA bubbles,” can be associated with binding affinity in the case of Fis, a well-studied nucleoid-associated protein in Escherichia coli. We found that a particular breathing profile is characteristic of high-affinity Fis binding sites and that DNA fragments known to bind Fis in vivo are statistically enriched for this profile. Furthermore, we used simulations of DNA breathing dynamics to guide design of gel-shift experiments aimed at testing the idea that local breathing influences Fis binding. As a result, we show that via nucleotide modifications but without modifying nucleotides that directly contact Fis, we were able to transform a low-affinity Fis binding site into a high-affinity site and vice versa. The nucleotide modifications were designed only based on DNA breathing simulations. Our study suggests that strong Fis-DNA binding depends on DNA breathing - a novel physicochemical characteristic that could be used for prediction and rational design of TF binding sites.

DnaA binding locus datA promotes DnaA-ATP hydrolysis to enable cell cycle-coordinated replication initiation

The initiation of chromosomal DNA replication is rigidly regulated to ensure that it occurs in a cell cycle-coordinated manner. To ensure this in Escherichia coli, multiple systems regulate the activity of the replication initiator ATP-DnaA. The level of ATP-DnaA increases before initiation after which it drops via DnaA-ATP hydrolysis, yielding initiation-inactive ADP-DnaA. DnaA-ATP hydrolysis is crucial to regulation of initiation and mainly occurs by a replication-coupled feedback mechanism named RIDA (regulatory inactivation of DnaA). Here, we report a second DnaA-ATP hydrolysis system that occurs at the chromosomal site datA. This locus has been annotated as a reservoir for DnaA that binds many DnaA molecules in a manner dependent upon the nucleoid-associated factor IHF (integration host factor), resulting in repression of untimely initiations; however, there is no direct evidence for the binding of many DnaA molecules at this locus. We reveal that a complex consisting of datA and IHF promotes DnaA-ATP hydrolysis in a manner dependent on specific inter-DnaA interactions. Deletion of datA or the ihf gene increased ATP-DnaA levels to the maximal attainable levels in RIDA-defective cells. Cell-cycle analysis suggested that IHF binds to datA just after replication initiation at a time when RIDA is activated. We propose a model in which cell cycle-coordinated ATP-DnaA inactivation is regulated in a concerted manner by RIDA and datA.

Improved predictions of transcription factor binding sites using physicochemical features of DNA

Typical approaches for predicting transcription factor binding sites (TFBSs) involve use of a position-specific weight matrix (PWM) to statistically characterize the sequences of the known sites. Recently, an alternative physicochemical approach, called SiteSleuth, was proposed. In this approach, a linear support vector machine (SVM) classifier is trained to distinguish TFBSs from background sequences based on local chemical and structural features of DNA. SiteSleuth appears to generally perform better than PWM-based methods. Here, we improve the SiteSleuth approach by considering both new physicochemical features and algorithmic modifications. New features are derived from Gibbs energies of amino acid-DNA interactions and hydroxyl radical cleavage profiles of DNA. Algorithmic modifications consist of inclusion of a feature selection step, use of a nonlinear kernel in the SVM classifier, and use of a consensus-based post-processing step for predictions. We also considered SVM classification based on letter features alone to distinguish performance gains from use of SVM-based models versus use of physicochemical features. The accuracy of each of the variant methods considered was assessed by cross validation using data available in the RegulonDB database for 54 Escherichia coli TFs, as well as by experimental validation using published ChIP-chip data available for Fis and Lrp.

Multiple-omic data analysis of Klebsiella pneumoniae MGH 78578 reveals its transcriptional architecture and regulatory features

Background

The increasing number of infections caused by strains of Klebsiella pneumoniae that are resistant to multiple antibiotics has developed into a major medical problem worldwide. The development of next-generation sequencing technologies now permits rapid sequencing of many K. pneumoniae isolates, but sequence information alone does not provide important structural and operational information for its genome.

Results

Here we take a systems biology approach to annotate the K. pneumoniae MGH 78578 genome at the structural and operational levels. Through the acquisition and simultaneous analysis of multiple sample-matched –omics data sets from two growth conditions, we detected 2677, 1227, and 1066 binding sites for RNA polymerase, RpoD, and RpoS, respectively, 3660 RNA polymerase-guided transcript segments, and 3585 transcription start sites throughout the genome. Moreover, analysis of the transcription start site data identified 83 probable leaderless mRNAs, while analysis of unannotated transcripts suggested the presence of 119 putative open reading frames, 15 small RNAs, and 185 antisense transcripts that are not currently annotated.

Conclusions

These findings highlight the strengths of systems biology approaches to the refinement of sequence-based annotations, and to provide new insight into fundamental genome-level biology for this important human pathogen.

Genome-wide PhoB binding and gene expression profiles reveal the hierarchical gene regulatory network of phosphate starvation in Escherichia coli

The phosphate starvation response in bacteria has been studied extensively for the past few decades and the phosphate-limiting signal is known to be mediated via the PhoBR two-component system. However, the global DNA binding profile of the response regulator PhoB and the PhoB downstream responses are currently unclear. In this study, chromatin immunoprecipitation for PhoB was combined with high-density tiling array (ChIP-chip) as well as gene expression microarray to reveal the first global down-stream responses of the responding regulator, PhoB in E. coli. Based on our ChIP-chip experimental data, forty-three binding sites were identified throughout the genome and the known PhoB binding pattern was updated by identifying the conserved pattern from these sites. From the gene expression microarray data analysis, 287 differentially expressed genes were identified in the presence of PhoB activity. By comparing the results obtained from our ChIP-chip and microarray experiments, we were also able to identify genes that were directly or indirectly affected through PhoB regulation. Nineteen out of these 287 differentially expressed genes were identified as the genes directly regulated by PhoB. Seven of the 19 directly regulated genes (including phoB) are transcriptional regulators. These transcriptional regulators then further pass the signal of phosphate starvation down to the remaining differentially expressed genes. Our results unveiled the genome-wide binding profile of PhoB and the downstream responses under phosphate starvation. We also present the hierarchical structure of the phosphate sensing regulatory network. The data suggest that PhoB plays protective roles in membrane integrity and oxidative stress reduction during phosphate starvation.

A flexible integrative approach based on random forest improves prediction of transcription factor binding sites

Transcription factor binding sites (TFBSs) are DNA sequences of 6–15 base pairs. Interaction of these TFBSs with transcription factors (TFs) is largely responsible for most spatiotemporal gene expression patterns. Here, we evaluate to what extent sequence-based prediction of TFBSs can be improved by taking into account the positional dependencies of nucleotides (NPDs) and the nucleotide sequence-dependent structure of DNA. We make use of the random forest algorithm to flexibly exploit both types of information. Results in this study show that both the structural method and the NPD method can be valuable for the prediction of TFBSs. Moreover, their predictive values seem to be complementary, even to the widely used position weight matrix (PWM) method. This led us to combine all three methods. Results obtained for five eukaryotic TFs with different DNA-binding domains show that our method improves classification accuracy for all five eukaryotic TFs compared with other approaches. Additionally, we contrast the results of seven smaller prokaryotic sets with high-quality data and show that with the use of high-quality data we can significantly improve prediction performance. Models developed in this study can be of great use for gaining insight into the mechanisms of TF binding.

Genomic analysis of DNA binding and gene regulation by homologous nucleoid-associated proteins IHF and HU in Escherichia coli K12

IHF and HU are two heterodimeric nucleoid-associated proteins (NAP) that belong to the same protein family but interact differently with the DNA. IHF is a sequence-specific DNA-binding protein that bends the DNA by over 160°. HU is the most conserved NAP, which binds non-specifically to duplex DNA with a particular preference for targeting nicked and bent DNA. Despite their importance, the in vivo interactions of the two proteins to the DNA remain to be described at a high resolution and on a genome-wide scale. Further, the effects of these proteins on gene expression on a global scale remain contentious. Finally, the contrast between the functions of the homo- and heterodimeric forms of proteins deserves the attention of further study. Here we present a genome-scale study of HU- and IHF binding to the Escherichia coli K12 chromosome using ChIP-seq. We also perform microarray analysis of gene expression in single- and double-deletion mutants of each protein to identify their regulons. The sequence-specific binding profile of IHF encompasses ∼30% of all operons, though the expression of <10% of these is affected by its deletion suggesting combinatorial control or a molecular backup. The binding profile for HU is reflective of relatively non-specific binding to the chromosome, however, with a preference for A/T-rich DNA. The HU regulon comprises highly conserved genes including those that are essential and possibly supercoiling sensitive. Finally, by performing ChIP-seq experiments, where possible, of each subunit of IHF and HU in the absence of the other subunit, we define genome-wide maps of DNA binding of the proteins in their hetero- and homodimeric forms.

A fundamental regulatory mechanism operating through OmpR and DNA topology controls expression of Salmonella pathogenicity islands SPI-1 and SPI-2

DNA topology has fundamental control over the ability of transcription factors to access their target DNA sites at gene promoters. However, the influence of DNA topology on protein–DNA and protein–protein interactions is poorly understood. For example, relaxation of DNA supercoiling strongly induces the well-studied pathogenicity gene ssrA (also called spiR) in Salmonella enterica, but neither the mechanism nor the proteins involved are known. We have found that relaxation of DNA supercoiling induces expression of the Salmonella pathogenicity island (SPI)-2 regulator ssrA as well as the SPI-1 regulator hilC through a mechanism that requires the two-component regulator OmpR-EnvZ. Additionally, the ompR promoter is autoregulated in the same fashion. Conversely, the SPI-1 regulator hilD is induced by DNA relaxation but is repressed by OmpR. Relaxation of DNA supercoiling caused an increase in OmpR binding to DNA and a concomitant decrease in binding by the nucleoid-associated protein FIS. The reciprocal occupancy of DNA by OmpR and FIS was not due to antagonism between these transcription factors, but was instead a more intrinsic response to altered DNA topology. Surprisingly, DNA relaxation had no detectable effect on the binding of the global repressor H-NS. These results reveal the underlying molecular mechanism that primes SPI genes for rapid induction at the onset of host invasion. Additionally, our results reveal novel features of the archetypal two-component regulator OmpR. OmpR binding to relaxed DNA appears to generate a locally supercoiled state, which may assist promoter activation by relocating supercoiling stress-induced destabilization of DNA strands. Much has been made of the mechanisms that have evolved to regulate horizontally-acquired genes such as SPIs, but parallels among the ssrA, hilC, and ompR promoters illustrate that a fundamental form of regulation based on DNA topology coordinates the expression of these genes regardless of their origins.

DNA is often considered to be a passive carrier of genetic information, but in fact DNA is an active participant in coordinating the expression of the genes it carries. This is because DNA is a dynamic molecule that can assume a wide range of topologies, and this has a direct impact on the formation of the protein–DNA complexes that drive gene expression. In a bacterium, the chromosome is supercoiled to variable levels according to environmental conditions, and supercoiling in turn governs the topology of gene promoters. Thus DNA supercoiling is able to transduce environmental signals to regulate promoter output. A previous study found that the intestinal pathogen Salmonella enterica may use changes in DNA supercoiling to detect when it has entered host immune cells, allowing the bacterium to induce the pathogenicity genes it requires to evade killing by macrophage. In dissecting the underlying molecular mechanisms, we have found that changes in DNA supercoiling also upregulate other key pathogenicity genes, and we have identified the proteins involved in this gene regulatory process. These findings indicate that a fundamental level of gene control arising from the interplay between protein transcription factors and DNA topology regulates Salmonella pathogenicity.

Architectural organization in E. coli nucleoid

In contrast to organized hierarchical structure of eukaryotic chromosome, bacterial chromosomes are believed not to have such structures. The genomes of bacteria are condensed into a compact structure called the nucleoid. Among many architectural, histone-like proteins which associate with the chromosomal DNA is HU which is implicated in folding DNA into a compact structure by bending and wrapping DNA. Unlike the majority of other histone-like proteins, HU is highly conserved in eubacteria and unique in its ability to bind RNA. Furthermore, an HU mutation profoundly alters the cellular transcription profile and consequently has global effects on physiology and the lifestyle of E. coli. Here we provide a short overview of the mechanisms by which the nucleoid is organized into different topological domains. We propose that HU is a major player in creating domain-specific superhelicities and thus influences the transcription profile from the constituent promoters. This article is part of a Special Issue entitled: Chromatin in time and space.