A comparative genomics approach to prediction of new members of regulons

Comprehensive Characterization of fucAO Operon Activation in Escherichia coli

Wildtype Escherichia coli cells cannot grow on L-1,2-propanediol, as the fucAO operon within the fucose (fuc) regulon is thought to be silent in the absence of L-fucose. Little information is available concerning the transcriptional regulation of this operon. Here, we first confirm that fucAO operon expression is highly inducible by fucose and is primarily attributable to the upstream operon promoter, while the fucO promoter within the 3'-end of fucA is weak and uninducible. Using 5'RACE, we identify the actual transcriptional start site (TSS) of the main fucAO operon promoter, refuting the originally proposed TSS. Several lines of evidence are provided showing that the fucAO locus is within a transcriptionally repressed region on the chromosome. Operon activation is dependent on FucR and Crp but not SrsR. Two Crp-cAMP binding sites previously found in the regulatory region are validated, where the upstream site plays a more critical role than the downstream site in operon activation. Furthermore, two FucR binding sites are identified, where the downstream site near the first Crp site is more important than the upstream site. Operon transcription relies on Crp-cAMP to a greater degree than on FucR. Our data strongly suggest that FucR mainly functions to facilitate the binding of Crp to its upstream site, which in turn activates the fucAO promoter by efficiently recruiting RNA polymerase.

Flexible comparative genomics of prokaryotic transcriptional regulatory networks

Background

Comparative genomics methods enable the reconstruction of bacterial regulatory networks using available experimental data. In spite of their potential for accelerating research into the composition and evolution of bacterial regulons, few comparative genomics suites have been developed for the automated analysis of these regulatory systems. Available solutions typically rely on precomputed databases for operon and ortholog predictions, limiting the scope of analyses to processed complete genomes, and several key issues such as the transfer of experimental information or the integration of regulatory information in a probabilistic setting remain largely unaddressed.

Results

Here we introduce CGB, a flexible platform for comparative genomics of prokaryotic regulons. CGB has few external dependencies and enables fully customized analyses of newly available genome data. The platform automates the merging of experimental information and uses a gene-centered, Bayesian framework to generate and integrate easily interpretable results. We demonstrate its flexibility and power by analyzing the evolution of type III secretion system regulation in pathogenic Proteobacteria and by characterizing the SOS regulon of a new bacterial phylum, the Balneolaeota.

Conclusions

Our results demonstrate the applicability of the CGB pipeline in multiple settings. CGB’s ability to automatically integrate experimental information from multiple sources and use complete and draft genomic data, coupled with its non-reliance on precomputed databases and its easily interpretable display of gene-centered posterior probabilities of regulation provide users with an unprecedented level of flexibility in launching comparative genomics analyses of prokaryotic transcriptional regulatory networks. The analyses of type III secretion and SOS response regulatory networks illustrate instances of convergent and divergent evolution of these regulatory systems, showcasing the power of formal ancestral state reconstruction at inferring the evolutionary history of regulatory networks.

Glycine Cleavage System and cAMP Receptor Protein Co-Regulate CRISPR/cas3 Expression to Resist Bacteriophage

The CRISPR/Cas system protects bacteria against bacteriophage and plasmids through a sophisticated mechanism where cas operon plays a crucial role consisting of cse1 and cas3. However, comprehensive studies on the regulation of cas3 operon of the Type I-E CRISPR/Cas system are scarce. Herein, we investigated the regulation of cas3 in Escherichia coli. The mutation in gcvP or crp reduced the CRISPR/Cas system interference ability and increased bacterial susceptibility to phage, when the casA operon of the CRISPR/Cas system was activated. The silence of the glycine cleavage system (GCS) encoded by gcvTHP operon reduced cas3 expression. Adding N⁵, N¹⁰-methylene tetrahydrofolate (N⁵, N¹⁰-mTHF), which is the product of GCS-catalyzed glycine, was able to activate cas3 expression. In addition, a cAMP receptor protein (CRP) encoded by crp activated cas3 expression via binding to the cas3 promoter in response to cAMP concentration. Since N⁵, N¹⁰-mTHF provides one-carbon unit for purine, we assumed GCS regulates cas3 through associating with CRP. It was evident that the mutation of gcvP failed to further reduce the cas3 expression with the crp deletion. These results illustrated a novel regulatory pathway which GCS and CRP co-regulate cas3 of the CRISPR/Cas system and contribute to the defence against invasive genetic elements, where CRP is indispensable for GCS regulation of cas3 expression.

Overproduction and purification of the Escherichia coli transcription factors "toxic" to a bacterial cell

Transcription factors play a crucial role in control of life of a bacterial cell, working as switchers to a different life style or pathogenicity. To reconstruct the network of regulatory events taking place in changing growth conditions, we need to know regulons of as many transcription factors as possible, and motifs recognized by them. Experimentally this can be attained via ChIP-seq in vivo, SELEX and DNAse I footprinting in vitro. All these approaches require large amounts of purified proteins. However, overproduction of transcription factors leading to their extensive binding to the regulatory elements on the DNA make them toxic to a bacterial cell thus significantly complicating production of a soluble protein. Here, on the example of three regulators from Escherichia coli, UxuR, ExuR, and LeuO, we show that stable production of toxic transcription factors in a soluble fraction can be significantly enhanced by holding the expression of a recombinant protein back at the early stages of bacterial growth. This can be achieved by cloning genes together with their regulatory regions containing repressor sites, with subsequent growth in a very rich media where activity of excessive regulators is not crucial, followed by induction with a very low concentration of an inducer. Schemes of further purification of these proteins were developed, and functional activity was confirmed.

RECTA: Regulon Identification Based on Comparative Genomics and Transcriptomics Analysis

Regulons, which serve as co-regulated gene groups contributing to the transcriptional regulation of microbial genomes, have the potential to aid in understanding of underlying regulatory mechanisms. In this study, we designed a novel computational pipeline, regulon identification based on comparative genomics and transcriptomics analysis (RECTA), for regulon prediction related to the gene regulatory network under certain conditions. To demonstrate the effectiveness of this tool, we implemented RECTA on Lactococcus lactis MG1363 data to elucidate acid-response regulons. A total of 51 regulons were identified, 14 of which have computational-verified significance. Among these 14 regulons, five of them were computationally predicted to be connected with acid stress response. Validated by literature, 33 genes in Lactococcus lactis MG1363 were found to have orthologous genes which were associated with six regulons. An acid response related regulatory network was constructed, involving two trans-membrane proteins, eight regulons (llrA, llrC, hllA, ccpA, NHP6A, rcfB, regulons #8 and #39), nine functional modules, and 33 genes with orthologous genes known to be associated with acid stress. The predicted response pathways could serve as promising candidates for better acid tolerance engineering in Lactococcus lactis. Our RECTA pipeline provides an effective way to construct a reliable gene regulatory network through regulon elucidation, and has strong application power and can be effectively applied to other bacterial genomes where the elucidation of the transcriptional regulation network is needed.

Phagocytosis and Epithelial Cell Invasion by Crohn's Disease-Associated Adherent-Invasive Escherichia coli Are Inhibited by the Anti-inflammatory Drug 6-Mercaptopurine

Adherent-invasive Escherichia coli (AIEC) strains are overrepresented in the dysbiotic microbiota of Crohn’s disease (CD) patients, and contribute to the onset of the chronic inflammation typical of the disease. However, the effects of anti-inflammatory drugs used for CD treatment on AIEC virulence have not yet been investigated. In this report, we show that exposure of AIEC LF82 strain to amino-6-mercaptopurine (6-MP) riboside, one of the most widely used anti-inflammatory drugs in CD, impairs its ability to adhere to, and consequently to invade, human epithelial cells. Notably, phagocytosis of LF82 treated with 6-MP by human macrophages is also reduced, suggesting that 6-MP affects AIEC cell surface determinants involved both in interaction with epithelial cells and in uptake by macrophages. Since a main target of 6-MP in bacterial cells is the inhibition of the important signal molecule c-di-GMP, we also tested whether perturbations in cAMP, another major signaling pathway in E. coli, might have similar effects on interactions with human cells. To this aim, we grew LF82 in the presence of glucose, which leads to inhibition of cAMP synthesis. Growth in glucose-supplemented medium resulted in a reduction in AIEC adhesion to epithelial cells and uptake by macrophages. Consistent with these results, both 6-MP and glucose can affect expression of cell adhesion-related genes, such as the csg genes, encoding thin aggregative fimbriae (curli). In addition, glucose strongly inhibits expression of the fim operon, encoding type 1 pili, a known AIEC determinant for adhesion to human cells. To further investigate whether 6-MP can indeed inhibit c-di-GMP signaling in AIEC, we performed biofilm and motility assays and determination of extracellular polysaccharides. 6-MP clearly affected biofilm formation and cellulose production, but also, unexpectedly, reduced cell motility, itself an important virulence factor for AIEC. Our results provide strong evidence that 6-MP can affect AIEC-host cell interaction by acting on the bacterial cell, thus strengthening the hypothesis that mercaptopurines might promote CD remission also by affecting gut microbiota composition and/or physiology, and suggesting that novel drugs targeting bacterial virulence and signaling might be effective in preventing chronic inflammation in CD.

Assessment of transfer methods for comparative genomics of regulatory networks in bacteria

Background

Comparative genomics can leverage the vast amount of available genomic sequences to reconstruct and analyze transcriptional regulatory networks in Bacteria, but the efficacy of this approach hinges on the ability to transfer regulatory network information from reference species to the genomes under analysis. Several methods have been proposed to transfer regulatory information between bacterial species, but the paucity and distributed nature of experimental information on bacterial transcriptional networks have prevented their systematic evaluation.

Results

We report the compilation of a large catalog of transcription factor-binding sites across Bacteria and its use to systematically benchmark proposed transfer methods across pairs of bacterial species. We evaluate motif- and accuracy-based metrics to assess the results of regulatory network transfer and we identify the precision-recall area-under-the-curve as the best metric for this purpose due to the large class-imbalanced nature of the problem. Methods assuming conservation of the transcription factor-binding motif (motif-based) are shown to substantially outperform those assuming conservation of regulon composition (network-based), even though their efficiency can decrease sharply with increasing phylogenetic distance. Variations of the basic motif-based transfer method do not yield significant improvements in transfer accuracy. Our results indicate that detection of a large enough number of regulated orthologs is critical for network-based transfer methods, but that relaxing orthology requirements does not improve results. Using the transcriptional regulators LexA and Fur as case examples, we also show how DNA-binding domain sequence similarity can yield confounding results as an indicator of transfer efficiency for motif-based methods.

Conclusions

Counter to standard practice, our evaluation of metrics to assess the efficiency of methods for regulatory network information transfer reveals that the area under precision-recall (PR) curves is a more precise and informative metric than that of receiver-operating-characteristic (ROC) curves, confirming similar findings in other class-imbalanced settings. Our systematic assessment of transfer methods reveals that simple approaches to both motif- and network-based transfer of regulatory information provide equal or better results than more elaborate methods. We also show that there are not effective predictors of transfer efficacy, substantiating the long-standing practice of manual curation in comparative genomics analyses.

Electronic supplementary material

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

Molecular and structural considerations of TF-DNA binding for the generation of biologically meaningful and accurate phylogenetic footprinting analysis: the LysR-type transcriptional regulator family as a study model

BACKGROUND

The goal of most programs developed to find transcription factor binding sites (TFBSs) is the identification of discrete sequence motifs that are significantly over-represented in a given set of sequences where a transcription factor (TF) is expected to bind. These programs assume that the nucleotide conservation of a specific motif is indicative of a selective pressure required for the recognition of a TF for its corresponding TFBS. Despite their extensive use, the accuracies reached with these programs remain low. In many cases, true TFBSs are excluded from the identification process, especially when they correspond to low-affinity but important binding sites of regulatory systems.

RESULTS

We developed a computational protocol based on molecular and structural criteria to perform biologically meaningful and accurate phylogenetic footprinting analyses. Our protocol considers fundamental aspects of the TF-DNA binding process, such as: i) the active homodimeric conformations of TFs that impose symmetric structures on the TFBSs, ii) the cooperative binding of TFs, iii) the effects of the presence or absence of co-inducers, iv) the proximity between two TFBSs or one TFBS and a promoter that leads to very long spurious motifs, v) the presence of AT-rich sequences not recognized by the TF but that are required for DNA flexibility, and vi) the dynamic order in which the different binding events take place to determine a regulatory response (i.e., activation or repression). In our protocol, the abovementioned criteria were used to analyze a profile of consensus motifs generated from canonical Phylogenetic Footprinting Analyses using a set of analysis windows of incremental sizes. To evaluate the performance of our protocol, we analyzed six members of the LysR-type TF family in Gammaproteobacteria.

CONCLUSIONS

The identification of TFBSs based exclusively on the significance of the over-representation of motifs in a set of sequences might lead to inaccurate results. The consideration of different molecular and structural properties of the regulatory systems benefits the identification of TFBSs and enables the development of elaborate, biologically meaningful and precise regulatory models that offer a more integrated view of the dynamics of the regulatory process of transcription.

A Bayesian inference method for the analysis of transcriptional regulatory networks in metagenomic data

Background

Metagenomics enables the analysis of bacterial population composition and the study of emergent population features, such as shared metabolic pathways. Recently, we have shown that metagenomics datasets can be leveraged to characterize population-wide transcriptional regulatory networks, or meta-regulons, providing insights into how bacterial populations respond collectively to specific triggers. Here we formalize a Bayesian inference framework to analyze the composition of transcriptional regulatory networks in metagenomes by determining the probability of regulation of orthologous gene sequences. We assess the performance of this approach on synthetic datasets and we validate it by analyzing the copper-homeostasis network of Firmicutes species in the human gut microbiome.

Results

Assessment on synthetic datasets shows that our method provides a robust and interpretable metric for assessing putative regulation by a transcription factor on sets of promoter sequences mapping to an orthologous gene cluster. The inference framework integrates the regulatory contribution of secondary sites and can discern false positives arising from multiple instances of a clonal sequence. Posterior probabilities for orthologous gene clusters decline sharply when less than 20 % of mapped promoters have binding sites, but we introduce a sensitivity adjustment procedure to speed up computation that enhances regulation assessment in heterogeneous ortholog clusters. Analysis of the copper-homeostasis regulon governed by CsoR in the human gut microbiome Firmicutes reveals that CsoR controls itself and copper-translocating P-type ATPases, but not CopZ-type copper chaperones. Our analysis also indicates that CsoR frequently targets promoters with dual CsoR-binding sites, suggesting that it exploits higher-order binding conformations to fine-tune its activity.

Conclusions

We introduce and validate a method for the analysis of transcriptional regulatory networks from metagenomic data that enables inference of meta-regulons in a systematic and interpretable way. Validation of this method on the CsoR meta-regulon of gut microbiome Firmicutes illustrates the usefulness of the approach, revealing novel properties of the copper-homeostasis network in poorly characterized bacterial species and putting forward evidence of new mechanisms of DNA binding for this transcriptional regulator. Our approach will enable the comparative analysis of regulatory networks across metagenomes, yielding novel insights into the evolution of transcriptional regulatory networks.

Electronic supplementary material

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

Bacterial regulon modeling and prediction based on systematic cis regulatory motif analyses

Regulons are the basic units of the response system in a bacterial cell, and each consists of a set of transcriptionally co-regulated operons. Regulon elucidation is the basis for studying the bacterial global transcriptional regulation network. In this study, we designed a novel co-regulation score between a pair of operons based on accurate operon identification and cis regulatory motif analyses, which can capture their co-regulation relationship much better than other scores. Taking full advantage of this discovery, we developed a new computational framework and built a novel graph model for regulon prediction. This model integrates the motif comparison and clustering and makes the regulon prediction problem substantially more solvable and accurate. To evaluate our prediction, a regulon coverage score was designed based on the documented regulons and their overlap with our prediction; and a modified Fisher Exact test was implemented to measure how well our predictions match the co-expressed modules derived from E. coli microarray gene-expression datasets collected under 466 conditions. The results indicate that our program consistently performed better than others in terms of the prediction accuracy. This suggests that our algorithms substantially improve the state-of-the-art, leading to a computational capability to reliably predict regulons for any bacteria.

Transcription factors with conserved binding sites near ATOH1 on the POU4F3 gene enhance the induction of cochlear hair cells

Overexpression of the transcription factor (TF) ATOH1 is known to induce the transformation of nonsensory cells in the organ of Corti into hair cells (HCs). Evaluating DNA 5Œ to the coding sequence of the pou4f3 gene, a target of ATOH1 in HCs, we identified in three regions containing clustered binding sites for ATOH1 and several other TFs that are expressed in developing inner ear sensory epithelia at the time of HC specification. These regions and sites are highly conserved across evolutionarily distant mammalian species. To test the hypothesis that the identified TFs act in combination to regulate the pou4f3 gene, we transfected by electroporation neonatal cochlear sensory epithelium from mice expressing green fluorescent protein (GFP) under the control of an 8.5-kb 5' pou4f3 genomic fragment. Plasmids encoding 21 TFs c-transfected with human ATOH1 (hATOH1). Cotransfection with hETV4, hNMYC, or hETS2 produced significantly more pou4f3/GFP and myosin 7A-positive nonsensory cells than hATOH1 alone. Co-transfection of hATOH1 with hHES1, hHES5, or hNEUROD1 reduced the effects of hATOH1. Chromatin immunoprecipitation (ChIP)of DNA from an inner ear cell line transfected with hNMYC,hETV4, or hETS2 revealed binding to a conserved region immediately proximal to the coding sequence. ChIP similarly revealed binding of hGATA3, hNMYC, and hTFE2 to a region several kilobases distal to the coding sequence, which we have previously shown to bind ATOH1. The results suggest that ATOH1 acts in concert with a subset of other TFs to directly regulate the pou4f3 gene and more broadly to regulate the HC phenotype.

A new insight into the role of intracellular nickel levels for the stress response, surface properties and twitching motility by Haemophilus influenzae

Nickel acts as a co-factor for a small number of enzymes in bacteria. Urease is one of the two nickel-dependent enzymes that have been identified in Haemophilus influenzae; glyoxalase I is the other. However, nickel has been suggested to have roles in H. influenzae that can not attributed to the function of these enzymes. We have previously shown that in the H. influenzae strain Rd KW20 the inability to acquire nickel led to alterations to the cell-type; an increased biofilm formation and changes in cell surface properties. Here we report the differences in the genome wide gene expression between Rd KW20 and a strain incapable of importing nickel (nikQ); revealing a link between intracellular nickel levels and genes involved in metabolic pathways, stress responses and genes associated with surface factors such as type IV pili. We have then taken a strain previously shown to use type IV pili both in biofilm formation and for twitching motility (86-028NP) and have shown its homologous genes (NTHI1417-1422; annotated as cobalt transporter, cbiKLMOQ) did import nickel and mutations in this locus had pleiotropic effects correlating to stress response and motility. Compared to wild type cells, the nickel depleted cells were more electronegativity charged, they aggregated and formed a biofilm. Correct intracellular nickel levels were also important for resistance to oxidative stress; the nickel depleted cells were more sensitive to oxidative stress. The nickel depleted cells were also non-motile, but the addition specifically of nickel returned these cells to a wild type motility state. We have also analysed the role of nickel uptake in a naturally, urease negative strain (the blood isolate R2866) and depleting intracellular nickel (a nikQ mutant) in this strain effected a similar range of cell functions. These data reveal a role for the capacity to acquire nickel from the environment and for the correct intracellular nickel levels as part of H. influenzae stress response and in signalling for a switch to a sessile bacterial lifestyle.

A discrete role for FNR in the transcriptional response to moderate changes in oxygen by Haemophilus influenzae Rd KW20

The survival by pathogenic bacteria within the specific conditions of an anatomical niche is critical for their persistence. These conditions include the combination of toxic chemicals, such as reactive oxygen (ROS) and reactive nitrogen species (RNS), with factors relevant to cell growth, such as oxygen. Haemophilus influenzae senses oxygen levels largely through the redox state of the intracellular fumarate-nitrate global regulator (FNR). H. influenzae certainly encounters oxygen levels that fluctuate, but in reality, these would rarely reach a state that results in FNR being fully reduced or oxidized. We were therefore interested in the response of H. influenzae to ROS and RNS at moderately high or low oxygen levels and the corresponding role of FNR. At these levels of oxygen, even though the growth rate of an H. influenzae fnr mutant was similar to wild type, its ROS and RNS tolerance was significantly different. Additionally, the subtle changes in oxygen did alter the whole cell transcriptional profile and this was different between the wild type and fnr mutant strains. It was the changed whole cell profile that impacted on ROS/RNS defence, but surprisingly, the FNR-regulated, anaerobic nitrite reductase (NrfA) continued to be expressed and had a role in this phenotype.

Effect of Global Regulators RpoS and Cyclic-AMP/CRP on the Catabolome and Transcriptome of Escherichia coli K12 during Carbon- and Energy-Limited Growth

For heterotrophic microbes, limited availability of carbon and energy sources is one of the major nutritional factors restricting the rate of growth in most ecosystems. Physiological adaptation to this hunger state requires metabolic versatility which usually involves expression of a wide range of different catabolic pathways and of high-affinity carbon transporters; together, this allows for simultaneous utilization of mixtures of carbonaceous compounds at low concentrations. In Escherichia coli the stationary phase sigma factor RpoS and the signal molecule cAMP are the major players in the regulation of transcription under such conditions; however, their interaction is still not fully understood. Therefore, during growth of E. coli in carbon-limited chemostat culture at different dilution rates, the transcriptomes, expression of periplasmic proteins and catabolomes of strains lacking one of these global regulators, either rpoS or adenylate cyclase (cya), were compared to those of the wild-type strain. The inability to synthesize cAMP exerted a strong negative influence on the expression of alternative carbon source uptake and degradation systems. In contrast, absence of RpoS increased the transcription of genes belonging to high-affinity uptake systems and central metabolism, presumably due to reduced competition of σ^D with σ^S. Phenotypical analysis confirmed this observation: The ability to respire alternative carbon substrates and to express periplasmic high-affinity binding proteins was eliminated in cya and crp mutants, while these properties were not affected in the rpoS mutant. As expected, transcription of numerous stress defence genes was negatively affected by the rpoS knock-out mutation. Interestingly, several genes of the RpoS stress response regulon were also down-regulated in the cAMP-negative strain indicating a coordinated global regulation. The results demonstrate that cAMP is crucial for catabolic flexibility during slow, carbon-limited growth, whereas RpoS is primarily involved in the regulation of stress response systems necessary for the survival of this bacterium under hunger conditions.

On the necessity and biological significance of threshold-free regulon prediction outputs

The in silico prediction of cis-acting elements in a genome is an efficient way to quickly obtain an overview of the biological processes controlled by a trans-acting factor, and connections between regulatory networks. Several regulon prediction web tools are available, designed to identify DNA motifs predicted to be bound by transcription factors using position weight matrix-based algorithms. In this paper we expose and discuss the conflicting objectives of software creators (bioinformaticians) and software users (biologists), who aim for reliable and exhaustive prediction outputs, respectively. Software makers, concerned with providing tools that minimise the number of false positive hits, often impose a stringent threshold score for a sequence to be included in the list of the putative cis-acting sites. This rigidity eventually results in the identification of strongly reliable but largely straightforward sites, i.e. those associated with genes already anticipated to be targeted by the studied transcription factor. Importantly, this biased identification of strongly bound sequences contrasts with the biological reality where, in many circumstances, a weak DNA-protein interaction is required for the appropriate gene's expression. We show here a series of transcriptionally controlled systems involving weakly bound cis-acting elements that could never have been discovered because of the policy of preventing software users from modifying the screening parameters. Proposing only trustworthy prediction outputs thus prevents biologists from fully utilising their knowledge background and deciding to analyse statistically irrelevant hits that could nonetheless be potentially involved in subtle, unexpected, though essential cis-trans relationships.

Comparative genomics and evolution of regulons of the LacI-family transcription factors

DNA-binding transcription factors (TFs) are essential components of transcriptional regulatory networks in bacteria. LacI-family TFs (LacI-TFs) are broadly distributed among certain lineages of bacteria. The majority of characterized LacI-TFs sense sugar effectors and regulate carbohydrate utilization genes. The comparative genomics approaches enable in silico identification of TF-binding sites and regulon reconstruction. To study the function and evolution of LacI-TFs, we performed genomics-based reconstruction and comparative analysis of their regulons. For over 1300 LacI-TFs from over 270 bacterial genomes, we predicted their cognate DNA-binding motifs and identified target genes. Using the genome context and metabolic subsystem analyses of reconstructed regulons, we tentatively assigned functional roles and predicted candidate effectors for 78 and 67% of the analyzed LacI-TFs, respectively. Nearly 90% of the studied LacI-TFs are local regulators of sugar utilization pathways, whereas the remaining 125 global regulators control large and diverse sets of metabolic genes. The global LacI-TFs include the previously known regulators CcpA in Firmicutes, FruR in Enterobacteria, and PurR in Gammaproteobacteria, as well as the three novel regulators—GluR, GapR, and PckR—that are predicted to control the central carbohydrate metabolism in three lineages of Alphaproteobacteria. Phylogenetic analysis of regulators combined with the reconstructed regulons provides a model of evolutionary diversification of the LacI protein family. The obtained genomic collection of in silico reconstructed LacI-TF regulons in bacteria is available in the RegPrecise database (http://regprecise.lbl.gov). It provides a framework for future structural and functional classification of the LacI protein family and identification of molecular determinants of the DNA and ligand specificity. The inferred regulons can be also used for functional gene annotation and reconstruction of sugar catabolic networks in diverse bacterial lineages.

There is a specific response to pH by isolates of Haemophilus influenzae and this has a direct influence on biofilm formation

Background

Haemophilus influenzae colonizes the nasopharynx as a commensal. Strain-specific factors allow some strains to migrate to particular anatomical niches, such as the middle ear, bronchi or blood, and induce disease by surviving within the conditions present at these sites in the body. It is established that H. influenzae colonization and in some cases survival is highly dependent on their ability to form a biofilm. Biofilm formation is a key trait in the development of chronic infection by certain isolates. This is exemplified by the contrast between the biofilm-forming strains found in middle ear infections and those isolates that survive within the blood and are rarely associated with biofilm development.

Results

Screening a group of H. influenzae strains revealed only slight variations in their growth across a range of pH conditions. However, some isolates responded to a pH of 8.0 by the formation of a biofilm. While the type b capsular blood isolate Eagan did not form a biofilm and grew at the same rate regardless of pH 6.8-8.0, transcriptomic analyses demonstrated that at pH 8.0 it uniquely induced a gluconate-uptake and metabolism pathway, which concurrently imports H⁺. A non-typeable H. influenzae, isolated from the middle ear, induced biofilm formation at pH 8.0, and at this pH it induced a series of iron acquisition genes, consistent with previous studies linking iron homeostasis to biofilm lifestyle.

Conclusions

Different strains of H. influenzae cope with changes in environmental factors using strain-specific mechanisms. These pathways define the scope and mode of niche-survival for an isolate. The pH is a property that is different from the middle ear (at least pH 8.0) compared to other sites that H. influenzae can colonize and infect. The transcriptional response to increasing pH by H. influenzae varies between strains, and pH is linked to pathways that allow strains to either continue free-living growth or induction of a biofilm. We showed that a biofilm-forming isolate induced iron metabolism pathways, whereas a strain that does not form biofilm at increasing pH induced mechanisms for growth and pH homeostasis based on sugar acid transport.

The cooperative binding energetics of CytR and cAMP receptor protein support a quantitative model of differential activation and repression of CytR-regulated class III Escherichia coli promoters

cAMP receptor protein (CRP) and CytR mediate positive and negative control of nine genes in Escherichia coli, most of which are involved in nucleoside catabolism and recycling. Five promoters share a common architecture in which tandem CRP sites flank an intervening CytR operator (CytO). CytR and CRP bind cooperatively to these promoters to form a three-protein, DNA-bound complex that controls activation and repression, the levels of which vary markedly among the promoters. To understand the specific combinatorial control mechanisms that are responsible for this outcome, we have used quantitative DNase I footprinting to generate individual site isotherms for each site of protein-DNA interaction. The intrinsic affinities of each transcription factor for its respective site and the specific patterns of cooperativity and competition underlying the molecular interactions at each promoter were determined by a global analysis of these titration data. Here we present results obtained for nupGP and tsxP2, adding to results published previously for deoP2, udpP, and cddP. These data allowed us to correlate the reported levels of activation, repression, and induction with the ligation states of these five promoters under physiologically relevant conditions. A general pattern of transcriptional regulation emerges that allows for complex patterns of regulation in this seemingly simple system.

RegPrecise 3.0--a resource for genome-scale exploration of transcriptional regulation in bacteria

Background

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in prokaryotes is one of the critical tasks of modern genomics. Bacteria from different taxonomic groups, whose lifestyles and natural environments are substantially different, possess highly diverged transcriptional regulatory networks. The comparative genomics approaches are useful for in silico reconstruction of bacterial regulons and networks operated by both transcription factors (TFs) and RNA regulatory elements (riboswitches).

Description

RegPrecise (http://regprecise.lbl.gov) is a web resource for collection, visualization and analysis of transcriptional regulons reconstructed by comparative genomics. We significantly expanded a reference collection of manually curated regulons we introduced earlier. RegPrecise 3.0 provides access to inferred regulatory interactions organized by phylogenetic, structural and functional properties. Taxonomy-specific collections include 781 TF regulogs inferred in more than 160 genomes representing 14 taxonomic groups of Bacteria. TF-specific collections include regulogs for a selected subset of 40 TFs reconstructed across more than 30 taxonomic lineages. Novel collections of regulons operated by RNA regulatory elements (riboswitches) include near 400 regulogs inferred in 24 bacterial lineages. RegPrecise 3.0 provides four classifications of the reference regulons implemented as controlled vocabularies: 55 TF protein families; 43 RNA motif families; ~150 biological processes or metabolic pathways; and ~200 effectors or environmental signals. Genome-wide visualization of regulatory networks and metabolic pathways covered by the reference regulons are available for all studied genomes. A separate section of RegPrecise 3.0 contains draft regulatory networks in 640 genomes obtained by an conservative propagation of the reference regulons to closely related genomes.

Conclusions

RegPrecise 3.0 gives access to the transcriptional regulons reconstructed in bacterial genomes. Analytical capabilities include exploration of: regulon content, structure and function; TF binding site motifs; conservation and variations in genome-wide regulatory networks across all taxonomic groups of Bacteria. RegPrecise 3.0 was selected as a core resource on transcriptional regulation of the Department of Energy Systems Biology Knowledgebase, an emerging software and data environment designed to enable researchers to collaboratively generate, test and share new hypotheses about gene and protein functions, perform large-scale analyses, and model interactions in microbes, plants, and their communities.

Comprehensive computational analysis of bacterial CRP/FNR superfamily and its target motifs reveals stepwise evolution of transcriptional networks

The cAMP receptor protein (CRP)/fumarate and nitrate reduction regulatory protein (FNR)-type transcription factors (TFs) are members of a well-characterized global TF family in bacteria and have two conserved domains: the N-terminal ligand-binding domain for small molecules (e.g., cAMP, NO, or O₂) and the C-terminal DNA-binding domain. Although the CRP/FNR-type TFs recognize very similar consensus DNA target sequences, they can regulate different sets of genes in response to environmental signals. To clarify the evolution of the CRP/FNR-type TFs throughout the bacterial kingdom, we undertook a comprehensive computational analysis of a large number of annotated CRP/FNR-type TFs and the corresponding bacterial genomes. Based on the amino acid sequence similarities among 1,455 annotated CRP/FNR-type TFs, spectral clustering classified the TFs into 12 representative groups, and stepwise clustering allowed us to propose a possible process of protein evolution. Although each cluster mainly consists of functionally distinct members (e.g., CRP, NTC, FNR-like protein, and FixK), FNR-related TFs are found in several groups and are distributed in a wide range of bacterial phyla in the sequence similarity network. This result suggests that the CRP/FNR-type TFs originated from an ancestral FNR protein, involved in nitrogen fixation. Furthermore, a phylogenetic profiling analysis showed that combinations of TFs and their target genes have fluctuated dynamically during bacterial evolution. A genome-wide analysis of TF-binding sites also suggested that the diversity of the transcriptional regulatory system was derived by the stepwise adaptation of TF-binding sites to the evolution of TFs.

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.

Hidden Markov models for evolution and comparative genomics analysis

The problem of reconstruction of ancestral states given a phylogeny and data from extant species arises in a wide range of biological studies. The continuous-time Markov model for the discrete states evolution is generally used for the reconstruction of ancestral states. We modify this model to account for a case when the states of the extant species are uncertain. This situation appears, for example, if the states for extant species are predicted by some program and thus are known only with some level of reliability; it is common for bioinformatics field. The main idea is formulation of the problem as a hidden Markov model on a tree (tree HMM, tHMM), where the basic continuous-time Markov model is expanded with the introduction of emission probabilities of observed data (e.g. prediction scores) for each underlying discrete state. Our tHMM decoding algorithm allows us to predict states at the ancestral nodes as well as to refine states at the leaves on the basis of quantitative comparative genomics. The test on the simulated data shows that the tHMM approach applied to the continuous variable reflecting the probabilities of the states (i.e. prediction score) appears to be more accurate then the reconstruction from the discrete states assignment defined by the best score threshold. We provide examples of applying our model to the evolutionary analysis of N-terminal signal peptides and transcription factor binding sites in bacteria. The program is freely available at http://bioinf.fbb.msu.ru/~nadya/tHMM and via web-service at http://bioinf.fbb.msu.ru/treehmmweb.

Genomic reconstruction of transcriptional regulatory networks in lactic acid bacteria

Background

Genome scale annotation of regulatory interactions and reconstruction of regulatory networks are the crucial problems in bacterial genomics. The Lactobacillales order of bacteria collates various microorganisms having a large economic impact, including both human and animal pathogens and strains used in the food industry. Nonetheless, no systematic genome-wide analysis of transcriptional regulation has been previously made for this taxonomic group.

Results

A comparative genomics approach was used for reconstruction of transcriptional regulatory networks in 30 selected genomes of lactic acid bacteria. The inferred networks comprise regulons for 102 orthologous transcription factors (TFs), including 47 novel regulons for previously uncharacterized TFs. Numerous differences between regulatory networks of the Streptococcaceae and Lactobacillaceae groups were described on several levels. The two groups are characterized by substantially different sets of TFs encoded in their genomes. Content of the inferred regulons and structure of their cognate TF binding motifs differ for many orthologous TFs between the two groups. Multiple cases of non-orthologous displacements of TFs that control specific metabolic pathways were reported.

Conclusions

The reconstructed regulatory networks substantially expand the existing knowledge of transcriptional regulation in lactic acid bacteria. In each of 30 studied genomes the obtained regulatory network contains on average 36 TFs and 250 target genes that are mostly involved in carbohydrate metabolism, stress response, metal homeostasis and amino acids biosynthesis. The inferred networks can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. All reconstructed regulons are captured within the Streptococcaceae and Lactobacillaceae collections in the RegPrecise database (http://regprecise.lbl.gov).

Inference of self-regulated transcriptional networks by comparative genomics

The assumption of basic properties, like self-regulation, in simple transcriptional regulatory networks can be exploited to infer regulatory motifs from the growing amounts of genomic and meta-genomic data. These motifs can in principle be used to elucidate the nature and scope of transcriptional networks through comparative genomics. Here we assess the feasibility of this approach using the SOS regulatory network of Gram-positive bacteria as a test case. Using experimentally validated data, we show that the known regulatory motif can be inferred through the assumption of self-regulation. Furthermore, the inferred motif provides a more robust search pattern for comparative genomics than the experimental motifs defined in reference organisms. We take advantage of this robustness to generate a functional map of the SOS response in Gram-positive bacteria. Our results reveal definite differences in the composition of the LexA regulon between Firmicutes and Actinobacteria, and confirm that regulation of cell-division inhibition is a widespread characteristic of this network among Gram-positive bacteria.

Binding Motifs in Bacterial Gene Promoters Modulate Transcriptional Effects of Global Regulators CRP and ArcA

Bacterial gene regulation involves transcription factors (TF) that bind to DNA recognition sequences in operon promoters. These recognition sequences, many of which are palindromic, are known as regulatory elements or transcription factor binding sites (TFBS). Some TFs are global regulators that can modulate the expression of hundreds of genes. In this study we examine global regulator half-sites, where a half-site, which we shall call a binding motif (BM), is one half of a palindromic TFBS. We explore the hypothesis that the number of BMs plays an important role in transcriptional regulation, examining empirical data from transcriptional profiling of the CRP and ArcA regulons. We compare the power of BM counts and of full TFBS characteristics to predict induced transcriptional activity. We find that CRP BM counts have a nonlinear effect on CRP-dependent transcriptional activity and predict this activity better than full TFBS quality or location.

Analysis of the SOS response of Vibrio and other bacteria with multiple chromosomes

Background

The SOS response is a well-known regulatory network present in most bacteria and aimed at addressing DNA damage. It has also been linked extensively to stress-induced mutagenesis, virulence and the emergence and dissemination of antibiotic resistance determinants. Recently, the SOS response has been shown to regulate the activity of integrases in the chromosomal superintegrons of the Vibrionaceae, which encompasses a wide range of pathogenic species harboring multiple chromosomes. Here we combine in silico and in vitro techniques to perform a comparative genomics analysis of the SOS regulon in the Vibrionaceae, and we extend the methodology to map this transcriptional network in other bacterial species harboring multiple chromosomes.

Results

Our analysis provides the first comprehensive description of the SOS response in a family (Vibrionaceae) that includes major human pathogens. It also identifies several previously unreported members of the SOS transcriptional network, including two proteins of unknown function. The analysis of the SOS response in other bacterial species with multiple chromosomes uncovers additional regulon members and reveals that there is a conserved core of SOS genes, and that specialized additions to this basic network take place in different phylogenetic groups. Our results also indicate that across all groups the main elements of the SOS response are always found in the large chromosome, whereas specialized additions are found in the smaller chromosomes and plasmids.

Conclusions

Our findings confirm that the SOS response of the Vibrionaceae is strongly linked with pathogenicity and dissemination of antibiotic resistance, and suggest that the characterization of the newly identified members of this regulon could provide key insights into the pathogenesis of Vibrio. The persistent location of key SOS genes in the large chromosome across several bacterial groups confirms that the SOS response plays an essential role in these organisms and sheds light into the mechanisms of evolution of global transcriptional networks involved in adaptability and rapid response to environmental changes, suggesting that small chromosomes may act as evolutionary test beds for the rewiring of transcriptional networks.

Prokaryotic genome regulation: a revolutionary paradigm

After determination of the whole genome sequence, the research frontier of bacterial molecular genetics has shifted to reveal the genome regulation under stressful conditions in nature. The gene selectivity of RNA polymerase is modulated after interaction with two groups of regulatory proteins, 7 sigma factors and 300 transcription factors. For identification of regulation targets of transcription factors in Escherichia coli, we have developed Genomic SELEX system and subjected to screening the binding sites of these factors on the genome. The number of regulation targets by a single transcription factor was more than those hitherto recognized, ranging up to hundreds of promoters. The number of transcription factors involved in regulation of a single promoter also increased to as many as 30 regulators. The multi-target transcription factors and the multi-factor promoters were assembled into complex networks of transcription regulation. The most complex network was identified in the regulation cascades of transcription of two master regulators for planktonic growth and biofilm formation.

BpsR modulates Bordetella biofilm formation by negatively regulating the expression of the Bps polysaccharide

Bordetella bacteria are Gram-negative respiratory pathogens of animals, birds, and humans. A hallmark feature of some Bordetella species is their ability to efficiently survive in the respiratory tract even after vaccination. Bordetella bronchiseptica and Bordetella pertussis form biofilms on abiotic surfaces and in the mouse respiratory tract. The Bps exopolysaccharide is one of the critical determinants for biofilm formation and the survival of Bordetella in the murine respiratory tract. In order to gain a better understanding of regulation of biofilm formation, we sought to study the mechanism by which Bps expression is controlled in Bordetella. Expression of bpsABCD (bpsA-D) is elevated in biofilms compared with levels in planktonically grown cells. We found that bpsA-D is expressed independently of BvgAS. Subsequently, we identified an open reading frame (ORF), BB1771 (designated here bpsR), that is located upstream of and in the opposite orientation to the bpsA-D locus. BpsR is homologous to the MarR family of transcriptional regulators. Measurement of bpsA and bpsD transcripts and the Bps polysaccharide levels from the wild-type and the ΔbpsR strains suggested that BpsR functions as a repressor. Consistent with enhanced production of Bps, the bpsR mutant displayed considerably more structured biofilms. We mapped the bpsA-D promoter region and showed that purified BpsR protein specifically bound to the bpsA-D promoter. Our results provide mechanistic insights into the regulatory strategy employed by Bordetella for control of the production of the Bps polysaccharide and biofilm formation.

Roles of hinge region, loops 3 and 4 in the activation of Escherichia coli cyclic AMP receptor protein

The cAMP receptor protein (CRP) requires cAMP for an allosteric change and regulates more than 150 genes in Escherichia coli. In this study, the modular half of cAMP receptor protein was used to investigate the allosteric signal transmission pathway induced by cAMP binding. The activation of CRP upon cAMP binding is indicated to be realignment of the two subunits within the CRP dimer. The interaction of loop 3 and Phe136 do not involve in signal transmission.

Genome-wide identification of in vivo binding sites of GlxR, a cyclic AMP receptor protein-type regulator in Corynebacterium glutamicum

Corynebacterium glutamicum GlxR is a cyclic AMP (cAMP) receptor protein-type regulator. Although over 200 GlxR-binding sites in the C. glutamicum genome are predicted in silico, studies on the physiological function of GlxR have been hindered by the severe growth defects of a glxR mutant. This study identified the GlxR regulon by chromatin immunoprecipitation in conjunction with microarray (ChIP-chip) analyses. In total, 209 regions were detected as in vivo GlxR-binding sites. In vitro binding assays and promoter-reporter assays demonstrated that GlxR directly activates expression of genes for aerobic respiration, ATP synthesis, and glycolysis and that it is required for expression of genes for cell separation and mechanosensitive channels. GlxR also directly represses a citrate uptake gene in the presence of citrate. Moreover, ChIP-chip analyses showed that GlxR was still able to interact with its target sites in a mutant with a deletion of cyaB, the sole adenylate cyclase gene in the genome, even though binding affinity was markedly decreased. Thus, GlxR is physiologically functional at the relatively low cAMP levels in the cyaB mutant, allowing the cyaB mutant to grow much better than the glxR mutant.

Escherichia coli exports cyclic AMP via TolC

In Escherichia coli more than 180 genes are regulated by the cyclic AMP (cAMP)-cAMP receptor protein (CRP) complex. However, more than 90% of cAMP that is made by intracellular adenylyl cyclases is found in the culture medium. How is cAMP exported from E. coli? In a tolC mutant, 0.03 mM IPTG (isopropyl-β-d-thiogalactopyranoside) was sufficient to induce β-galactosidase compared to 0.1 mM IPTG in the parent strain. In a cya mutant unable to produce cAMP about 1 mM extracellular cAMP was required to induce β-galactosidase, whereas in a cya tolC mutant 0.1 mM cAMP was sufficient. When cAMP in E. coli cya was generated intracellularly by a recombinant, weakly active adenylyl cyclase from Corynebacterium glutamicum, the critical level of cAMP necessary for induction of maltose degradation was only achieved in a tolC mutant and not in the parent strain. Deletion of a putative cAMP phosphodiesterase of E. coli, CpdA, resulted in a slightly similar, yet more diffuse phenotype. The data demonstrate that export of cAMP via TolC is a most efficient way of E. coli to lower high concentrations of cAMP in the cell and maintain its sensitivity in changing metabolic environments.

RegPredict: an integrated system for regulon inference in prokaryotes by comparative genomics approach

RegPredict web server is designed to provide comparative genomics tools for reconstruction and analysis of microbial regulons using comparative genomics approach. The server allows the user to rapidly generate reference sets of regulons and regulatory motif profiles in a group of prokaryotic genomes. The new concept of a cluster of co-regulated orthologous operons allows the user to distribute the analysis of large regulons and to perform the comparative analysis of multiple clusters independently. Two major workflows currently implemented in RegPredict are: (i) regulon reconstruction for a known regulatory motif and (ii) ab initio inference of a novel regulon using several scenarios for the generation of starting gene sets. RegPredict provides a comprehensive collection of manually curated positional weight matrices of regulatory motifs. It is based on genomic sequences, ortholog and operon predictions from the MicrobesOnline. An interactive web interface of RegPredict integrates and presents diverse genomic and functional information about the candidate regulon members from several web resources. RegPredict is freely accessible at http://regpredict.lbl.gov.

Activation of sigma 28-dependent transcription in Escherichia coli by the cyclic AMP receptor protein requires an unusual promoter organization

The Escherichia coli aer regulatory region contains a single promoter that is recognized by RNA polymerase containing the flagellar sigma factor, sigma(28). Expression from this promoter is dependent on direct activation by the cyclic AMP receptor protein, which binds to a target centred 49.5 base pairs upstream from the transcript start. Activator-dependent transcription from the aer promoter was reconstituted in vitro, and a tethered inorganic nuclease was used to find the position of the C-terminal domains of the RNA polymerase alpha subunits in transcriptionally competent open complexes. We report that the ternary activator--RNA polymerase--aer promoter open complex is organized differently from complexes at previously characterized promoters. Among other E. coli promoters recognized by RNA polymerase containing sigma(28), only the trg promoter is activated directly by the cyclic AMP receptor protein. The organization of the different promoter elements and the activator binding site at the trg promoter is the same as at the aer promoter, suggesting a common activation mechanism.

Genome-wide de novo prediction of cis-regulatory binding sites in prokaryotes

Although cis-regulatory binding sites (CRBSs) are at least as important as the coding sequences in a genome, our general understanding of them in most sequenced genomes is very limited due to the lack of efficient and accurate experimental and computational methods for their characterization, which has largely hindered our understanding of many important biological processes. In this article, we describe a novel algorithm for genome-wide de novo prediction of CRBSs with high accuracy. We designed our algorithm to circumvent three identified difficulties for CRBS prediction using comparative genomics principles based on a new method for the selection of reference genomes, a new metric for measuring the similarity of CRBSs, and a new graph clustering procedure. When operon structures are correctly predicted, our algorithm can predict 81% of known individual binding sites belonging to 94% of known cis-regulatory motifs in the Escherichia coli K12 genome, while achieving high prediction specificity. Our algorithm has also achieved similar prediction accuracy in the Bacillus subtilis genome, suggesting that it is very robust, and thus can be applied to any other sequenced prokaryotic genome. When compared with the prior state-of-the-art algorithms, our algorithm outperforms them in both prediction sensitivity and specificity.

Resistance of Haemophilus influenzae to reactive nitrogen donors and gamma interferon-stimulated macrophages requires the formate-dependent nitrite reductase regulator-activated ytfE gene

Haemophilus influenzae efficiently colonizes and persists at the human nasopharyngeal mucosa, causing disease when it spreads to other sites. Nitric oxide (NO) represents a major antimicrobial defense deployed by host cells in locations colonized by H. influenzae during pathogenesis that are likely to vary in oxygen levels. Formate-dependent nitrite reductase regulator (FNR) is an oxygen-sensitive regulator in several bacterial pathogens. We report that fnr of H. influenzae is required for anaerobic defense against exposure to NO donors and to resist NO-dependent effects of gamma interferon (IFN-gamma)-activated murine bone marrow-derived macrophages. To understand the mechanism of resistance, we investigated the role of FNR-regulated genes in defense against NO sources. Expression analysis revealed FNR-dependent activation of nrfA, dmsA, napA, and ytfE. Nonpolar deletion mutants of nrfA and ytfE exhibited sensitivity to NO donors, and the ytfE gene was more critical for survival. Compared to the wild-type strain, the ytfE mutant exhibited decreased survival when exposed to macrophages, a defect that was more pronounced after prior stimulation of macrophages with IFN-gamma or lipopolysaccharide. Complementation restored survival of the mutant to the level in the parental strain. Increased sensitivity of the ytfE mutant relative to that of the parent was abrogated by treatment of macrophages with a NO synthase inhibitor, implicating YtfE in resistance to a NO-dependent pathway. These results identify a requirement for FNR in positive control of ytfE and indicate a critical role for ytfE in resistance of H. influenzae to reactive nitrogen species and the antibacterial effects of macrophages.

The iron/heme regulated genes of Haemophilus influenzae: comparative transcriptional profiling as a tool to define the species core modulon

BACKGROUND

Haemophilus influenzae requires heme for aerobic growth and possesses multiple mechanisms to obtain this essential nutrient. Although an understanding of the heme acquisition mechanisms of H. influenzae is emerging, significant gaps in our knowledge remain. Unresolved issues include the identities of all genes exhibiting altered transcription in response to iron and heme availability, the fraction of such genes functioning in iron/heme acquisition, and the heterogeneity of this gene set among clinical isolates. Previously we utilized H. influenzae strain Rd KW20 to demonstrate the utility of transcriptional profiling in defining the genes exhibiting altered transcription in response to environmental iron and heme levels. The current study expands upon those observations by determining the iron/heme modulons of two clinical isolates, the type b isolate 10810 and the nontypeable isolate R2866. These data are used to begin to define the core iron/heme modulon of the species.

RESULTS

Microarray studies were performed to compare gene expression on transition from iron/heme-restricted to iron/heme-replete conditions for each isolate. Of 1820 ORFs on the array corresponding to R2866 genes, 363 were significantly differentially expressed: 233 were maximally transcribed under iron/heme-replete conditions and 130 under iron/heme-restricted conditions. Of the 1883 ORFs representing genes of strain 10810, 353 were significantly differentially transcribed: 150 were preferentially transcribed under iron/heme-replete conditions and 203 under iron/heme-restricted conditions. Comparison of the data sets indicated that 163 genes exhibited similar regulation in both isolates and that 74 of these exhibited similar patterns of regulation in Rd KW20. These comprise the putative core iron/heme modulon.

CONCLUSION

This study provides evidence for a conserved core of H. influenzae genes the transcription of which is altered by the availability of iron and/or heme in the growth environment. Elucidation of this modulon provides a means to identify genes with unrecognized roles in iron/heme acquisition or homeostasis, unanticipated responsiveness to environmental levels of the micronutrients or potential roles in virulence. Defining these core genes is also of potential importance in identifying targets for therapeutic and vaccine designs since products of these genes are likely to be preferentially expressed during growth in iron/heme restricted sites of the human body.

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

Predicting sigma28 promoters in eleven Shewanella genomes

An iterative position-specific score matrix (PSSM)-based approach was used to predict sigma(28) promoters in 11 Shewanella genomes. The Shewanella Correlation Browser was used to distinguish true-positive predictions from false-positive predictions in Shewanella oneidensis MR-1 by generating a sigma(28)-regulated transcriptional network from transcriptional profiling data. This dual-pronged approach identified several genes that have sigma(28) promoters and that may be involved with motility or chemotaxis in Shewanella.

Impact of Transcription Units rearrangement on the evolution of the regulatory network of gamma-proteobacteria

BACKGROUND

In the past years, several studies begun to unravel the structure, dynamical properties, and evolution of transcriptional regulatory networks. However, even those comparative studies that focus on a group of closely related organisms are limited by the rather scarce knowledge on regulatory interactions outside a few model organisms, such as E. coli among the prokaryotes.

RESULTS

In this paper we used the information annotated in Tractor_DB (a database of regulatory networks in gamma-proteobacteria) to calculate a normalized Site Orthology Score (SOS) that quantifies the conservation of a regulatory link across thirty genomes of this subclass. Then we used this SOS to assess how regulatory connections have evolved in this group, and how the variation of basic regulatory connection is reflected on the structure of the chromosome. We found that individual regulatory interactions shift between different organisms, a process that may be described as rewiring the network. At this evolutionary scale (the gamma-proteobacteria subclass) this rewiring process may be an important source of variation of regulatory incoming interactions for individual networks. We also noticed that the regulatory links that form feed forward motifs are conserved in a better correlated manner than triads of random regulatory interactions or pairs of co-regulated genes. Furthermore, the rewiring process that takes place at the most basic level of the regulatory network may be linked to rearrangements of genetic material within bacterial chromosomes, which change the structure of Transcription Units and therefore the regulatory connections between Transcription Factors and structural genes.

CONCLUSION

The rearrangements that occur in bacterial chromosomes-mostly inversion or horizontal gene transfer events - are important sources of variation of gene regulation at this evolutionary scale.

PREDetector: a new tool to identify regulatory elements in bacterial genomes

In the post-genomic area, the prediction of transcription factor regulons by position weight matrix-based programmes is a powerful approach to decipher biological pathways and to modelize regulatory networks in bacteria. The main difficulty once a regulon prediction is available is to estimate its reliability prior to start expensive experimental validations and therefore trying to find a way how to identify true positive hits from an endless list of potential target genes of a regulatory protein. Here we introduce PREDetector (Prokaryotic Regulatory Elements Detector), a tool developed for predicting regulons of DNA-binding proteins in bacterial genomes that, beside the automatic prediction, scoring and positioning of potential binding sites and their respective target genes in annotated bacterial genomes, it also provides an easy way to estimate the thresholds where to find reliable possible new target genes. PREDetector can be downloaded freely at http://www.montefiore.ulg.ac.be/~hiard/PreDetector/PreDetector.php.

Comparative genomics - a perspective

The rapidly emerging field of comparative genomics has yielded dramatic results. Comparative genome analysis has become feasible with the availability of a number of completely sequenced genomes. Comparison of complete genomes between organisms allow for global views on genome evolution and the availability of many completely sequenced genomes increases the predictive power in deciphering the hidden information in genome design, function and evolution. Thus, comparison of human genes with genes from other genomes in a genomic landscape could help assign novel functions for un-annotated genes. Here, we discuss the recently used techniques for comparative genomics and their derived inferences in genome biology.