Identifying global regulators in transcriptional regulatory networks in bacteria

The logic layout of the TOL network of Pseudomonas putida pWW0 plasmid stems from a metabolic amplifier motif (MAM) that optimizes biodegradation of m-xylene

Background

The genetic network of the TOL plasmid pWW0 of the soil bacterium Pseudomonas putida mt-2 for catabolism of m-xylene is an archetypal model for environmental biodegradation of aromatic pollutants. Although nearly every metabolic and transcriptional component of this regulatory system is known to an extraordinary molecular detail, the complexity of its architecture is still perplexing. To gain an insight into the inner layout of this network a logic model of the TOL system was implemented, simulated and experimentally validated. This analysis made sense of the specific regulatory topology out on the basis of an unprecedented network motif around which the entire genetic circuit for m-xylene catabolism gravitates.

Results

The most salient feature of the whole TOL regulatory network is the control exerted by two distinct but still intertwined regulators (XylR and XylS) on expression of two separated catabolic operons (upper and lower) for catabolism of m-xylene. Following model reduction, a minimal modular circuit composed by five basic variables appeared to suffice for fully describing the operation of the entire system. In silico simulation of the effect of various perturbations were compared with experimental data in which specific portions of the network were activated with selected inducers: m-xylene, o-xylene, 3-methylbenzylalcohol and 3-methylbenzoate. The results accredited the ability of the model to faithfully describe network dynamics. This analysis revealed that the entire regulatory structure of the TOL system enables the action an unprecedented metabolic amplifier motif (MAM). This motif synchronizes expression of the upper and lower portions of a very long metabolic system when cells face the head pathway substrate, m-xylene.

Conclusion

Logic modeling of the TOL circuit accounted for the intricate regulatory topology of this otherwise simple metabolic device. The found MAM appears to ensure a simultaneous expression of the upper and lower segments of the m-xylene catabolic route that would be difficult to bring about with a standard substrate-responsive single promoter. Furthermore, it is plausible that the MAM helps to avoid biochemical conflicts between competing plasmid-encoded and chromosomally-encoded pathways in this bacterium.

Within-species lateral genetic transfer and the evolution of transcriptional regulation in Escherichia coli and Shigella

BACKGROUND

Changes in transcriptional regulation underlie many of the phenotypic differences observed within and between species of bacteria. Lateral genetic transfer (LGT) can significantly impact the transcription factor (TF) genes which drive these transcriptional changes. Although much emphasis has been placed on LGT of intact genes, the units of transfer and recombination do not necessarily correspond to regions delineated by exact gene boundaries. Here we apply phylogenetic and network-based methods to investigate the relationship between units of lateral transfer and recombination within the Escherichia coli - Shigella clade and the topological properties of genes in the E. coli transcriptional regulatory network (TRN).

RESULTS

We carried out a systematic phylogenetic study of genetic transfer among 5282 sets of putatively orthologous genes from 27 strains belonging to the E. coli - Shigella clade. We then used these results to examine the evolutionary histories of TF genes, as well as the transcriptional regulation of lateral genes. We found evidence of LGT in 2655 (50.3%) gene sets: 678 (12.8%) show evidence of recombination breakpoints within the gene boundaries. Thus, within- and whole- gene lateral transfer is widespread among strains of E. coli and Shigella. We found that unlike global regulators, which have mostly evolved vertically, neighbour regulators (genes which regulate adjacent genes on the chromosome) have frequently been subject to transfer within the E. coli - Shigella clade. At least 56 (62%) of the 90 neighbour regulator gene sets examined show evidence of LGT, 19 (34%) of which have internal recombination breakpoints. Neighbour regulators show no evidence of co-transfer with their nearby target genes. Rather, the frequency of recombination breakpoints, and conflicting evolutionary histories among neighbour regulators and their target genes, suggest that the genomic regions encoding these genes have been constructed through successive layering of LGT events within the clade. We find no difference in the relative complexity of regulation (i.e. the number of regulators) of lateral versus vertical genes.

CONCLUSIONS

Neighbour regulators show higher frequencies of transfer than other types of regulatory genes. This implicates the topological properties of regulatory genes in the TRN, and their physical proximity to targets on the chromosome, as contributing to successful LGT. The prevalence of recombination breakpoints within regulatory and target gene sets indicates that within-gene transfer has had a significant cumulative effect on the evolution of regulatory interactions in E. coli and Shigella.

TFRank: network-based prioritization of regulatory associations underlying transcriptional responses

MOTIVATION

Uncovering mechanisms underlying gene expression control is crucial to understand complex cellular responses. Studies in gene regulation often aim to identify regulatory players involved in a biological process of interest, either transcription factors coregulating a set of target genes or genes eventually controlled by a set of regulators. These are frequently prioritized with respect to a context-specific relevance score. Current approaches rely on relevance measures accounting exclusively for direct transcription factor-target interactions, namely overrepresentation of binding sites or target ratios. Gene regulation has, however, intricate behavior with overlapping, indirect effect that should not be neglected. In addition, the rapid accumulation of regulatory data already enables the prediction of large-scale networks suitable for higher level exploration by methods based on graph theory. A paradigm shift is thus emerging, where isolated and constrained analyses will likely be replaced by whole-network, systemic-aware strategies.

RESULTS

We present TFRank, a graph-based framework to prioritize regulatory players involved in transcriptional responses within the regulatory network of an organism, whereby every regulatory path containing genes of interest is explored and incorporated into the analysis. TFRank selected important regulators of yeast adaptation to stress induced by quinine and acetic acid, which were missed by a direct effect approach. Notably, they reportedly confer resistance toward the chemicals. In a preliminary study in human, TFRank unveiled regulators involved in breast tumor growth and metastasis when applied to genes whose expression signatures correlated with short interval to metastasis.

Quantifying the effect of environment stability on the transcription factor repertoire of marine microbes

Background

DNA-binding transcription factors (TFs) regulate cellular functions in prokaryotes, often in response to environmental stimuli. Thus, the environment exerts constant selective pressure on the TF gene content of microbial communities. Recently a study on marine Synechococcus strains detected differences in their genomic TF content related to environmental adaptation, but so far the effect of environmental parameters on the content of TFs in bacterial communities has not been systematically investigated.

Results

We quantified the effect of environment stability on the transcription factor repertoire of marine pelagic microbes from the Global Ocean Sampling (GOS) metagenome using interpolated physico-chemical parameters and multivariate statistics. Thirty-five percent of the difference in relative TF abundances between samples could be explained by environment stability. Six percent was attributable to spatial distance but none to a combination of both spatial distance and stability. Some individual TFs showed a stronger relationship to environment stability and space than the total TF pool.

Conclusions

Environmental stability appears to have a clearly detectable effect on TF gene content in bacterioplanktonic communities described by the GOS metagenome. Interpolated environmental parameters were shown to compare well to in situ measurements and were essential for quantifying the effect of the environment on the TF content. It is demonstrated that comprehensive and well-structured contextual data will strongly enhance our ability to interpret the functional potential of microbes from metagenomic data.

A composite feed-forward loop I4-FFL involving IHF and Crc stabilizes expression of the XylR regulator of Pseudomonas putida mt-2 from growth phase perturbations

Genetic networks are typically composed of a series of connected motifs that confer specific logic and dynamic properties to the resulting circuits. While some feed forward loop (FFL) variants abound in such networks, others (e.g. the type-4 incoherent FFL or I4-FFL) are virtually absent from the known regulatory devices. We report here that the key node that rules the expression of the m-xylene biodegradation pathway of the soil bacterium Pseudomonas putida mt-2 merges opposite physiological effects of the growth phase by means of a regulatory device based on the rarely found I4-FFL motif. Specifically, the FFL includes the integration host factor (IHF), which both co-activates the master P(u) promoter and represses transcription of its cognate regulatory gene xylR at the onset of the stationary phase. On the other hand, the catabolite repression control (Crc) protein inhibits translation of XylR during exponential growth. By computing these two conflicting regulatory actions within a composite I4-FFL gate, cells shield the expression of XylR from perturbations caused by the growth phase, thereby ensuring a steady supply of the regulator regardless of physiological conditions. This device thus endows xylR expression with a degree of robustness in respect to the growth phase that could hardly be achieved with e.g. a simple constitutive promoter.

Deciphering the regulatory codes in bacterial genomes

Interactions between cis-regulatory elements and trans-acting factors are fundamental for cellular functions such as transcription. With the revolution in microarrays and sequencing technologies, genome-wide binding locations of trans-acting factors are being determined in large numbers. The richness of the genome-scale information has revealed that the nature of the bacterial transcriptome and regulome are considerably more complex than previously expected. In addition, the emerging view of the bacterial transcriptome is revising the concept of the operon organization of the genome. This review describes current advances in the genome-scale analysis of the interaction between cis-regulatory elements and trans-acting factors in microorganisms.

The Regulatory Network of Pseudomonas aeruginosa

Background

Pseudomonas aeruginosa is an important bacterial model due to its metabolic and pathogenic abilities, which allow it to interact and colonize a wide range of hosts, including plants and animals. In this work we compile and analyze the structure and organization of an experimentally supported regulatory network in this bacterium.

Results

The regulatory network consists of 690 genes and 1020 regulatory interactions between their products (12% of total genes: 54% sigma and 16% of transcription factors). This complex interplay makes the third largest regulatory network of those reported in bacteria. The entire network is enriched for activating interactions and, peculiarly, self-activation seems to occur more prominent for transcription factors (TFs), which contrasts with other biological networks where self-repression is dominant. The network contains a giant component of 650 genes organized into 11 hierarchies, encompassing important biological processes, such as, biofilms formation, production of exopolysaccharide alginate and several virulence factors, and of the so-called quorum sensing regulons.

Conclusions

The study of gene regulation in P. aeruginosa is biased towards pathogenesis and virulence processes, all of which are interconnected. The network shows power-law distribution -input degree -, and we identified the top ten global regulators, six two-element cycles, the longest paths have ten steps, six biological modules and the main motifs containing three and four elements. We think this work can provide insights for the design of further studies to cover the many gaps in knowledge of this important bacterial model, and for the design of systems strategies to combat this bacterium.

The regulation of the secondary metabolism of Streptomyces: new links and experimental advances

Streptomycetes and other actinobacteria are renowned as a rich source of natural products of clinical, agricultural and biotechnological value. They are being mined with renewed vigour, supported by genome sequencing efforts, which have revealed a coding capacity for secondary metabolites in vast excess of expectations that were based on the detection of antibiotic activities under standard laboratory conditions. Here we review what is known about the control of production of so-called secondary metabolites in streptomycetes, with an emphasis on examples where details of the underlying regulatory mechanisms are known. Intriguing links between nutritional regulators, primary and secondary metabolism and morphological development are discussed, and new data are included on the carbon control of development and antibiotic production, and on aspects of the regulation of the biosynthesis of microbial hormones. Given the tide of antibiotic resistance emerging in pathogens, this review is peppered with approaches that may expand the screening of streptomycetes for new antibiotics by awakening expression of cryptic antibiotic biosynthetic genes. New technologies are also described that have potential to greatly further our understanding of gene regulation in what is an area fertile for discovery and exploitation

Manipulation of pH shift to enhance the growth and antibiotic activity of Xenorhabdus nematophila

To evaluate the effects of pH control strategy on cell growth and the production of antibiotic (cyclo(2-Me-BABA-Gly)) by Xenorhabdus nematophila and enhance the antibiotic activity. The effects of uncontrolled- (different initial pH) and controlled-pH (different constant pH and pH-shift) operations on cell growth and antibiotic activity of X. nematophila YL00I were examined. Experiments showed that the optimal initial pH for cell growth and antibiotic production of X. nematophila YL001 occurred at 7.0. Under different constant pH, a pH level of 7.5 was found to be optimal for biomass and antibiotic activity at 23.71 g/L and 100.0 U/mL, respectively. Based on the kinetic information relating to the different constant pH effects on the fermentation of X. nematophila YL001, a two-stage pH control strategy in which pH 6.5 was maintained for the first 24 h, and then switched to 7.5 after 24 h, was established to improve biomass production and antibiotic activity. By applying this pH-shift strategy, the maximal antibiotic activity and productivity were significantly improved and reaching 185.0 U/mL and 4.41 U/mL/h, respectively, compared to values obtained from constant pH operation (100.0 U/mL and 1.39 U/mL/h).

Novel two-component systems implied in antibiotic production in Streptomyces coelicolor

The abundance of two-component systems (TCSs) in Streptomyces coelicolor A3(2) genome indicates their importance in the physiology of this soil bacteria. Currently, several TCSs have been related to antibiotic regulation, and the purpose in this study was the characterization of five TCSs, selected by sequence homology with the well-known absA1A2 system, that could also be associated with this important process. Null mutants of the five TCSs were obtained and two mutants (ΔSCO1744/1745 and ΔSCO4596/4597/4598) showed significant differences in both antibiotic production and morphological differentiation, and have been renamed as abr (antibiotic regulator). No detectable changes in antibiotic production were found in the mutants in the systems that include the ORFs SCO3638/3639, SCO3640/3641 and SCO2165/2166 in any of the culture conditions assayed. The system SCO1744/1745 (AbrA1/A2) was involved in negative regulation of antibiotic production, and acted also as a negative regulator of the morphological differentiation. By contrast, the system SCO4596/4597/4598 (AbrC1/C2/C3), composed of two histidine kinases and one response regulator, had positive effects on both morphological development and antibiotic production. Microarray analyses of the ΔabrC1/C2/C3 and wild-type transcriptomes revealed downregulation of actII-ORF4 and cdaR genes, the actinorhodin and calcium-dependent antibiotic pathway-specific regulators respectively. These results demonstrated the involvement of these new two-component systems in antibiotic production and morphological differentiation by different approaches. One is a pleiotropic negative regulator: abrA1/A2. The other one is a positive regulator composed of three elements, two histidine kinases and one response regulator: abrC1/C2/C3.

Implementing an OR-NOT (ORN) logic gate with components of the SOS regulatory network of Escherichia coli

Whether biological or electronic, man-engineered computation is based on logic circuits assembled with binary gates that are interconnected to perform Boolean operations. We report here the rewiring of the SOS system of Escherichia in a fashion that makes the output of both the recA and lexA promoters to faithfully follow the pattern of a binary composite OR-NOT gate (ORN) in which the inputs are DNA damage (e.g. nalidixic acid addition) and IPTG as an exogenous signal. Unlike other non-natural gates whose implementation requires changes in genes and promoters of the genome of the host cells, this ORN was brought about by the sole addition of wild-type bacteria with a plasmid encoding a module for LacI(q)-dependent expression of lexA. Specifically, we demonstrate that the interplay between native, chromosomally-encoded components of the SOS system and the extra parts engineered in such a plasmid made the desired performance to happen without any modification of the core DNA-damage response network. It is thus possible to artificially interface autonomous cell networks with a predetermined logic by means of Boolean gates built with regulatory elements already functioning in the recipient organism.

The PurR regulon in Escherichia coli K-12 MG1655

The PurR transcription factor plays a critical role in transcriptional regulation of purine metabolism in enterobacteria. Here, we elucidate the role of PurR under exogenous adenine stimulation at the genome-scale using high-resolution chromatin immunoprecipitation (ChIP)–chip and gene expression data obtained under in vivo conditions. Analysis of microarray data revealed that adenine stimulation led to changes in transcript level of about 10% of Escherichia coli genes, including the purine biosynthesis pathway. The E. coli strain lacking the purR gene showed that a total of 56 genes are affected by the deletion. From the ChIP–chip analysis, we determined that over 73% of genes directly regulated by PurR were enriched in the biosynthesis, utilization and transport of purine and pyrimidine nucleotides, and 20% of them were functionally unknown. Compared to the functional diversity of the regulon of the other general transcription factors in E. coli, the functions and size of the PurR regulon are limited.

Isolation and characterization of Ehrlichia chaffeensis RNA polymerase and its use in evaluating p28 outer membrane protein gene promoters

Background

Ehrlichia chaffeensis is a tick-transmitted rickettsial pathogen responsible for an important emerging disease, human monocytic ehrlichiosis. To date how E. chaffeensis and many related tick-borne rickettsial pathogens adapt and persist in vertebrate and tick hosts remain largely unknown. In recent studies, we demonstrated significant host-specific differences in protein expression in E. chaffeensis originating from its tick and vertebrate host cells. The adaptive response of the pathogen to different host environments entails switch of gene expression regulated at the level of transcription, possibly by altering RNA polymerase activity.

Results

In an effort to understand the molecular basis of pathogen gene expression differences, we isolated native E. chaffeensis RNA polymerase using a heparin-agarose purification method and developed an in vitro transcription system to map promoter regions of two differentially expressed genes of the p28 outer membrane protein locus, p28-Omp14 and p28-Omp19. We also prepared a recombinant protein of E. chaffeensis σ70 homologue and used it for in vitro promoter analysis studies. The possible role of one or more proteins presents in E. chaffeensis lysates in binding to the promoter segments and on the modulation of in vitro transcription was also assessed.

Conclusions

Our experiments demonstrated that both the native and recombinant proteins are functional and have similar enzyme properties in driving the transcription from E. chaffeensis promoters. This is the first report of the functional characterization of E. chaffeensis RNA polymerase and in vitro mapping of the pathogen promoters using the enzyme. This study marks the beginning to broadly characterize the mechanisms controlling the transcription by Anaplasmataceae pathogens.

Effect of iclR and arcA knockouts on biomass formation and metabolic fluxes in Escherichia coli K12 and its implications on understanding the metabolism of Escherichia coli BL21 (DE3)

Background

Gene expression is regulated through a complex interplay of different transcription factors (TFs) which can enhance or inhibit gene transcription. ArcA is a global regulator that regulates genes involved in different metabolic pathways, while IclR as a local regulator, controls the transcription of the glyoxylate pathway genes of the aceBAK operon. This study investigates the physiological and metabolic consequences of arcA and iclR deletions on E. coli K12 MG1655 under glucose abundant and limiting conditions and compares the results with the metabolic characteristics of E. coli BL21 (DE3).

Results

The deletion of arcA and iclR results in an increase in the biomass yield both under glucose abundant and limiting conditions, approaching the maximum theoretical yield of 0.65 c-mole/c-mole glucose under glucose abundant conditions. This can be explained by the lower flux through several CO₂ producing pathways in the E. coli K12 ΔarcAΔiclR double knockout strain. Due to iclR gene deletion, the glyoxylate pathway is activated resulting in a redirection of 30% of the isocitrate molecules directly to succinate and malate without CO₂ production. Furthermore, a higher flux at the entrance of the TCA was noticed due to arcA gene deletion, resulting in a reduced production of acetate and less carbon loss. Under glucose limiting conditions the flux through the glyoxylate pathway is further increased in the ΔiclR knockout strain, but this effect was not observed in the double knockout strain. Also a striking correlation between the glyoxylate flux data and the isocitrate lyase activity was observed for almost all strains and under both growth conditions, illustrating the transcriptional control of this pathway. Finally, similar central metabolic fluxes were observed in E. coli K12 ΔarcA ΔiclR compared to the industrially relevant E. coli BL21 (DE3), especially with respect to the pentose pathway, the glyoxylate pathway, and the TCA fluxes. In addition, a comparison of the genome sequences of the two strains showed that BL21 possesses two mutations in the promoter region of iclR and rare codons are present in arcA implying a lower tRNA acceptance. Both phenomena presumably result in a reduced ArcA and IclR synthesis in BL21, which contributes to the similar physiology as observed in E. coli K12 ΔarcAΔiclR.

Conclusions

The deletion of arcA results in a decrease of repression on transcription of TCA cycle genes under glucose abundant conditions, without significantly affecting the glyoxylate pathway activity. IclR clearly represses transcription of glyoxylate pathway genes under glucose abundance, a condition in which Crp activation is absent. Under glucose limitation, Crp is responsible for the high glyoxylate flux, but IclR still represses transcription. Finally, in E. coli BL21 (DE3), ArcA and IclR are poorly expressed, explaining the similar fluxes observed compared to the ΔarcAΔiclR strain.

Stress response regulators identified through genome-wide transcriptome analysis of the (p)ppGpp-dependent response in Rhizobium etli

BACKGROUND

The alarmone (p)ppGpp mediates a global reprogramming of gene expression upon nutrient limitation and other stresses to cope with these unfavorable conditions. Synthesis of (p)ppGpp is, in most bacteria, controlled by RelA/SpoT (Rsh) proteins. The role of (p)ppGpp has been characterized primarily in Escherichia coli and several Gram-positive bacteria. Here, we report the first in-depth analysis of the (p)ppGpp-regulon in an α-proteobacterium using a high-resolution tiling array to better understand the pleiotropic stress phenotype of a relA/rsh mutant.

RESULTS

We compared gene expression of the Rhizobium etli wild type and rsh (previously rel) mutant during exponential and stationary phase, identifying numerous (p)ppGpp targets, including small non-coding RNAs. The majority of the 834 (p)ppGpp-dependent genes were detected during stationary phase. Unexpectedly, 223 genes were expressed (p)ppGpp-dependently during early exponential phase, indicating the hitherto unrecognized importance of (p)ppGpp during active growth. Furthermore, we identified two (p)ppGpp-dependent key regulators for survival during heat and oxidative stress and one regulator putatively involved in metabolic adaptation, namely extracytoplasmic function sigma factor EcfG2/PF00052, transcription factor CH00371, and serine protein kinase PrkA.

CONCLUSIONS

The regulatory role of (p)ppGpp in R. etli stress adaptation is far-reaching in redirecting gene expression during all growth phases. Genome-wide transcriptome analysis of a strain deficient in a global regulator, and exhibiting a pleiotropic phenotype, enables the identification of more specific regulators that control genes associated with a subset of stress phenotypes. This work is an important step toward a full understanding of the regulatory network underlying stress responses in α-proteobacteria.

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Physiological roles of ArcA, Crp, and EtrA and their interactive control on aerobic and anaerobic respiration in Shewanella oneidensis

In the genome of Shewanella oneidensis, genes encoding the global regulators ArcA, Crp, and EtrA have been identified. All these proteins deviate from their counterparts in E. coli significantly in terms of functionality and regulon. It is worth investigating the involvement and relationship of these global regulators in aerobic and anaerobic respiration in S. oneidensis. In this study, the impact of the transcriptional factors ArcA, Crp, and EtrA on aerobic and anaerobic respiration in S. oneidensis were assessed. While all these proteins appeared to be functional in vivo, the importance of individual proteins in these two major biological processes differed. The ArcA transcriptional factor was critical in aerobic respiration while the Crp protein was indispensible in anaerobic respiration. Using a newly developed reporter system, it was found that expression of arcA and etrA was not influenced by growth conditions but transcription of crp was induced by removal of oxygen. An analysis of the impact of each protein on transcription of the others revealed that Crp expression was independent of the other factors whereas ArcA repressed both etrA and its own transcription while EtrA also repressed arcA transcription. Transcriptional levels of arcA in the wild type, crp, and etrA strains under either aerobic or anaerobic conditions were further validated by quantitative immunoblotting with a polyclonal antibody against ArcA. This extensive survey demonstrated that all these three global regulators are functional in S. oneidensis. In addition, the reporter system constructed in this study will facilitate in vivo transcriptional analysis of targeted promoters.

Unexpected coregulator range for the global regulator Lrp of Escherichia coli and Proteus mirabilis

The Lrp/AsnC family of transcription factors links gene regulation to metabolism in bacteria and archaea. Members of this family, collectively, respond to a wide range of amino acids as coregulators. In Escherichia coli, Lrp regulates over 200 genes directly and is well known to respond to leucine and, to a somewhat lesser extent, alanine. We focused on Lrp from Proteus mirabilis and E. coli, orthologs with 98% identity overall and identical helix-turn-helix motifs, for which a previous study nevertheless found functional differences. Sequence differences between these orthologs, within and adjacent to the amino acid-responsive RAM domain, led us to test for differential sensitivity to coregulatory amino acids. In the course of this investigation, we found, via in vivo reporter fusion assays and in vitro electrophoretic mobility shift experiments, that E. coli Lrp itself responded to a broader range of amino acids than was previously appreciated. In particular, for both the E. coli and P. mirabilis orthologs, Lrp responsiveness to methionine was similar in magnitude to that to leucine. Both Lrp orthologs are also fairly sensitive to Ile, His, and Thr. These observations suggest that Lrp ties gene expression in the Enterobacteriaceae rather extensively to physiological status, as reflected in amino acid pools. These findings also have substantial implications for attempts to model regulatory architecture from transcriptome measurements or to infer such architecture from genome sequences, and they suggest that even well-studied regulators deserve ongoing exploration.

A new framework for identifying cis-regulatory motifs in prokaryotes

We present a new algorithm, BOBRO, for prediction of cis-regulatory motifs in a given set of promoter sequences. The algorithm substantially improves the prediction accuracy and extends the scope of applicability of the existing programs based on two key new ideas: (i) we developed a highly effective method for reliably assessing the possibility for each position in a given promoter to be the (approximate) start of a conserved sequence motif; and (ii) we developed a highly reliable way for recognition of actual motifs from the accidental ones based on the concept of ‘motif closure’. These two key ideas are embedded in a classical framework for motif finding through finding cliques in a graph but have made this framework substantially more sensitive as well as more selective in motif finding in a very noisy background. A comparative analysis shows that the performance coefficient was improved from 29% to 41% by our program compared to the best among other six state-of-the-art prediction tools on a large-scale data sets of promoters from one genome, and also consistently improved by substantial margins on another kind of large-scale data sets of orthologous promoters across multiple genomes. The power of BOBRO in dealing with noisy data was further demonstrated through identification of the motifs of the global transcriptional regulators by running it over 2390 promoter sequences of Escherichia coli K12.

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

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

Reconstruction of the regulatory network of Lactobacillus plantarum WCFS1 on basis of correlated gene expression and conserved regulatory motifs

Gene regulatory networks can be reconstructed by combining transcriptome data from many different experiments to elucidate relations between the activity of certain transcription factors and the genes they control. To obtain insight in the regulatory network of Lactobacillus plantarum, microarray transcriptome data from more than 70 different experimental conditions were combined and the expression profiles of the transcriptional units (TUs) were compared. The TUs that displayed correlated expression were used to identify putative cis‐regulatory elements by searching the upstream regions of the TUs for conserved motifs. Predicted motifs were extended and refined by searching for motifs in the upstream regions of additional TUs with correlated expression. In this way, cis‐acting elements were identified for 41 regulons consisting of at least four TUs (correlation > 0.7). This set of regulons included the known regulons of CtsR and LexA, but also several novel ones encompassing genes with coherent biological functions. Visualization of the regulons and their connections revealed a highly interconnected regulatory network. This network contains several subnetworks that encompass genes of correlated biological function, such as sugar and energy metabolism, nitrogen metabolism and stress response.

The methanogen-specific transcription factor MsvR regulates the fpaA-rlp-rub oxidative stress operon adjacent to msvR in Methanothermobacter thermautotrophicus

Methanogens represent some of the most oxygen-sensitive organisms in laboratory culture. Recent studies indicate that they have developed mechanisms to deal with brief oxygen exposure. MsvR is a transcriptional regulator that has a domain architecture unique to a select group of methanogens. Here, runoff in vitro transcription assays were used to demonstrate that MsvR regulates transcription of the divergently transcribed fpaA-rlp-rub operon in Methanothermobacter thermautotrophicus in addition to transcription from its own promoter. The protein products of the fpaA-rlp-rub operon have previously been implicated in oxidative stress responses in M. thermautotrophicus. Additionally, electrophoretic mobility shift assays (EMSAs) and DNase I footprinting were used to confirm a binding site inferred by bioinformatic analysis. Sequence mutations within these binding sites did not significantly alter EMSA shifting patterns on longer templates but did on shorter 50-bp fragments encompassing only the region containing the binding sites. Footprinting confirmed that the regions protected for the longer mutant templates are at different positions within the intergenic region compared to those seen in the intact intergenic region. Oxidized and reduced preparations of MsvR demonstrated different EMSA binding patterns and regions of protection on the intergenic sequence, suggesting that MsvR may play a role in detecting the redox state of the cell.

Regulatory exaptation of the catabolite repression protein (Crp)-cAMP system in Pseudomonas putida

The genome of the soil bacterium Pseudomonas putida KT2440 encodes singular orthologues of genes crp (encoding the catabolite repression protein, Crp) and cyaA (adenylate cyclase) of Escherichia coli. The levels of cAMP formed by P. putida cells were below detection with a Dictyostelium biosensor in vivo. The cyaA(P. putida) gene was transcribed in vivo but failed to complement the lack of maltose consumption of a cyaA mutant of E. coli, thereby indicating that cyaA(P. putida) was poorly translated or rendered non-functional in the heterologous host. Yet, generation of cAMP by CyaA(P. putida) could be verified by expressing the cyaA(P. putida) gene in a hypersensitive E. coli strain. On the other hand, the crp(P. putida) gene restored the metabolic capacities of an equivalent crp mutant of E. coli, but not in a double crp/cyaA strain, suggesting that the ability to regulate such functions required cAMP. In order to clarify the breadth of the Crp/cAMP system in P. putida, crp and cyaA mutants were generated and passed through a battery of phenotypic tests for recognition of gross metabolic properties and stress-endurance abilities. These assays revealed that the loss of each gene led in most (but not all) cases to the same phenotypic behaviour, indicating a concerted functionality. Unexpectedly, none of the mutations affected the panel of carbon compounds that can be used by P. putida as growth substrates, the mutants being impaired only in the use of various dipeptides as N sources. Furthermore, the lack of crp or cyaA had little influence on the gross growth fingerprinting of the cells. The poor physiological profile of the Crp-cAMP system of P. putida when compared with E. coli exposes a case of regulatory exaptation, i.e. the process through which a property evolved for a particular function is co-opted for a new use.

Novel sequence-based method for identifying transcription factor binding sites in prokaryotic genomes

MOTIVATION

Computational techniques for microbial genomic sequence analysis are becoming increasingly important. With next-generation sequencing technology and the human microbiome project underway, current sequencing capacity is significantly greater than the speed at which organisms of interest can be studied experimentally. Most related computational work has been focused on sequence assembly, gene annotation and metabolic network reconstruction. We have developed a method that will primarily use available sequence data in order to determine prokaryotic transcription factor (TF) binding specificities.

RESULTS

Specificity determining residues (critical residues) were identified from crystal structures of DNA-protein complexes and TFs with the same critical residues were grouped into specificity classes. The putative binding regions for each class were defined as the set of promoters for each TF itself (autoregulatory) and the immediately upstream and downstream operons. MEME was used to find putative motifs within each separate class. Tests on the LacI and TetR TF families, using RegulonDB annotated sites, showed the sensitivity of prediction 86% and 80%, respectively.

AVAILABILITY

http://ural.wustl.edu/∼gsahota/HTHmotif/

Identifying modules of coexpressed transcript units and their organization of Saccharopolyspora erythraea from time series gene expression profiles

Background

The Saccharopolyspora erythraea genome sequence was released in 2007. In order to look at the gene regulations at whole transcriptome level, an expression microarray was specifically designed on the S. erythraea strain NRRL 2338 genome sequence. Based on these data, we set out to investigate the potential transcriptional regulatory networks and their organization.

Methodology/Principal Findings

In view of the hierarchical structure of bacterial transcriptional regulation, we constructed a hierarchical coexpression network at whole transcriptome level. A total of 27 modules were identified from 1255 differentially expressed transcript units (TUs) across time course, which were further classified in to four groups. Functional enrichment analysis indicated the biological significance of our hierarchical network. It was indicated that primary metabolism is activated in the first rapid growth phase (phase A), and secondary metabolism is induced when the growth is slowed down (phase B). Among the 27 modules, two are highly correlated to erythromycin production. One contains all genes in the erythromycin-biosynthetic (ery) gene cluster and the other seems to be associated with erythromycin production by sharing common intermediate metabolites. Non-concomitant correlation between production and expression regulation was observed. Especially, by calculating the partial correlation coefficients and building the network based on Gaussian graphical model, intrinsic associations between modules were found, and the association between those two erythromycin production-correlated modules was included as expected.

Conclusions

This work created a hierarchical model clustering transcriptome data into coordinated modules, and modules into groups across the time course, giving insight into the concerted transcriptional regulations especially the regulation corresponding to erythromycin production of S. erythraea. This strategy may be extendable to studies on other prokaryotic microorganisms.

De novo evolution of complex, global and hierarchical gene regulatory mechanisms

Gene regulatory networks exhibit complex, hierarchical features such as global regulation and network motifs. There is much debate about whether the evolutionary origins of such features are the results of adaptation, or the by-products of non-adaptive processes of DNA replication. The lack of availability of gene regulatory networks of ancestor species on evolutionary timescales makes this a particularly difficult problem to resolve. Digital organisms, however, can be used to provide a complete evolutionary record of lineages. We use a biologically realistic evolutionary model that includes gene expression, regulation, metabolism and biosynthesis, to investigate the evolution of complex function in gene regulatory networks. We discover that: (i) network architecture and complexity evolve in response to environmental complexity, (ii) global gene regulation is selected for in complex environments, (iii) complex, inter-connected, hierarchical structures evolve in stages, with energy regulation preceding stress responses, and stress responses preceding growth rate adaptations and (iv) robustness of evolved models to mutations depends on hierarchical level: energy regulation and stress responses tend not to be robust to mutations, whereas growth rate adaptations are more robust and non-lethal when mutated. These results highlight the adaptive and incremental evolution of complex biological networks, and the value and potential of studying realistic in silico evolutionary systems as a way of understanding living systems.

Defining the plasticity of transcription factor binding sites by Deconstructing DNA consensus sequences: the PhoP-binding sites among gamma/enterobacteria

Transcriptional regulators recognize specific DNA sequences. Because these sequences are embedded in the background of genomic DNA, it is hard to identify the key cis-regulatory elements that determine disparate patterns of gene expression. The detection of the intra- and inter-species differences among these sequences is crucial for understanding the molecular basis of both differential gene expression and evolution. Here, we address this problem by investigating the target promoters controlled by the DNA-binding PhoP protein, which governs virulence and Mg(2+) homeostasis in several bacterial species. PhoP is particularly interesting; it is highly conserved in different gamma/enterobacteria, regulating not only ancestral genes but also governing the expression of dozens of horizontally acquired genes that differ from species to species. Our approach consists of decomposing the DNA binding site sequences for a given regulator into families of motifs (i.e., termed submotifs) using a machine learning method inspired by the "Divide & Conquer" strategy. By partitioning a motif into sub-patterns, computational advantages for classification were produced, resulting in the discovery of new members of a regulon, and alleviating the problem of distinguishing functional sites in chromatin immunoprecipitation and DNA microarray genome-wide analysis. Moreover, we found that certain partitions were useful in revealing biological properties of binding site sequences, including modular gains and losses of PhoP binding sites through evolutionary turnover events, as well as conservation in distant species. The high conservation of PhoP submotifs within gamma/enterobacteria, as well as the regulatory protein that recognizes them, suggests that the major cause of divergence between related species is not due to the binding sites, as was previously suggested for other regulators. Instead, the divergence may be attributed to the fast evolution of orthologous target genes and/or the promoter architectures resulting from the interaction of those binding sites with the RNA polymerase.

Exploiting a global regulator for small molecule discovery in Photorhabdus luminescens

Bacterially produced small molecules demonstrate a remarkable range of structural and functional diversity and include some of our most useful biological probes and therapeutic agents. Annotations of bacterial genomes reveal a large gap between the number of known small molecules and the number of biosynthetic genes/loci that could produce such small molecules, a gap that most likely originates from tight regulatory control by the producing organism. This study coupled a global transcriptional regulator, HexA, to secondary metabolite production in Photorhabdus luminescens, a member of the Gammaproteobacteria that participates in a complex symbiosis with nematode worms and insect larvae. HexA is a LysR-type transcriptional repressor, and knocking it out to create a P. luminescens DeltahexA mutant led to dramatic upregulation of biosynthesized small molecules. Use of this mutant expanded a family of stilbene-derived small molecules, which were known to play important roles in the symbiosis, from three members to at least nine members.

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

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

The regulatory action of the myxobacterial CarD/CarG complex: a bacterial enhanceosome?

A global regulatory complex made up of two unconventional transcriptional factors, CarD and CarG, is implicated in the control of various processes in Myxococcus xanthus, a Gram-negative bacterium that serves as a prokaryotic model system for multicellular development and the response to blue light. CarD has a unique two-domain architecture composed of: (1) a C-terminal DNA-binding domain that resembles eukaryotic high mobility group A (HMGA) proteins, which are relatively abundant, nonhistone components of chromatin that remodel DNA and prime it for the assembly of multiprotein-DNA complexes essential for various DNA transactions, and (2) an N-terminal domain involved in interactions with CarG and RNA polymerase, which is also the founding member of the large CarD_TRCF family of bacterial proteins. CarG, which does not bind DNA directly, has a zinc-binding motif of the type found in the archaemetzincin class of metalloproteases that, in CarG, appears to play a purely structural role. This review aims to provide an overview of the known molecular details and insights emerging from the study of the singular CarD-CarG prokaryotic regulatory complex and its parallels with enhanceosomes, the higher order, nucleoprotein transcription complexes in eukaryotes.

Inclusion of neighboring base interdependencies substantially improves genome-wide prokaryotic transcription factor binding site prediction

Prediction of transcription factor binding sites is an important challenge in genome analysis. The advent of next generation genome sequencing technologies makes the development of effective computational approaches particularly imperative. We have developed a novel training-based methodology intended for prokaryotic transcription factor binding site prediction. Our methodology extends existing models by taking into account base interdependencies between neighbouring positions using conditional probabilities and includes genomic background weighting. This has been tested against other existing and novel methodologies including position-specific weight matrices, first-order Hidden Markov Models and joint probability models. We have also tested the use of gapped and ungapped alignments and the inclusion or exclusion of background weighting. We show that our best method enhances binding site prediction for all of the 22 Escherichia coli transcription factors with at least 20 known binding sites, with many showing substantial improvements. We highlight the advantage of using block alignments of binding sites over gapped alignments to capture neighbouring position interdependencies. We also show that combining these methods with ChIP-on-chip data has the potential to further improve binding site prediction. Finally we have developed the ungapped likelihood under positional background platform: a user friendly website that gives access to the prediction method devised in this work.

Many chromosomal genes modulate MarA-mediated multidrug resistance in Escherichia coli

Multidrug resistance (MDR) in clinical isolates of Escherichia coli can be associated with overexpression of marA, a transcription factor that upregulates multidrug efflux and downregulates membrane permeability. Using random transposome mutagenesis, we found that many chromosomal genes and environmental stimuli affected MarA-mediated antibiotic resistance. Seven genes affected resistance mediated by MarA in an antibiotic-specific way; these were mostly genes encoding unrelated enzymes, transporters, and unknown proteins. Other genes affected MarA-mediated resistance to all antibiotics tested. These genes were acrA, acrB, and tolC (which encode the major MarA-regulated multidrug efflux pump AcrAB-TolC), crp, cyaA, hns, and pcnB (four genes involved in global regulation of gene expression), and the unknown gene damX. The last five genes affected MarA-mediated MDR by altering marA expression or MarA function specifically on acrA. These findings demonstrate that MarA-mediated MDR is regulated at multiple levels by different genes and stimuli, which makes it both complex and fine-tuned and interconnects it with global cell regulation and metabolism. Such a regulation could contribute to the adaptation and spread of MDR strains and may be targeted to treat antibiotic-resistant E. coli and related pathogens.

Stochasticity versus determinism: consequences for realistic gene regulatory network modelling and evolution

Gene regulation is one important mechanism in producing observed phenotypes and heterogeneity. Consequently, the study of gene regulatory network (GRN) architecture, function and evolution now forms a major part of modern biology. However, it is impossible to experimentally observe the evolution of GRNs on the timescales on which living species evolve. In silico evolution provides an approach to studying the long-term evolution of GRNs, but many models have either considered network architecture from non-adaptive evolution, or evolution to non-biological objectives. Here, we address a number of important modelling and biological questions about the evolution of GRNs to the realistic goal of biomass production. Can different commonly used simulation paradigms, in particular deterministic and stochastic Boolean networks, with and without basal gene expression, be used to compare adaptive with non-adaptive evolution of GRNs? Are these paradigms together with this goal sufficient to generate a range of solutions? Will the interaction between a biological goal and evolutionary dynamics produce trade-offs between growth and mutational robustness? We show that stochastic basal gene expression forces shrinkage of genomes due to energetic constraints and is a prerequisite for some solutions. In systems that are able to evolve rates of basal expression, two optima, one with and one without basal expression, are observed. Simulation paradigms without basal expression generate bloated networks with non-functional elements. Further, a range of functional solutions was observed under identical conditions only in stochastic networks. Moreover, there are trade-offs between efficiency and yield, indicating an inherent intertwining of fitness and evolutionary dynamics.

An unusual pattern of CytR and CRP binding energetics at Escherichia coli cddP suggests a unique blend of class I and class II mediated activation

Two transcription factors, CRP and CytR, mediate positive and negative control of nine cistrons involved in nucleoside catabolism and recycling in Escherichia coli. The ability of multiple transcription factors to combine in different ways to confer differential gene regulation is of significant interest in both prokaryotic and eukaryotic gene regulation. Analysis of cooperative interactions between CytR and CRP at the deoP2 and udpP promoters has implicated the importance of promoter architecture in controlling repression and induction. These studies have also identified competition between CytR and CRP as an additional contributor to differential regulation. The pattern and energetics of CytR and CRP interactions at the cdd promoter, the most strongly activated of the CytR-regulated promoters, have been delineated using DNase I footprinting. Surprisingly, CRP has greater affinity for the promoter proximal site at cddP, CRP1, than for the distal site, CRP2, in contrast to promoters studied previously. This difference is a major contributor to unusually high CRP-mediated activation of cddP. Additionally, while cytidine binding to CytR nearly eliminates the pairwise interactions between CytR and CRP bound at CRP1, it has little effect on pairwise cooperativity between CytR and CRP bound at CRP2 or as a consequence on the overall cooperativity of the three-protein complex in which CRP is bound to both sites. The effect of cytidine binding on cooperativity differs between the three promoters studied thus far. We propose that the different patterns of interaction reflect the spacing between CytR half-sites and the location of the CytR operator in relation to the two CRP sites.

Model of transcriptional activation by MarA in Escherichia coli

The AraC family transcription factor MarA activates ∼40 genes (the marA/soxS/rob regulon) of the Escherichia coli chromosome resulting in different levels of resistance to a wide array of antibiotics and to superoxides. Activation of marA/soxS/rob regulon promoters occurs in a well-defined order with respect to the level of MarA; however, the order of activation does not parallel the strength of MarA binding to promoter sequences. To understand this lack of correspondence, we developed a computational model of transcriptional activation in which a transcription factor either increases or decreases RNA polymerase binding, and either accelerates or retards post-binding events associated with transcription initiation. We used the model to analyze data characterizing MarA regulation of promoter activity. The model clearly explains the lack of correspondence between the order of activation and the MarA-DNA affinity and indicates that the order of activation can only be predicted using information about the strength of the full MarA-polymerase-DNA interaction. The analysis further suggests that MarA can activate without increasing polymerase binding and that activation can even involve a decrease in polymerase binding, which is opposite to the textbook model of activation by recruitment. These findings are consistent with published chromatin immunoprecipitation assays of interactions between polymerase and the E. coli chromosome. We find that activation involving decreased polymerase binding yields lower latency in gene regulation and therefore might confer a competitive advantage to cells. Our model yields insights into requirements for predicting the order of activation of a regulon and enables us to suggest that activation might involve a decrease in polymerase binding which we expect to be an important theme of gene regulation in E. coli and beyond.

When environmental conditions change, cell survival can depend on sudden production of proteins that are normally in low demand. Protein production is controlled by transcription factors which bind to DNA near genes and either increase or decrease RNA production. Many puzzles remain concerning the ways transcription factors do this. Recently we collected data relating the intracellular level of a single transcription factor, MarA, to the increase in expression of several genes related to antibiotic and superoxide resistance in Escherichia coli. These data indicated that target genes are turned on in a well-defined order with respect to the level of MarA, enabling cells to mount a response that is commensurate to the level of threat detected in the environment. Here we develop a computational model to yield insight into how MarA turns on its target genes. The modeling suggests that MarA can increase the frequency with which a transcript is made while decreasing the overall presence of the transcription machinery at the start of a gene. This mechanism is opposite to the textbook model of transcriptional activation; nevertheless it enables cells to respond quickly to environmental challenges and is likely of general importance for gene regulation in E. coli and beyond.

New insights into the regulatory networks of paralogous genes in bacteria

Extensive genomic studies on gene duplication in model organisms such as Escherichia coli and Saccharomyces cerevisiae have recently been undertaken. In these models, it is commonly considered that a duplication event may include a transcription factor (TF), a target gene, or both. Following a gene duplication episode, varying scenarios have been postulated to describe the evolution of the regulatory network. However, in most of these, the TFs have emerged as the most important and in some cases the only factor shaping the regulatory network as the organism responds to a natural selection process, in order to fulfil its metabolic needs. Recent findings concerning the regulatory role played by elements other than TFs have indicated the need to reassess these early models. Thus, we performed an exhaustive review of paralogous gene regulation in E. coli and Bacillus subtilis based on published information, available in the NCBI PubMed database and in well-established regulatory databases. Our survey reinforces the notion that despite TFs being the most prominent components shaping the regulatory networks, other elements are also important. These include small RNAs, riboswitches, RNA-binding proteins, sigma factors, protein-protein interactions and DNA supercoiling, which modulate the expression of genes involved in particular metabolic processes or induce a more complex response in terms of the regulatory networks of paralogous genes in an integrated interplay with TFs.

Inferring the transcriptional network of Bacillus subtilis

The adaptation of bacteria to the vigorous environmental changes they undergo is crucial to their survival. They achieve this adaptation partly via intricate regulation of the transcription of their genes. In this study, we infer the transcriptional network of the Gram-positive model organism, Bacillus subtilis. We use a data integration workflow, exploiting both motif and expression data, towards the generation of condition-dependent transcriptional modules. In building the motif data, we rely on both known and predicted information. Known motifs were derived from DBTBS, while predicted motifs were generated by a de novo motif detection method that utilizes comparative genomics. The expression data consists of a compendium of microarrays across different platforms. Our results indicate that a considerable part of the B. subtilis network is yet undiscovered; we could predict 417 new regulatory interactions for known regulators and 453 interactions for yet uncharacterized regulators. The regulators in our network showed a preference for regulating modules in certain environmental conditions. Also, substantial condition-dependent intra-operonic regulation seems to take place. Global regulators seem to require functional flexibility to attain their roles by acting as both activators and repressors.

Metabolic flux analysis of Escherichia coli creB and arcA mutants reveals shared control of carbon catabolism under microaerobic growth conditions

Escherichia coli has several elaborate sensing mechanisms for response to availability of oxygen and other electron acceptors, as well as the carbon source in the surrounding environment. Among them, the CreBC and ArcAB two-component signal transduction systems are responsible for regulation of carbon source utilization and redox control in response to oxygen availability, respectively. We assessed the role of CreBC and ArcAB in regulating the central carbon metabolism of E. coli under microaerobic conditions by means of (13)C-labeling experiments in chemostat cultures of a wild-type strain, DeltacreB and DeltaarcA single mutants, and a DeltacreB DeltaarcA double mutant. Continuous cultures were conducted at D = 0.1 h(-1) under carbon-limited conditions with restricted oxygen supply. Although all experimental strains metabolized glucose mainly through the Embden-Meyerhof-Parnas pathway, mutant strains had significantly lower fluxes in both the oxidative and the nonoxidative pentose phosphate pathways. Significant differences were also found at the pyruvate branching point. Both pyruvate-formate lyase and the pyruvate dehydrogenase complex contributed to acetyl-coenzyme A synthesis from pyruvate, and their activity seemed to be modulated by both ArcAB and CreBC. Strains carrying the creB deletion showed a higher biomass yield on glucose compared to the wild-type strain and its DeltaarcA derivative, which also correlated with higher fluxes from building blocks to biomass. Glyoxylate shunt and lactate dehydrogenase were active mainly in the DeltaarcA strain. Finally, it was observed that the tricarboxylic acid cycle reactions operated in a rather cyclic fashion under our experimental conditions, with reduced activity in the mutant strains.

Plasticity of transcriptional machinery in bacteria is increased by the repertoire of regulatory families

Escherichia coli K12 and Bacillus subtilis 168 are two of the best characterized bacterial organisms with a long history in molecular biology for understanding various mechanisms in prokaryotic species. However, at the level of transcriptional regulation little is known on a comparative scale. Here we address the question of the degree to which transcription factors (TFs) and their evolutionary families are shared between them. We found that 59 proteins and 28 families are shared between these two bacteria, whereas different subsets were lineage specific. We demonstrate that majority of the common families expand in a lineage-specific manner. More specifically, we found that AraC, ColD, Ebp, LuxR and LysR families are over-represented in E. coli, while ArsR, AsnC, MarR, MerR and TetR families have significantly expanded in B. subtilis. We introduce the notion of regulatory superfamilies based on an empirical number of functional categories regulated by them and show that these families are essentially different in the two bacteria. We further show that global regulators seem to be constrained to smaller regulatory families and generally originate from lineage-specific families. We find that although TF families may be conserved across genomes their functional roles might evolve in a lineage-specific manner and need not be conserved, indicating convergence to be an important phenomenon involved in the functional evolution of TFs of the same family. Although topologically the networks of transcriptional interactions among TF families are similar in both the genomes, we found that the players are different, suggesting different evolutionary origins for the transcriptional regulatory machinery in both bacteria. This study provides evidence from complete repertoires that not only novel families originate in different lineages but conserved TF families expand/contrast in a lineage-specific manner, and suggests that part of the global regulatory mechanisms might originate independently in different lineages.

Global regulators orchestrate group II intron retromobility

Group II introns are hypothesized to share common ancestry with both nuclear spliceosomal introns and retrotransposons, which collectively occupy the majority of genome space in higher eukaryotes. These phylogenetically diverse introns are mobile retroelements that move through an RNA intermediate. Disruption of Escherichia coli genes encoding enzymes that catalyze synthesis of global regulators cAMP and ppGpp inhibits group II intron retromobility. These small molecules program genetic transitions between nutrient excess and starvation. Accordingly, we demonstrated that glucose depletion of wild-type cells and cAMP supplementation of mutants stimulated retromobility. Likewise, amino acid starvation, which induces the alarmone ppGpp, activated retromobility. In both cases, retrotransposition to ectopic sites was favored over retrohoming. Interestingly, these stimulatory effects are mediated at the level of the DNA target, rather than of expression of the retroelement. Thereby, during metabolic stress, cAMP and ppGpp control group II intron movement in concert with the cell's global genetic circuitry, stimulating genetic diversity.

Cell envelope perturbation induces oxidative stress and changes in iron homeostasis in Vibrio cholerae

The Vibrio cholerae type II secretion (T2S) machinery is a multiprotein complex that spans the cell envelope. When the T2S system is inactivated, cholera toxin and other exoproteins accumulate in the periplasmic compartment. Additionally, loss of secretion via the T2S system leads to a reduced growth rate, compromised outer membrane integrity, and induction of the extracytoplasmic stress factor RpoE (A. E. Sikora, S. R. Lybarger, and M. Sandkvist, J. Bacteriol. 189:8484-8495, 2007). In this study, gene expression profiling reveals that inactivation of the T2S system alters the expression of genes encoding cell envelope components and proteins involved in central metabolism, chemotaxis, motility, oxidative stress, and iron storage and acquisition. Consistent with the gene expression data, molecular and biochemical analyses indicate that the T2S mutants suffer from internal oxidative stress and increased levels of intracellular ferrous iron. By using a tolA mutant of V. cholerae that shares a similar compromised membrane phenotype but maintains a functional T2S machinery, we show that the formation of radical oxygen species, induction of oxidative stress, and changes in iron physiology are likely general responses to cell envelope damage and are not unique to T2S mutants. Finally, we demonstrate that disruption of the V. cholerae cell envelope by chemical treatment with polymyxin B similarly results in induction of the RpoE-mediated stress response, increased sensitivity to oxidants, and a change in iron metabolism. We propose that many types of extracytoplasmic stresses, caused either by genetic alterations of outer membrane constituents or by chemical or physical damage to the cell envelope, induce common signaling pathways that ultimately lead to internal oxidative stress and misregulation of iron homeostasis.

Deciphering the connectivity structure of biological networks using MixNet

Background

As biological networks often show complex topological features, mathematical methods are required to extract meaningful information. Clustering methods are useful in this setting, as they allow the summary of the network's topology into a small number of relevant classes. Different strategies are possible for clustering, and in this article we focus on a model-based strategy that aims at clustering nodes based on their connectivity profiles.

Results

We present MixNet, the first publicly available computer software that analyzes biological networks using mixture models. We apply this method to various networks such as the E. coli transcriptional regulatory network, the macaque cortex network, a foodweb network and the Buchnera aphidicola metabolic network. This method is also compared with other approaches such as module identification or hierarchical clustering.

Conclusion

We show how MixNet can be used to extract meaningful biological information, and to give a summary of the networks topology that highlights important biological features. This approach is powerful as MixNet is adaptive to the network under study, and finds structural information without any a priori on the structure that is investigated. This makes MixNet a very powerful tool to summarize and decipher the connectivity structure of biological networks.