Comparative genomic reconstruction of transcriptional regulatory networks in bacteria

Genomic reconstruction of transcriptional regulatory networks in lactic acid bacteria

Background

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

Results

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

Conclusions

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

Functional diversification of ROK-family transcriptional regulators of sugar catabolism in the Thermotogae phylum

Large and functionally heterogeneous families of transcription factors have complex evolutionary histories. What shapes specificities toward effectors and DNA sites in paralogous regulators is a fundamental question in biology. Bacteria from the deep-branching lineage Thermotogae possess multiple paralogs of the repressor, open reading frame, kinase (ROK) family regulators that are characterized by carbohydrate-sensing domains shared with sugar kinases. We applied an integrated genomic approach to study functions and specificities of regulators from this family. A comparative analysis of 11 Thermotogae genomes revealed novel mechanisms of transcriptional regulation of the sugar utilization networks, DNA-binding motifs and specific functions. Reconstructed regulons for seven groups of ROK regulators were validated by DNA-binding assays using purified recombinant proteins from the model bacterium Thermotoga maritima. All tested regulators demonstrated specific binding to their predicted cognate DNA sites, and this binding was inhibited by specific effectors, mono- or disaccharides from their respective sugar catabolic pathways. By comparing ligand-binding domains of regulators with structurally characterized kinases from the ROK family, we elucidated signature amino acid residues determining sugar-ligand regulator specificity. Observed correlations between signature residues and the sugar-ligand specificities provide the framework for structure functional classification of the entire ROK family.

Transcription factor family-based reconstruction of singleton regulons and study of the Crp/Fnr, ArsR, and GntR families in Desulfovibrionales genomes

Accurate detection of transcriptional regulatory elements is essential for high-quality genome annotation, metabolic reconstruction, and modeling of regulatory networks. We developed a computational approach for reconstruction of regulons operated by transcription factors (TFs) from large protein families and applied this novel approach to three TF families in 10 Desulfovibrionales genomes. Phylogenetic analyses of 125 regulators from the ArsR, Crp/Fnr, and GntR families revealed that 65% of these regulators (termed reference TFs) are well conserved in Desulfovibrionales, while the remaining 35% of regulators (termed singleton TFs) are species specific and show a mosaic distribution. For regulon reconstruction in the group of singleton TFs, the standard orthology-based approach was inefficient, and thus, we developed a novel approach based on the simultaneous study of all homologous TFs from the same family in a group of genomes. As a result, we identified binding for 21 singleton TFs and for all reference TFs in all three analyzed families. Within each TF family we observed structural similarities between DNA-binding motifs of different reference and singleton TFs. The collection of reconstructed regulons is available at the RegPrecise database (http://regprecise.lbl.gov/RegPrecise/Desulfovibrionales.jsp).

Evolution of transcriptional regulation in closely related bacteria

BACKGROUND

The exponential growth of the number of fully sequenced genomes at varying taxonomic closeness allows one to characterize transcriptional regulation using comparative-genomics analysis instead of time-consuming experimental methods. A transcriptional regulatory unit consists of a transcription factor, its binding site and a regulated gene. These units constitute a graph which contains so-called "network motifs", subgraphs of a given structure. Here we consider genomes of closely related Enterobacteriales and estimate the fraction of conserved network motifs and sites as well as positions under selection in various types of non-coding regions.

RESULTS

Using a newly developed technique, we found that the highest fraction of positions under selection, approximately 50%, was observed in synvergon spacers (between consecutive genes from the same strand), followed by ~45% in divergon spacers (common 5'-regions), and ~10% in convergon spacers (common 3'-regions). The fraction of selected positions in functional regions was higher, 60% in transcription factor-binding sites and ~45% in terminators and promoters. Small, but significant differences were observed between Escherichia coli and Salmonella enterica. This fraction is similar to the one observed in eukaryotes.The conservation of binding sites demonstrated some differences between types of regulatory units. In E. coli, strains the interactions of the type "local transcriptional factor gene" turned out to be more conserved in feed-forward loops (FFLs) compared to non-motif interactions. The coherent FFLs tend to be less conserved than the incoherent FFLs. A natural explanation is that the former imply functional redundancy.

CONCLUSIONS

A naïve hypothesis that FFL would be highly conserved turned out to be not entirely true: its conservation depends on its status in the transcriptional network and also from its usage. The fraction of positions under selection in intergenic regions of bacterial genomes is roughly similar to that of eukaryotes. Known regulatory sites explain 20±5% of selected positions.

Shrinking the metabolic solution space using experimental datasets

Constraint-based models of metabolism have been used in a variety of studies on drug discovery, metabolic engineering, evolution, and multi-species interactions. These genome-scale models can be generated for any sequenced organism since their main parameters (i.e., reaction stoichiometry) are highly conserved. Their relatively low parameter requirement makes these models easy to develop; however, these models often result in a solution space with multiple possible flux distributions, making it difficult to determine the precise flux state in the cell. Recent research efforts in this modeling field have investigated how additional experimental data, including gene expression, protein expression, metabolite concentrations, and kinetic parameters, can be used to reduce the solution space. This mini-review provides a summary of the data-driven computational approaches that are available for reducing the solution space and thereby improve predictions of intracellular fluxes by constraint-based models.

Tagaturonate-fructuronate epimerase UxaE, a novel enzyme in the hexuronate catabolic network in Thermotoga maritima

Thermotoga maritima is a marine hyperthermophilic microorganism that degrades a wide range of simple and complex carbohydrates including pectin and produces fermentative hydrogen at high yield. Galacturonate and glucuronate, two abundant hexuronic acids in pectin and xylan, respectively, are catabolized via committed metabolic pathways to supply carbon and energy for a variety of microorganisms. By a combination of bioinformatics and experimental techniques we identified a novel enzyme family (named UxaE) catalysing a previously unknown reaction in the hexuronic acid catabolic pathway, epimerization of tagaturonate to fructuronate. The enzymatic activity of the purified recombinant tagaturonate epimerase from T. maritima was directly confirmed and kinetically characterized. Its function was also confirmed by genetic complementation of the growth of the Escherichia coli uxaB knockout mutant strain on galacturonate. An inferred novel galacturonate to mannonate catabolic pathway in T. maritima was reconstituted in vitro using a mixture of recombinant purified enzymes UxaE, UxaC and UxuB. Members of the newly identified UxaE family were identified in ~50 phylogenetically diverse heterotrophic bacteria from aquatic and soil environments. The genomic context of respective genes and reconstruction of associated pathways suggest that UxaE enzymatic and biological function remains conserved in all of these species.

Transcriptional repression mediated by a TetR family protein, PfmR, from Thermus thermophilus HB8

PfmR is one of four TetR family transcriptional regulators found in the extremely thermophilic bacterium, Thermus thermophilus HB8. We identified three promoters with strong negative regulation by PfmR, both in vivo and in vitro. PfmR binds pseudopalindromic sequences, with the consensus sequence of 5'-TACCGACCGNTNGGTN-3' surrounding the promoters. According to the amino acid sequence and three-dimensional structure analyses of the PfmR-regulated gene products, they are predicted to be involved in phenylacetic acid and fatty acid metabolism. In vitro analyses revealed that PfmR weakly cross-regulated with the TetR family repressor T. thermophilus PaaR, which controls the expression of the paa gene cluster putatively involved in phenylacetic acid degradation but not with another functionally identified TetR family repressor, T. thermophilus FadR, which is involved in fatty acid degradation. The X-ray crystal structure of the N-terminal DNA-binding domain of PfmR and the nucleotide sequence of the predicted PfmR-binding site are quite similar to those of the TetR family repressor QacR from Staphylococcus aureus. Similar to QacR, two PfmR dimers bound per target DNA. The bases recognized by QacR within the QacR-binding site are conserved in the predicted PfmR-binding site, and they were important for PfmR to recognize the binding site and properly assemble on it. The center of the PfmR molecule contains a tunnel-like pocket, which may be the ligand-binding site of this regulator.

Differential regulation of antibiotic biosynthesis by DraR-K, a novel two-component system in Streptomyces coelicolor

A novel two-component system (TCS) designated as DraR-K (sco3063/sco3062) was identified to be involved in differential regulation of antibiotic biosynthesis in Streptomyces coelicolor. The S. coelicolor mutants with deletion of either or both of draR and draK exhibited significantly reduced actinorhodin (ACT) but increased undecylprodigiosin (RED) production on minimal medium (MM) supplemented separately with high concentration of different nitrogen sources. These mutants also overproduced a yellow-pigmented type I polyketide (yCPK) on MM with glutamate (Glu). It was confirmed that DraR-K activates ACT but represses yCPK production directly through the pathway-specific activator genes actII-ORF4 and kasO, respectively, while its role on RED biosynthesis was independent of pathway-specific activator genes redD/redZ. DNase I footprinting assays revealed that the DNA binding sites for DraR were at -124 to -98 nt and -24 to -1 nt relative to the respective transcription start point of actII-ORF4 and kasO. Comparison of the binding sites allowed the identification of a consensus DraR-binding sequence, 5'-AMAAWYMAKCA-3' (M: A or C; W: A or T; Y: C or T; K: G or T). By genome screening and gel-retardation assay, 11 new targets of DraR were further identified in the genome of S. coelicolor. Functional analysis of these tentative targets revealed the involvement of DraR-K in primary metabolism. DraR-K homologues are widely spread in different streptomycetes. Interestingly, deletion of draR-Ksav (sav_3481/sav_3480, homologue of draR-K) in the industrial model strain S. avermitilis NRRL-8165 led to similar abnormal antibiotic biosynthesis, showing higher avermectin while slightly decreased oligomycin A production, suggesting that DraR-K-mediated regulation system might be conserved in streptomycetes. This study further reveals the complexity of TCS in regulation of antibiotic biosynthesis in Streptomyces.

N-acetylgalactosamine utilization pathway and regulon in proteobacteria: genomic reconstruction and experimental characterization in Shewanella

We used a comparative genomics approach to reconstruct the N-acetyl-d-galactosamine (GalNAc) and galactosamine (GalN) utilization pathways and transcriptional regulons in Proteobacteria. The reconstructed GalNAc/GalN utilization pathways include multiple novel genes with specific functional roles. Most of the pathway variations were attributed to the amino sugar transport, phosphorylation, and deacetylation steps, whereas the downstream catabolic enzymes in the pathway were largely conserved. The predicted GalNAc kinase AgaK, the novel variant of GalNAc-6-phosphate deacetylase AgaA(II) and the GalN-6-phosphate deaminase AgaS from Shewanella sp. ANA-3 were validated in vitro using individual enzymatic assays and reconstitution of the three-step pathway. By using genetic techniques, we confirmed that AgaS but not AgaI functions as the main GalN-6-P deaminase in the GalNAc/GalN utilization pathway in Escherichia coli. Regulons controlled by AgaR repressors were reconstructed by bioinformatics in most proteobacterial genomes encoding GalNAc pathways. Candidate AgaR-binding motifs share a common sequence with consensus CTTTC that was found in multiple copies and arrangements in regulatory regions of aga genes. This study provides comprehensive insights into the common and distinctive features of the GalNAc/GalN catabolism and its regulation in diverse Proteobacteria.

RegPrecise web services interface: programmatic access to the transcriptional regulatory interactions in bacteria reconstructed by comparative genomics

Web services application programming interface (API) was developed to provide a programmatic access to the regulatory interactions accumulated in the RegPrecise database (http://regprecise.lbl.gov), a core resource on transcriptional regulation for the microbial domain of the Department of Energy (DOE) Systems Biology Knowledgebase. RegPrecise captures and visualize regulogs, sets of genes controlled by orthologous regulators in several closely related bacterial genomes, that were reconstructed by comparative genomics. The current release of RegPrecise 2.0 includes >1400 regulogs controlled either by protein transcription factors or by conserved ribonucleic acid regulatory motifs in >250 genomes from 24 taxonomic groups of bacteria. The reference regulons accumulated in RegPrecise can serve as a basis for automatic annotation of regulatory interactions in newly sequenced genomes. The developed API provides an efficient access to the RegPrecise data by a comprehensive set of 14 web service resources. The RegPrecise web services API is freely accessible at http://regprecise.lbl.gov/RegPrecise/services.jsp with no login requirements.

Phenylacetic acid catabolism and its transcriptional regulation in Corynebacterium glutamicum

The industrially important organism Corynebacterium glutamicum has been characterized in recent years for its robust ability to assimilate aromatic compounds. In this study, C. glutamicum strain AS 1.542 was investigated for its ability to catabolize phenylacetic acid (PAA). The paa genes were identified; they are organized as a continuous paa gene cluster. The type strain of C. glutamicum, ATCC 13032, is not able to catabolize PAA, but the recombinant strain ATCC 13032/pEC-K18mob2::paa gained the ability to grow on PAA. The paaR gene, encoding a TetR family transcription regulator, was studied in detail. Disruption of paaR in strain AS 1.542 resulted in transcriptional increases of all paa genes. Transcription start sites and putative promoter regions were determined. An imperfect palindromic motif (5'-ACTNACCGNNCGNNCGGTNAGT-3'; 22 bp) was identified in the upstream regions of paa genes. Electrophoretic mobility shift assays (EMSA) demonstrated specific binding of PaaR to this motif, and phenylacetyl coenzyme A (PA-CoA) blocked binding. It was concluded that PaaR is the negative regulator of PAA degradation and that PA-CoA is the PaaR effector. In addition, GlxR binding sites were found, and binding to GlxR was confirmed. Therefore, PAA catabolism in C. glutamicum is regulated by the pathway-specific repressor PaaR, and also likely by the global transcription regulator GlxR. By comparative genomic analysis, we reconstructed orthologous PaaR regulons in 57 species, including species of Actinobacteria, Proteobacteria, and Flavobacteria, that carry PAA utilization genes and operate by conserved binding motifs, suggesting that PaaR-like regulation might commonly exist in these bacteria.

Regulation and evolution of malonate and propionate catabolism in proteobacteria

Bacteria catabolize malonate via two pathways, encoded by the mdc and mat genes. In various bacteria, transcription of these genes is controlled by the GntR family transcription factors (TFs) MatR/MdcY and/or the LysR family transcription factor MdcR. Propionate is metabolized via the methylcitrate pathway, comprising enzymes encoded by the prp and acn genes. PrpR, the Fis family sigma 54-dependent transcription factor, is known to be a transcriptional activator of the prp genes. Here, we report a detailed comparative genomic analysis of malonate and propionate metabolism and its regulation in proteobacteria. We characterize genomic loci and gene regulation and identify binding motifs for four new TFs and also new regulon members, in particular, tripartite ATP-independent periplasmic (TRAP) transporters. We describe restructuring of the genomic loci and regulatory interactions during the evolution of proteobacteria.

Characterization of the biotin uptake system encoded by the biotin-inducible bioYMN operon of Corynebacterium glutamicum

Background

The amino acid-producing Gram-positive Corynebacterium glutamicum is auxotrophic for biotin although biotin ring assembly starting from the precursor pimeloyl-CoA is still functional. It possesses AccBC, the α-subunit of the acyl-carboxylases involved in fatty acid and mycolic acid synthesis, and pyruvate carboxylase as the only biotin-containing proteins. Comparative genome analyses suggested that the putative transport system BioYMN encoded by cg2147, cg2148 and cg2149 might be involved in biotin uptake by C. glutamicum.

Results

By comparison of global gene expression patterns of cells grown with limiting or excess supply of biotin or with dethiobiotin as supplement replacing biotin revealed that expression of genes coding for enzymes of biotin ring assembly and for the putative uptake system was regulated according to biotin availability. RT-PCR and 5'-RACE experiments demonstrated that the genes bioY, bioM, and bioN are transcribed from one promoter as a single transcript. Biochemical analyses revealed that BioYMN catalyzes the effective uptake of biotin with a concentration of 60 nM biotin supporting a half-maximal transport rate. Maximal biotin uptake rates were at least five fold higher in biotin-limited cells as compared to cells grown with excess biotin. Overexpression of bioYMN led to an at least 50 fold higher biotin uptake rate as compared to the empty vector control. Overproduction of BioYMN alleviated biotin limitation and interfered with triggering L-glutamate production by biotin limitation.

Conclusions

The operon bioYMN from C. glutamicum was shown to be induced by biotin limitation. Transport assays with radio-labeled biotin revealed that BioYMN functions as a biotin uptake system. Overexpression of bioYMN affected L-glutamate production triggered by biotin limitation.

Genome-wide mapping of the binding sites of proteins that interact with DNA

The coordinated regulation of the expression of a group of genes by a specific transcription factor frequently lies at the heart of the ability of a bacterium to respond to an environmental signal, or to progress through a developmental program. Thus, in many situations, it is of interest to identify all of the genes that are under the control of a particular regulatory protein. This chapter begins with a brief overview of some of the methods that have been used in attempts to identify some or all of the members of a regulon (i.e., those genes that are the targets for a transcriptional activator or repressor). Thereafter, the chapter will focus on one technique, chromatin immunoprecipitation and microarray analysis (ChIP-chip) and some of its variants. Design considerations and some protocols for ChIP-chip experiments are provided, along with some considerations related to downstream data analysis. ChIP-chip is a method for the genome-wide localization of protein-binding sites. In a typical ChIP-chip protocol, proteins are cross-linked nonspecifically to DNA in vivo. Chromatin is extracted and sheared, and specific protein-DNA complexes are immunoprecipitated with a suitable antibody. After purification, the DNA is hybridized to a microarray (after an amplification step in some protocols), together with a differentially labeled reference sample. Features on the microarray that show an elevated fluorescence ratio reveal DNA sequences that were enriched by immunoprecipitation. The corresponding genomic locations are those that were enriched, and are therefore close to sites of binding. The use of high-density tiled microarrays allows for binding site localization with quite high resolution. It is likely that ChIP-chip will soon be superseded by ChIP-seq, in which the immunoprecipitated DNA is analyzed directly by next-generation sequencing technologies. ChIP-chip and ChIP-seq applications are not confined to regulatory proteins, since they can be used with any protein that binds to DNA, either directly, or indirectly via an interaction with another protein. Thus, ChIP-chip has been used successfully to map binding sites for nucleoid proteins, and proteins involved in DNA replication.

Transcriptional regulation of central carbon and energy metabolism in bacteria by redox-responsive repressor Rex

Redox-sensing repressor Rex was previously implicated in the control of anaerobic respiration in response to the cellular NADH/NAD(+) levels in gram-positive bacteria. We utilized the comparative genomics approach to infer candidate Rex-binding DNA motifs and assess the Rex regulon content in 119 genomes from 11 taxonomic groups. Both DNA-binding and NAD-sensing domains are broadly conserved in Rex orthologs identified in the phyla Firmicutes, Thermotogales, Actinobacteria, Chloroflexi, Deinococcus-Thermus, and Proteobacteria. The identified DNA-binding motifs showed significant conservation in these species, with the only exception detected in Clostridia, where the Rex motif deviates in two positions from the generalized consensus, TTGTGAANNNNTTCACAA. Comparative analysis of candidate Rex sites revealed remarkable variations in functional repertoires of candidate Rex-regulated genes in various microorganisms. Most of the reconstructed regulatory interactions are lineage specific, suggesting frequent events of gain and loss of regulator binding sites in the evolution of Rex regulons. We identified more than 50 novel Rex-regulated operons encoding functions that are essential for resumption of the NADH:NAD(+) balance. The novel functional role of Rex in the control of the central carbon metabolism and hydrogen production genes was validated by in vitro DNA binding assays using the TM0169 protein in the hydrogen-producing bacterium Thermotoga maritima.

Ribulokinase and transcriptional regulation of arabinose metabolism in Clostridium acetobutylicum

The transcription factor AraR controls utilization of L-arabinose in Bacillus subtilis. In this study, we combined a comparative genomic reconstruction of AraR regulons in nine Clostridium species with detailed experimental characterization of AraR-mediated regulation in Clostridium acetobutylicum. Based on the reconstructed AraR regulons, a novel ribulokinase, AraK, present in all analyzed Clostridium species was identified, which was a nonorthologous replacement of previously characterized ribulokinases. The predicted function of the araK gene was confirmed by inactivation of the araK gene in C. acetobutylicum and biochemical assays using purified recombinant AraK. In addition to the genes involved in arabinose utilization and arabinoside degradation, extension of the AraR regulon to the pentose phosphate pathway genes in several Clostridium species was revealed. The predicted AraR-binding sites in the C. acetobutylicum genome and the negative effect of L-arabinose on DNA-regulator complex formation were verified by in vitro binding assays. The predicted AraR-controlled genes in C. acetobutylicum were experimentally validated by testing gene expression patterns in both wild-type and araR-inactivated mutant strains during growth in the absence or presence of L-arabinose.

Negative transcriptional control of biotin metabolism genes by the TetR-type regulator BioQ in biotin-auxotrophic Corynebacterium glutamicum ATCC 13032

Genomic context analysis in actinobacteria revealed that biotin biosynthesis and transport (bio) genes are co-localized in several genomes with a gene encoding a transcription regulator of the TetR protein family, now named BioQ. Comparative analysis of the upstream regions of bio genes identified the common 13-bp palindromic motif TGAAC-N3-GTTAC as candidate BioQ-binding site. To verify the role of BioQ in controlling the transcription of bio genes, a deletion in the bioQ coding region (cg2309) was constructed in Corynebacterium glutamicum ATCC 13032, resulting in the mutant strain C. glutamicum IB2309. Comparative whole-genome DNA microarray hybridizations and subsequent expression analyses by real-time reverse transcriptase PCR revealed enhanced transcript levels of all bio genes in C. glutamicum IB2309, when compared with the wild-type strain ATCC 13032. Accordingly, the BioQ protein of C. glutamicum acts as a repressor of ten genes that are organized in four transcription units: bioA-bioD, cg2884-cg2883, bioB-cg0096-cg0097, and bioY-bioM-bioN. DNA band shift assays with an intein-tagged BioQ protein demonstrated the specific binding of the purified protein to DNA fragments containing the candidate BioQ-binding sites, which were located within the mapped promoter regions of bioA, cg2884, bioB, and bioY. These data confirmed the direct regulatory role of BioQ in the control of biotin biosynthesis and transport genes in C. glutamicum. Differential expression of bio genes in C. glutamicum IB2309 was moreover complemented by bioQ genes cloned from other corynebacterial genomes.

Complete genome sequence of Corynebacterium variabile DSM 44702 isolated from the surface of smear-ripened cheeses and insights into cheese ripening and flavor generation

Background

Corynebacterium variabile is part of the complex microflora on the surface of smear-ripened cheeses and contributes to the development of flavor and textural properties during cheese ripening. Still little is known about the metabolic processes and microbial interactions during the production of smear-ripened cheeses. Therefore, the gene repertoire contributing to the lifestyle of the cheese isolate C. variabile DSM 44702 was deduced from the complete genome sequence to get a better understanding of this industrial process.

Results

The chromosome of C. variabile DSM 44702 is composed of 3, 433, 007 bp and contains 3, 071 protein-coding regions. A comparative analysis of this gene repertoire with that of other corynebacteria detected 1, 534 predicted genes to be specific for the cheese isolate. These genes might contribute to distinct metabolic capabilities of C. variabile, as several of them are associated with metabolic functions in cheese habitats by playing roles in the utilization of alternative carbon and sulphur sources, in amino acid metabolism, and fatty acid degradation. Relevant C. variabile genes confer the capability to catabolize gluconate, lactate, propionate, taurine, and gamma-aminobutyric acid and to utilize external caseins. In addition, C. variabile is equipped with several siderophore biosynthesis gene clusters for iron acquisition and an exceptional repertoire of AraC-regulated iron uptake systems. Moreover, C. variabile can produce acetoin, butanediol, and methanethiol, which are important flavor compounds in smear-ripened cheeses.

Conclusions

The genome sequence of C. variabile provides detailed insights into the distinct metabolic features of this bacterium, implying a strong adaption to the iron-depleted cheese surface habitat. By combining in silico data obtained from the genome annotation with previous experimental knowledge, occasional observations on genes that are involved in the complex metabolic capacity of C. variabile were integrated into a global view on the lifestyle of this species.

Control of proteobacterial central carbon metabolism by the HexR transcriptional regulator: a case study in Shewanella oneidensis

Bacteria exploit multiple mechanisms for controlling central carbon metabolism (CCM). Thus, a bioinformatic analysis together with some experimental data implicated the HexR transcriptional factor as a global CCM regulator in some lineages of Gammaproteobacteria operating as a functional replacement of the Cra regulator characteristic of Enterobacteriales. In this study, we combined a large scale comparative genomic reconstruction of HexR-controlled regulons in 87 species of Proteobacteria with the detailed experimental analysis of the HexR regulatory network in the Shewanella oneidensis model system. Although nearly all of the HexR-controlled genes are associated with CCM, remarkable variations were revealed in the scale (from 1 to 2 target operons in Enterobacteriales up to 20 operons in Aeromonadales) and gene content of HexR regulons between 11 compared lineages. A predicted 17-bp pseudo-palindrome with a consensus tTGTAATwwwATTACa was confirmed as a HexR-binding motif for 15 target operons (comprising 30 genes) by in vitro binding assays. The negative effect of the key CCM intermediate, 2-keto-3-deoxy-6-phosphogluconate, on the DNA-regulator complex formation was verified. A dual mode of HexR action on various target promoters, repression of genes involved in catabolic pathways and activation of gluconeogenic genes, was for the first time predicted by the bioinformatic analysis and experimentally verified by changed gene expression pattern in S. oneidensis ΔhexR mutant. Phenotypic profiling revealed the inability of this mutant to grow on lactate or pyruvate as a single carbon source. A comparative metabolic flux analysis of wild-type and mutant strains of S. oneidensis using [(13)C]lactate labeling and GC-MS analysis confirmed the hypothesized HexR role as a master regulator of gluconeogenic flux from pyruvate via the transcriptional activation of phosphoenolpyruvate synthase (PpsA).

Comparative analyses imply that the enigmatic Sigma factor 54 is a central controller of the bacterial exterior

Background

Sigma-54 is a central regulator in many pathogenic bacteria and has been linked to a multitude of cellular processes like nitrogen assimilation and important functional traits such as motility, virulence, and biofilm formation. Until now it has remained obscure whether these phenomena and the control by Sigma-54 share an underlying theme.

Results

We have uncovered the commonality by performing a range of comparative genome analyses. A) The presence of Sigma-54 and its associated activators was determined for all sequenced prokaryotes. We observed a phylum-dependent distribution that is suggestive of an evolutionary relationship between Sigma-54 and lipopolysaccharide and flagellar biosynthesis. B) All Sigma-54 activators were identified and annotated. The relation with phosphotransfer-mediated signaling (TCS and PTS) and the transport and assimilation of carboxylates and nitrogen containing metabolites was substantiated. C) The function annotations, that were represented within the genomic context of all genes encoding Sigma-54, its activators and its promoters, were analyzed for intra-phylum representation and inter-phylum conservation. Promoters were localized using a straightforward scoring strategy that was formulated to identify similar motifs. We found clear highly-represented and conserved genetic associations with genes that concern the transport and biosynthesis of the metabolic intermediates of exopolysaccharides, flagella, lipids, lipopolysaccharides, lipoproteins and peptidoglycan.

Conclusion

Our analyses directly implicate Sigma-54 as a central player in the control over the processes that involve the physical interaction of an organism with its environment like in the colonization of a host (virulence) or the formation of biofilm.

Comparative genomic reconstruction of transcriptional networks controlling central metabolism in the Shewanella genus

Background

Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in bacteria is one of the critical tasks of modern genomics. The Shewanella genus is comprised of metabolically versatile gamma-proteobacteria, whose lifestyles and natural environments are substantially different from Escherichia coli and other model bacterial species. The comparative genomics approaches and computational identification of regulatory sites are useful for the in silico reconstruction of transcriptional regulatory networks in bacteria.

Results

To explore conservation and variations in the Shewanella transcriptional networks we analyzed the repertoire of transcription factors and performed genomics-based reconstruction and comparative analysis of regulons in 16 Shewanella genomes. The inferred regulatory network includes 82 transcription factors and their DNA binding sites, 8 riboswitches and 6 translational attenuators. Forty five regulons were newly inferred from the genome context analysis, whereas others were propagated from previously characterized regulons in the Enterobacteria and Pseudomonas spp.. Multiple variations in regulatory strategies between the Shewanella spp. and E. coli include regulon contraction and expansion (as in the case of PdhR, HexR, FadR), numerous cases of recruiting non-orthologous regulators to control equivalent pathways (e.g. PsrA for fatty acid degradation) and, conversely, orthologous regulators to control distinct pathways (e.g. TyrR, ArgR, Crp).

Conclusions

We tentatively defined the first reference collection of ~100 transcriptional regulons in 16 Shewanella genomes. The resulting regulatory network contains ~600 regulated genes per genome that are mostly involved in metabolism of carbohydrates, amino acids, fatty acids, vitamins, metals, and stress responses. Several reconstructed regulons including NagR for N-acetylglucosamine catabolism were experimentally validated in S. oneidensis MR-1. Analysis of correlations in gene expression patterns helps to interpret the reconstructed regulatory network. The inferred regulatory interactions will provide an additional regulatory constrains for an integrated model of metabolism and regulation in S. oneidensis MR-1.

The transcriptional response of Listeria monocytogenes during adaptation to growth on lactate and diacetate includes synergistic changes that increase fermentative acetoin production

The organic acids lactate and diacetate are commonly used in combination in ready-to-eat foods because they show synergistic ability to inhibit the growth of Listeria monocytogenes. Full-genome microarrays were used to investigate the synergistic transcriptomic responses of two L. monocytogenes strains, H7858 (serotype 4b) and F6854 (serotype 1/2a), to these two organic acids under conditions representing osmotic and cold stress encountered in foods. Strains were exposed to brain heart infusion (BHI) broth at 7°C with 4.65% water-phase (w.p.) NaCl at pH 6.1 with (i) 2% w.p. potassium lactate, (ii) 0.14% w.p. sodium diacetate, (iii) the combination of both at the same levels, or (iv) no organic acids as a control. RNA was extracted 8 h after exposure, during lag phase, to capture gene transcription changes during adaptation to the organic acid stress. Significant differential transcription of 1,041 genes in H7858 and 640 genes in F6854 was observed in at least one pair of the 4 different treatments. The effects of combined treatment with lactate and diacetate included (i) synergistic transcription differences for 474 and 209 genes in H7858 and F6854, respectively, (ii) differential transcription of genes encoding cation transporters and ABC transporters of metals, and (iii) altered metabolism, including induction of a nutrient-limiting stress response, reduction of menaquinone biosynthesis, and a shift from fermentative production of acetate and lactate to energetically less favorable, neutral acetoin. These data suggest that additional treatments that interfere with cellular energy generation processes could more efficiently inhibit the growth of L. monocytogenes.

Comparative genomic analysis of the hexuronate metabolism genes and their regulation in gammaproteobacteria

The hexuronate metabolism in Escherichia coli is regulated by two related transcription factors from the FadR subfamily of the GntR family, UxuR and ExuR. UxuR controls the d-glucuronate metabolism, while ExuR represses genes involved in the metabolism of all hexuronates. We use a comparative genomics approach to reconstruct the hexuronate metabolic pathways and transcriptional regulons in gammaproteobacteria. We demonstrate differences in the binding motifs of UxuR and ExuR, identify new candidate members of the UxuR/ExuR regulons, and describe the links between the UxuR/ExuR regulons and the adjacent regulons UidR, KdgR, and YjjM. We provide experimental evidence that two predicted members of the UxuR regulon, yjjM and yjjN, are the subject of complex regulation by this transcription factor in E. coli.

Inference of the transcriptional regulatory network in Staphylococcus aureus by integration of experimental and genomics-based evidence

Transcriptional regulatory networks are fine-tuned systems that help microorganisms respond to changes in the environment and cell physiological state. We applied the comparative genomics approach implemented in the RegPredict Web server combined with SEED subsystem analysis and available information on known regulatory interactions for regulatory network reconstruction for the human pathogen Staphylococcus aureus and six related species from the family Staphylococcaceae. The resulting reference set of 46 transcription factor regulons contains more than 1,900 binding sites and 2,800 target genes involved in the central metabolism of carbohydrates, amino acids, and fatty acids; respiration; the stress response; metal homeostasis; drug and metal resistance; and virulence. The inferred regulatory network in S. aureus includes ∼320 regulatory interactions between 46 transcription factors and ∼550 candidate target genes comprising 20% of its genome. We predicted ∼170 novel interactions and 24 novel regulons for the control of the central metabolic pathways in S. aureus. The reconstructed regulons are largely variable in the Staphylococcaceae: only 20% of S. aureus regulatory interactions are conserved across all studied genomes. We used a large-scale gene expression data set for S. aureus to assess relationships between the inferred regulons and gene expression patterns. The predicted reference set of regulons is captured within the Staphylococcus collection in the RegPrecise database (http://regprecise.lbl.gov).

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background

Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor σ^S (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain.

Results

Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms.

Conclusions

The ppGpp and σ^S variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

Carbohydrate metabolism in Bifidobacteria

Members of the genus Bifidobacterium can be found as components of the gastrointestinal microbiota, and are believed to play an important role in maintaining and promoting human health by eliciting a number of beneficial properties. Bifidobacteria can utilize a diverse range of dietary carbohydrates that escape degradation in the upper parts of the intestine, many of which are plant-derived oligo- and polysaccharides. The gene content of a bifidobacterial genome reflects this apparent metabolic adaptation to a complex carbohydrate-rich gastrointestinal tract environment as it encodes a large number of predicted carbohydrate-modifying enzymes. Different bifidobacterial strains may possess different carbohydrate utilizing abilities, as established by a number of studies reviewed here. Carbohydrate-degrading activities described for bifidobacteria and their relevance to the deliberate enhancement of number and/or activity of bifidobacteria in the gut are also discussed in this review.

The complete genome sequence of Corynebacterium pseudotuberculosis FRC41 isolated from a 12-year-old girl with necrotizing lymphadenitis reveals insights into gene-regulatory networks contributing to virulence

BACKGROUND

Corynebacterium pseudotuberculosis is generally regarded as an important animal pathogen that rarely infects humans. Clinical strains are occasionally recovered from human cases of lymphadenitis, such as C. pseudotuberculosis FRC41 that was isolated from the inguinal lymph node of a 12-year-old girl with necrotizing lymphadenitis. To detect potential virulence factors and corresponding gene-regulatory networks in this human isolate, the genome sequence of C. pseudotuberculosis FCR41 was determined by pyrosequencing and functionally annotated.

RESULTS

Sequencing and assembly of the C. pseudotuberculosis FRC41 genome yielded a circular chromosome with a size of 2,337,913 bp and a mean G+C content of 52.2%. Specific gene sets associated with iron and zinc homeostasis were detected among the 2,110 predicted protein-coding regions and integrated into a gene-regulatory network that is linked with both the central metabolism and the oxidative stress response of FRC41. Two gene clusters encode proteins involved in the sortase-mediated polymerization of adhesive pili that can probably mediate the adherence to host tissue to facilitate additional ligand-receptor interactions and the delivery of virulence factors. The prominent virulence factors phospholipase D (Pld) and corynebacterial protease CP40 are encoded in the genome of this human isolate. The genome annotation revealed additional serine proteases, neuraminidase H, nitric oxide reductase, an invasion-associated protein, and acyl-CoA carboxylase subunits involved in mycolic acid biosynthesis as potential virulence factors. The cAMP-sensing transcription regulator GlxR plays a key role in controlling the expression of several genes contributing to virulence.

CONCLUSION

The functional data deduced from the genome sequencing and the extended knowledge of virulence factors indicate that the human isolate C. pseudotuberculosis FRC41 is equipped with a distinct gene set promoting its survival under unfavorable environmental conditions encountered in the mammalian host.

Predicting DNA-binding specificities of eukaryotic transcription factors

Today, annotated amino acid sequences of more and more transcription factors (TFs) are readily available. Quantitative information about their DNA-binding specificities, however, are hard to obtain. Position frequency matrices (PFMs), the most widely used models to represent binding specificities, are experimentally characterized only for a small fraction of all TFs. Even for some of the most intensively studied eukaryotic organisms (i.e., human, rat and mouse), roughly one-sixth of all proteins with annotated DNA-binding domain have been characterized experimentally. Here, we present a new method based on support vector regression for predicting quantitative DNA-binding specificities of TFs in different eukaryotic species. This approach estimates a quantitative measure for the PFM similarity of two proteins, based on various features derived from their protein sequences. The method is trained and tested on a dataset containing 1 239 TFs with known DNA-binding specificity, and used to predict specific DNA target motifs for 645 TFs with high accuracy.

Evolution of prokaryotic genes by shift of stop codons

De novo origin of coding sequence remains an obscure issue in molecular evolution. One of the possible paths for addition (subtraction) of DNA segments to (from) a gene is stop codon shift. Single nucleotide substitutions can destroy the existing stop codon, leading to uninterrupted translation up to the next stop codon in the gene's reading frame, or create a premature stop codon via a nonsense mutation. Furthermore, short indels-caused frameshifts near gene's end may lead to premature stop codons or to translation past the existing stop codon. Here, we describe the evolution of the length of coding sequence of prokaryotic genes by change of positions of stop codons. We observed cases of addition of regions of 3'UTR to genes due to mutations at the existing stop codon, and cases of subtraction of C-terminal coding segments due to nonsense mutations upstream of the stop codon. Many of the observed stop codon shifts cannot be attributed to sequencing errors or rare deleterious variants segregating within bacterial populations. The additions of regions of 3'UTR tend to occur in those genes in which they are facilitated by nearby downstream in-frame triplets which may serve as new stop codons. Conversely, subtractions of coding sequence often give rise to in-frame stop codons located nearby. The amino acid composition of the added region is significantly biased, compared to the overall amino acid composition of the genes. Our results show that in prokaryotes, shift of stop codon is an underappreciated contributor to functional evolution of gene length.

Premature terminator analysis sheds light on a hidden world of bacterial transcriptional attenuation

BACKGROUND

Bacterial transcription attenuation occurs through a variety of cis-regulatory elements that control gene expression in response to a wide range of signals. The signal-sensing structures in attenuators are so diverse and rapidly evolving that only a small fraction have been properly annotated and characterized to date. Here we apply a broad-spectrum detection tool in order to achieve a more complete view of the transcriptional attenuation complement of key bacterial species.

RESULTS

Our protocol seeks gene families with an unusual frequency of 5' terminators found across multiple species. Many of the detected attenuators are part of annotated elements, such as riboswitches or T-boxes, which often operate through transcriptional attenuation. However, a significant fraction of candidates were not previously characterized in spite of their unmistakable footprint. We further characterized some of these new elements using sequence and secondary structure analysis. We also present elements that may control the expression of several non-homologous genes, suggesting co-transcription and response to common signals. An important class of such elements, which we called mobile attenuators, is provided by 3' terminators of insertion sequences or prophages that may be exapted as 5' regulators when inserted directly upstream of a cellular gene.

CONCLUSIONS

We show here that attenuators involve a complex landscape of signal-detection structures spanning the entire bacterial domain. We discuss possible scenarios through which these diverse 5' regulatory structures may arise or evolve.

Positive transcriptional control of the pyridoxal phosphate biosynthesis genes pdxST by the MocR-type regulator PdxR of Corynebacterium glutamicum ATCC 13032

The pdxR (cg0897) gene of Corynebacterium glutamicum ATCC 13032 encodes a regulatory protein belonging to the MocR subfamily of GntR-type transcription regulators and consisting of an amino-terminal winged helix-turn-helix DNA-binding domain and a carboxy-terminal aminotransferase-like domain. A defined deletion in the pdxR gene resulted in the decreased expression of the divergently orientated pdxST genes coding for the subunits of pyridoxal 5'-phosphate synthase. The pdxST mutant C. glutamicum NJ0898 and the pdxR mutant C. glutamicum AMH17 showed vitamin B(6) auxotrophy that was restored by supplementing the growth medium with either pyridoxal, pyridoxal 5'-phosphate or pyridoxamine. The genetic organization of the 89 bp pdxR-pdxST intergenic region was elucidated by mapping the 5' ends of the respective transcripts, followed by detection of typical promoter sequences. Bioinformatic pattern searches and comparative genomics revealed three DNA motifs with the consensus sequence AAAGTGGW(-/T)CTA, overlapping the deduced promoter sequences and serving as candidate DNA-binding sites for PdxR. DNA band shift assays with the purified PdxR protein demonstrated the specific binding of the transcription regulator to double-stranded 40-mer sequences containing the detected motifs, thereby confirming the direct regulatory role of PdxR in activating the expression of the pdxST genes.

Genomic encyclopedia of sugar utilization pathways in the Shewanella genus

Background

Carbohydrates are a primary source of carbon and energy for many bacteria. Accurate projection of known carbohydrate catabolic pathways across diverse bacteria with complete genomes constitutes a substantial challenge due to frequent variations in components of these pathways. To address a practically and fundamentally important challenge of reconstruction of carbohydrate utilization machinery in any microorganism directly from its genomic sequence, we combined a subsystems-based comparative genomic approach with experimental validation of selected bioinformatic predictions by a combination of biochemical, genetic and physiological experiments.

Results

We applied this integrated approach to systematically map carbohydrate utilization pathways in 19 genomes from the Shewanella genus. The obtained genomic encyclopedia of sugar utilization includes ~170 protein families (mostly metabolic enzymes, transporters and transcriptional regulators) spanning 17 distinct pathways with a mosaic distribution across Shewanella species providing insights into their ecophysiology and adaptive evolution. Phenotypic assays revealed a remarkable consistency between predicted and observed phenotype, an ability to utilize an individual sugar as a sole source of carbon and energy, over the entire matrix of tested strains and sugars.

Comparison of the reconstructed catabolic pathways with E. coli identified multiple differences that are manifested at various levels, from the presence or absence of certain sugar catabolic pathways, nonorthologous gene replacements and alternative biochemical routes to a different organization of transcription regulatory networks.

Conclusions

The reconstructed sugar catabolome in Shewanella spp includes 62 novel isofunctional families of enzymes, transporters, and regulators. In addition to improving our knowledge of genomics and functional organization of carbohydrate utilization in Shewanella, this study led to a substantial expansion of our current version of the Genomic Encyclopedia of Carbohydrate Utilization. A systematic and iterative application of this approach to multiple taxonomic groups of bacteria will further enhance it, creating a knowledge base adequate for the efficient analysis of any newly sequenced genome as well as of the emerging metagenomic data.

A eubacterial riboswitch class that senses the coenzyme tetrahydrofolate

Comparative sequence analyses of bacterial genomes are revealing many structured RNA motifs that function as metabolite-binding riboswitches. We have identified an RNA motif frequently positioned in the 5' UTRs of folate transport and biosynthesis genes in Firmicute genomes. Biochemical experiments confirm that representatives of this new-found RNA class selectively bind derivatives of the vitamin folate, including di- and tetrahydrofolate coenzymes. In addition, representatives of this aptamer class occasionally reside upstream of RNA structures that are predicted to control translation initiation in response to ligand binding. These findings expand the number of coenzymes that are directly sensed by RNA and reveal possible riboswitch-controlled regulons that respond to changes in single-carbon metabolism.

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

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

OptORF: Optimal metabolic and regulatory perturbations for metabolic engineering of microbial strains

Background

Computational modeling and analysis of metabolic networks has been successful in metabolic engineering of microbial strains for valuable biochemical production. Limitations of currently available computational methods for metabolic engineering are that they are often based on reaction deletions rather than gene deletions and do not consider the regulatory networks that control metabolism. Due to the presence of multi-functional enzymes and isozymes, computational designs based on reaction deletions can sometimes result in strategies that are genetically complicated or infeasible. Additionally, strains might not be able to grow initially due to regulatory restrictions. To overcome these limitations, we have developed a new approach (OptORF) for identifying metabolic engineering strategies based on gene deletion and overexpression.

Results

Here we propose an effective method to systematically integrate transcriptional regulatory networks and metabolic networks. This allows for the formulation of linear optimization problems that search for metabolic and/or regulatory perturbations that couple biomass and biochemical production, thus proposing adaptive evolutionary strain designs. Using genome-scale models of Escherichia coli, we have implemented the OptORF algorithm (which considers gene deletions and transcriptional regulation) and compared its metabolic engineering strategies for ethanol production to those found using OptKnock (which considers reaction deletions). Our results found that the reaction-based strategies often require more gene deletions to remove the identified reactions (2 more genes than reactions), and result in lethal growth phenotypes when transcriptional regulation is considered (162 out of 200 cases). Finally, we present metabolic engineering strategies for producing ethanol and higher alcohols (e.g. isobutanol) in E. coli using our OptORF approach. We have found common genetic modifications such as deletion of pgi and overexpression of edd, as well as chemical specific strategies for producing different alcohols.

Conclusions

By taking regulatory effects into account, OptORF can propose changes such as the overexpression of metabolic genes or deletion of transcriptional factors, in addition to the deletion of metabolic genes, that may lead to faster evolutionary trajectories. While biofuel production in E. coli is evaluated here, the developed OptORF approach is general and can be applied to optimize the production of different compounds in other biological systems.

Reconstruction of xylose utilization pathway and regulons in Firmicutes

Background

Many Firmicutes bacteria, including solvent-producing clostridia such as Clostridium acetobutylicum, are able to utilize xylose, an abundant carbon source in nature. Nevertheless, homology searches failed to recognize all the genes for the complete xylose and xyloside utilization pathway in most of them. Moreover, the regulatory mechanisms of xylose catabolism in many Firmicutes except Bacillus spp. still remained unclear.

Results

A comparative genomic approach was used to reconstruct the xylose and xyloside utilization pathway and analyze its regulatory mechanisms in 24 genomes of the Firmicutes. A novel xylose isomerase that is not homologous to previously characterized xylose isomerase, was identified in C. acetobutylicum and several other Clostridia species. The candidate genes for the xylulokinase, xylose transporters, and the transcriptional regulator of xylose metabolism (XylR), were unambiguously assigned in all of the analyzed species based on the analysis of conserved chromosomal gene clustering and regulons. The predicted functions of these genes in C. acetobutylicum were experimentally confirmed through a combination of genetic and biochemical techniques. XylR regulons were reconstructed by identification and comparative analysis of XylR-binding sites upstream of xylose and xyloside utilization genes. A novel XylR-binding DNA motif, which is exceptionally distinct from the DNA motif known for Bacillus XylR, was identified in three Clostridiales species and experimentally validated in C. acetobutylicum by an electrophoretic mobility shift assay.

Conclusions

This study provided comprehensive insights to the xylose catabolism and its regulation in diverse Firmicutes bacteria especially Clostridia species, and paved ways for improving xylose utilization capability in C. acetobutylicum by genetic engineering in the future.

Network inference and network response identification: moving genome-scale data to the next level of biological discovery

The escalating amount of genome-scale data demands a pragmatic stance from the research community. How can we utilize this deluge of information to better understand biology, cure diseases, or engage cells in bioremediation or biomaterial production for various purposes? A research pipeline moving new sequence, expression and binding data towards practical end goals seems to be necessary. While most individual researchers are not motivated by such well-articulated pragmatic end goals, the scientific community has already self-organized itself to successfully convert genomic data into fundamentally new biological knowledge and practical applications. Here we review two important steps in this workflow: network inference and network response identification, applied to transcriptional regulatory networks. Among network inference methods, we concentrate on relevance networks due to their conceptual simplicity. We classify and discuss network response identification approaches as either data-centric or network-centric. Finally, we conclude with an outlook on what is still missing from these approaches and what may be ahead on the road to biological discovery.

The Zur regulon of Corynebacterium glutamicum ATCC 13032

Background

Zinc is considered as an essential element for all living organisms, but it can be toxic at large concentrations. Bacteria therefore tightly regulate zinc metabolism. The Cg2502 protein of Corynebacterium glutamicum was a candidate to control zinc metabolism in this species, since it was classified as metalloregulator of the zinc uptake regulator (Zur) subgroup of the ferric uptake regulator (Fur) family of DNA-binding transcription regulators.

Results

The cg2502 (zur) gene was deleted in the chromosome of C. glutamicum ATCC 13032 by an allelic exchange procedure to generate the zur-deficient mutant C. glutamicum JS2502. Whole-genome DNA microarray hybridizations and real-time RT-PCR assays comparing the gene expression in C. glutamicum JS2502 with that of the wild-type strain detected 18 genes with enhanced expression in the zur mutant. The expression data were combined with results from cross-genome comparisons of shared regulatory sites, revealing the presence of candidate Zur-binding sites in the mapped promoter regions of five transcription units encoding components of potential zinc ABC-type transporters (cg0041-cg0042/cg0043; cg2911-cg2912-cg2913), a putative secreted protein (cg0040), a putative oxidoreductase (cg0795), and a putative P-loop GTPase of the COG0523 protein family (cg0794). Enhanced transcript levels of the respective genes in C. glutamicum JS2502 were verified by real-time RT-PCR, and complementation of the mutant with a wild-type zur gene reversed the effect of differential gene expression. The zinc-dependent expression of the putative cg0042 and cg2911 operons was detected in vivo with a gfp reporter system. Moreover, the zinc-dependent binding of purified Zur protein to double-stranded 40-mer oligonucleotides containing candidate Zur-binding sites was demonstrated in vitro by DNA band shift assays.

Conclusion

Whole-genome expression profiling and DNA band shift assays demonstrated that Zur directly represses in a zinc-dependent manner the expression of nine genes organized in five transcription units. Accordingly, the Zur (Cg2502) protein is the key transcription regulator for genes involved in zinc homeostasis in C. glutamicum.

Zinc-independent folate biosynthesis: genetic, biochemical, and structural investigations reveal new metal dependence for GTP cyclohydrolase IB

GTP cyclohydrolase I (GCYH-I) is an essential Zn(2+)-dependent enzyme that catalyzes the first step of the de novo folate biosynthetic pathway in bacteria and plants, the 7-deazapurine biosynthetic pathway in Bacteria and Archaea, and the biopterin pathway in mammals. We recently reported the discovery of a new prokaryotic-specific GCYH-I (GCYH-IB) that displays no sequence identity to the canonical enzyme and is present in approximately 25% of bacteria, the majority of which lack the canonical GCYH-I (renamed GCYH-IA). Genomic and genetic analyses indicate that in those organisms possessing both enzymes, e.g., Bacillus subtilis, GCYH-IA and -IB are functionally redundant, but differentially expressed. Whereas GCYH-IA is constitutively expressed, GCYH-IB is expressed only under Zn(2+)-limiting conditions. These observations are consistent with the hypothesis that GCYH-IB functions to allow folate biosynthesis during Zn(2+) starvation. Here, we present biochemical and structural data showing that bacterial GCYH-IB, like GCYH-IA, belongs to the tunneling-fold (T-fold) superfamily. However, the GCYH-IA and -IB enzymes exhibit significant differences in global structure and active-site architecture. While GCYH-IA is a unimodular, homodecameric, Zn(2+)-dependent enzyme, GCYH-IB is a bimodular, homotetrameric enzyme activated by a variety of divalent cations. The structure of GCYH-IB and the broad metal dependence exhibited by this enzyme further underscore the mechanistic plasticity that is emerging for the T-fold superfamily. Notably, while humans possess the canonical GCYH-IA enzyme, many clinically important human pathogens possess only the GCYH-IB enzyme, suggesting that this enzyme is a potential new molecular target for antibacterial development.

RegPrecise: a database of curated genomic inferences of transcriptional regulatory interactions in prokaryotes

The RegPrecise database (http://regprecise.lbl.gov) was developed for capturing, visualization and analysis of predicted transcription factor regulons in prokaryotes that were reconstructed and manually curated by utilizing the comparative genomic approach. A significant number of high-quality inferences of transcriptional regulatory interactions have been already accumulated for diverse taxonomic groups of bacteria. The reconstructed regulons include transcription factors, their cognate DNA motifs and regulated genes/operons linked to the candidate transcription factor binding sites. The RegPrecise allows for browsing the regulon collections for: (i) conservation of DNA binding sites and regulated genes for a particular regulon across diverse taxonomic lineages; (ii) sets of regulons for a family of transcription factors; (iii) repertoire of regulons in a particular taxonomic group of species; (iv) regulons associated with a metabolic pathway or a biological process in various genomes. The initial release of the database includes ∼11 500 candidate binding sites for ∼400 orthologous groups of transcription factors from over 350 prokaryotic genomes. Majority of these data are represented by genome-wide regulon reconstructions in Shewanella and Streptococcus genera and a large-scale prediction of regulons for the LacI family of transcription factors. Another section in the database represents the results of accurate regulon propagation to the closely related genomes.

A subset of the diverse COG0523 family of putative metal chaperones is linked to zinc homeostasis in all kingdoms of life

Background

COG0523 proteins are, like the nickel chaperones of the UreG family, part of the G3E family of GTPases linking them to metallocenter biosynthesis. Even though the first COG0523-encoding gene, cobW, was identified almost 20 years ago, little is known concerning the function of other members belonging to this ubiquitous family.

Results

Based on a combination of comparative genomics, literature and phylogenetic analyses and experimental validations, the COG0523 family can be separated into at least fifteen subgroups. The CobW subgroup involved in cobalamin synthesis represents only one small sub-fraction of the family. Another, larger subgroup, is suggested to play a predominant role in the response to zinc limitation based on the presence of the corresponding COG0523-encoding genes downstream from putative Zur binding sites in many bacterial genomes. Zur binding sites in these genomes are also associated with candidate zinc-independent paralogs of zinc-dependent enzymes. Finally, the potential role of COG0523 in zinc homeostasis is not limited to Bacteria. We have predicted a link between COG0523 and regulation by zinc in Archaea and show that two COG0523 genes are induced upon zinc depletion in a eukaryotic reference organism, Chlamydomonas reinhardtii.

Conclusion

This work lays the foundation for the pursuit by experimental methods of the specific role of COG0523 members in metal trafficking. Based on phylogeny and comparative genomics, both the metal specificity and the protein target(s) might vary from one COG0523 subgroup to another. Additionally, Zur-dependent expression of COG0523 and putative paralogs of zinc-dependent proteins may represent a mechanism for hierarchal zinc distribution and zinc sparing in the face of inadequate zinc nutrition.

Mechanisms and evolution of control logic in prokaryotic transcriptional regulation

A major part of organismal complexity and versatility of prokaryotes resides in their ability to fine-tune gene expression to adequately respond to internal and external stimuli. Evolution has been very innovative in creating intricate mechanisms by which different regulatory signals operate and interact at promoters to drive gene expression. The regulation of target gene expression by transcription factors (TFs) is governed by control logic brought about by the interaction of regulators with TF binding sites (TFBSs) in cis-regulatory regions. A factor that in large part determines the strength of the response of a target to a given TF is motif stringency, the extent to which the TFBS fits the optimal TFBS sequence for a given TF. Advances in high-throughput technologies and computational genomics allow reconstruction of transcriptional regulatory networks in silico. To optimize the prediction of transcriptional regulatory networks, i.e., to separate direct regulation from indirect regulation, a thorough understanding of the control logic underlying the regulation of gene expression is required. This review summarizes the state of the art of the elements that determine the functionality of TFBSs by focusing on the molecular biological mechanisms and evolutionary origins of cis-regulatory regions.

Positive and negative feedback regulatory loops of thiol-oxidative stress response mediated by an unstable isoform of sigmaR in actinomycetes

Alternate sigma factors provide an effective way of diversifying bacterial gene expression in response to environmental changes. In Streptomyces coelicolor where more than 65 sigma factors are predicted, sigma(R) is the major regulator for response to thiol-oxidative stresses. sigma(R) becomes available when its bound anti-sigma factor RsrA is oxidized at sensitive cysteine thiols to form disulphide bonds. sigma(R) regulon includes genes for itself and multiple thiol-reducing systems, which constitute positive and negative feedback loops respectively. We found that the positive amplification loop involves an isoform of sigma(R) (sigma(R')) with an N-terminal extension of 55 amino acids, produced from an upstream start codon. A major difference between constitutive sigma(R) and inducible sigma(R') is that the latter is markedly unstable (t(1/2) approximately 10 min) compared with the former (> 70 min). The rapid turnover of sigma(R') is partly due to induced ClpP1/P2 proteases from the sigma(R) regulon. This represents a novel way of elaborating positive and negative feedback loops in a control circuit. Similar phenomenon may occur in other actinomycetes that harbour multiple start codons in the sigR homologous gene. We observed that sigH gene, the sigR orthologue in Mycobacterium smegmatis, produces an unstable larger isoform of sigma(H) upon induction by thiol-oxidative stress.

Fur controls iron homeostasis and oxidative stress defense in the oligotrophic alpha-proteobacterium Caulobacter crescentus

In most bacteria, the ferric uptake regulator (Fur) is a global regulator that controls iron homeostasis and other cellular processes, such as oxidative stress defense. In this work, we apply a combination of bioinformatics, in vitro and in vivo assays to identify the Caulobacter crescentus Fur regulon. A C. crescentus fur deletion mutant showed a slow growth phenotype, and was hypersensitive to H(2)O(2) and organic peroxide. Using a position weight matrix approach, several predicted Fur-binding sites were detected in the genome of C. crescentus, located in regulatory regions of genes not only involved in iron uptake and usage but also in other functions. Selected Fur-binding sites were validated using electrophoretic mobility shift assay and DNAse I footprinting analysis. Gene expression assays revealed that genes involved in iron uptake were repressed by iron-Fur and induced under conditions of iron limitation, whereas genes encoding iron-using proteins were activated by Fur under conditions of iron sufficiency. Furthermore, several genes that are regulated via small RNAs in other bacteria were found to be directly regulated by Fur in C. crescentus. In conclusion, Fur functions as an activator and as a repressor, integrating iron metabolism and oxidative stress response in C. crescentus.