Genomic analysis of protein-DNA interactions in bacteria: insights into transcription and chromosome organization

Combined transcriptomic and ChIPseq analyses of the Bordetella pertussis RisA regulon

The regulation of Bordetella pertussis virulence is mediated by the two-component system BvgA/S, which activates the transcription of virulence-activated genes (vags). In the avirulent phase, the vags are not expressed, but instead, virulence-repressed genes (vrgs) are expressed, under the control of another two-component system, RisA/K. Here, we combined transcriptomic and chromatin immunoprecipitation sequencing (ChIPseq) data to examine the RisA/K regulon. We performed RNAseq analyses of RisA-deficient and RisA-phosphoablative B. pertussis mutants cultivated in virulent and avirulent conditions. We confirmed that the expression of most vrgs is regulated by phosphorylated RisA. However, the expression of some, including those involved in flagellum biosynthesis and chemotaxis, requires RisA independently of phosphorylation. Many RisA-regulated genes encode proteins with regulatory functions, suggesting multiple RisA regulation cascades. By ChIPseq analyses, we identified 430 RisA-binding sites, 208 within promoter regions, 201 within open reading frames, and 21 in non-coding regions. RisA binding was demonstrated in the promoter regions of most vrgs and, surprisingly, of some vags, as well as for other genes not identified as vags or vrgs. Unexpectedly, many genes, including some vags, like prn, brpL, bipA, and cyaA, contain a BvgA-binding site and a RisA-binding site, which increases the complexity of the RisAK/BvgAS network in B. pertussis virulence regulation.IMPORTANCEThe expression of virulence-activated genes (vags) of Bordetella pertussis, the etiological agent of whooping cough, is under the transcriptional control of the two-component system BvgA/S, which allows the bacterium to switch between virulent and avirulent phases. In addition, the more recently identified two-component system RisA/K is required for the expression of B. pertussis genes, collectively named vrgs, that are repressed during the virulent phase but activated during the avirulent phase. We have characterized the RisA/K regulon by combined transcriptomic and chromatin immunoprecipitation sequencing analyses. We identified more than 400 RisA-binding sites. Many of them are localized in promoter regions, especially vrgs, but some were found within open reading frames and in non-coding regions. Surprisingly, RisA-binding sites were also found in promoter regions of some vags, illustrating the previously underappreciated complexity of virulence regulation in B. pertussis.

Genomic SELEX Screening of Regulatory Targets of Transcription Factors

The genome of Escherichia coli K-12 is transcribed by a single species of RNA polymerase. The selectivity of transcriptional targets is determined via interaction with one of seven species of the sigma subunit and a total of approximately 300 species of transcription factor (TFs). For comprehensive identification of the regulatory targets of these two groups of regulatory proteins on the genome, we developed an in vitro approach, "Genomic SELEX" (gSELEX) screening. Here we describe a detailed protocol of the gSELEX screening system, which uses purified regulatory proteins and fragments of genomic DNA from E. coli. Moreover, we describe methods and examples of results using cell-free synthetic proteins.

ChIP-qPCR of FLAG-Tagged Proteins in Bacteria

DNA-protein interactions occur in biological processes such as genome replication, gene transcription, DNA repair, and chromatin compaction and organization. Mapping the distribution of the DNA-bound proteins on the chromosome is essential for understanding their associated biological process. Chromatin immunoprecipitation (ChIP) involves the antibody-mediated enrichment of DNA fragments bound by a target protein and has become one of the most powerful techniques for exploring the distribution of proteins on the chromosome. By incorporating quantitative polymerase chain reaction (qPCR) downstream of the ChIP assay, ChIP-qPCR was developed to describe binding profiles of DNA-associated proteins at a candidate locus. In this chapter, we describe ChIP-qPCR. We provide a step-by-step protocol for the preparation of a ChIP library of a 3× FLAG-tagged protein in bacteria, describe how downstream qPCR experiments can be performed with the appropriate controls, and explain how the data is analyzed. This chapter provides reliable technical guidance for ChIP-qPCR studies in bacteria.

Functional Characterization of TetR-like Transcriptional Regulator PA3973 from Pseudomonas aeruginosa

Pseudomonas aeruginosa, a human opportunistic pathogen, is a common cause of nosocomial infections. Its ability to survive under different conditions relies on a complex regulatory network engaging transcriptional regulators controlling metabolic pathways and capabilities to efficiently use the available resources. P. aeruginosa PA3973 encodes an uncharacterized TetR family transcriptional regulator. In this study, we applied a transcriptome profiling (RNA-seq), genome-wide identification of binding sites using ChIP-seq, as well as the phenotype analyses to unravel the biological role of PA3973. Transcriptional profiling of P. aeruginosa PAO1161 overexpressing PA3973 showed changes in the mRNA level of 648 genes. Concomitantly, ChIP-seq analysis identified more than 300 PA3973 binding sites in the P. aeruginosa genome. A 13 bp sequence motif was indicated as the binding site of PA3973. The PA3973 regulon encompasses the PA3972-PA3971 genes encoding a probable acyl-CoA dehydrogenase and a thioesterase. In vitro analysis showed PA3973 binding to PA3973p. Accordingly, the lack of PA3973 triggered increased expression of PA3972 and PA3971. The ∆PA3972-71 PAO1161 strain demonstrated impaired growth in the presence of stress-inducing agents hydroxylamine or hydroxyurea, thus suggesting the role of PA3972-71 in pathogen survival upon stress. Overall our results showed that TetR-type transcriptional regulator PA3973 has multiple binding sites in the P. aeruginosa genome and influences the expression of diverse genes, including PA3972-PA3971, encoding proteins with a proposed role in stress response.

Dissecting psa Locus Regulation in Yersinia pestis

The pH 6 antigen (PsaA) of Yersinia pestis is a virulence factor that is expressed in response to high temperature (37°C) and low pH (6.0). Previous studies have implicated the PsaE and PsaF regulators in the temperature- and pH-dependent regulation of psaA. Here, we show that PsaE levels are themselves controlled by pH and temperature, explaining the regulation of psaA. We identify hundreds of binding sites for PsaE across the Y. pestis genome, with the majority of binding sites located in intergenic regions bound by the nucleoid-associated protein H-NS. However, we detect direct regulation of only two transcripts by PsaE, likely due to displacement of H-NS from the corresponding promoter regions; our data suggest that most PsaE binding sites are nonregulatory or that they require additional environmental cues. We also identify the precise binding sites for PsaE that are required for temperature- and pH-dependent regulation of psaA and psaE. Thus, our data reveal the critical role that PsaE plays in the regulation of psaA and suggest that PsaE may have many additional regulatory targets. IMPORTANCE Y. pestis, the etiologic agent of plague, has been responsible for high mortality in several epidemics throughout human history. The plague bacillus has been used as a biological weapon during human history and is currently one of the most likely biological threats. PsaA and PsaE appear to play important roles during Y. pestis infection. Understanding their regulation by environmental cues would facilitate a solution to impede Y. pestis infection.

Hierarchy of transcription factor network in Escherichia coli K-12: H-NS-mediated silencing and Anti-silencing by global regulators

Transcriptional regulation for genome expression determines growth and adaptation of single-cell bacteria that are directly exposed to environment. The transcriptional apparatus in Escherichia coli K-12 is composed of RNA polymerase core enzyme and two groups of its regulatory proteins, seven species of promoter-recognition subunit sigma and about 300 species of transcription factors. The identification of regulatory targets for all these regulatory proteins is critical toward understanding the genome regulation as a whole. For this purpose, we performed a systematic search in vitro of the whole set of binding sites for each factor by gSELEX system. This review summarizes the accumulated knowledge of regulatory targets for more than 150 TFs from E. coli K-12. Overall TFs could be classified into four families: nucleoid-associated bifunctional TFs; global regulators; local regulators; and single-target regulators, in which the regulatory functions remain uncharacterized for the nucleoid-associated TFs. Here we overview the regulatory targets of two nucleoid-associated TFs, H-NS and its paralog StpA, both together playing the silencing role of a set of non-essential genes. Participation of LeuO and other global regulators have been indicated for the anti-silencing. Finally, we propose the hierarchy of TF network as a key framework of the bacterial genome regulation.

Next-generation microbiology: from comparative genomics to gene function

Microbiology is at a turning point in its 120-year history. Widespread next-generation sequencing has revealed genetic complexity among bacteria that could hardly have been imagined by pioneers such as Pasteur, Escherich and Koch. This data cascade brings enormous potential to improve our understanding of individual bacterial cells and the genetic basis of phenotype variation. However, this revolution in data science cannot replace established microbiology practices, presenting the challenge of how to integrate these new techniques. Contrasting comparative and functional genomic approaches, we evoke molecular microbiology theory and established practice to present a conceptual framework and practical roadmap for next-generation microbiology.

Combined RNAseq and ChIPseq Analyses of the BvgA Virulence Regulator of Bordetella pertussis

Bordetella pertussis, the etiological agent of whooping cough, remains a major global health problem. Despite the global usage of whole-cell vaccines since the 1950s and of acellular vaccines in the 1990s, it still is one of the most prevalent vaccine-preventable diseases in industrialized countries. Virulence of B. pertussis is controlled by BvgA/S, a two-component system responsible for upregulation of virulence-activated genes (vags) and downregulation of virulence-repressed genes (vrgs). By transcriptome sequencing (RNAseq) analyses, we identified more than 270 vags or vrgs, and chromatin immunoprecipitation sequencing (ChIPseq) analyses revealed 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Some vags, such as dnt and fhaL, do not contain a BvgA-binding site, suggesting indirect regulation. In contrast, several vrgs and some genes not identified by RNAseq analyses under laboratory conditions contain strong BvgA-binding sites, indicating previously unappreciated complexities of BvgA/S biology.

Bordetella pertussis regulates the production of its virulence factors by the two-component system BvgAS. In the virulence phase, BvgS phosphorylates BvgA, which then activates the transcription of virulence-activated genes (vags). In the avirulence phase, such as during growth in the presence of MgSO₄, BvgA is not phosphorylated and the vags are not expressed. Instead, a set of virulence-repressed genes (vrgs) is expressed. Here, we performed transcriptome sequencing (RNAseq) analyses on B. pertussis cultivated with or without MgSO₄ and on a BvgA-deficient Tohama I derivative. We observed that 146 genes were less expressed under modulating conditions or in the BvgA-deficient strain than under the nonmodulating condition, while 130 genes were more expressed. Some of the genes code for proteins with regulatory functions, suggesting a BvgA/S regulation cascade. To determine which genes are directly regulated by BvgA, we performed chromatin immunoprecipitation sequencing (ChIPseq) analyses. We identified 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Among the former, 32 are in BvgA-regulated putative promoter regions. Some vags, such as dnt and fhaL, contain no BvgA-binding site, suggesting indirect BvgA regulation. Unexpectedly, BvgA also bound to some vrg putative promoter regions. Together, these observations indicate an unrecognized complexity of BvgA/S biology.

IMPORTANCEBordetella pertussis, the etiological agent of whooping cough, remains a major global health problem. Despite the global usage of whole-cell vaccines since the 1950s and of acellular vaccines in the 1990s, it still is one of the most prevalent vaccine-preventable diseases in industrialized countries. Virulence of B. pertussis is controlled by BvgA/S, a two-component system responsible for upregulation of virulence-activated genes (vags) and downregulation of virulence-repressed genes (vrgs). By transcriptome sequencing (RNAseq) analyses, we identified more than 270 vags or vrgs, and chromatin immunoprecipitation sequencing (ChIPseq) analyses revealed 148 BvgA-binding sites, 91 within putative promoter regions, 52 within open reading frames, and 5 in noncoding regions. Some vags, such as dnt and fhaL, do not contain a BvgA-binding site, suggesting indirect regulation. In contrast, several vrgs and some genes not identified by RNAseq analyses under laboratory conditions contain strong BvgA-binding sites, indicating previously unappreciated complexities of BvgA/S biology.

Global Transcriptional Programs in Archaea Share Features with the Eukaryotic Environmental Stress Response

The environmental stress response (ESR), a global transcriptional program originally identified in yeast, is characterized by a rapid and transient transcriptional response composed of large, oppositely regulated gene clusters. Genes induced during the ESR encode core components of stress tolerance, macromolecular repair, and maintenance of homeostasis. In this review, we investigate the possibility for conservation of the ESR across the eukaryotic and archaeal domains of life. We first re-analyze existing transcriptomics data sets to illustrate that a similar transcriptional response is identifiable in Halobacterium salinarum, an archaeal model organism. To substantiate the archaeal ESR, we calculated gene-by-gene correlations, gene function enrichment, and comparison of temporal dynamics. We note reported examples of variation in the ESR across fungi, then synthesize high-level trends present in expression data of other archaeal species. In particular, we emphasize the need for additional high-throughput time series expression data to further characterize stress-responsive transcriptional programs in the Archaea. Together, this review explores an open question regarding features of global transcriptional stress response programs shared across domains of life.

Single-target regulators form a minor group of transcription factors in Escherichia coli K-12

The identification of regulatory targets of all TFs is critical for understanding the entire network of the genome regulation. The lac regulon of Escherichia coli K-12 W3110 is composed of the lacZYA operon and its repressor lacI gene, and has long been recognized as the seminal model of transcription regulation in bacteria with only one highly preferred target. After the Genomic SELEX screening in vitro of more than 200 transcription factors (TFs) from E. coli K-12, however, we found that most TFs regulate multiple target genes. With respect to the number of regulatory targets, a total of these 200 E. coli TFs form a hierarchy ranging from a single target to as many as 1000 targets. Here we focus a total of 13 single-target TFs, 9 known TFs (BetI, KdpE, LacI, MarR, NanR, RpiR, TorR, UlaR and UxuR) and 4 uncharacterized TFs (YagI, YbaO, YbiH and YeaM), altogether forming only a minor group of TFs in E. coli. These single-target TFs were classified into three groups based on their functional regulation.

From biophysics to 'omics and systems biology

Recent decades brought a revolution to biology, driven mainly by exponentially increasing amounts of data coming from "'omics" sciences. To handle these data, bioinformatics often has to combine biologically heterogeneous signals, for which methods from statistics and engineering (e.g. machine learning) are often used. While such an approach is sometimes necessary, it effectively treats the underlying biological processes as a black box. Similarly, systems biology deals with inherently complex systems, characterized by a large number of degrees of freedom, and interactions that are highly non-linear. To deal with this complexity, the underlying physical interactions are often (over)simplified, such as in Boolean modelling of network dynamics. In this review, we argue for the utility of applying a biophysical approach in bioinformatics and systems biology, including discussion of two examples from our research which address sequence analysis and understanding intracellular gene expression dynamics.

ArgR of Streptomyces coelicolor Is a Pleiotropic Transcriptional Regulator: Effect on the Transcriptome, Antibiotic Production, and Differentiation in Liquid Cultures

ArgR is a well-characterized transcriptional repressor controlling the expression of arginine and pyrimidine biosynthetic genes in bacteria. In this work, the biological role of Streptomyces coelicolor ArgR was analyzed by comparing the transcriptomes of S. coelicolor ΔargR and its parental strain, S. coelicolor M145, at five different times over a 66-h period. The effect of S. coelicolor ArgR was more widespread than that of the orthologous protein of Escherichia coli, affecting the expression of 1544 genes along the microarray time series. This S. coelicolor regulator repressed the expression of arginine and pyrimidine biosynthetic genes, but it also modulated the expression of genes not previously described to be regulated by ArgR: genes involved in nitrogen metabolism and nitrate utilization; the act, red, and cpk genes for antibiotic production; genes for the synthesis of the osmotic stress protector ectoine; genes related to hydrophobic cover formation and sporulation (chaplins, rodlins, ramR, and whi genes); all the cwg genes encoding proteins for glycan cell wall biosynthesis; and genes involved in gas vesicle formation. Many of these genes contain ARG boxes for ArgR binding. ArgR binding to seven new ARG boxes, located upstream or near the ectA-ectB, afsS, afsR, glnR, and redH genes, was tested by DNA band-shift assays. These data and those of previously assayed fragments permitted the construction of an improved model of the ArgR binding site. Interestingly, the overexpression of sporulation genes observed in the ΔargR mutant in our culture conditions correlated with a sporulation-like process, an uncommon phenotype.

Genomic SELEX Screening of Regulatory Targets of Escherichia coli Transcription Factors

The genome of Escherichia coli K-12 is transcribed by a single species of RNA polymerase. The selectivity of its transcriptional targets is modulated via two-steps of protein-protein interaction: at the first step, seven species of the sigma subunit are involved, at the second step, a total of approximately 300 species of transcription factor (TFs). For the identification of the regulatory targets of these two groups of regulatory proteins, we developed two in vitro approaches, "Genomic SELEX" (currently designated as gSELEX) and "PS (promoter-specific)-TF" screenings. Here, we describe a detailed protocol of the genomic SELEX screening system which uses purified regulatory proteins and fragments of genomic DNA from E. coli.

RNA polymerase supply and flux through the lac operon in Escherichia coli

Chromatin immunoprecipitation, followed by quantification of immunoprecipitated DNA, can be used to measure RNA polymerase binding to any DNA segment in Escherichia coli. By calibrating measurements against the signal from a single RNA polymerase bound at a single promoter, we can calculate both promoter occupancy levels and the flux of transcribing RNA polymerase through transcription units. Here, we have applied the methodology to the E. coli lactose operon promoter. We confirm that promoter occupancy is limited by recruitment and that the supply of RNA polymerase to the lactose operon promoter depends on its location in the E. coli chromosome. Measurements of RNA polymerase binding to DNA segments within the lactose operon show that flux of RNA polymerase through the operon is low, with, on average, over 18 s elapsing between the passage of transcribing polymerases. Similar low levels of flux were found when semi-synthetic promoters were used to drive transcript initiation, even when the promoter elements were changed to ensure full occupancy of the promoter by RNA polymerase.

This article is part of the themed issue ‘The new bacteriology’.

Scoring Targets of Transcription in Bacteria Rather than Focusing on Individual Binding Sites

Reliable identification of targets of bacterial regulators is necessary to understand bacterial gene expression regulation. These targets are commonly predicted by searching for high-scoring binding sites in the upstream genomic regions, which typically leads to a large number of false positives. In contrast to the common approach, here we propose a novel concept, where overrepresentation of the scoring distribution that corresponds to the entire searched region is assessed, as opposed to predicting individual binding sites. We explore two implementations of this concept, based on Kolmogorov-Smirnov (KS) and Anderson-Darling (AD) tests, which both provide straightforward P-value estimates for predicted targets. This approach is implemented for pleiotropic bacterial regulators, including σ70 (bacterial housekeeping σ factor) target predictions, which is a classical bioinformatics problem characterized by low specificity. We show that KS based approach is both faster and more accurate, departing from the current paradigm of AD being slower, but more accurate. Moreover, KS approach leads to a significant increase in the search accuracy compared to the standard approach, while at the same time straightforwardly assigning well established P-values to each potential target. Consequently, the new KS based method proposed here, which assigns P-values to fixed length upstream regions, provides a fast and accurate approach for predicting bacterial transcription targets.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

Comprehensive mapping of the Helicobacter pylori NikR regulon provides new insights in bacterial nickel responses

Nickel homeostasis is important for pathogenic and ureolytic bacteria, which use this metal ion as enzymatic cofactor. For example, in the human pathogen Helicobacter pylori an optimal balance between nickel uptake and incorporation in metallo-enzymes is fundamental for colonization of the host. Nickel is also used as cofactor to modulate DNA binding of the NikR regulator, which controls transcription of genes involved in nickel trafficking or infection in many bacteria. Accordingly, there is much interest in a systematic characterization of NikR regulation. Herein we use H. pylori as a model to integrate RNA-seq and ChIP-seq data demonstrating that NikR not only regulates metal-ion transporters but also virulence factors, non-coding RNAs, as well as toxin-antitoxin systems in response to nickel stimulation. Altogether, results provide new insights into the pathobiology of H. pylori and contribute to understand the responses to nickel in other bacteria.

The genome-scale DNA-binding profile of BarR, a β-alanine responsive transcription factor in the archaeon Sulfolobus acidocaldarius

Background

The Leucine-responsive Regulatory Protein (Lrp) family is a widespread family of regulatory transcription factors in prokaryotes. BarR is an Lrp-like transcription factor in the model archaeon Sulfolobus acidocaldarius that activates the expression of a β-alanine aminotransferase gene, which is involved in β-alanine degradation. In contrast to classical Lrp-like transcription factors, BarR is not responsive to any of the α-amino acids but interacts specifically with β-alanine. Besides the juxtaposed β-alanine aminotransferase gene, other regulatory targets of BarR have not yet been identified although β-alanine is the precursor of coenzyme A and thus an important central metabolite. The aim of this study is to extend the knowledge of the DNA-binding characteristics of BarR and of its corresponding regulon from a local to a genome-wide perspective.

Results

We characterized the genome-wide binding profile of BarR using chromatin immunoprecipation combined with high-throughput sequencing (ChIP-seq). This revealed 21 genomic binding loci. High-enrichment binding regions were validated to interact with purified BarR protein in vitro using electrophoretic mobility shift assays and almost all targets were also shown to harbour a conserved semi-palindromic binding motif. Only a small subset of enriched genomic sites are located in intergenic regions at a relative short distance to a promoter, and qRT-PCR analysis demonstrated that only one additional operon is under activation of BarR, namely the glutamine synthase operon. The latter is also a target of other Lrp-like transcription factors. Detailed inspection of the BarR ChIP-seq profile at the β-alanine aminotransferase promoter region in combination with binding motif predictions indicate that the operator structure is more complicated than previously anticipated, consisting of multiple (major and auxiliary) operators.

Conclusions

BarR has a limited regulon, and includes also glutamine synthase genes besides the previously characterized β-alanine aminotransferase. Regulation of glutamine synthase is suggestive of a link between β-alanine and α-amino acid metabolism in S. acidocaldarius. Furthermore, this work reveals that the BarR regulon overlaps with that of other Lrp-like regulators.

Electronic supplementary material

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

The unexpected complexity of bacterial genomes

Gene organization and control are described by models conceived in the 1960s. These models explain basic gene regulatory mechanisms and underpin current genome annotation. However, such models struggle to explain recent genome-scale observations. For example, accounts of RNA synthesis initiating within genes, widespread antisense transcription and non-canonical DNA binding by gene regulatory proteins are difficult to reconcile with traditional thinking. As a result, unexpected observations have often been dismissed and downstream consequences ignored. In this paper I will argue that, to fully understand the biology of bacterial chromosomes, we must embrace their hidden layers of complexity.

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.

Transcription profile of Escherichia coli: genomic SELEX search for regulatory targets of transcription factors

Bacterial genomes are transcribed by DNA-dependent RNA polymerase (RNAP), which achieves gene selectivity through interaction with sigma factors that recognize promoters, and transcription factors (TFs) that control the activity and specificity of RNAP holoenzyme. To understand the molecular mechanisms of transcriptional regulation, the identification of regulatory targets is needed for all these factors. We then performed genomic SELEX screenings of targets under the control of each sigma factor and each TF. Here we describe the assembly of 156 SELEX patterns of a total of 116 TFs performed in the presence and absence of effector ligands. The results reveal several novel concepts: (i) each TF regulates more targets than hitherto recognized; (ii) each promoter is regulated by more TFs than hitherto recognized; and (iii) the binding sites of some TFs are located within operons and even inside open reading frames. The binding sites of a set of global regulators, including cAMP receptor protein, LeuO and Lrp, overlap with those of the silencer H-NS, suggesting that certain global regulators play an anti-silencing role. To facilitate sharing of these accumulated SELEX datasets with the research community, we compiled a database, ‘Transcription Profile of Escherichia coli’ (www.shigen.nig.ac.jp/ecoli/tec/).

Genomic Targets and Features of BarA-UvrY (-SirA) Signal Transduction Systems

The two-component signal transduction system BarA-UvrY of Escherichia coli and its orthologs globally regulate metabolism, motility, biofilm formation, stress resistance, virulence of pathogens and quorum sensing by activating the transcription of genes for regulatory sRNAs, e.g. CsrB and CsrC in E. coli. These sRNAs act by sequestering the RNA binding protein CsrA (RsmA) away from lower affinity mRNA targets. In this study, we used ChIP-exo to identify, at single nucleotide resolution, genomic sites for UvrY (SirA) binding in E. coli and Salmonella enterica. The csrB and csrC genes were the strongest targets of crosslinking, which required UvrY phosphorylation by the BarA sensor kinase. Crosslinking occurred at two sites, an inverted repeat sequence far upstream of the promoter and a site near the -35 sequence. DNAse I footprinting revealed specific binding of UvrY in vitro only to the upstream site, indicative of additional binding requirements and/or indirect binding to the downstream site. Additional genes, including cspA, encoding the cold-shock RNA-binding protein CspA, showed weaker crosslinking and modest or negligible regulation by UvrY. We conclude that the global effects of UvrY/SirA on gene expression are primarily mediated by activating csrB and csrC transcription. We also used in vivo crosslinking and other experimental approaches to reveal new features of csrB/csrC regulation by the DeaD and SrmB RNA helicases, IHF, ppGpp and DksA. Finally, the phylogenetic distribution of BarA-UvrY was analyzed and found to be uniquely characteristic of γ-Proteobacteria and strongly anti-correlated with fliW, which encodes a protein that binds to CsrA and antagonizes its activity in Bacillus subtilis. We propose that BarA-UvrY and orthologous TCS transcribe sRNA antagonists of CsrA throughout the γ-Proteobacteria, but rarely or never perform this function in other species.

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.

Systems biology approaches to defining transcription regulatory networks in halophilic archaea

To survive complex and changing environmental conditions, microorganisms use gene regulatory networks (GRNs) composed of interacting regulatory transcription factors (TFs) to control the timing and magnitude of gene expression. Genome-wide datasets; such as transcriptomics and protein-DNA interactions; and experiments such as high throughput growth curves; facilitate the construction of GRNs and provide insight into TF interactions occurring under stress. Systems biology approaches integrate these datasets into models of GRN architecture as well as statistical and/or dynamical models to understand the function of networks occurring in cells. Previously, these types of studies have focused on traditional model organisms (e.g. Escherichia coli, yeast). However, recent advances in archaeal genetics and other tools have enabled a systems approach to understanding GRNs in these relatively less studied archaeal model organisms. In this report, we outline a systems biology workflow for generating and integrating data focusing on the TF regulator. We discuss experimental design, outline the process of data collection, and provide the tools required to produce high confidence regulons for the TFs of interest. We provide a case study as an example of this workflow, describing the construction of a GRN centered on multi-TF coordinate control of gene expression governing the oxidative stress response in the hypersaline-adapted archaeon Halobacterium salinarum.

Single cell FRET analysis for the identification of optimal FRET-pairs in Bacillus subtilis using a prototype MEM-FLIM system

Protein-protein interactions can be studied in vitro, e.g. with bacterial or yeast two-hybrid systems or surface plasmon resonance. In contrast to in vitro techniques, in vivo studies of protein-protein interactions allow examination of spatial and temporal behavior of such interactions in their native environment. One approach to study protein-protein interactions in vivo is via Förster Resonance Energy Transfer (FRET). Here, FRET efficiency of selected FRET-pairs was studied at the single cell level using sensitized emission and Frequency Domain-Fluorescence Lifetime Imaging Microscopy (FD-FLIM). For FRET-FLIM, a prototype Modulated Electron-Multiplied FLIM system was used, which is, to the best of our knowledge, the first account of Frequency Domain FLIM to analyze FRET in single bacterial cells. To perform FRET-FLIM, we first determined and benchmarked the best fluorescent protein-pair for FRET in Bacillus subtilis using a novel BglBrick-compatible integration vector. We show that GFP-tagRFP is an excellent donor-acceptor pair for B. subtilis in vivo FRET studies. As a proof of concept, selected donor and acceptor fluorescent proteins were fused using a linker that contained a tobacco etch virus (TEV)-protease recognition sequence. Induction of TEV-protease results in loss of FRET efficiency and increase in fluorescence lifetime. The loss of FRET efficiency after TEV induction can be followed in time in single cells via time-lapse microscopy. This work will facilitate future studies of in vivo dynamics of protein complexes in single B. subtilis cells.

Transcriptional analysis of the MrpJ network: modulation of diverse virulence-associated genes and direct regulation of mrp fimbrial and flhDC flagellar operons in Proteus mirabilis

The enteric bacterium Proteus mirabilis is associated with a significant number of catheter-associated urinary tract infections (UTIs). Strict regulation of the antagonistic processes of adhesion and motility, mediated by fimbriae and flagella, respectively, is essential for disease progression. Previously, the transcriptional regulator MrpJ, which is encoded by the mrp fimbrial operon, has been shown to repress both swimming and swarming motility. Here we show that MrpJ affects an array of cellular processes beyond adherence and motility. Microarray analysis found that expression of mrpJ mimicking levels observed during UTIs leads to differential expression of 217 genes related to, among other functions, bacterial virulence, type VI secretion, and metabolism. We probed the molecular mechanism of transcriptional regulation by MrpJ using transcriptional reporters and chromatin immunoprecipitation (ChIP). Binding of MrpJ to two virulence-associated target gene promoters, the promoters of the flagellar master regulator flhDC and mrp itself, appears to be affected by the condensation state of the native chromosome, although both targets share a direct MrpJ binding site proximal to the transcriptional start. Furthermore, an mrpJ deletion mutant colonized the bladders of mice at significantly lower levels in a transurethral model of infection. Additionally, we observed that mrpJ is widely conserved in a collection of recent clinical isolates. Altogether, these findings support a role of MrpJ as a global regulator of P. mirabilis virulence.

The TyrR transcription factor regulates the divergent akr-ipdC operons of Enterobacter cloacae UW5

The TyrR transcription factor regulates genes involved in the uptake and biosynthesis of aromatic amino acids in Enterobacteriaceae. Genes may be positively or negatively regulated depending on the presence or absence of each aromatic amino acid, all three of which function as cofactors for TyrR. In this report we detail the transcriptional control of two divergently transcribed genes, akr and ipdC, by TyrR, elucidated by promoter fusion expression assays and electrophoretic mobility shift assays to assess protein-DNA interactions. Expression of both genes was shown to be controlled by TyrR via interactions with two TyrR boxes located within the akr-ipdC intergenic region. Expression of ipdC required TyrR bound to the proximal strong box, and is strongly induced by phenylalanine, and to a lesser extent by tryptophan and tyrosine. Down-regulation of akr was reliant on interactions with the weak box, and may also require a second, as yet unidentified protein for further repression. Tyrosine enhanced repression of akr. Electrophoretic mobility shift assays demonstrated that TyrR interacts with both the strong and weak boxes, and that binding of the weak box in vitro requires an intact adjacent strong box. While the strong box shows a high degree of conservation with the TyrR binding site consensus sequence, the weak box has atypical spacing of the two half sites comprising the palindromic arms. Site-directed mutagenesis demonstrated sequence-specific interaction between TyrR and the weak box. This is the first report of TyrR-controlled expression of two divergent protein-coding genes, transcribed from independent promoters. Moreover, the identification of a predicted aldo-keto reductase as a member of the TyrR regulon further extends the function of the TyrR regulon.

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.

The DNA-binding network of Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) infects 30% of all humans and kills someone every 20-30 s. Here we report genome-wide binding for ~80% of all predicted MTB transcription factors (TFs), and assayed global expression following induction of each TF. The MTB DNA-binding network consists of ~16,000 binding events from 154 TFs. We identify >50 TF-DNA consensus motifs and >1,150 promoter-binding events directly associated with proximal gene regulation. An additional ~4,200 binding events are in promoter windows and represent strong candidates for direct transcriptional regulation under appropriate environmental conditions. However, we also identify >10,000 'dormant' DNA-binding events that cannot be linked directly with proximal transcriptional control, suggesting that widespread DNA binding may be a common feature that should be considered when developing global models of coordinated gene expression.

ChIP-seq for genome-scale analysis of bacterial DNA-binding proteins

Protein-DNA interactions are central to many basic biological processes, including transcription regulation, DNA replication, and DNA repair. Chromatin Immunoprecipitation (ChIP) is used to determine the position and strength of protein-DNA interactions in vivo. Coupling ChIP with microarrays (ChIP-chip), and more recently with deep sequencing (ChIP-seq), has allowed genome-wide profiling of DNA binding events in vivo. In this chapter we outline the steps to generate ChIP-seq libraries from bacterial samples and briefly discuss basic analysis of the data.

Mapping Transcription Regulatory Networks with ChIP-seq and RNA-seq

Bacterial genomes encode numerous transcription factors, DNA-binding proteins that regulate transcription initiation. Identifying the regulatory targets of transcription factors is a major challenge of systems biology. Here I describe the use of two genome-scale approaches, ChIP-seq and RNA-seq, that are used to map transcription factor regulons. ChIP-seq maps the association of transcription factors with DNA, and RNA-seq determines changes in RNA levels associated with transcription factor perturbation. I discuss the strengths and weaknesses of these and related approaches, and I describe how ChIP-seq and RNA-seq can be combined to map individual transcription factor regulons and entire regulatory networks.

High-resolution digital profiling of the epigenome

The widespread adoption of short-read DNA sequencing as a digital epigenomic readout platform has motivated the development of genome-wide tools that achieve base-pair resolution. New methods for footprinting and affinity purification of nucleosomes, RNA polymerases, chromatin remodellers and transcription factors have increased the resolution of epigenomic profiling by two orders of magnitude, leading to new insights into how the chromatin landscape affects gene regulation. These digital epigenomic tools have also been applied to directly profile both turnover kinetics and transcription in situ. In this Review, we describe how these new genome-wide tools allow interrogation of diverse aspects of the epigenome.

Comprehensive mapping of the Escherichia coli flagellar regulatory network

Flagellar synthesis is a highly regulated process in all motile bacteria. In Escherichia coli and related species, the transcription factor FlhDC is the master regulator of a multi-tiered transcription network. FlhDC activates transcription of a number of genes, including some flagellar genes and the gene encoding the alternative Sigma factor FliA. Genes whose expression is required late in flagellar assembly are primarily transcribed by FliA, imparting temporal regulation of transcription and coupling expression to flagellar assembly. In this study, we use ChIP-seq and RNA-seq to comprehensively map the E. coli FlhDC and FliA regulons. We define a surprisingly restricted FlhDC regulon, including two novel regulated targets and two binding sites not associated with detectable regulation of surrounding genes. In contrast, we greatly expand the known FliA regulon. Surprisingly, 30 of the 52 FliA binding sites are located inside genes. Two of these intragenic promoters are associated with detectable noncoding RNAs, while the others either produce highly unstable RNAs or are inactive under these conditions. Together, our data redefine the E. coli flagellar regulatory network, and provide new insight into the temporal orchestration of gene expression that coordinates the flagellar assembly process.

Flagella are surface-associated appendages that propel bacteria and are involved in diverse functions such as chemotaxis, surface attachment, and host cell invasion. Flagella are incredibly complex macromolecular machines that are energetically costly to produce, assemble, and power. Flagellar production is tightly regulated and flagellar components are only synthesized when flagellar motility is advantageous. Regulation also ensures that flagellar components are produced in roughly the same order in which they are needed, increasing efficiency of the assembly process. The transcriptional regulation of flagellar genes has been studied extensively in the model organism Escherichia coli; however, many questions remain. We have used an unbiased, genome-wide approach to comprehensively identify all of the binding sites and regulatory targets for the two key regulators of flagellar synthesis, FlhDC and FliA. Our results redefine the flagellar regulatory network, and suggest that FliA binds many sites that are not associated with productive transcription. This work is important because it suggests possible new functions for FliA outside of the transcription of canonical mRNAs, and it provides new insight into the temporal orchestration of gene expression that coordinates the flagellar assembly process.

Lies and deception in bacterial gene regulation: the roles of nucleic acid decoys

Bacteria use intricately interconnected mechanisms acting at the transcriptional and post-transcriptional level to adjust gene expression to their needs. An intriguing example found in the chitosugar utilization systems of Escherichia coli and Salmonella is uncovered in a study by Plumbridge and colleagues. Three transcription factors (TFs), a small regulatory RNA (sRNA) and a sRNA trap cooperate to set thresholds and dynamics in regulation of chitosugar utilization. Specifically, under inducing conditions a decoy site on the polycistronic chitobiose (chbBCARFG) mRNA sequesters sRNA ChiX, which represses synthesis of the separately encoded chitoporin ChiP. Base-pairing of ChiX with its decoy has no role for the chb genes themselves when the mRNA is in excess. In the absence of substrate, however, this base-pairing tightly represses chbC encoding a subunit of the chitosugar transporter. Thus, one and the same sRNA/mRNA interaction serves different regulatory functions under different environmental conditions. The employment of RNA decoys to control the activities of post-transcriptional regulators themselves is an increasingly recognized mechanism in gene regulation. Another observation in the current study highlights the possibility that decoy sites might even exist on the DNA controlling the availability of TFs for their target promoters.

Chromatin immunoprecipitation for ChIP-chip and ChIP-seq

Bacterial adaptation to given environmental conditions is largely achieved by complex gene regulatory processes. To address the question how and to what extend single transcriptional regulators modulate gene expression, chromatin immunoprecipitation (ChIP) coupled to DNA microarrays (ChIP-chip) or to next-generation sequencing (ChIP-seq) is one of the preferred methods. Both ChIP-chip and ChIP-seq can generate genome-wide maps of protein-DNA interactions and thus identify primary regulons of transcription factors. In combination with transcriptome analyses, the obtained data can be used to compile complex regulatory networks which in terms will advance our understanding of bacterial adaptation processes to specific environmental conditions.

Identification of HilD-regulated genes in Salmonella enterica serovar Typhimurium

Salmonella enterica serovar Typhimurium (S. Typhimurium) pathogenicity island 1 (SPI-1) encodes a type III secretion system required for invasion of host gut epithelial cells. Expression of SPI-1 virulence genes is controlled by a complex hierarchy of transcription factors encoded within and outside SPI-1. The master regulator of SPI-1, HilA, is itself regulated by three homologous transcription factors, HilD, HilC, and RtsA. HilD activates transcription of hilA and other target genes in response to environmental conditions associated with the intestinal microenvironment of the host. We have mapped the binding of HilD across the S. Typhimurium genome using chromatin immunoprecipitation-sequencing (ChIP-seq). Thus, we have identified 17 regions bound by HilD, including 11 novel targets. The majority of HilD targets are located outside SPI-1. We demonstrate transcription activation of 8 genes by HilD; four of these genes have not been previously described as being regulated by HilD, including lpxR, which encodes a lipid A deacylase important for immune evasion. We also show that HilD-activated genes are frequently activated by HilC and RtsA, indicating extensive overlap of the HilD, HilC, and RtsA regulons.

Light-dependent, plastome-wide association of the plastid-encoded RNA polymerase with chloroplast DNA

Plastid genes are transcribed by two types of RNA polymerases: a plastid-encoded eubacterial-type RNA polymerase (PEP) and nuclear-encoded phage-type RNA polymerases (NEPs). To investigate the spatio-temporal expression of PEP, we tagged its α-subunit with a hemagglutinin epitope (HA). Transplastomic tobacco plants were generated and analyzed for the distribution of the tagged polymerase in plastid sub-fractions, and associated genes were identified under various light conditions. RpoA:HA was detected as early as the 3rd day after imbibition, and was constitutively expressed in green tissue over 60 days of plant development. We found that the tagged polymerase subunit preferentially associated with the plastid membranes, and was less abundant in the soluble stroma fraction. Attachment of RpoA:HA to the membrane fraction during early seedling development was independent of DNA, but at later stages of development, DNA appears to facilitate attachment of the polymerase to membranes. To survey PEP-dependent transcription units, we probed for nucleic acids enriched in RpoA:HA precipitates using a tobacco chloroplast whole-genome tiling array. The most strongly co-enriched DNA fragments represent photosynthesis genes (e.g. psbA, psbC, psbD and rbcL), whose expression is known to be driven by PEP promoters, while NEP-dependent genes were less abundant in RpoA:HA precipitates. Additionally, we demonstrate that the association of PEP with photosynthesis-related genes was reduced during the dark period, indicating that plastome-wide PEP-DNA association is a light-dependent process.

Genome-scale analyses of Escherichia coli and Salmonella enterica AraC reveal noncanonical targets and an expanded core regulon

Escherichia coli AraC is a well-described transcription activator of genes involved in arabinose metabolism. Using complementary genomic approaches, chromatin immunoprecipitation (ChIP)-chip, and transcription profiling, we identify direct regulatory targets of AraC, including five novel target genes: ytfQ, ydeN, ydeM, ygeA, and polB. Strikingly, only ytfQ has an established connection to arabinose metabolism, suggesting that AraC has a broader function than previously described. We demonstrate arabinose-dependent repression of ydeNM by AraC, in contrast to the well-described arabinose-dependent activation of other target genes. We also demonstrate unexpected read-through of transcription at the Rho-independent terminators downstream of araD and araE, leading to significant increases in the expression of polB and ygeA, respectively. AraC is highly conserved in the related species Salmonella enterica. We use ChIP sequencing (ChIP-seq) and RNA sequencing (RNA-seq) to map the AraC regulon in S. enterica. A comparison of the E. coli and S. enterica AraC regulons, coupled with a bioinformatic analysis of other related species, reveals a conserved regulatory network across the family Enterobacteriaceae comprised of 10 genes associated with arabinose transport and metabolism.

Identification of the alternative sigma factor SigX regulon and its implications for Pseudomonas aeruginosa pathogenicity

Pseudomonas aeruginosa is distinguished by its broad metabolic diversity and its remarkable capability for adaptation, which relies on a large collection of transcriptional regulators and alternative sigma (σ) factors. The largest group of alternative σ factors is that of the extracytoplasmic function (ECF) σ factors, which control key transduction pathways for maintenance of envelope homeostasis in response to external stress and cell growth. In addition, there are specific roles of alternative σ factors in regulating the expression of virulence and virulence-associated genes. Here, we analyzed a deletion mutant of the ECF σ factor SigX and applied mRNA profiling to define the SigX-dependent regulon in P. aeruginosa in response to low-osmolarity-medium conditions. Furthermore, the combination of transcriptional data with chromatin immunoprecipitation (ChIP) followed by high-throughput sequencing (ChIP-seq) led to the identification of the DNA binding motif of SigX. Genome-wide mapping of SigX-binding regions revealed enrichment of downstream genes involved in fatty acid biosynthesis, type III secretion, swarming and cyclic di-GMP (c-di-GMP) signaling. In accordance, a sigX deletion mutant exhibited altered fatty acid composition of the cell membrane, reduced cytotoxicity, impaired swarming activity, elevated c-di-GMP levels, and increased biofilm formation. In conclusion, a combination of ChIP-seq with transcriptional profiling and bioinformatic approaches to define consensus DNA binding sequences proved to be effective for the elucidation of the regulon of the alternative σ factor SigX, revealing its role in complex virulence-associated phenotypes in P. aeruginosa.

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.

Integrated stress response of Escherichia coli to methylglyoxal: transcriptional readthrough from the nemRA operon enhances protection through increased expression of glyoxalase I

Methylglyoxal (MG) elicits activation of K⁺ efflux systems to protect cells against the toxicity of the electrophile. ChIP-chip targeting RNA polymerase, supported by a range of other biochemical measurements and mutant creation, was used to identify genes transcribed in response to MG and which complement this rapid response. The SOS DNA repair regulon is induced at cytotoxic levels of MG, even when exposure to MG is transient. Glyoxalase I alone among the core MG protective systems is induced in response to MG exposure. Increased expression is an indirect consequence of induction of the upstream nemRA operon, encoding an enzyme system that itself does not contribute to MG detoxification. Moreover, this induction, via nemRA only occurs when cells are exposed to growth inhibitory concentrations of MG. We show that the kdpFABCDE genes are induced and that this expression occurs as a result of depletion of cytoplasmic K⁺ consequent upon activation of the KefGB K⁺ efflux system. Finally, our analysis suggests that the transcriptional changes in response to MG are a culmination of the damage to DNA and proteins, but that some integrate specific functions, such as DNA repair, to augment the allosteric activation of the main protective system, KefGB.

Non-canonical protein-DNA interactions identified by ChIP are not artifacts

Background

ChIP-chip and ChIP-seq are widely used methods to map protein-DNA interactions on a genomic scale in vivo. Waldminghaus and Skarstad recently reported, in this journal, a modified method for ChIP-chip. Based on a comparison of our previously-published ChIP-chip data for Escherichia coli σ³² with their own data, Waldminghaus and Skarstad concluded that many of the σ³² targets identified in our earlier work are false positives. In particular, we identified many non-canonical σ³² targets that are located inside genes or are associated with genes that show no detectable regulation by σ³². Waldminghaus and Skarstad propose that such non-canonical sites are artifacts, identified due to flaws in the standard ChIP methodology. Waldminghaus and Skarstad suggest specific changes to the standard ChIP procedure that reportedly eliminate the claimed artifacts.

Results

We reanalyzed our published ChIP-chip datasets for σ³² and the datasets generated by Waldminghaus and Skarstad to assess data quality and reproducibility. We also performed targeted ChIP/qPCR for σ³² and an unrelated transcription factor, AraC, using the standard ChIP method and the modified ChIP method proposed by Waldminghaus and Skarstad. Furthermore, we determined the association of core RNA polymerase with disputed σ³² promoters, with and without overexpression of σ³². We show that (i) our published σ³² ChIP-chip datasets have a consistently higher dynamic range than those of Waldminghaus and Skarstad, (ii) our published σ³² ChIP-chip datasets are highly reproducible, whereas those of Waldminghaus and Skarstad are not, (iii) non-canonical σ³² target regions are enriched in a σ³² ChIP in a heat shock-dependent manner, regardless of the ChIP method used, (iv) association of core RNA polymerase with some disputed σ³² target genes is induced by overexpression of σ³², (v) σ³² targets disputed by Waldminghaus and Skarstad are predominantly those that are most weakly bound, and (vi) the modifications to the ChIP method proposed by Waldminghaus and Skarstad reduce enrichment of all protein-bound genomic regions.

Conclusions

The modifications to the ChIP-chip method suggested by Waldminghaus and Skarstad reduce rather than increase the quality of ChIP data. Hence, the non-canonical σ³² targets identified in our previous study are likely to be genuine. We propose that the failure of Waldminghaus and Skarstad to identify many of these σ³² targets is due predominantly to the lower data quality in their study. We conclude that surprising ChIP-chip results are not artifacts to be ignored, but rather indications that our understanding of DNA-binding proteins is incomplete.

Evidence classification of high-throughput protocols and confidence integration in RegulonDB

RegulonDB provides curated information on the transcriptional regulatory network of Escherichia coli and contains both experimental data and computationally predicted objects. To account for the heterogeneity of these data, we introduced in version 6.0, a two-tier rating system for the strength of evidence, classifying evidence as either ‘weak’ or ‘strong’ (Gama-Castro,S., Jimenez-Jacinto,V., Peralta-Gil,M. et al. RegulonDB (Version 6.0): gene regulation model of Escherichia Coli K-12 beyond transcription, active (experimental) annotated promoters and textpresso navigation. Nucleic Acids Res., 2008;36:D120–D124.). We now add to our classification scheme the classification of high-throughput evidence, including chromatin immunoprecipitation (ChIP) and RNA-seq technologies. To integrate these data into RegulonDB, we present two strategies for the evaluation of confidence, statistical validation and independent cross-validation. Statistical validation involves verification of ChIP data for transcription factor-binding sites, using tools for motif discovery and quality assessment of the discovered matrices. Independent cross-validation combines independent evidence with the intention to mutually exclude false positives. Both statistical validation and cross-validation allow to upgrade subsets of data that are supported by weak evidence to a higher confidence level. Likewise, cross-validation of strong confidence data extends our two-tier rating system to a three-tier system by introducing a third confidence score ‘confirmed’.

Database URL:http://regulondb.ccg.unam.mx/

Nanobody(R)-based chromatin immunoprecipitation/micro-array analysis for genome-wide identification of transcription factor DNA binding sites

Nanobodies® are single-domain antibody fragments derived from camelid heavy-chain antibodies. Because of their small size, straightforward production in Escherichia coli, easy tailoring, high affinity, specificity, stability and solubility, nanobodies® have been exploited in various biotechnological applications. A major challenge in the post-genomics and post-proteomics era is the identification of regulatory networks involving nucleic acid-protein and protein-protein interactions. Here, we apply a nanobody® in chromatin immunoprecipitation followed by DNA microarray hybridization (ChIP-chip) for genome-wide identification of DNA-protein interactions. The Lrp-like regulator Ss-LrpB, arguably one of the best-studied specific transcription factors of the hyperthermophilic archaeon Sulfolobus solfataricus, was chosen for this proof-of-principle nanobody®-assisted ChIP. Three distinct Ss-LrpB-specific nanobodies®, each interacting with a different epitope, were generated for ChIP. Genome-wide ChIP-chip with one of these nanobodies® identified the well-established Ss-LrpB binding sites and revealed several unknown target sequences. Furthermore, these ChIP-chip profiles revealed auxiliary operator sites in the open reading frame of Ss-lrpB. Our work introduces nanobodies® as a novel class of affinity reagents for ChIP. Taking into account the unique characteristics of nanobodies®, in particular, their short generation time, nanobody®-based ChIP is expected to further streamline ChIP-chip and ChIP-Seq experiments, especially in organisms with no (or limited) possibility of genetic manipulation.

Diverse control of metabolism and other cellular processes in Streptomyces coelicolor by the PhoP transcription factor: genome-wide identification of in vivo targets

Streptomycetes sense and respond to the stress of phosphate starvation via the two-component PhoR–PhoP signal transduction system. To identify the in vivo targets of PhoP we have undertaken a chromatin-immunoprecipitation-on-microarray analysis of wild-type and phoP mutant cultures and, in parallel, have quantified their transcriptomes. Most (ca. 80%) of the previously in vitro characterized PhoP targets were identified in this study among several hundred other putative novel PhoP targets. In addition to activating genes for phosphate scavenging systems PhoP was shown to target two gene clusters for cell wall/extracellular polymer biosynthesis. Furthermore PhoP was found to repress an unprecedented range of pathways upon entering phosphate limitation including nitrogen assimilation, oxidative phosphorylation, nucleotide biosynthesis and glycogen catabolism. Moreover, PhoP was shown to target many key genes involved in antibiotic production and morphological differentiation, including afsS, atrA, bldA, bldC, bldD, bldK, bldM, cdaR, cdgA, cdgB and scbR-scbA. Intriguingly, in the PhoP-dependent cpk polyketide gene cluster, PhoP accumulates substantially at three specific sites within the giant polyketide synthase-encoding genes. This study suggests that, following phosphate limitation, Streptomyces coelicolor PhoP functions as a ‘master’ regulator, suppressing central metabolism, secondary metabolism and developmental pathways until sufficient phosphate is salvaged to support further growth and, ultimately, morphological development.