Transcriptional regulation of the carbohydrate utilization network in Thermotoga maritima

Sugar transport in thermophiles: Bridging lignocellulose deconstruction and bioconversion

Biomass degrading thermophiles play an indispensable role in building lignocellulose-based supply chains. They operate at high temperatures to improve process efficiencies and minimize mesophilic contamination, can overcome lignocellulose recalcitrance through their native carbohydrate-active enzyme (CAZyme) inventory, and can utilize a wide range of sugar substrates. However, sugar transport in thermophiles is poorly understood and investigated, as compared to enzymatic lignocellulose deconstruction and metabolic conversion of sugars to value-added chemicals. Here, we review the general modes of sugar transport in thermophilic bacteria and archaea, covering the structural, molecular, and biophysical basis of their high-affinity sugar uptake. We also discuss recent genetic studies on sugar transporter function. With this understanding of sugar transport, we discuss strategies for how sugar transport can be engineered in thermophiles, with the potential to enhance the conversion of lignocellulosic biomass into renewable products.

ONE-SENTENCE SUMMARY

Sugar transport is the understudied link between extracellular biomass deconstruction and intracellular sugar metabolism in thermophilic lignocellulose bioprocessing.

A Cytoplasmic NAD(P)H-Dependent Polysulfide Reductase with Thiosulfate Reductase Activity from the Hyperthermophilic Bacterium Thermotoga maritima

Thermotoga maritima is an anaerobic hyperthermophilic bacterium that efficiently produces H2 by fermenting carbohydrates. High concentration of H2 inhibits the growth of T. maritima, and S0 could eliminate the inhibition and stimulate the growth through its reduction. The mechanism of T. maritima sulfur reduction, however, has not been fully understood. Herein, based on its similarity with archaeal NAD(P)H-dependent sulfur reductases (NSR), the ORF THEMA_RS02810 was identified and expressed in Escherichia coli, and the recombinant protein was characterized. The purified flavoprotein possessed NAD(P)H-dependent S0 reductase activity (1.3 U/mg for NADH and 0.8 U/mg for NADPH), polysulfide reductase activity (0.32 U/mg for NADH and 0.35 U/mg for NADPH), and thiosulfate reductase activity (2.3 U/mg for NADH and 2.5 U/mg for NADPH), which increased 3~4-folds by coenzyme A stimulation. Quantitative RT-PCR analysis showed that nsr was upregulated together with the mbx, yeeE, and rnf genes when the strain grew in S0- or thiosulfate-containing medium. The mechanism for sulfur reduction in T. maritima was discussed, which may affect the redox balance and energy metabolism of T. maritima. Genome search revealed that NSR homolog is widely distributed in thermophilic bacteria and archaea, implying its important role in the sulfur cycle of geothermal environments. IMPORTANCE The reduction of S0 and thiosulfate is essential in the sulfur cycle of geothermal environments, in which thermophiles play an important role. Despite previous research on some sulfur reductases of thermophilic archaea, the mechanism of sulfur reduction in thermophilic bacteria is still not clearly understood. Herein, we confirmed the presence of a cytoplasmic NAD(P)H-dependent polysulfide reductase (NSR) from the hyperthermophile T. maritima, with S0, polysulfide, and thiosulfate reduction activities, in contrast to other sulfur reductases. When grown in S0- or thiosulfate-containing medium, its expression was upregulated. And the putative membrane-bound MBX and Rnf may also play a role in the metabolism, which might influence the redox balance and energy metabolism of T. maritima. This is distinct from the mechanism of sulfur reduction in mesophiles such as Wolinella succinogenes. NSR homologs are widely distributed among heterotrophic thermophiles, suggesting that they may be vital in the sulfur cycle in geothermal environments.

Genus Thermotoga: A valuable home of multifunctional glycoside hydrolases (GHs) for industrial sustainability

Nature is a dexterous and prolific chemist for cataloging a number of hostile niches that are the ideal residence of various thermophiles. Apart from having other species, these subsurface environments are considered a throne of bacterial genus Thermotoga. The genome sequence of Thermotogales encodes complex and incongruent clusters of glycoside hydrolases (GHs), which are superior to their mesophilic counterparts and play a prominent role in various applications due to their extreme intrinsic stability. They have a tremendous capacity to use a wide variety of simple and multifaceted carbohydrates through GHs, formulate fermentative hydrogen and bioethanol at extraordinary yield, and catalyze high-temperature reactions for various biotechnological applications. Nevertheless, no stringent rules exist for the thermo-stabilization of biocatalysts present in the genus Thermotoga. These enzymes endure immense attraction in fundamental aspects of how these polypeptides attain and stabilize their distinctive three-dimensional (3D) structures to accomplish their physiological roles. Moreover, numerous genome sequences from Thermotoga species have revealed a significant fraction of genes most closely related to those of archaeal species, thus firming a staunch belief of lateral gene transfer mechanism. However, the question of its magnitude is still in its infancy. In addition to GHs, this genus is a paragon of encapsulins which carry pharmacological and industrial significance in the field of life sciences. This review highlights an intricate balance between the genomic organizations, factors inducing the thermostability, and pharmacological and industrial applications of GHs isolated from genus Thermotoga.

Novel functional insights into a modified sugar-binding protein from Synechococcus MITS9220

Paradigms of metabolic strategies employed by photoautotrophic marine picocyanobacteria have been challenged in recent years. Based on genomic annotations, picocyanobacteria are predicted to assimilate organic nutrients via ATP-binding cassette importers, a process mediated by substrate-binding proteins. We report the functional characterisation of a modified sugar-binding protein, MsBP, from a marine Synechococcus strain, MITS9220. Ligand screening of MsBP shows a specific affinity for zinc (K_D ~ 1.3 μM) and a preference for phosphate-modified sugars, such as fructose-1,6-biphosphate, in the presence of zinc (K_D ~ 5.8 μM). Our crystal structures of apo MsBP (no zinc or substrate-bound) and Zn-MsBP (with zinc-bound) show that the presence of zinc induces structural differences, leading to a partially-closed substrate-binding cavity. The Zn-MsBP structure also sequesters several sulphate ions from the crystallisation condition, including two in the binding cleft, appropriately placed to mimic the orientation of adducts of a biphosphate hexose. Combined with a previously unseen positively charged binding cleft in our two structures and our binding affinity data, these observations highlight novel molecular variations on the sugar-binding SBP scaffold. Our findings lend further evidence to a proposed sugar acquisition mechanism in picocyanobacteria alluding to a mixotrophic strategy within these ubiquitous photosynthetic bacteria.

Characterisation of a hyperthermophilic transketolase from Thermotoga maritima DSM3109 as a biocatalyst for 7-keto-octuronic acid synthesis

Transketolase (TK) is a fundamentally important enzyme in industrial biocatalysis which carries out a stereospecific carbon-carbon bond formation, and is widely used in the synthesis of prochiral ketones. This study describes the biochemical and molecular characterisation of a novel and unusual hyperthermophilic TK from Thermotoga maritima DSM3109 (TKtmar). TKtmar has a low protein sequence homology compared to the already described TKs, with key amino acid residues in the active site highly conserved. TKtmar has a very high optimum temperature (>90 °C) and shows pronounced stability at high temperature (e.g. t1/2 99 and 9.3 h at 50 and 80 °C, respectively) and in presence of organic solvents commonly used in industry (DMSO, acetonitrile and methanol). Substrate screening showed activity towards several monosaccharides and aliphatic aldehydes. In addition, for the first time, TK specificity towards uronic acids was achieved with TKtmar catalysing the efficient conversion of d-galacturonic acid and lithium hydroxypyruvate into 7-keto-octuronic acid, a very rare C8 uronic acid, in high yields (98%, 49 mM).

Transcriptional Regulation of Plant Biomass Degradation and Carbohydrate Utilization Genes in the Extreme Thermophile Caldicellulosiruptor bescii

Extremely thermophilic bacteria from the genus Caldicellulosiruptor can degrade polysaccharide components of plant cell walls and subsequently utilize the constituting mono- and oligosaccharides. Through metabolic engineering, ethanol and other industrially important end products can be produced. Previous experimental studies identified a variety of carbohydrate-active enzymes in model species Caldicellulosiruptor saccharolyticus and Caldicellulosiruptor bescii, while prior transcriptomic experiments identified their putative carbohydrate uptake transporters. We investigated the mechanisms of transcriptional regulation of carbohydrate utilization genes using a comparative genomics approach applied to 14 Caldicellulosiruptor species. The reconstruction of carbohydrate utilization regulatory network includes the predicted binding sites for 34 mostly local regulators and point to the regulatory mechanisms controlling expression of genes involved in degradation of plant biomass. The Rex and CggR regulons control the central glycolytic and primary redox reactions. The identified transcription factor binding sites and regulons were validated with transcriptomic and transcription start site experimental data for C. bescii grown on cellulose, cellobiose, glucose, xylan, and xylose. The XylR and XynR regulons control xylan-induced transcriptional response of genes involved in degradation of xylan and xylose utilization. The reconstructed regulons informed the carbohydrate utilization reconstruction analysis and improved functional annotations of 51 transporters and 11 catabolic enzymes. Using gene deletion, we confirmed that the shared ATPase component MsmK is essential for growth on oligo- and polysaccharides but not for the utilization of monosaccharides. By elucidating the carbohydrate utilization framework in C. bescii, strategies for metabolic engineering can be pursued to optimize yields of bio-based fuels and chemicals from lignocellulose. IMPORTANCE To develop functional metabolic engineering platforms for nonmodel microorganisms, a comprehensive understanding of the physiological and metabolic characteristics is critical. Caldicellulosiruptor bescii and other species in this genus have untapped potential for conversion of unpretreated plant biomass into industrial fuels and chemicals. The highly interactive and complex machinery used by C. bescii to acquire and process complex carbohydrates contained in lignocellulose was elucidated here to complement related efforts to develop a metabolic engineering platform with this bacterium. Guided by the findings here, a clearer picture of how C. bescii natively drives carbohydrate utilization is provided and strategies to engineer this bacterium for optimal conversion of lignocellulose to commercial products emerge.

The role of solute binding proteins in signal transduction

The solute binding proteins (SBPs) of prokaryotes are present in the extracytosolic space. Although their primary function is providing substrates to transporters, SBPs also stimulate different signaling proteins, including chemoreceptors, sensor kinases, diguanylate cyclases/phosphodiesterases and Ser/Thr kinases, thereby causing a wide range of responses. While relatively few such systems have been identified, several pieces of evidence suggest that SBP-mediated receptor activation is a widespread mechanism. (1) These systems have been identified in Gram-positive and Gram-negative bacteria and archaea. (2) There is a structural diversity in the receptor domains that bind SBPs. (3) SBPs belonging to thirteen different families interact with receptor ligand binding domains (LBDs). (4) For the two most abundant receptor LBD families, dCache and four-helix-bundle, there are different modes of interaction with SBPs. (5) SBP-stimulated receptors carry out many different functions. The advantage of SBP-mediated receptor stimulation is attributed to a strict control of SBP levels, which allows a precise adjustment of the systeḿs sensitivity. We have compiled information on the effect of ligands on the transcript/protein levels of their cognate SBPs. In 87 % of the cases analysed, ligands altered SBP expression levels. The nature of the regulatory effect depended on the ligand family. Whereas inorganic ligands typically downregulate SBP expression, an upregulation was observed in response to most sugars and organic acids. A major unknown is the role that SBPs play in signaling and in receptor stimulation. This review attempts to summarize what is known and to present new information to narrow this gap in knowledge.

Metabolic Phenotypes as Potential Biomarkers for Linking Gut Microbiome With Inflammatory Bowel Diseases

The gut microbiome is of utmost importance to human health. While a healthy microbiome can be represented by a variety of structures, its functional capacity appears to be more important. Gene content of the community can be assessed by “shotgun” metagenomics, but this approach is still too expensive. High-throughput amplicon-based surveys are a method of choice for large-scale surveys of links between microbiome, diseases, and diet, but the algorithms for predicting functional composition need to be improved to achieve good precision. Here we show how feature engineering based on microbial phenotypes, an advanced method for functional prediction from 16S rRNA sequencing data, improves identification of alterations of the gut microbiome linked to the disease. We processed a large collection of published gut microbial datasets of inflammatory bowel disease (IBD) patients to derive their community phenotype indices (CPI)—high-precision semiquantitative profiles aggregating metabolic potential of the community members based on genome-wide metabolic reconstructions. The list of selected metabolic functions included metabolism of short-chain fatty acids, vitamins, and carbohydrates. The machine-learning approach based on microbial phenotypes allows us to distinguish the microbiome profiles of healthy controls from patients with Crohn's disease and from ones with ulcerative colitis. The classifiers were comparable in quality to conventional taxonomy-based classifiers but provided new findings giving insights into possible mechanisms of pathogenesis. Feature-wise partial dependence plot (PDP) analysis of contribution to the classification result revealed a diversity of patterns. These observations suggest a constructive basis for defining functional homeostasis of the healthy human gut microbiome. The developed features are promising interpretable candidate biomarkers for assessing microbiome contribution to disease risk for the purposes of personalized medicine and clinical trials.

Effect of Cultivation Parameters on Fermentation and Hydrogen Production in the Phylum Thermotogae

The phylum Thermotogae is composed of a single class (Thermotogae), 4 orders (Thermotogales, Kosmotogales, Petrotogales, Mesoaciditogales), 5 families (Thermatogaceae, Fervidobacteriaceae, Kosmotogaceae, Petrotogaceae, Mesoaciditogaceae), and 13 genera. They have been isolated from extremely hot environments whose characteristics are reflected in the metabolic and phenotypic properties of the Thermotogae species. The metabolic versatility of Thermotogae members leads to a pool of high value-added products with application potentials in many industry fields. The low risk of contamination associated with their extreme culture conditions has made most species of the phylum attractive candidates in biotechnological processes. Almost all members of the phylum, especially those in the order Thermotogales, can produce bio-hydrogen from a variety of simple and complex sugars with yields close to the theoretical Thauer limit of 4 mol H₂/mol consumed glucose. Acetate, lactate, and L-alanine are the major organic end products. Thermotagae fermentation processes are influenced by various factors, such as hydrogen partial pressure, agitation, gas sparging, culture/headspace ratio, inoculum, pH, temperature, nitrogen sources, sulfur sources, inorganic compounds, metal ions, etc. Optimization of these parameters will help to fully unleash the biotechnological potentials of Thermotogae and promote their applications in industry. This article gives an overview of how these operational parameters could impact Thermotogae fermentation in terms of sugar consumption, hydrogen yields, and organic acids production.

Construction and Validation of a Genome-Scale Metabolic Network of Thermotoga sp. Strain RQ7

Thermotoga are anaerobic hyperthermophiles that have a deep lineage to the last universal ancestor and produce biological hydrogen gas accompanying cell growth. In recent years, systems-level approaches have been used to elucidate their metabolic capacities, by integrating mathematical modeling and experimental results. To assist biochemical engineering studies of T. sp. strain RQ7, this work aims at building a metabolic model of the bacterium that quantitatively simulates its metabolism at the genome scale. The constructed model, RQ7_iJG408, consists of 408 genes, 692 reactions, and 538 metabolites. Constraint-based flux balance analyses were used to simulate cell growth in both the complex and defined media. Quantitative comparison of the predicted and measured growth rates resulted in good agreements. This model serves as a foundation for an integrated biochemical description of T. sp. strain RQ7. It is a useful tool in designing growth media, identifying metabolic engineering strategies, and exploiting the physiological potentials of this biotechnologically significant organism.

A review of methods for the reconstruction and analysis of integrated genome-scale models of metabolism and regulation

The current survey aims to describe the main methodologies for extending the reconstruction and analysis of genome-scale metabolic models and phenotype simulation with Flux Balance Analysis mathematical frameworks, via the integration of Transcriptional Regulatory Networks and/or gene expression data. Although the surveyed methods are aimed at improving phenotype simulations obtained from these models, the perspective of reconstructing integrated genome-scale models of metabolism and gene expression for diverse prokaryotes is still an open challenge.

Gene targets for engineering osmotolerance in Caldicellulosiruptor bescii

BACKGROUND

Caldicellulosiruptor bescii, a promising biocatalyst being developed for use in consolidated bioprocessing of lignocellulosic materials to ethanol, grows poorly and has reduced conversion at elevated medium osmolarities. Increasing tolerance to elevated fermentation osmolarities is desired to enable performance necessary of a consolidated bioprocessing (CBP) biocatalyst.

RESULTS

Two strains of C. bescii showing growth phenotypes in elevated osmolarity conditions were identified. The first strain, ORCB001, carried a deletion of the FapR fatty acid biosynthesis and malonyl-CoA metabolism repressor and had a severe growth defect when grown in high-osmolarity conditions-introduced as the addition of either ethanol, NaCl, glycerol, or glucose to growth media. The second strain, ORCB002, displayed a growth rate over three times higher than its genetic parent when grown in high-osmolarity medium. Unexpectedly, a genetic complement ORCB002 exhibited improved growth, failing to revert the observed phenotype, and suggesting that mutations other than the deleted transcription factor (the fruR/cra gene) are responsible for the growth phenotype observed in ORCB002. Genome resequencing identified several other genomic alterations (three deleted regions, three substitution mutations, one silent mutation, and one frameshift mutation), which may be responsible for the observed increase in osmolarity tolerance in the fruR/cra-deficient strain, including a substitution mutation in dnaK, a gene previously implicated in osmoresistance in bacteria. Differential expression analysis and transcription factor binding site inference indicates that FapR negatively regulates malonyl-CoA and fatty acid biosynthesis, as it does in many other bacteria. FruR/Cra regulates neighboring fructose metabolism genes, as well as other genes in global manner.

CONCLUSIONS

Two systems able to effect tolerance to elevated osmolarities in C. bescii are identified. The first is fatty acid biosynthesis. The other is likely the result of one or more unintended, secondary mutations present in another transcription factor deletion strain. Though the locus/loci and mechanism(s) responsible remain unknown, candidate mutations are identified, including a mutation in the dnaK chaperone coding sequence. These results illustrate both the promise of targeted regulatory manipulation for osmotolerance (in the case of fapR) and the challenges (in the case of fruR/cra).

CO2-Induced Transcriptional Reorganization: Molecular Basis of Capnophillic Lactic Fermentation in Thermotoga neapolitana

Capnophilic lactic fermentation (CLF) is a novel anaplerotic pathway able to convert sugars to lactic acid (LA) and hydrogen using CO₂ as carbon enhancer in the hyperthermophilic bacterium Thermotoga neapolitana. In order to give further insights into CLF metabolic networks, we investigated the transcriptional modification induced by CO₂ using a RNA-seq approach. Transcriptomic analysis revealed 1601 differentially expressed genes (DEGs) in an enriched CO₂ atmosphere over a total of 1938 genes of the T. neapolitana genome. Transcription of PFOR and LDH genes belonging to the CLF pathway was up-regulated by CO₂ together with 6-phosphogluconolactonase (6PGL) and 6-phosphogluconate dehydratase (EDD) of the Entner–Doudoroff (ED) pathway. The transcriptomic study also revealed up-regulation of genes coding for the flavin-based enzymes NADH-dependent reduced ferredoxin:NADP oxidoreductase (NFN) and NAD-ferredoxin oxidoreductase (RNF) that control supply of reduced ferredoxin and NADH and allow energy conservation-based sodium translocation through the cell membrane. These results support the hypothesis that CO₂ induces rearrangement of the central carbon metabolism together with activation of mechanisms that increase availability of the reducing equivalents that are necessary to sustain CLF. In this view, this study reports a first rationale of the molecular basis of CLF in T. neapolitana and provides a list of target genes for the biotechnological implementation of this process.

Biochemical characterisation of four rhamnosidases from thermophilic bacteria of the genera Thermotoga, Caldicellulosiruptor and Thermoclostridium

Carbohydrate active enzymes are classified in databases based on sequence and structural similarity. However, their function can vary considerably within a similarity-based enzyme family, which makes biochemical characterisation indispensable to unravel their physiological role and to arrive at a meaningful annotation of the corresponding genes. In this study, we biochemically characterised the four related enzymes Tm_Ram106B, Tn_Ram106B, Cb_Ram106B and Ts_Ram106B from the thermophilic bacteria Thermotoga maritima MSB8, Thermotoga neapolitana Z2706-MC24, Caldicellulosiruptor bescii DSM 6725 and Thermoclostridium stercorarium DSM 8532, respectively, as α-l-rhamnosidases. Cobalt, nickel, manganese and magnesium ions stimulated while EDTA and EGTA inhibited all four enzymes. The kinetic parameters such as K_m, V_max and k_cat were about average compared to other rhamnosidases. The enzymes were inhibited by rhamnose, with half-maximal inhibitory concentrations (IC₅₀) between 5 mM and 8 mM. The α-l-rhamnosidases removed the terminal rhamnose moiety from the rutinoside in naringin, a natural flavonone glycoside. The Thermotoga sp. enzymes displayed the highest optimum temperatures and thermostabilities of all rhamnosidases reported to date. The four thermophilic and divalent ion-dependent rhamnosidases are the first biochemically characterised orthologous enzymes recently assigned to glycoside hydrolase family 106.

Novel Metabolic Pathways and Regulons for Hexuronate Utilization in Proteobacteria

We used comparative genomics to reconstruct d-galacturonic and d-glucuronic acid catabolic pathways and associated transcriptional regulons involving the tripartite ATP-independent periplasmic (TRAP) family transporters that bind hexuronates in proteobacteria. The reconstructed catabolic network involves novel transcription factors, catabolic enzymes, and transporters for utilization of both hexuronates and aldarates (d-glucarate and meso-galactarate). The reconstructed regulons for a novel GntR family transcription factor, GguR, include the majority of hexuronate/aldarate utilization genes in 47 species from the Burkholderiaceae, Comamonadaceae, Halomonadaceae, and Pseudomonadaceae families. GudR, GulR, and UdhR are additional local regulators of some hexuronate/aldarate utilization genes in some of the above-mentioned organisms. The predicted DNA binding motifs of GguR and GudR regulators from Ralstonia pickettii and Polaromonas were validated by in vitro binding assays. Genes from the GulR- and GguR-controlled loci were differentially expressed in R. pickettii grown on hexuronates and aldarates. By a combination of bioinformatics and experimental techniques we identified a novel variant of the oxidative pathway for hexuronate utilization, including two previously uncharacterized subfamilies of lactone hydrolases (UxuL and UxuF). The genomic context of respective genes and reconstruction of associated pathways suggest that both enzymes catalyze the conversion of d-galactaro- and d-glucaro-1,5-lactones to the ring-opened aldarates. The activities of the purified recombinant enzymes, UxuL and UxuF, from four proteobacterial species were directly confirmed and kinetically characterized. The inferred novel aldarate-specific transporter from the tripartite tricarboxylate transporter (TTT) family transporter TctC was confirmed to bind d-glucarate in vitro This study expands our knowledge of bacterial carbohydrate catabolic pathways by identifying novel families of catabolic enzymes, transcriptional regulators, and transporters.IMPORTANCE Hexuronate catabolic pathways and their transcriptional networks are highly variable among different bacteria. We identified novel transcriptional regulators that control the hexuronate and aldarate utilization genes in four families of proteobacteria. By regulon reconstruction and genome context analysis we identified several novel components of the common hexuronate/aldarate utilization pathways, including novel uptake transporters and catabolic enzymes. Two novel families of lactonases involved in the oxidative pathway of hexuronate catabolism were characterized. Novel transcriptional regulons were validated via in vitro binding assays and gene expression studies with Polaromonas and Ralstonia species. The reconstructed catabolic pathways are interconnected with each other metabolically and coregulated via the GguR regulons in proteobacteria.

Comparative Genomics Reveals the Regulatory Complexity of Bifidobacterial Arabinose and Arabino-Oligosaccharide Utilization

Members of the genus Bifidobacterium are common inhabitants of the human gastrointestinal tract. Previously it was shown that arabino-oligosaccharides (AOS) might act as prebiotics and stimulate the bifidobacterial growth in the gut. However, despite the rapid accumulation of genomic data, the precise mechanisms by which these sugars are utilized and associated transcription control still remain unclear. In the current study, we used a comparative genomic approach to reconstruct arabinose and AOS utilization pathways in over 40 bacterial species belonging to the Bifidobacteriaceae family. The results indicate that the gene repertoire involved in the catabolism of these sugars is highly diverse, and even phylogenetically close species may differ in their utilization capabilities. Using bioinformatics analysis we identified potential DNA-binding motifs and reconstructed putative regulons for the arabinose and AOS utilization genes in the Bifidobacteriaceae genomes. Six LacI-family transcriptional factors (named AbfR, AauR, AauU1, AauU2, BauR1 and BauR2) and a TetR-family regulator (XsaR) presumably act as local repressors for AOS utilization genes encoding various α- or β-L-arabinofuranosidases and predicted AOS transporters. The ROK-family regulator AraU and the LacI-family regulator AraQ control adjacent operons encoding putative arabinose transporters and catabolic enzymes, respectively. However, the AraQ regulator is universally present in all Bifidobacterium species including those lacking the arabinose catabolic genes araBDA, suggesting its control of other genes. Comparative genomic analyses of prospective AraQ-binding sites allowed the reconstruction of AraQ regulons and a proposed binary repression/activation mechanism. The conserved core of reconstructed AraQ regulons in bifidobacteria includes araBDA, as well as genes from the central glycolytic and fermentation pathways (pyk, eno, gap, tkt, tal, galM, ldh). The current study expands the range of genes involved in bifidobacterial arabinose/AOS utilization and demonstrates considerable variations in associated metabolic pathways and regulons. Detailed comparative and phylogenetic analyses allowed us to hypothesize how the identified reconstructed regulons evolved in bifidobacteria. Our findings may help to improve carbohydrate catabolic phenotype prediction and metabolic modeling, while it may also facilitate rational development of novel prebiotics.

A novel bifunctional transcriptional regulator of riboflavin metabolism in Archaea

Riboflavin (vitamin B2) is the precursor of flavin mononucleotide (FMN) and flavin adenine dinucleotide, which are essential coenzymes in all free-living organisms. Riboflavin biosynthesis in many Bacteria but not in Archaea is controlled by FMN-responsive riboswitches. We identified a novel bifunctional riboflavin kinase/regulator (RbkR), which controls riboflavin biosynthesis and transport genes in major lineages of Crenarchaeota, Euryarchaeota and Thaumarchaeota. RbkR proteins are composed of the riboflavin kinase domain and a DNA-binding winged helix-turn-helix-like domain. Using comparative genomics, we predicted RbkR operator sites and reconstructed RbkR regulons in 94 archaeal genomes. While the identified RbkR operators showed significant variability between archaeal lineages, the conserved core of RbkR regulons includes riboflavin biosynthesis genes, known/predicted vitamin uptake transporters and the rbkR gene. The DNA motifs and CTP-dependent riboflavin kinase activity of two RbkR proteins were experimentally validated in vitro. The DNA binding activity of RbkR was stimulated by CTP and suppressed by FMN, a product of riboflavin kinase. The crystallographic structure of RbkR from Thermoplasma acidophilum was determined in complex with CTP and its DNA operator revealing key residues for operator and ligand recognition. Overall, this study contributes to our understanding of metabolic and regulatory networks for vitamin homeostasis in Archaea.

Comparative Genomic Analysis of the Human Gut Microbiome Reveals a Broad Distribution of Metabolic Pathways for the Degradation of Host-Synthetized Mucin Glycans and Utilization of Mucin-Derived Monosaccharides

The colonic mucus layer is a dynamic and complex structure formed by secreted and transmembrane mucins, which are high-molecular-weight and heavily glycosylated proteins. Colonic mucus consists of a loose outer layer and a dense epithelium-attached layer. The outer layer is inhabited by various representatives of the human gut microbiota (HGM). Glycans of the colonic mucus can be used by the HGM as a source of carbon and energy when dietary fibers are not sufficiently available. Both commensals and pathogens can utilize mucin glycans. Commensals are mostly involved in the cleavage of glycans, while pathogens mostly utilize monosaccharides released by commensals. This HGM-derived degradation of the mucus layer increases pathogen susceptibility and causes many other health disorders. Here, we analyzed 397 individual HGM genomes to identify pathways for the cleavage of host-synthetized mucin glycans to monosaccharides as well as for the catabolism of the derived monosaccharides. Our key results are as follows: (i) Genes for the cleavage of mucin glycans were found in 86% of the analyzed genomes, which significantly higher than a previous estimation. (ii) Genes for the catabolism of derived monosaccharides were found in 89% of the analyzed genomes. (iii) Comparative genomic analysis identified four alternative forms of the monosaccharide-catabolizing enzymes and four alternative forms of monosaccharide transporters. (iv) Eighty-five percent of the analyzed genomes may be involved in potential feeding pathways for the monosaccharides derived from cleaved mucin glycans. (v) The analyzed genomes demonstrated different abilities to degrade known mucin glycans. Generally, the ability to degrade at least one type of mucin glycan was predicted for 81% of the analyzed genomes. (vi) Eighty-two percent of the analyzed genomes can form mutualistic pairs that are able to degrade mucin glycans and are not degradable by any of the paired organisms alone. Taken together, these findings provide further insight into the inter-microbial communications of the HGM as well as into host-HGM interactions.

Electrochemically applied potentials induce growth and metabolic shift changes in the hyperthermophilic bacterium Thermotoga maritima MSB8

Bioelectrochemical systems are an attractive technology for regulating microbial activity. The effect of an applied potential on hydrolysis of starch in Thermotoga maritima as a model bacterium was investigated in this study. A cathodic potential (−0.6 and −0.8 V) induced 5-h earlier growth initiation of T. maritima with starch as the polymeric substrate than that without electrochemical regulation. Moreover, metabolic patterns of starch consumption were altered by the cathodic potential. While acetate, H2, and CO2 were the major products of starch consumption in the control experiment without electrolysis, lactate accumulation was detected rather than decreased acetate and H2 levels in the bioelectrochemical system experiments with the cathodic potential. These results indicate that the applied potential could control microbial activities related to the hydrolysis of polymeric organic substances and shift carbon and electron flux to a lactate-producing reaction in T. maritima.

Genomic Reconstruction of Carbohydrate Utilization Capacities in Microbial-Mat Derived Consortia

Two nearly identical unicyanobacterial consortia (UCC) were previously isolated from benthic microbial mats that occur in a heliothermal saline lake in northern Washington State. Carbohydrates are a primary source of carbon and energy for most heterotrophic bacteria. Since CO₂ is the only carbon source provided, the cyanobacterium must provide a source of carbon to the heterotrophs. Available genomic sequences for all members of the UCC provide opportunity to investigate the metabolic routes of carbon transfer between autotroph and heterotrophs. Here, we applied a subsystem-based comparative genomics approach to reconstruct carbohydrate utilization pathways and identify glycohydrolytic enzymes, carbohydrate transporters and pathway-specific transcriptional regulators in 17 heterotrophic members of the UCC. The reconstructed metabolic pathways include 800 genes, near a one-fourth of which encode enzymes, transporters and regulators with newly assigned metabolic functions resulting in discovery of novel functional variants of carbohydrate utilization pathways. The in silico analysis revealed the utilization capabilities for 40 carbohydrates and their derivatives. Two Halomonas species demonstrated the largest number of sugar catabolic pathways. Trehalose, sucrose, maltose, glucose, and beta-glucosides are the most commonly utilized saccharides in this community. Reconstructed regulons for global regulators HexR and CceR include central carbohydrate metabolism genes in the members of Gammaproteobacteria and Alphaproteobacteria, respectively. Genomics analyses were supplemented by experimental characterization of metabolic phenotypes in four isolates derived from the consortia. Measurements of isolate growth on the defined medium supplied with individual carbohydrates confirmed most of the predicted catabolic phenotypes. Not all consortia members use carbohydrates and only a few use complex polysaccharides suggesting a hierarchical carbon flow from cyanobacteria to each heterotroph. In summary, the genomics-based identification of carbohydrate utilization capabilities provides a basis for future experimental studies of carbon flow in UCC.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Comparative genomics and evolution of transcriptional regulons in Proteobacteria

Comparative genomics approaches are broadly used for analysis of transcriptional regulation in bacterial genomes. In this work, we identified binding sites and reconstructed regulons for 33 orthologous groups of transcription factors (TFs) in 196 reference genomes from 21 taxonomic groups of Proteobacteria. Overall, we predict over 10 600 TF binding sites and identified more than 15 600 target genes for 1896 TFs constituting the studied orthologous groups of regulators. These include a set of orthologues for 21 metabolism-associated TFs from Escherichia coli and/or Shewanella that are conserved in five or more taxonomic groups and several additional TFs that represent non-orthologous substitutions of the metabolic regulators in some lineages of Proteobacteria. By comparing gene contents of the reconstructed regulons, we identified the core, taxonomy-specific and genome-specific TF regulon members and classified them by their metabolic functions. Detailed analysis of ArgR, TyrR, TrpR, HutC, HypR and other amino-acid-specific regulons demonstrated remarkable differences in regulatory strategies used by various lineages of Proteobacteria. The obtained genomic collection of in silico reconstructed TF regulons contains a large number of new regulatory interactions that await future experimental validation. The collection provides a framework for future evolutionary studies of transcriptional regulatory networks in Bacteria. It can be also used for functional annotation of putative metabolic transporters and enzymes that are abundant in the reconstructed regulons.

Genomic characterization of Defluviitoga tunisiensis L3, a key hydrolytic bacterium in a thermophilic biogas plant and its abundance as determined by metagenome fragment recruitment

The genome sequence of Defluviitoga tunisiensis L3 originating from a thermophilic biogas-production plant was established and recently published as Genome Announcement by our group. The circular chromosome of D. tunisiensis L3 has a size of 2,053,097bp and a mean GC content of 31.38%. To analyze the D. tunisiensis L3 genome sequence in more detail, a phylogenetic analysis of completely sequenced Thermotogae strains based on shared core genes was performed. It appeared that Petrotoga mobilis DSM 10674(T), originally isolated from a North Sea oil-production well, is the closest relative of D. tunisiensis L3. Comparative genome analyses of P. mobilis DSM 10674(T) and D. tunisiensis L3 showed moderate similarities regarding occurrence of orthologous genes. Both genomes share a common set of 1351 core genes. Reconstruction of metabolic pathways important for the biogas production process revealed that the D. tunisiensis L3 genome encodes a large set of genes predicted to facilitate utilization of a variety of complex polysaccharides including cellulose, chitin and xylan. Ethanol, acetate, hydrogen (H2) and carbon dioxide (CO2) were found as possible end-products of the fermentation process. The latter three metabolites are considered to represent substrates for methanogenic Archaea, the key organisms in the final step of the anaerobic digestion process. To determine the degree of relatedness between D. tunisiensis L3 and dominant biogas community members within the thermophilic biogas-production plant, metagenome sequences obtained from the corresponding microbial community were mapped onto the L3 genome sequence. This fragment recruitment revealed that the D. tunisiensis L3 genome is almost completely covered with metagenome sequences featuring high matching accuracy. This result indicates that strains highly related or even identical to the reference strain D. tunisiensis L3 play a dominant role within the community of the thermophilic biogas-production plant.

Gut bacteria that prevent growth impairments transmitted by microbiota from malnourished children

Undernourished children exhibit impaired development of their gut microbiota. Transplanting microbiota from 6- and 18-month-old healthy or undernourished Malawian donors into young germ-free mice that were fed a Malawian diet revealed that immature microbiota from undernourished infants and children transmit impaired growth phenotypes. The representation of several age-discriminatory taxa in recipient animals correlated with lean body mass gain; liver, muscle, and brain metabolism; and bone morphology. Mice were cohoused shortly after receiving microbiota from healthy or severely stunted and underweight infants; age- and growth-discriminatory taxa from the microbiota of the former were able to invade that of the latter, which prevented growth impairments in recipient animals. Adding two invasive species, Ruminococcus gnavus and Clostridium symbiosum, to the microbiota from undernourished donors also ameliorated growth and metabolic abnormalities in recipient animals. These results provide evidence that microbiota immaturity is causally related to undernutrition and reveal potential therapeutic targets and agents.

Transcriptional Regulation of Carbohydrate Utilization Pathways in the Bifidobacterium Genus

Bifidobacteria, which represent common commensals of mammalian gut, are believed to have positive effects on human health. The influence of certain non-digestible carbohydrates (and their use as so-called prebiotics) on growth and metabolic activity of bifidobacteria is of increasing interest; however, mechanisms of transcriptional control of carbohydrate metabolism are poorly understood in these species. We used a comparative genomics approach to reconstruct carbohydrate utilization pathways and transcriptional regulons in 10 Bifidobacterium genomes. Analysis of regulatory gene regions revealed candidate DNA motifs and reconstructed regulons for 268 transcription factors from the LacI, ROK, DeoR, AraC, GntR, and TetR families that form 64 orthologous groups of regulators. Most of the reconstructed regulons are local and control specific catabolic pathways for host- and diet-derived glycans and monosaccharides. Mosaic distributions of many of these local regulators across Bifidobacterium species correlate with distribution of corresponding catabolic pathways. In contrast, the maltose, galactose, sucrose, and fructose regulons, as well as a novel global LacI-family regulator that is predicted to control the central carbohydrate metabolism and arabinose catabolism genes, are universally present in all 10 studied bifidobacteria. A novel group of TetR-family regulators presumably controls the glucoside and galactoside utilization pathways. Paralogs of the ribose repressor RbsR control the pyrimidine nucleoside utilization genes. Multiple paralogs of the maltose regulator MalR co-regulate large sets of genes involved in maltodextrin utilization. The inferred metabolic regulons provide new insights on diverse carbohydrate utilization networks in bifidobacteria that can be employed in metabolic modeling, phenotype prediction and the rational development of novel prebiotics.

Comparative genomics of pyridoxal 5'-phosphate-dependent transcription factor regulons in Bacteria

The MocR-subfamily transcription factors (MocR-TFs) characterized by the GntR-family DNA-binding domain and aminotransferase-like sensory domain are broadly distributed among certain lineages of Bacteria. Characterized MocR-TFs bind pyridoxal 5'-phosphate (PLP) and control transcription of genes involved in PLP, gamma aminobutyric acid (GABA) and taurine metabolism via binding specific DNA operator sites. To identify putative target genes and DNA binding motifs of MocR-TFs, we performed comparative genomics analysis of over 250 bacterial genomes. The reconstructed regulons for 825 MocR-TFs comprise structural genes from over 200 protein families involved in diverse biological processes. Using the genome context and metabolic subsystem analysis we tentatively assigned functional roles for 38 out of 86 orthologous groups of studied regulators. Most of these MocR-TF regulons are involved in PLP metabolism, as well as utilization of GABA, taurine and ectoine. The remaining studied MocR-TF regulators presumably control genes encoding enzymes involved in reduction/oxidation processes, various transporters and PLP-dependent enzymes, for example aminotransferases. Predicted DNA binding motifs of MocR-TFs are generally similar in each orthologous group and are characterized by two to four repeated sequences. Identified motifs were classified according to their structures. Motifs with direct and/or inverted repeat symmetry constitute the majority of inferred DNA motifs, suggesting preferable TF dimerization in head-to-tail or head-to-head configuration. The obtained genomic collection of in silico reconstructed MocR-TF motifs and regulons in Bacteria provides a basis for future experimental characterization of molecular mechanisms for various regulators in this family.

Regulation of the rhaEWRBMA Operon Involved in l-Rhamnose Catabolism through Two Transcriptional Factors, RhaR and CcpA, in Bacillus subtilis

The Bacillus subtilis rhaEWRBMA (formerly yuxG-yulBCDE) operon consists of four genes encoding enzymes for l-rhamnose catabolism and the rhaR gene encoding a DeoR-type transcriptional regulator. DNase I footprinting analysis showed that the RhaR protein specifically binds to the regulatory region upstream of the rhaEW gene, in which two imperfect direct repeats are included. Gel retardation analysis revealed that the direct repeat farther upstream is essential for the high-affinity binding of RhaR and that the DNA binding of RhaR was effectively inhibited by L-rhamnulose-1-phosphate, an intermediate of L-rhamnose catabolism. Moreover, it was demonstrated that the CcpA/P-Ser-HPr complex, primarily governing the carbon catabolite control in B. subtilis, binds to the catabolite-responsive element, which overlaps the RhaR binding site. In vivo analysis of the rhaEW promoter-lacZ fusion in the background of ccpA deletion showed that the L-rhamnose-responsive induction of the rhaEW promoter was negated by the disruption of rhaA or rhaB but not rhaEW or rhaM, whereas rhaR disruption resulted in constitutive rhaEW promoter activity. These in vitro and in vivo results clearly indicate that RhaR represses the operon by binding to the operator site, which is detached by L-rhamnulose-1-phosphate formed from L-rhamnose through a sequence of isomerization by RhaA and phosphorylation by RhaB, leading to the derepression of the operon. In addition, the lacZ reporter analysis using the strains with or without the ccpA deletion under the background of rhaR disruption supported the involvement of CcpA in the carbon catabolite repression of the operon.

IMPORTANCE

Since L-rhamnose is a component of various plant-derived compounds, it is a potential carbon source for plant-associating bacteria. Moreover, it is suggested that L-rhamnose catabolism plays a significant role in some bacteria-plant interactions, e.g., invasion of plant pathogens and nodulation of rhizobia. Despite the physiological importance of L-rhamnose catabolism for various bacterial species, the transcriptional regulation of the relevant genes has been poorly understood, except for the regulatory system of Escherichia coli. In this study, we show that, in Bacillus subtilis, one of the plant growth-promoting rhizobacteria, the rhaEWRBMA operon for L-rhamnose catabolism is controlled by RhaR and CcpA. This regulatory system can be another standard model for better understanding the regulatory mechanisms of L-rhamnose catabolism in other bacterial species.

Insights into thermoadaptation and the evolution of mesophily from the bacterial phylum Thermotogae

Thermophiles are extremophiles that grow optimally at temperatures >45 °C. To survive and maintain function of their biological molecules, they have a suite of characteristics not found in organisms that grow at moderate temperature (mesophiles). At the cellular level, thermophiles have mechanisms for maintaining their membranes, nucleic acids, and other cellular structures. At the protein level, each of their proteins remains stable and retains activity at temperatures that would denature their mesophilic homologs. Conversely, cellular structures and proteins from thermophiles may not function optimally at moderate temperatures. These differences between thermophiles and mesophiles presumably present a barrier for evolutionary transitioning between the 2 lifestyles. Therefore, studying closely related thermophiles and mesophiles can help us determine how such lifestyle transitions may happen. The bacterial phylum Thermotogae contains hyperthermophiles, thermophiles, mesophiles, and organisms with temperature ranges wide enough to span both thermophilic and mesophilic temperatures. Genomic, proteomic, and physiological differences noted between other bacterial thermophiles and mesophiles are evident within the Thermotogae. We argue that the Thermotogae is an ideal group of organisms for understanding of the response to fluctuating temperature and of long-term evolutionary adaptation to a different growth temperature range.

Adaptive Evolution of Thermotoga maritima Reveals Plasticity of the ABC Transporter Network

Thermotoga maritima is a hyperthermophilic anaerobe that utilizes a vast network of ABC transporters to efficiently metabolize a variety of carbon sources to produce hydrogen. For unknown reasons, this organism does not metabolize glucose as readily as it does glucose di- and polysaccharides. The leading hypothesis implicates the thermolability of glucose at the physiological temperatures at which T. maritima lives. After a 25-day laboratory evolution, phenotypes were observed with growth rates up to 1.4 times higher than and glucose utilization rates exceeding 50% those of the wild type. Genome resequencing revealed mutations in evolved cultures related to glucose-responsive ABC transporters. The native glucose ABC transporter, GluEFK, has more abundant transcripts either as a result of gene duplication-amplification or through mutations to the operator sequence regulating this operon. Conversely, BglEFGKL, a transporter of beta-glucosides, is substantially downregulated due to a nonsense mutation to the solute binding protein or due to a deletion of the upstream promoter. Analysis of the ABC2 uptake porter families for carbohydrate and peptide transport revealed that the solute binding protein, often among the transcripts detected at the highest levels, is predominantly downregulated in the evolved cultures, while the membrane-spanning domain and nucleotide binding components are less varied. Similar trends were observed in evolved strains grown on glycerol, a substrate that is not dependent on ABC transporters. Therefore, improved growth on glucose is achieved through mutations favoring GluEFK expression over BglEFGKL, and in lieu of carbon catabolite repression, the ABC transporter network is modulated to achieve improved growth fitness.

Every Site Counts: Submitting Transcription Factor-Binding Site Information through the CollecTF Portal

Experimentally verified transcription factor-binding sites represent an information-rich and highly applicable data type that aptly summarizes the results of time-consuming experiments and inference processes. Currently, there is no centralized repository for this type of data, which is routinely embedded in articles and extremely hard to mine. CollecTF provides the first standardized resource for submission and deposition of these data into the NCBI RefSeq database, maximizing its accessibility and prompting the community to adopt direct submission policies.

Prediction and characterization of enzymatic activities guided by sequence similarity and genome neighborhood networks

Metabolic pathways in eubacteria and archaea often are encoded by operons and/or gene clusters (genome neighborhoods) that provide important clues for assignment of both enzyme functions and metabolic pathways. We describe a bioinformatic approach (genome neighborhood network; GNN) that enables large scale prediction of the in vitro enzymatic activities and in vivo physiological functions (metabolic pathways) of uncharacterized enzymes in protein families. We demonstrate the utility of the GNN approach by predicting in vitro activities and in vivo functions in the proline racemase superfamily (PRS; InterPro IPR008794). The predictions were verified by measuring in vitro activities for 51 proteins in 12 families in the PRS that represent ∼85% of the sequences; in vitro activities of pathway enzymes, carbon/nitrogen source phenotypes, and/or transcriptomic studies confirmed the predicted pathways. The synergistic use of sequence similarity networks³ and GNNs will facilitate the discovery of the components of novel, uncharacterized metabolic pathways in sequenced genomes.

DOI:http://dx.doi.org/10.7554/eLife.03275.001

DNA molecules are polymers in which four nucleotides—guanine, adenine, thymine, and cytosine—are arranged along a sugar backbone. The sequence of these four nucleotides along the DNA strand determines the genetic code of the organism, and can be deciphered using various genome sequencing techniques. Microbial genomes are particularly easy to sequence as they contain fewer than several million nucleotides, compared with the 3 billion or so nucleotides that are present in the human genome.

Reading a genome sequence is straight forward, but predicting the physiological functions of the proteins encoded by the genes in the sequence can be challenging. In a process called genome annotation, the function of protein is predicted by comparing the relevant gene to the genes of proteins with known functions. However, microbial genomes and proteins are hugely diverse and over 50% of the microbial genomes that have been sequenced have not yet been related to any physiological function. With thousands of microbial genomes waiting to be deciphered, large scale approaches are needed.

Zhao et al. take advantage of a particular characteristic of microbial genomes. DNA sequences that code for two proteins required for the same task tend to be closer to each other in the genome than two sequences that code for unrelated functions. Operons are an extreme example; an operon is a unit of DNA that contains several genes that are expressed as proteins at the same time.

Zhao et al. have developed a bioinformatic method called the genome neighbourhood network approach to work out the function of proteins based on their position relative to other proteins in the genome. When applied to the proline racemase superfamily (PRS), which contains enzymes with similar sequences that can catalyze three distinct chemical reactions, the new approach was able to assign a function to the majority of proteins in a public database of PRS enzymes, and also revealed new members of the PRS family. Experiments confirmed that the proteins behaved as predicted. The next challenge is to develop the genome neighbourhood network approach so that it can be applied to more complex systems.

DOI:http://dx.doi.org/10.7554/eLife.03275.002

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

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

Comparative genomics and functional analysis of rhamnose catabolic pathways and regulons in bacteria

L-rhamnose (L-Rha) is a deoxy-hexose sugar commonly found in nature. L-Rha catabolic pathways were previously characterized in various bacteria including Escherichia coli. Nevertheless, homology searches failed to recognize all the genes for the complete L-Rha utilization pathways in diverse microbial species involved in biomass decomposition. Moreover, the regulatory mechanisms of L-Rha catabolism have remained unclear in most species. A comparative genomics approach was used to reconstruct the L-Rha catabolic pathways and transcriptional regulons in the phyla Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria, and Thermotogae. The reconstructed pathways include multiple novel enzymes and transporters involved in the utilization of L-Rha and L-Rha-containing polymers. Large-scale regulon inference using bioinformatics revealed remarkable variations in transcriptional regulators for L-Rha utilization genes among bacteria. A novel bifunctional enzyme, L-rhamnulose-phosphate aldolase (RhaE) fused to L-lactaldehyde dehydrogenase (RhaW), which is not homologous to previously characterized L-Rha catabolic enzymes, was identified in diverse bacteria including Chloroflexi, Bacilli, and Alphaproteobacteria. By using in vitro biochemical assays we validated both enzymatic activities of the purified recombinant RhaEW proteins from Chloroflexus aurantiacus and Bacillus subtilis. Another novel enzyme of the L-Rha catabolism, L-lactaldehyde reductase (RhaZ), was identified in Gammaproteobacteria and experimentally validated by in vitro enzymatic assays using the recombinant protein from Salmonella typhimurium. C. aurantiacus induced transcription of the predicted L-Rha utilization genes when L-Rha was present in the growth medium and consumed L-Rha from the medium. This study provided comprehensive insights to L-Rha catabolism and its regulation in diverse Bacteria.

Polysaccharides utilization in human gut bacterium Bacteroides thetaiotaomicron: comparative genomics reconstruction of metabolic and regulatory networks

Background

Bacteroides thetaiotaomicron, a predominant member of the human gut microbiota, is characterized by its ability to utilize a wide variety of polysaccharides using the extensive saccharolytic machinery that is controlled by an expanded repertoire of transcription factors (TFs). The availability of genomic sequences for multiple Bacteroides species opens an opportunity for their comparative analysis to enable characterization of their metabolic and regulatory networks.

Results

A comparative genomics approach was applied for the reconstruction and functional annotation of the carbohydrate utilization regulatory networks in 11 Bacteroides genomes. Bioinformatics analysis of promoter regions revealed putative DNA-binding motifs and regulons for 31 orthologous TFs in the Bacteroides. Among the analyzed TFs there are 4 SusR-like regulators, 16 AraC-like hybrid two-component systems (HTCSs), and 11 regulators from other families. Novel DNA motifs of HTCSs and SusR-like regulators in the Bacteroides have the common structure of direct repeats with a long spacer between two conserved sites.

Conclusions

The inferred regulatory network in B. thetaiotaomicron contains 308 genes encoding polysaccharide and sugar catabolic enzymes, carbohydrate-binding and transport systems, and TFs. The analyzed TFs control pathways for utilization of host and dietary glycans to monosaccharides and their further interconversions to intermediates of the central metabolism. The reconstructed regulatory network allowed us to suggest and refine specific functional assignments for sugar catabolic enzymes and transporters, providing a substantial improvement to the existing metabolic models for B. thetaiotaomicron. The obtained collection of reconstructed TF regulons is available in the RegPrecise database (http://regprecise.lbl.gov).

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

Background

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

Description

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

Conclusions

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