PredictRegulon: a web server for the prediction of the regulatory protein binding sites and operons in prokaryote genomes

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.

Mining for Microbial Gems: Integrating Proteomics in the Postgenomic Natural Product Discovery Pipeline

Natural products (NPs) are a major source of compounds for medical, agricultural, and biotechnological industries. Many of these compounds are of microbial origin, and, in particular, from Actinobacteria or filamentous fungi. To successfully identify novel compounds that correlate to a bioactivity of interest, or discover new enzymes with desired functions, systematic multiomics approaches have been developed over the years. Bioinformatics tools harness the rapidly expanding wealth of genome sequence information, revealing previously unsuspected biosynthetic diversity. Varying growth conditions or application of elicitors are applied to activate cryptic biosynthetic gene clusters, and metabolomics provide detailed insights into the NPs they specify. Combining these technologies with proteomics-based approaches to profile the biosynthetic enzymes provides scientists with insights into the full biosynthetic potential of microorganisms. The proteomics approaches include enrichment strategies such as employing activity-based probes designed by chemical biology, as well as unbiased (quantitative) proteomics methods. In this review, the opportunities and challenges in microbial NP research are discussed, and, in particular, the application of proteomics to link biosynthetic enzymes to the molecules they produce, and vice versa.

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.

Phylogenetic footprinting: a boost for microbial regulatory genomics

Phylogenetic footprinting is a method for the discovery of regulatory elements in a set of homologous regulatory regions, usually collected from multiple species. It does so by identifying the best conserved motifs in those homologous regions. There are two popular sets of methods-alignment-based and motif-based, which are generally employed for phylogenetic methods. However, serious efforts have lacked to develop a tool exclusively for phylogenetic footprinting, based on either of these methods. Nevertheless, a number of software and tools exist that can be applied for prediction of phylogenetic footprinting with variable degree of success. The output from these tools may get affected by a number of factors associated with current state of knowledge, techniques and other resources available. We here present a critical apprehension of various phylogenetic approaches with reference to prokaryotes outlining the available resources and also discussing various factors affecting footprinting in order to make a clear idea about the proper use of this approach on prokaryotes.

Peroxide stimulon and role of PerR in group A Streptococcus

We have characterized group A Streptococcus (GAS) genome-wide responses to hydrogen peroxide and assessed the role of the peroxide response regulator (PerR) in GAS under oxidative stress. Comparison of transcriptome changes elicited by peroxide in wild-type bacteria with those in a perR deletion mutant showed that 76 out of 237 peroxide-regulated genes are PerR dependent. Unlike the PerR-mediated upregulation of peroxidases and other peroxide stress defense mechanisms previously reported in gram-positive species, PerR-dependent genes in GAS were almost exclusively downregulated and encoded proteins involved in purine and deoxyribonucleotide biosynthesis, heme uptake, and amino acid/peptide transport, but they also included a strongly activated putative transcriptional regulator (SPy1198). Of the 161 PerR-independent loci, repressed genes (86 of 161) encoded proteins with functions similar to those coordinated by PerR, in contrast to upregulated loci that encoded proteins that function in DNA damage repair, cofactor metabolism, reactive oxygen species detoxification, pilus biosynthesis, and hypothetical proteins. Complementation of the perR deletion mutant with wild-type PerR restored PerR-dependent regulation, whereas complementation with either one of two PerR variants carrying single mutations in two predicted metal-binding sites did not rescue the mutant phenotype. Metal content analyses of the recombinant wild type and respective PerR mutants, in addition to regulation studies in metal-supplemented and iron-depleted media, showed binding of zinc and iron by PerR and an iron requirement for optimal responses to peroxide. Our findings reveal a novel physiological contribution of PerR in coordinating DNA and protein metabolic functions in peroxide and identify GAS adaptive responses that may serve to enhance oxidative stress resistance and virulence in the host.

A novel metal transporter mediating manganese export (MntX) regulates the Mn to Fe intracellular ratio and Neisseria meningitidis virulence

Neisseria meningitidis (Nm) and N. gonorrhoeae (Ng) are adapted to different environments within their human host. If the basis of this difference has not yet been fully understood, previous studies (including our own data) have reported that, unlike Ng, Nm tolerates high manganese concentrations. As transition metals are essential regulators of cell growth and host pathogen interactions, we aimed to address mechanisms of Nm Mn²⁺ tolerance and its pathogenic consequences. Using bioinformatics, gene deletion and heterologous expression we identified a conserved bacterial manganese resistance factor MntX (formerly YebN). The predicted structure suggests that MntX represents a new family of transporters exporting Mn. In the Neisseria genus, this exporter is present and functional in all Nm isolates but it is mutated in a majority of Ng strains and commonly absent in nonpathogenic species. In Nm, Mn²⁺ export via MntX regulates the intracellular Mn/Fe ratio and protects against manganese toxicity that is exacerbated in low iron conditions. MntX is also important for N. meningitidis to resist killing by human serum and for survival in mice blood during septicemia. The present work thus points to new clues about Mn homeostasis, its interplay with Fe metabolism and the influence on N. meningitidis physiology and pathogenicity.

PromBase: a web resource for various genomic features and predicted promoters in prokaryotic genomes

Background

As more and more genomes are being sequenced, an overview of their genomic features and annotation of their functional elements, which control the expression of each gene or transcription unit of the genome, is a fundamental challenge in genomics and bioinformatics.

Findings

Relative stability of DNA sequence has been used to predict promoter regions in 913 microbial genomic sequences with GC-content ranging from 16.6% to 74.9%. Irrespective of the genome GC-content the relative stability based promoter prediction method has already been proven to be robust in terms of recall and precision. The predicted promoter regions for the 913 microbial genomes have been accumulated in a database called PromBase. Promoter search can be carried out in PromBase either by specifying the gene name or the genomic position. Each predicted promoter region has been assigned to a reliability class (low, medium, high, very high and highest) based on the difference between its average free energy and the downstream region. The recall and precision values for each class are shown graphically in PromBase. In addition, PromBase provides detailed information about base composition, CDS and CG/TA skews for each genome and various DNA sequence dependent structural properties (average free energy, curvature and bendability) in the vicinity of all annotated translation start sites (TLS).

Conclusion

PromBase is a database, which contains predicted promoter regions and detailed analysis of various genomic features for 913 microbial genomes. PromBase can serve as a valuable resource for comparative genomics study and help the experimentalist to rapidly access detailed information on various genomic features and putative promoter regions in any given genome. This database is freely accessible for academic and non- academic users via the worldwide web http://nucleix.mbu.iisc.ernet.in/prombase/.

Genetic and biochemical analysis of PadR-padC promoter interactions during the phenolic acid stress response in Bacillus subtilis 168

Bacillus subtilis 168 is resistant to phenolic acids by expression of an inducible enzyme, the phenolic acid decarboxylase (PadC), that decarboxylates these acids into less toxic vinyl derivatives. In the phenolic acid stress response (PASR), the repressor of padC, PadR, is inactivated by these acids. Inactivation of PadR is followed by a strong expression of padC. To elucidate the functional interaction between PadR and the padC promoter, we performed (i) footprinting assays to identify the region protected by PadR, (ii) electrophoretic mobility shift assays (EMSAs) with a modified padC promoter protected region to determine the interacting sequences, and (iii) random mutagenesis of padR to identify amino acid residues essential for the function of PadR. We identified an important consensus dyad sequence called IR1-2 (ATGT-8N-ACAT) overlapping a second dyad element (GTGT-8N-ACAT) that we named dIR1-2bis. The entire dIR1-2bis/IR1-2 sequence permits binding of two PadR dimers in EMSAs, which may be observed for bacteria grown under noninduced conditions where the padC promoter is completely repressed. Three groups of modified PadRs giving a PASR phenotype were characterized in vivo. The DNA sequences of certain mutant padR alleles indicate that important residues are all located in the region containing the coiled-coil leucine zipper domain that is involved in dimerization. These substitutions reduce the affinity of PadR binding to the padC promoter. Of particular interest are residue L128, located at the center of the putative coiled-coil leucine zipper domain, and residue E97, which is conserved among all PadRs.

Unraveling networks of co-regulated genes on the sole basis of genome sequences

With the growing number of available microbial genome sequences, regulatory signals can now be revealed as conserved motifs in promoters of orthologous genes (phylogenetic footprints). A next challenge is to unravel genome-scale regulatory networks. Using as sole input genome sequences, we predicted cis-regulatory elements for each gene of the yeast Saccharomyces cerevisiae by discovering over-represented motifs in the promoters of their orthologs in 19 Saccharomycetes species. We then linked all genes displaying similar motifs in their promoter regions and inferred a co-regulation network including 56,919 links between 3171 genes. Comparison with annotated regulons highlights the high predictive value of the method: a majority of the top-scoring predictions correspond to already known co-regulations. We also show that this inferred network is as accurate as a co-expression network built from hundreds of transcriptome microarray experiments. Furthermore, we experimentally validated 14 among 16 new functional links between orphan genes and known regulons. This approach can be readily applied to unravel gene regulatory networks from hundreds of microbial genomes for which no other information is available except the sequence. Long-term benefits can easily be perceived when considering the exponential increase of new genome sequences.

FITBAR: a web tool for the robust prediction of prokaryotic regulons

Background

The binding of regulatory proteins to their specific DNA targets determines the accurate expression of the neighboring genes. The in silico prediction of new binding sites in completely sequenced genomes is a key aspect in the deeper understanding of gene regulatory networks. Several algorithms have been described to discriminate against false-positives in the prediction of new binding targets; however none of them has been implemented so far to assist the detection of binding sites at the genomic scale.

Results

FITBAR (Fast Investigation Tool for Bacterial and Archaeal Regulons) is a web service designed to identify new protein binding sites on fully sequenced prokaryotic genomes. This tool consists in a workbench where the significance of the predictions can be compared using different statistical methods, a feature not found in existing resources. The Local Markov Model and the Compound Importance Sampling algorithms have been implemented to compute the P-value of newly discovered binding sites. In addition, FITBAR provides two optimized genomic scanning algorithms using either log-odds or entropy-weighted position-specific scoring matrices. Other significant features include the production of a detailed genomic context map for each detected binding site and the export of the search results in spreadsheet and portable document formats. FITBAR discovery of a high affinity Escherichia coli NagC binding site was validated experimentally in vitro as well as in vivo and published.

Conclusions

FITBAR was developed in order to allow fast, accurate and statistically robust predictions of prokaryotic regulons. This feature constitutes the main advantage of this web tool over other matrix search programs and does not impair its performance. The web service is available at http://archaea.u-psud.fr/fitbar.

Phyloscan: locating transcription-regulating binding sites in mixed aligned and unaligned sequence data

The transcription of a gene from its DNA template into an mRNA molecule is the first, and most heavily regulated, step in gene expression. Especially in bacteria, regulation is typically achieved via the binding of a transcription factor (protein) or small RNA molecule to the chromosomal region upstream of a regulated gene. The protein or RNA molecule recognizes a short, approximately conserved sequence within a gene's promoter region and, by binding to it, either enhances or represses expression of the nearby gene. Since the sought-for motif (pattern) is short and accommodating to variation, computational approaches that scan for binding sites have trouble distinguishing functional sites from look-alikes. Many computational approaches are unable to find the majority of experimentally verified binding sites without also finding many false positives. Phyloscan overcomes this difficulty by exploiting two key features of functional binding sites: (i) these sites are typically more conserved evolutionarily than are non-functional DNA sequences; and (ii) these sites often occur two or more times in the promoter region of a regulated gene. The website is free and open to all users, and there is no login requirement. Address: (http://bayesweb.wadsworth.org/phyloscan/).

Identification of ColR binding consensus and prediction of regulon of ColRS two-component system

BACKGROUND

Conserved two-component system ColRS of Pseudomonas genus has been implicated in several unrelated phenotypes. For instance, deficiency of P. putida ColRS system results in lowered phenol tolerance, hindrance of transposition of Tn4652 and lysis of a subpopulation of glucose-grown bacteria. In order to discover molecular mechanisms behind these phenotypes, we focused here on identification of downstream components of ColRS signal transduction pathway.

RESULTS

First, highly similar ColR binding sites were mapped upstream of outer membrane protein-encoding oprQ and a putative methyltransferase-encoding PP0903. These two ColR binding sequences were used as an input in computational genome-wide screening for new potential ColR recognition boxes upstream of different genes in P. putida. Biological relevance of a set of in silico predicted ColR-binding sites was analysed in vivo by studying the effect of ColR on transcription from promoters carrying these sites. This analysis disclosed seven novel genes of which six were positively and one negatively regulated by ColR. Interestingly, all promoters tested responded more significantly to the over-expression than to the absence of ColR suggesting that either ColR is limiting or ColS-activating signal is low under the conditions applied. The binding sites of ColR in the promoters analysed were validated by gel mobility shift and/or DNase I footprinting assays. ColR binding consensus was defined according to seven ColR binding motifs mapped by DNase I protection assay and this consensus was used to predict minimal regulon of ColRS system.

CONCLUSION

Combined usage of experimental and computational approach enabled us to define the binding consensus for response regulator ColR and to discover several new ColR-regulated genes. For instance, genes of outer membrane lipid A 3-O-deacylase PagL and cytoplasmic membrane diacylglycerol kinase DgkA are the members of ColR regulon. Furthermore, over 40 genes were predicted to be putatively controlled by ColRS two-component system in P. putida. It is notable that many of ColR-regulated genes encode membrane-related products thus confirming the previously proposed role of ColRS system in regulation of membrane functionality.

IdeR in Mycobacteria: From Target Recognition to Physiological Function

In mycobacteria, iron dependent transcription regulator (IdeR) regulates transcription of genes in response to iron levels. The IdeR regulated genes have been investigated mostly in M. tuberculosis, M. smegmatis, and in few of the other related species. Recent advances in crystal structure solution and computational as well as experimental identification of IdeR targets has provided insight into IdeR structure and function. Here in this review we take stock of current state of knowledge on IdeR and its targets to understand the underlying design of the IdeR regulon and its role in mycobacterial physiology.

PerR confers phagocytic killing resistance and allows pharyngeal colonization by group A Streptococcus

The peroxide response transcriptional regulator, PerR, is thought to contribute to virulence of group A Streptococcus (GAS); however, the specific mechanism through which it enhances adaptation for survival in the human host remains unknown. Here, we identify a critical role of PerR-regulated gene expression in GAS phagocytosis resistance and in virulence during pharyngeal infection. Deletion of perR in M-type 3 strain 003Sm was associated with reduced resistance to phagocytic killing in human blood and by murine macrophages in vitro. The increased phagocytic killing of the perR mutant was abrogated in the presence of the general oxidative burst inhibitor diphenyleneiodonium chloride (DPI), a result that suggests PerR-dependent gene expression counteracts the phagocyte oxidative burst. Moreover, an isogenic perR mutant was severely attenuated in a baboon model of GAS pharyngitis. In competitive infection experiments, the perR mutant was cleared from two animals at 24 h and from four of five animals by day 14, in sharp contrast to wild-type bacteria that persisted in the same five animals for 28 to 42 d. GAS genomic microarrays were used to compare wild-type and perR mutant transcriptomes in order to characterize the PerR regulon of GAS. These studies identified 42 PerR-dependent loci, the majority of which had not been previously recognized. Surprisingly, a large proportion of these loci are involved in sugar utilization and transport, in addition to oxidative stress adaptive responses and virulence. This finding suggests a novel role for PerR in mediating sugar uptake and utilization that, together with phagocytic killing resistance, may contribute to GAS fitness in the infected host. We conclude that PerR controls expression of a diverse regulon that enhances GAS resistance to phagocytic killing and allows adaptation for survival in the pharynx.

Detecting cis-regulatory binding sites for cooperatively binding proteins

Several methods are available to predict cis-regulatory modules in DNA based on position weight matrices. However, the performance of these methods generally depends on a number of additional parameters that cannot be derived from sequences and are difficult to estimate because they have no physical meaning. As the best way to detect cis-regulatory modules is the way in which the proteins recognize them, we developed a new scoring method that utilizes the underlying physical binding model. This method requires no additional parameter to account for multiple binding sites; and the only necessary parameters to model homotypic cooperative interactions are the distances between adjacent protein binding sites in basepairs, and the corresponding cooperative binding constants. The heterotypic cooperative binding model requires one more parameter per cooperatively binding protein, which is the concentration multiplied by the partition function of this protein. In a case study on the bacterial ferric uptake regulator, we show that our scoring method for homotypic cooperatively binding proteins significantly outperforms other PWM-based methods where biophysical cooperativity is not taken into account.

A generic approach to identify Transcription Factor-specific operator motifs; Inferences for LacI-family mediated regulation in Lactobacillus plantarum WCFS1

Background

A key problem in the sequence-based reconstruction of regulatory networks in bacteria is the lack of specificity in operator predictions. The problem is especially prominent in the identification of transcription factor (TF) specific binding sites. More in particular, homologous TFs are abundant and, as they are structurally very similar, it proves difficult to distinguish the related operators by automated means. This also holds for the LacI-family, a family of TFs that is well-studied and has many members that fulfill crucial roles in the control of carbohydrate catabolism in bacteria including catabolite repression. To overcome the specificity problem, a comprehensive footprinting approach was formulated to identify TF-specific operator motifs and was applied to the LacI-family of TFs in the model gram positive organism, Lactobacillus plantarum WCFS1. The main premise behind the approach is that only orthologous sequences that share orthologous genomic context will share equivalent regulatory sites.

Results

When the approach was applied to the 12 LacI-family TFs of the model species, a specific operator motif was identified for each of them. With the TF-specific operator motifs, potential binding sites were found on the genome and putative minimal regulons could be defined. Moreover, specific inducers could in most cases be linked to the TFs through phylogeny, thereby unveiling the biological role of these regulons. The operator predictions indicated that the LacI-family TFs can be separated into two subfamilies with clearly distinct operator motifs. They also established that the operator related to the 'global' regulator CcpA is not inherently distinct from that of other LacI-family members, only more degenerate. Analysis of the chromosomal position of the identified putative binding sites confirmed that the LacI-family TFs are mostly auto-regulatory and relate mainly to carbohydrate uptake and catabolism.

Conclusion

Our approach to identify specific operator motifs for different TF-family members is specific and in essence generic. The data infer that, although the specific operator motifs can be used to identify minimal regulons, experimental knowledge on TF activity especially is essential to determine complete regulons as well as to estimate the overlap between TF affinities.

Evolution of the mycobacterial SigK regulon

Previous studies have established that members of the Mycobacterium tuberculosis complex exhibit variable production of the antigenic proteins MPT70 and MPT83 due to mutations in their positive regulator, SigK (sigma factor K), and their negative regulator, RskA (regulator of sigma K). To further understand this highly specific SigK-controlled regulon, we have undertaken evolutionary studies to determine the presence of homologues of SigK-regulated genes in other organisms and to predict its transcriptional network. Evolutionary analysis indicates that the positive and negative regulators are conserved across many organisms, but that the genes under their control are variable. Moreover, the addition, loss, and movement of various genes in the mpt70/83 locus suggest that these genes are unlikely to be cotranscribed. To test predictions from sequence analysis, we have used promoter luciferase fusions and Northern blots to show that the majority of genes in this locus have their own promoters, of which a subset are SigK regulated (mpt83, dipZ, mpt70, and Rv0449c). Next, we have shown that the intracellular inducibility of mpt70 and mpt83 is a conserved property, shared between M. tuberculosis and Mycobacterium marinum. In addition, we have shown that SigK and RskA from an environmental mycobacterium isolate (M. gilvum PYR-GCK) complemented the regulatory activity of M. tuberculosis delta sigK rskA. Together, our data indicate that the regulatory system SigK/RskA is conserved across the Mycobacterium genus, whereas the regulon under its control varies considerably across species.

The Sulfolobus solfataricus radA paralogue sso0777 is DNA damage inducible and positively regulated by the Sta1 protein

Little is known about the regulation of the DNA damage-mediated gene expression in archaea. Here we report that the addition of actinomycin D to Sulfolobus solfataricus cultures triggers the expression of the radA paralogue sso0777. Furthermore, a specific retarded band is observed when electrophoretic mobility shift assays (EMSAs) with crude S. solfataricus cell extracts and the sso0777 promoter were carried out. The protein that binds to this promoter was isolated and identified as Sta1. Footprinting experiments have shown that the Sta1 DNA-binding site is included in the ATTTTTTATTTTCACATGTAAGATGTTTATT sequence, which is located upstream the putative TTG translation starting codon of the sso0777 gene. Additionally, gel electrophoretic mobility retardation experiments using mutant sso0777 promoter derivatives show the presence of three essential motifs (TTATT, CANGNA and TTATT) that are absolutely required for Sta1 DNA binding. Finally, in vitro transcription experiments confirm that Sta1 functions as an activator for sso0777 gene expression being the first identified archaeal regulatory protein associated with the DNA damage-mediated induction of gene expression.

Identification of the genes directly controlled by the response regulator CiaR in Streptococcus pneumoniae: five out of 15 promoters drive expression of small non-coding RNAs

The two-component regulatory system CiaRH of Streptococcus pneumoniae has been implicated in beta-lactam resistance, maintenance of cell integrity, competence and virulence, but the genes that are regulated directly by the system have not been defined. Using transcriptional mapping, in vitro CiaR binding, and in vivo analysis of CiaR-mediated regulation, 15 promoters were identified to be directly controlled by the response regulator CiaR. A direct repeat, TTTAAG-N5-TTTAAG, was found to be essential for CiaR binding and regulation. It is present, either completely or with subtle changes, in all promoter regions. Fourteen promoters of the regulon are activated by CiaR, and one was found to be controlled negatively. The genes that are transcribed from these promoters included ciaRH, loci that are predicted to be involved in the modification of teichoic acids (lic), in sugar metabolism (mal, man), stress response (htrA), chromosome segregation (parB), protease maturation (ppmA) and unknown functions. Remarkably, the five strongest promoters of the CiaR regulon drive expression of small RNAs. These small RNAs, designated csRNAs for cia-dependent small RNAs, are non-coding, between 87 and 151 nt in size, and show a high degree of similarity to each other. The analysis of deletion mutants in the csRNA genes revealed that csRNA4 and csRNA5 affect stationary-phase autolysis. The identification of five small non-coding regulatory RNAs opens new perspectives to approach the physiological role of the CiaRH two-component regulatory system.

GntR family of regulators in Mycobacterium smegmatis: a sequence and structure based characterization

Background

Mycobacterium smegmatis is fast growing non-pathogenic mycobacteria. This organism has been widely used as a model organism to study the biology of other virulent and extremely slow growing species like Mycobacterium tuberculosis. Based on the homology of the N-terminal DNA binding domain, the recently sequenced genome of M. smegmatis has been shown to possess several putative GntR regulators. A striking characteristic feature of this family of regulators is that they possess a conserved N-terminal DNA binding domain and a diverse C-terminal domain involved in the effector binding and/or oligomerization. Since the physiological role of these regulators is critically dependent upon effector binding and operator sites, we have analysed and classified these regulators into their specific subfamilies and identified their potential binding sites.

Results

The sequence analysis of M. smegmatis putative GntRs has revealed that FadR, HutC, MocR and the YtrA-like regulators are encoded by 45, 8, 8 and 1 genes respectively. Further out of 45 FadR-like regulators, 19 were classified into the FadR group and 26 into the VanR group. All these proteins showed similar secondary structural elements specific to their respective subfamilies except MSMEG_3959, which showed additional secondary structural elements. Using the reciprocal BLAST searches, we further identified the orthologs of these regulators in Bacillus subtilis and other mycobacteria. Since the expression of many regulators is auto-regulatory, we have identified potential operator sites for a number of these GntR regulators by analyzing the upstream sequences.

Conclusion

This study helps in extending the annotation of M. smegmatis GntR proteins. It identifies the GntR regulators of M. smegmatis that could serve as a model for studying orthologous regulators from virulent as well as other saprophytic mycobacteria. This study also sheds some light on the nucleotide preferences in the target-motifs of GntRs thus providing important leads for initiating the experimental characterization of these proteins, construction of the gene regulatory network for these regulators and an understanding of the influence of these proteins on the physiology of the mycobacteria.

Glutathione-dependent alcohol dehydrogenase AdhC is required for defense against nitrosative stress in Haemophilus influenzae

In Haemophilus influenzae Rd KW20, we identified a gene, adhC, which encodes a class III alcohol dehydrogenase (AdhC) and has S-nitrosoglutathione reductase activity. adhC exists on an operon with estD, which encodes an esterase. Divergent to the adhC-estD operon is the Haemophilus influenzae nmlR gene (nmlR(HI)), which encodes a MerR family regulator that is homologous to the Neisseria MerR-like regulator (NmlR). Analysis of an nmlR(HI) mutant indicated that expression of the adhC-estD operon is regulated by NmlR(HI) in strain Rd KW20. Chromosomal inactivation of either adhC or nmlR(HI) resulted in sensitivity to S-nitrosoglutathione and decreased S-nitrosoglutathione reductase activity. Examination of the NmlR(HI)-AdhC system in the genome sequences of nontypeable H. influenzae strains R2846, R2866, and 86-028NP identified significant variations. The adhC gene of 86-028NP was predicted to be nonfunctional due to a premature stop codon. Polymorphisms in the operator/promoter region of R2866 resulted in reduced enzyme activity. This correlated with an increased sensitivity to S-nitrosoglutathione. The adhC-nmlR(HI) system was examined in thirty-three clinical isolates (both capsular and nontypeable strains). Nucleic acid sequence data showed that only strain 86-028NP contained a premature stop codon. There were some variations in the DNA sequence of the operator/promoter region which altered the nmlR(HI) promoter. However, the clinical isolates still possessed S-nitrosoglutathione reductase activity and showed at least the equivalent ability to grow in the presence of S-nitrosoglutathione as Rd KW20. These data suggest that the nmlR(HI)-adhC system has a role in the defense against nitrosative stress in Haemophilus influenzae.

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.

The IclR-type transcriptional repressor LtbR regulates the expression of leucine and tryptophan biosynthesis genes in the amino acid producer Corynebacterium glutamicum

The transcriptional regulator Cg1486 of Corynebacterium glutamicum ATCC 13032 is a member of the IclR protein family and belongs to the conserved set of regulatory proteins in corynebacteria. A defined deletion in the cg1486 gene, now designated ltbR (leucine and tryptophan biosynthesis regulator), led to the mutant strain C. glutamicum IB1486. According to whole-genome expression analysis by DNA microarray hybridizations, transcription of the leuB and leuCD genes encoding enzymes of the leucine biosynthesis pathway was enhanced in C. glutamicum IB1486 compared with the wild-type strain. Moreover, the genes of the trpEGDCFBA operon involved in tryptophan biosynthesis of C. glutamicum showed an enhanced expression in the cg1486 mutant strain. Bioinformatics pattern searches in the upstream regions of the differentially expressed genes revealed the common 12-bp motif CA(T/C)ATAGTG(A/G)GA that is located downstream of the -10 region of the mapped promoter sequences. DNA band shift assays with a streptavidin-tagged LtbR protein demonstrated the specific binding of the purified protein to 40-mers containing the 12-bp motif localized in front of leuB, leuC, and trpE, thereby confirming the direct regulatory role of LtbR in the expression of the leucine and tryptophan biosynthesis pathway genes of C. glutamicum. Genes homologous with ltbR were detected upstream of the leuCD genes in almost all sequenced genomes of bacteria belonging to the taxonomic class Actinobacteria. The ltbR-like genes of Corynebacterium diphtheriae, Corynebacterium jeikeium, Mycobacterium bovis, and Bifidobacterium longum were cloned and shown to complement the deregulation of leuB, leuCD, and trpE gene expression in C. glutamicum IB1486.

PhyloScan: identification of transcription factor binding sites using cross-species evidence

Background

When transcription factor binding sites are known for a particular transcription factor, it is possible to construct a motif model that can be used to scan sequences for additional sites. However, few statistically significant sites are revealed when a transcription factor binding site motif model is used to scan a genome-scale database.

Methods

We have developed a scanning algorithm, PhyloScan, which combines evidence from matching sites found in orthologous data from several related species with evidence from multiple sites within an intergenic region, to better detect regulons. The orthologous sequence data may be multiply aligned, unaligned, or a combination of aligned and unaligned. In aligned data, PhyloScan statistically accounts for the phylogenetic dependence of the species contributing data to the alignment and, in unaligned data, the evidence for sites is combined assuming phylogenetic independence of the species. The statistical significance of the gene predictions is calculated directly, without employing training sets.

Results

In a test of our methodology on synthetic data modeled on seven Enterobacteriales, four Vibrionales, and three Pasteurellales species, PhyloScan produces better sensitivity and specificity than MONKEY, an advanced scanning approach that also searches a genome for transcription factor binding sites using phylogenetic information. The application of the algorithm to real sequence data from seven Enterobacteriales species identifies novel Crp and PurR transcription factor binding sites, thus providing several new potential sites for these transcription factors. These sites enable targeted experimental validation and thus further delineation of the Crp and PurR regulons in E. coli.

Conclusion

Better sensitivity and specificity can be achieved through a combination of (1) using mixed alignable and non-alignable sequence data and (2) combining evidence from multiple sites within an intergenic region.

Nitrogen Metabolism and Nitrogen Control in Corynebacteria: Variations of a Common Theme

The published genome sequences of Corynebacterium diphtheriae, Corynebacterium efficiens, Corynebacterium glutamicum and Corynebacterium jeikeium were screened for genes encoding central components of nitrogen source uptake, nitrogen assimilation and nitrogen control systems. Interestingly, the soil-living species C. efficiens and C. glutamicum exhibit a broader spectrum of genes for nitrogen transport and metabolism than the pathogenic species C. diphtheriae and C. jeikeium. The latter are characterized by gene decay and loss of functions like urea metabolism and nitrogen-dependent transcription control. The global regulator of nitrogen regulation AmtR and its DNA-binding motif are conserved in C. diphtheriae, C. efficiens and C. glutamicum, while in C. jeikeium, an AmtR-encoding gene as well as putative AmtR-binding motifs are missing.

MycoperonDB: a database of computationally identified operons and transcriptional units in Mycobacteria

Background

A key post genomics challenge is to identify how genes in an organism come together and perform physiological functions. An important first step in this direction is to identify transcriptional units, operons and regulons in a genome. Here we implement and report a strategy to computationally identify transcriptional units and operons of mycobacteria and construct a database-MycoperonDB.

Description

We have predicted transcriptional units and operons in mycobacteria and organized these predictions in the form of relational database called MycoperonDB. MycoperonDB database at present consists of 18053 genes organized as 8256 predicted operons and transcriptional units from five closely related species of mycobacteria. The database further provides literature links for experimentally characterized operons. All known promoters and related information is collected, analysed and stored. It provides a user friendly interface to allow a web based navigation of transcription units and operons. The web interface provides search tools to locate transcription factor binding DNA motif upstream to various genes. The reliability of operon prediction has been assessed by comparing the predicted operons with a set of known operons.

Conclusion

MycoperonDB is a publicly available structured relational database which has information about mycobacterial genes, transcriptional units and operons. We expect this database to assist molecular biologists/microbiologists in general, to hypothesize functional linkages between operonic genes of mycobacteria, their experimental characterization and validation. The database is freely available from our website .

The transcriptional regulator RovS controls the attachment of Streptococcus agalactiae to human epithelial cells and the expression of virulence genes

Streptococcus agalactiae is part of the normal flora of the human gastrointestinal tract and also the leading cause of bacterial infections in human newborns and immunocompromised adults. The colonization and infection of different regions within the human host require a regulatory network in S. agalactiae that senses environmental stimuli and controls the formation of specific virulence factors. In the present study, we characterized an Rgg-like transcriptional regulator, designated RovS (regulator of virulence in Streptococcus agalactiae). Deletion of the rovS gene in the genome of S. agalactiae resulted in strain 6313 DeltarovS, which exhibited an increased attachment to immobilized fibrinogen and a significant increase in adherence to the eukaryotic lung epithelial cell line A549. Quantification of expression levels of known and putative S. agalactiae virulence genes by real-time PCR revealed that RovS influences the expression of fbsA, gbs0230, sodA, rogB, and the cyl operon. The altered gene expression in mutant 6313 DeltarovS was restored by plasmid-mediated expression of rovS, confirming the RovS deficiency as the cause for the observed changes in virulence gene expression in S. agalactiae. DNA electrophoretic mobility shift assays showed that RovS specifically binds to the promoter regions of fbsA, gbs0230, sodA, and the cyl operon, indicating that RovS directly regulates their expression. Deletion and mutation studies in the promoter region of fbsA, encoding the main fibrinogen receptor in S. agalactiae, identified a RovS DNA motif. Similar motifs were also found in the promoter regions of gbs0230, sodA, and the cyl operon, and alignments allowed us to propose a consensus sequence for the DNA-binding site of RovS.

RegTransBase--a database of regulatory sequences and interactions in a wide range of prokaryotic genomes

RegTransBase is a manually curated database of regulatory interactions in prokaryotes that captures the knowledge in public scientific literature using a controlled vocabulary. Although several databases describing interactions between regulatory proteins and their binding sites are already being maintained, they either focus mostly on the model organisms Escherichia coli and Bacillus subtilis or are entirely computationally derived. RegTransBase describes a large number of regulatory interactions reported in many organisms and contains the following types of experimental data: the activation or repression of transcription by an identified direct regulator, determining the transcriptional regulatory function of a protein (or RNA) directly binding to DNA (RNA), mapping or prediction of a binding site for a regulatory protein and characterization of regulatory mutations. Currently, RegTransBase content is derived from about 3000 relevant articles describing over 7000 experiments in relation to 128 microbes. It contains data on the regulation of about 7500 genes and evidence for 6500 interactions with 650 regulators. RegTransBase also contains manually created position weight matrices (PWM) that can be used to identify candidate regulatory sites in over 60 species. RegTransBase is available at .

Pip, a novel activator of phenazine biosynthesis in Pseudomonas chlororaphis PCL1391

Secondary metabolites are important factors for interactions between bacteria and other organisms. Pseudomonas chlororaphis PCL1391 produces the antifungal secondary metabolite phenazine-1-carboxamide (PCN) that inhibits growth of Fusarium oxysporum f. sp. radius lycopersici the causative agent of tomato foot and root rot. Our previous work unraveled a cascade of genes regulating the PCN biosynthesis operon, phzABCDEFGH. Via a genetic screen, we identify in this study a novel TetR/AcrR regulator, named Pip (phenazine inducing protein), which is essential for PCN biosynthesis. A combination of a phenotypical characterization of a pip mutant, in trans complementation assays of various mutant strains, and electrophoretic mobility shift assays identified Pip as the fifth DNA-binding protein so far involved in regulation of PCN biosynthesis. In this regulatory pathway, Pip is positioned downstream of PsrA (Pseudomonas sigma factor regulator) and the stationary-phase sigma factor RpoS, while it is upstream of the quorum-sensing system PhzI/PhzR. These findings provide further evidence that the path leading to the expression of secondary metabolism gene clusters in Pseudomonas species is highly complex.

iCR: a web tool to identify conserved targets of a regulatory protein across the multiple related prokaryotic species

Gene regulatory circuits are often commonly shared between two closely related organisms. Our web tool iCR (identify Conserved target of a Regulon) makes use of this fact and identify conserved targets of a regulatory protein. iCR is a special refined extension of our previous tool PredictRegulon- that predicts genome wide, the potential binding sites and target operons of a regulatory protein in a single user selected genome. Like PredictRegulon, the iCR accepts known binding sites of a regulatory protein as ungapped multiple sequence alignment and provides the potential binding sites. However important differences are that the user can select more than one genome at a time and the output reports the genes that are common in two or more species. In order to achieve this, iCR makes use of Cluster of Orthologous Group (COG) indices for the genes. This tool analyses the upstream region of all user-selected prokaryote genome and gives the output based on conservation target orthologs. iCR also reports the Functional class codes based on COG classification for the encoded proteins of downstream genes which helps user understand the nature of the co-regulated genes at the result page itself. iCR is freely accessible at .

The SloR/Dlg metalloregulator modulates Streptococcus mutans virulence gene expression

Metal ion availability in the human oral cavity plays a putative role in Streptococcus mutans virulence gene expression and in appropriate formation of the plaque biofilm. In this report, we present evidence that supports such a role for the DtxR-like SloR metalloregulator (called Dlg in our previous publications) in this oral pathogen. Specifically, the results of gel mobility shift assays revealed the sloABC, sloR, comDE, ropA, sod, and spaP promoters as targets of SloR binding. We confirmed differential expression of these genes in a GMS584 SloR-deficient mutant versus the UA159 wild-type progenitor by real-time semiquantitative reverse transcriptase PCR experiments. The results of additional expression studies support a role for SloR in S. mutans control of glucosyltransferases, glucan binding proteins, and genes relevant to antibiotic resistance. Phenotypic analysis of GMS584 revealed that it forms aberrant biofilms on an abiotic surface, is compromised for genetic competence, and demonstrates heightened incorporation of iron and manganese as well as resistance to oxidative stress compared to the wild type. Taken together, these findings support a role for SloR in S. mutans adherence, biofilm formation, genetic competence, metal ion homeostasis, oxidative stress tolerance, and antibiotic gene regulation, all of which contribute to S. mutans-induced disease.

Comparative analysis of iron regulated genes in mycobacteria

Iron dependent regulator, IdeR, regulates the expression of genes in response to intracellular iron levels in M. tuberculosis. Orthologs of IdeR are present in all the sequenced genomes of mycobacteria. We have used a computational approach to identify conserved IdeR regulated genes across the mycobacteria and the genes that are specific to each of the mycobacteria. Novel iron regulated genes that code for a predicted 4-hydroxy benzoyl coA hydrolase (Rv1847) and a protease dependent antibiotic regulatory system (Rv1846c, Rv0185c) are conserved across the mycobacteria. Although Mycobacterium natural-resistance-associated macrophage protein (Mramp) is present in all mycobacteria, it is, as predicted, an iron-regulated gene in only one species, M. avium subsp. paratuberculosis. We also observed an additional iron-regulated exochelin biosynthetic operon, which is present only in non-pathogenic Mycobacterium, M. smegmatis.

Evolution of the NadR regulon in Enterobacteriaceae

The NAD biosynthetic pathway and NAD transformations in E. coli and S. typhi are well characterized. Using comparative genomics methods we describe the NadR regulon in other Enterobacteriaceae, identity new candidate regulon members and demonstrate that even a very simple regulon covering an essential methabolic pathway could be different in closely related genomes.

Computational prediction and experimental verification of novel IdeR binding sites in the upstream sequences of Mycobacterium tuberculosis open reading frames

IdeR (iron-dependent regulator) is a key regulator of virulence factors and iron acquisition systems in Mycobacterium tuberculosis. Despite the wealth of information available on IdeR-regulated genes of M.tuberculosis, there is still an underlying possibility that there are novel genes/pathways that have gone undetected, the identification of which could give new insights into understanding the pathogenesis of M.tuberculosis. We describe an in silico approach employing the positional relative entropy method to identify potential IdeR binding sites in the upstream sequences of all the 3919 ORFs of M.tuberculosis. While many of the predictions made by this approach overlapped with the ones already identified by microarray experiments and binding assays, pointing to the accuracy of our method, a few genes for which there has been no evidence for IdeR regulation were additionally identified. Our results have implications on the iron-dependent regulatory mechanism of M.tuberculosis vis-a-vis the activity of urease operon and novel transcription regulators and transporters.

Differences in LexA regulon structure among Proteobacteria through in vivo assisted comparative genomics

The LexA regulon encompasses an ensemble of genes involved in preserving cell viability under massive DNA damage and is present in most bacterial phyla. Up to date, however, the scope of this network had only been assessed in the Gamma Proteobacteria. Here, we report the structure of the LexA regulon in the Alpha Proteobacteria, using a combined approach that makes use of in vitro and in vivo techniques to assist and validate the comparative genomics in silico methodology. This leads to the first experimentally validated description of the LexA regulon in the Alpha Proteobacteria, and comparison of regulon core structures in both classes suggests that a least common multiple set of genes (recA, ssb, uvrA and ruvCAB) might be a defining property of the Proteobacteria LexA network.

Prediction of DtxR regulon: identification of binding sites and operons controlled by Diphtheria toxin repressor in Corynebacterium diphtheriae

BACKGROUND

The diphtheria toxin repressor, DtxR, of Corynebacterium diphtheriae has been shown to be an iron-activated transcription regulator that controls not only the expression of diphtheria toxin but also of iron uptake genes. This study aims to identify putative binding sites and operons controlled by DtxR to understand the role of DtxR in patho-physiology of Corynebacterium diphtheriae.

RESULT

Positional Shannon relative entropy method was used to build the DtxR-binding site recognition profile and the later was used to identify putative regulatory sites of DtxR within C. diphtheriae genome. In addition, DtxR-regulated operons were also identified taking into account the predicted DtxR regulatory sites and genome annotation. Few of the predicted motifs were experimentally validated by electrophoretic mobility shift assay. The analysis identifies motifs upstream to the novel iron-regulated genes that code for Formamidopyrimidine-DNA glycosylase (FpG), an enzyme involved in DNA-repair and starvation inducible DNA-binding protein (Dps) which is involved in iron storage and oxidative stress defense. In addition, we have found the DtxR motifs upstream to the genes that code for sortase which catalyzes anchoring of host-interacting proteins to the cell wall of pathogenic bacteria and the proteins of secretory system which could be involved in translocation of various iron-regulated virulence factors including diphtheria toxin.

CONCLUSIONS

We have used an in silico approach to identify the putative binding sites and genes controlled by DtxR in Corynebacterium diphtheriae. Our analysis shows that DtxR could provide a molecular link between Fe+2-induced Fenton's reaction and protection of DNA from oxidative damage. DtxR-regulated Dps prevents lethal combination of Fe+2 and H2O2 and also protects DNA by nonspecific DNA-binding. In addition DtxR could play an important role in host interaction and virulence by regulating the levels of sortase, a potential vaccine candidate and proteins of secretory system.