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

Crp and Arc system directly regulate the transcription of NADH dehydrogenase genes in Shewanella oneidensis nitrate and nitrite respiration

NADH oxidation by NADH dehydrogenases (NDHs) is crucial for feeding respiratory quinone pool and maintaining the balance of NADH/NAD+. In the respiratory model organism Shewanella oneidensis, which possesses four NDHs, the longstanding notion had been that NDHs were not required under anoxic conditions until recent studies demonstrated their role in extracellular electron transfer. However, the role of each NDH, particularly under anoxic conditions, has not been characterized. Here, we systematically investigated the role of each NDH in aerobic and anaerobic nitrate and nitrite respiration using NDH triple mutants. We corroborated the involvement of NDHs in anaerobic nitrate/nitrite respiration, revealing different repertoires of NDHs employed by S. oneidensis in response to electron acceptor (EA) availability. The transcript levels of two nqrs were modulated by the EA conversion from nitrate to nitrite. Furthermore, we demonstrated that the global regulators Crp and the Arc system both directly controlled the transcription of four NDHs during nitrate/nitrite respiration. This study confirms the requirement of NDHs for anaerobic nitrate and nitrite respiration and sheds light on the respiratory remodeling mechanism whereby global regulators coordinate NDH genes transcription to adapt to redox-stratified environments.IMPORTANCENADH is an important electron source for the respiratory quinone pool. Multiple NADH dehydrogenases (NDHs) are widely present in prokaryotes, indicating the flexibility in NADH oxidation. As a renowned respiratory versatile model strain, Shewanella oneidensis possesses four NDHs, encompassing all three types of NDHs, with varying ion-translocating efficiencies. The redundancy of NDHs may confer advantages for S. oneidensis to survive and thrive in redox-stratified environments. However, the roles of each NDH, especially in anaerobic respiration, are less understood. Here, we evaluated the role of each NDH in aerobic and anaerobic nitrate/nitrite respiration. We found that the conversion of electron acceptor from nitrate to nitrite triggered the changes in the transcriptional levels of NDH genes, and global regulators Crp and the Arc system were involved in these processes. These findings elucidate the mechanism of the respiratory chain remodeling at the NADH oxidation step in response to different electron acceptors.

Functional genomics of chitin degradation by Vibrio parahaemolyticus reveals finely integrated metabolic contributions to support environmental fitness

Vibrio species are marine prokaryotes that inhabit diverse ecological niches, colonizing abiotic and biotic surfaces. These bacteria are vital players in the global carbon cycle, assimilating billions of tonnes of chitin for carbon (and nitrogen) metabolites. Many bacterial proteins involved in the process-including chitinases, sugar transporters, and modifying enzymes-have been well studied. However, the genetic functional interplay and key drivers of Vibrio competitive survival in the presence of chitin as the dominant carbon source is not understood. To address this question, we carried out transposon sequencing (Tn-seq) to determine the genetic fitness of Vibrio parahaemolyticus mutants grown on chitin as a sole carbon source. Along with validating known Vibrio genes associated with chitin metabolism, our data newly identified vital roles for an unclassified OprD-like import chitoporin and a HexR family transcriptional regulator. Furthermore, we functionally implicated HexR in regulating multiple physiological processes involved in V. parahaemolyticus environmental survival including carbon assimilation and cell growth, biofilm formation, and cell motility. Under nutrient limiting conditions, our data revealed a requirement for HexR in filamentous cell morphology, a critical trait for V. parahaemolyticus environmental fitness. Therefore, a vital import porin and genomic regulation mediated by HexR support multiple physiological processes for Vibrio chitinolytic growth and environmental fitness.

Bioinformatics in Russia: history and present-day landscape

Bioinformatics has become an interdisciplinary subject due to its universal role in molecular biology research. The current status of Russia's bioinformatics research in Russia is not known. Here, we review the history of bioinformatics in Russia, present the current landscape, and highlight future directions and challenges. Bioinformatics research in Russia is driven by four major industries: information technology, pharmaceuticals, biotechnology, and agriculture. Over the past three decades, despite a delayed start, the field has gained momentum, especially in protein and nucleic acid research. Dedicated and shared centers for genomics, proteomics, and bioinformatics are active in different regions of Russia. Present-day bioinformatics in Russia is characterized by research issues related to genetics, metagenomics, OMICs, medical informatics, computational biology, environmental informatics, and structural bioinformatics. Notable developments are in the fields of software (tools, algorithms, and pipelines), use of high computation power (e.g. by the Siberian Supercomputer Center), and large-scale sequencing projects (the sequencing of 100 000 human genomes). Government funding is increasing, policies are being changed, and a National Genomic Information Database is being established. An increased focus on eukaryotic genome sequencing, the development of a common place for developers and researchers to share tools and data, and the use of biological modeling, machine learning, and biostatistics are key areas for future focus. Universities and research institutes have started to implement bioinformatics modules. A critical mass of bioinformaticians is essential to catch up with the global pace in the discipline.

Molecular mechanism of the one-component regulator RccR on bacterial metabolism and virulence

The regulation of carbon metabolism and virulence is critical for the rapid adaptation of pathogenic bacteria to host conditions. In Pseudomonas aeruginosa, RccR is a transcriptional regulator of genes involved in primary carbon metabolism and is associated with bacterial resistance and virulence, although the exact mechanism is unclear. Our study demonstrates that PaRccR is a direct repressor of the transcriptional regulator genes mvaU and algU. Biochemical and structural analyses reveal that PaRccR can switch its DNA recognition mode through conformational changes triggered by KDPG binding or release. Mutagenesis and functional analysis underscore the significance of allosteric communication between the SIS domain and the DBD domain. Our findings suggest that, despite its overall structural similarity to other bacterial RpiR-type regulators, RccR displays a more complex regulatory element binding mode induced by ligands and a unique regulatory mechanism.

Interplay between autotrophic and heterotrophic prokaryotic metabolism in the bathypelagic realm revealed by metatranscriptomic analyses

BACKGROUND

Heterotrophic microbes inhabiting the dark ocean largely depend on the settling of organic matter from the sunlit ocean. However, this sinking of organic materials is insufficient to cover their demand for energy and alternative sources such as chemoautotrophy have been proposed. Reduced sulfur compounds, such as thiosulfate, are a potential energy source for both auto- and heterotrophic marine prokaryotes.

METHODS

Seawater samples were collected from Labrador Sea Water (LSW, ~ 2000 m depth) in the North Atlantic and incubated in the dark at in situ temperature unamended, amended with 1 µM thiosulfate, or with 1 µM thiosulfate plus 10 µM glucose and 10 µM acetate (thiosulfate plus dissolved organic matter, DOM). Inorganic carbon fixation was measured in the different treatments and samples for metatranscriptomic analyses were collected after 1 h and 72 h of incubation.

RESULTS

Amendment of LSW with thiosulfate and thiosulfate plus DOM enhanced prokaryotic inorganic carbon fixation. The energy generated via chemoautotrophy and heterotrophy in the amended prokaryotic communities was used for the biosynthesis of glycogen and phospholipids as storage molecules. The addition of thiosulfate stimulated unclassified bacteria, sulfur-oxidizing Deltaproteobacteria (SAR324 cluster bacteria), Epsilonproteobacteria (Sulfurimonas sp.), and Gammaproteobacteria (SUP05 cluster bacteria), whereas, the amendment with thiosulfate plus DOM stimulated typically copiotrophic Gammaproteobacteria (closely related to Vibrio sp. and Pseudoalteromonas sp.).

CONCLUSIONS

The gene expression pattern of thiosulfate utilizing microbes specifically of genes involved in energy production via sulfur oxidation and coupled to CO2 fixation pathways coincided with the change in the transcriptional profile of the heterotrophic prokaryotic community (genes involved in promoting energy storage), suggesting a fine-tuned metabolic interplay between chemoautotrophic and heterotrophic microbes in the dark ocean. Video Abstract.

Glucose-6-Phosphate Dehydrogenase, ZwfA, a Dual Cofactor-Specific Isozyme Is Predominantly Involved in the Glucose Metabolism of Pseudomonas bharatica CSV86T

Glucose-6-phosphate dehydrogenase (Zwf) is an important enzyme in glucose metabolism via the Entner-Doudoroff pathway and the first enzyme in the oxidative pentose-phosphate pathway. It generates NAD(P)H during the conversion of glucose-6-phosphate (G6P) to 6-phosphogluconolactone, thus aiding in anabolic processes, energy yield, and oxidative stress responses. Pseudomonas bharatica CSV86T preferentially utilized aromatic compounds over glucose and exhibited a significantly lower growth rate on glucose (0.24 h-1) with a prolonged lag phase (~10 h). In strain CSV86T, glucose was metabolized via the intracellular phosphorylative route only because it lacked an oxidative (gluconate and 2-ketogluconate) route. The genome harbored three genes zwfA, zwfB, and zwfC encoding three Zwf isozymes. The present study aimed to understand gene arrangement, gene expression profiling, and molecular and kinetic properties of the purified enzymes to unveil their physiological significance in the strain CSV86T. The zwfA was found to be a part of the zwfA-pgl-eda operon, which was proximal to other glucose transport and metabolic clusters. The zwfB was found to be arranged as a gnd-zwfB operon, while zwfC was present independently. Among the three, zwfA was transcribed maximally, and the purified ZwfA displayed the highest catalytic efficiency, cooperativity with respect to G6P, and dual cofactor specificity. Isozymes ZwfB and ZwfC were NADP+-preferring and NADP+-specific, respectively. Among other functionally characterized Zwfs, ZwfA from strain CSV86T displayed poor catalytic efficiency and the further absence of oxidative routes of glucose metabolism reflected its lower growth rate on glucose compared to P. putida KT2440 and could be probable reasons for the unique carbon source utilization hierarchy. IMPORTANCE Pseudomonas bharatica CSV86T metabolizes glucose exclusively via the intracellular phosphorylative Entner-Doudoroff pathway leading the entire glucose flux through Zwf as the strain lacks oxidative routes. This may lead to limiting the concentration of downstream metabolic intermediates. The strain CSV86T possesses three isoforms of glucose-6-phosphate dehydrogenase, ZwfA, ZwfB, and ZwfC. The expression profile and kinetic properties of purified enzymes will help to understand glucose metabolism. Isozyme ZwfA dominated in terms of expression and displayed cooperativity with dual cofactor specificity. ZwfB preferred NADP+, and ZwfC was NADP+ specific, which may aid in redox cofactor balance. Such beneficial metabolic flexibility facilitated the regulation of metabolic pathways giving survival/fitness advantages in dynamic environments. Additionally, multiple genes allowed the distribution of function among these isoforms where the primary function was allocated to one of the isoforms.

Human Milk Oligosaccharide Utilization in Intestinal Bifidobacteria Is Governed by Global Transcriptional Regulator NagR

Bifidobacterium longum subsp. infantis is a prevalent beneficial bacterium that colonizes the human neonatal gut and is uniquely adapted to efficiently use human milk oligosaccharides (HMOs) as a carbon and energy source. Multiple studies have focused on characterizing the elements of HMO utilization machinery in B. longum subsp. infantis; however, the regulatory mechanisms governing the expression of these catabolic pathways remain poorly understood. A bioinformatic regulon reconstruction approach used in this study implicated NagR, a transcription factor from the ROK family, as a negative global regulator of gene clusters encoding lacto-N-biose/galacto-N-biose (LNB/GNB), lacto-N-tetraose (LNT), and lacto-N-neotetraose (LNnT) utilization pathways in B. longum subsp. infantis. This conjecture was corroborated by transcriptome profiling upon nagR genetic inactivation and experimental assessment of binding of recombinant NagR to predicted DNA operators. The latter approach also implicated N-acetylglucosamine (GlcNAc), a universal intermediate of LNT and LNnT catabolism, and its phosphorylated derivatives as plausible NagR transcriptional effectors. Reconstruction of NagR regulons in various Bifidobacterium lineages revealed multiple potential regulon expansion events, suggesting evolution from a local regulator of GlcNAc catabolism in ancestral bifidobacteria to a global regulator controlling the utilization of mixtures of GlcNAc-containing host glycans in B. longum subsp. infantis and Bifidobacterium bifidum. IMPORTANCE The predominance of bifidobacteria in the gut of breastfed infants is attributed to the ability of these bacteria to metabolize human milk oligosaccharides (HMOs). Thus, individual HMOs such as lacto-N-tetraose (LNT) and lacto-N-neotetraose (LNnT) are considered promising prebiotics that would stimulate the growth of bifidobacteria and confer multiple health benefits to preterm and malnourished children suffering from impaired (stunted) gut microbiota development. However, the rational selection of HMO-based prebiotics is hampered by the incomplete knowledge of regulatory mechanisms governing HMO utilization in target bifidobacteria. This study describes NagR-mediated transcriptional regulation of LNT and LNnT utilization in Bifidobacterium longum subsp. infantis. The elucidated regulatory network appears optimally adapted to simultaneous utilization of multiple HMOs, providing a rationale to add HMO mixtures (rather than individual components) to infant formulas. The study also provides insights into the evolutionary trajectories of complex regulatory networks controlling carbohydrate metabolism in bifidobacteria.

CRISPR/dCas9-RpoD-Mediated Simultaneous Transcriptional Activation and Repression in Shewanella oneidensis MR-1

Extracellular electron transfer (EET) of electroactive microorganisms (EAMs) is the dominating factor for versatile applications of bio-electrochemical systems. Shewanella oneidensis MR-1 is one of the model EAMs for the study of EET, which is associated with a variety of cellular activities. However, due to the lack of a transcriptional activation tool, regulation of multiple genes is labor-intensive and time-consuming, which hampers the advancement of improving the EET efficiency in S. oneidensis. In this study, we developed an easily operated and multifunctional regulatory tool, that is, a simultaneous clustered regularly interspaced short palindromic repeats (CRISPR)-mediated transcriptional activation (CRISPRa) and interference (CRISPRi) system, for application in S. oneidensis. First, a large number of activators were screened, and RpoD (σ⁷⁰) was determined as the optimal activator. Second, the effective activation range was identified to be 190-216 base upstream of the transcriptional start site. Third, up- and downregulation was achieved in concert by two orthogonal single guide RNAs targeting different positions. The activation of the cell division gene (minCDE) and repression of the cytotoxic gene (SO_3166) were concurrently implemented, increasing the power density by 2.5-fold and enhancing the degradation rate of azo dyes by 2.9-fold. The simultaneous CRISPRa and CRISPRi system enables simultaneous multiplex genetic regulation, offering the potential to further advance studies of the EET mechanism and application in S. oneidensis.

Differences in gene expression and endophytic bacterial diversity in Atractylodes macrocephala Koidz. rhizomes from different growth years

Atractylodes macrocephala Koidz. (AMK) is widely used owing to its pharmacological activity in traditional Chinese medicine (TCM). Here, we aimed to characterize the differentially expressed genes (DEGs) of one- and three-year growth (OYG and TYG) rhizomes of AMK combined with the endophytic bacterial diversity analysis using high-throughput RNA-sequencing. 114,572 unigenes were annotated in six public databases. 3570 DEGs revealed a clear difference, of which 936 and 2634 genes were up- and down-regulated, respectively. The results of KEGG pathway analysis indicated that DEGs corresponding to the terpenoid synthesis gene were downregulated in TYG rhizomes. 414,424 sequences corresponding to the 16S rRNA gene were divided into 1267 operational taxonomic units (OTUs). Moreover, the diversity of endophytic bacteria changed with species in OYG (773) and TYG (1201) rhizomes at OTU level, and Proteobacteria, Actinobacteria, and Bacteroidetes were the dominant phyla. Comparison of species differences among different growth years revealed that some species were significantly different, such as Actinomycetes, Variovorax, Cloacibacterium, etc. Interestingly, the decrease in the function-related metabolism of terpenoids and polyketides was found to be correlated the low expression of terpene synthesis genes in TYG rhizomes assessed using PICRUSt2. These data provide a scientific basis for elucidating the mechanism underlying metabolite accumulation and endophytic bacterial diversity in relation to the growth years in AMK.

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.

Engineering a Synthetic Pathway for Gentisate in Pseudomonas Chlororaphis P3

Pseudomonas chlororaphis P3 has been well-engineered as a platform organism for biologicals production due to enhanced shikimate pathway and excellent physiological and genetic characteristics. Gentisate displays high antiradical and antioxidant activities and is an important intermediate that can be used as a precursor for drugs. Herein, a plasmid-free biosynthetic pathway of gentisate was constructed by connecting the endogenous degradation pathway from 3-hydroxybenzoate in Pseudomonas for the first time. As a result, the production of gentisate reached 365 mg/L from 3-HBA via blocking gentisate conversion and enhancing the gentisate precursors supply through the overexpression of the rate-limiting step. With a close-up at the future perspectives, a series of bioactive compounds could be achieved by constructing synthetic pathways in conventional Pseudomonas to establish a cell factory.

Pectobacterium atrosepticum KDPG aldolase, Eda, participates in the Entner-Doudoroff pathway and independently inhibits expression of virulence determinants

Pectobacterium carotovorum has an incomplete Entner-Doudoroff (ED) pathway, including enzyme 2-keto-3-deoxy-6-phosphogluconate aldolase (Eda) but lacking phosphogluconate dehydratase (Edd), while P. atrosepticum (Pba) has a complete pathway. To understand the role of the ED pathway in Pectobacterium infection, mutants of these two key enzymes, Δeda and Δedd, were constructed in Pba SCRI1039. Δeda exhibited significant decreased virulence on potato tubers and colonization in planta and was greatly attenuated in pectinase activity and the ability to use pectin breakdown products, including polygalacturonic acid (PGA) and galacturonic acid. These reduced phenotypes were restored following complementation with an external vector expressing eda. Quantitative reverse transcription PCR analysis revealed that expression of the pectinase genes pelA, pelC, pehN, pelW, and pmeB in Δeda cultured in pyruvate, with or without PGA, was significantly reduced compared to the wild type, while genes for virulence regulators (kdgR, hexR, hexA, and rsmA) remained unchanged. However, Δedd showed similar phenotypes to the wild type. To our knowledge, this is the first demonstration that disruption of eda has a feedback effect on inhibiting pectin degradation and that Eda is involved in building the arsenal of pectinases needed during infection by Pectobacterium.

Lignin induced iron reduction by novel sp., Tolumonas lignolytic BRL6-1

Lignin is the second most abundant carbon polymer on earth and despite having more fuel value than cellulose, it currently is considered a waste byproduct in many industrial lignocellulose applications. Valorization of lignin relies on effective and green methods of de-lignification, with a growing interest in the use of microbes. Here we investigate the physiology and molecular response of the novel facultative anaerobic bacterium, Tolumonas lignolytica BRL6-1, to lignin under anoxic conditions. Physiological and biochemical changes were compared between cells grown anaerobically in either lignin-amended or unamended conditions. In the presence of lignin, BRL6-1 accumulates higher biomass and has a shorter lag phase compared to unamended conditions, and 14% of the proteins determined to be significantly higher in abundance by log2 fold-change of 2 or greater were related to Fe(II) transport in late logarithmic phase. Ferrozine assays of the supernatant confirmed that Fe(III) was bound to lignin and reduced to Fe(II) only in the presence of BRL6-1, suggesting redox activity by the cells. LC-MS/MS analysis of the secretome showed an extra band at 20 kDa in lignin-amended conditions. Protein sequencing of this band identified a protein of unknown function with homology to enzymes in the radical SAM superfamily. Expression of this protein in lignin-amended conditions suggests its role in radical formation. From our findings, we suggest that BRL6-1 is using a protein in the radical SAM superfamily to interact with the Fe(III) bound to lignin and reducing it to Fe(II) for cellular use, increasing BRL6-1 yield under lignin-amended conditions. This interaction potentially generates organic free radicals and causes a radical cascade which could modify and depolymerize lignin. Further research should clarify the extent to which this mechanism is similar to previously described aerobic chelator-mediated Fenton chemistry or radical producing lignolytic enzymes, such as lignin peroxidases, but under anoxic conditions.

Transcriptome of Pectobacterium carotovorum subsp. carotovorum PccS1 infected in calla plants in vivo highlights a spatiotemporal expression pattern of genes related to virulence, adaptation, and host response

Bacterial pathogens from the genus Pectobacterium cause soft rot in various plants, and result in important economic losses worldwide. We understand much about how these pathogens digest their hosts and protect themselves against plant defences, as well as some regulatory networks in these processes. However, the spatiotemporal expression of genome-wide infection of Pectobacterium remains unclear, although researchers analysed this in some phytopathogens. In the present work, comparing the transcriptome profiles from cellular infection with growth in minimal and rich media, RNA-Seq analyses revealed that the differentially expressed genes (log₂ -fold ratio ≥ 1.0) in the cells of Pectobacterium carotovorum subsp. carotovorum PccS1 recovered at a series of time points after inoculation in the host in vivo covered approximately 50% of genes in the genome. Based on the dynamic expression changes in infection, the significantly differentially expressed genes (log₂ -fold ratio ≥ 2.0) were classified into five types, and the main expression pattern of the genes for carbohydrate metabolism underlying the processes of infection was identified. The results are helpful to our understanding of the inducement of host plant and environmental adaption of Pectobacterium. In addition, our results demonstrate that maceration caused by PccS1 is due to the depression of callose deposition in the plant for resistance by the pathogenesis-related genes and the superlytic ability of pectinolytic enzymes produced in PccS1, rather than the promotion of plant cell death elicited by the T3SS of bacteria as described in previous work.

Transcriptional regulation of a gonococcal gene encoding a virulence factor (L-lactate permease)

GdhR is a GntR-type regulator of Neisseria gonorrhoeae encoded by a gene (gdhR) belonging to the MtrR regulon, which comprises multiple genes required for antibiotic resistance such as the mtrCDE efflux pump genes. In previous work we showed that loss of gdhR results in enhanced gonococcal fitness in a female mouse model of lower genital tract infection. Here, we used RNA-Seq to perform a transcriptional profiling study to determine the GdhR regulon. GdhR was found to regulate the expression of 2.3% of all the genes in gonococcal strain FA19, of which 39 were activated and 11 were repressed. Within the GdhR regulon we found that lctP, which encodes a unique L-lactate transporter and has been associated with gonococcal pathogenesis, was the highest of GdhR-repressed genes. By using in vitro transcription and DNase I footpriting assays we mapped the lctP transcriptional start site (TSS) and determined that GdhR directly inhibits transcription by binding to an inverted repeat sequence located 9 bases downstream of the lctP TSS. Epistasis analysis revealed that, while loss of lctP increased susceptibility of gonococci to hydrogen peroxide (H2O2) the loss of gdhR enhanced resistance; however, this GdhR-endowed property was reversed in a double gdhR lctP null mutant. We assessed the effect of different carbon sources on lctP expression and found that D-glucose, but not L-lactate or pyruvate, repressed lctP expression within a physiological concentration range but in a GdhR-independent manner. Moreover, we found that adding glucose to the medium enhanced susceptibility of gonococci to hydrogen peroxide. We propose a model for the role of lctP regulation via GdhR and glucose in the pathogenesis of N. gonorrhoeae.

The Expression of Virulence Factors in Vibrio anguillarum Is Dually Regulated by Iron Levels and Temperature

Vibrio anguillarum causes a hemorrhagic septicemia that affects cold- and warm-water adapted fish species. The main goal of this work was to determine the temperature-dependent changes in the virulence factors that could explain the virulence properties of V. anguillarum for fish cultivated at different temperatures. We have found that although the optimal growth temperature is around 25°C, the degree of virulence of V. anguillarum RV22 is higher at 15°C. To explain this result, an RNA-Seq analysis was performed to compare the whole transcriptome profile of V. anguillarum RV22 cultured under low-iron availability at either 25 or 15°C, which would mimic the conditions that V. anguillarum finds during colonization of fish cultivated at warm- or cold-water temperatures. The comparative analysis of transcriptomes at high- and low-iron conditions showed profound metabolic adaptations to grow under low iron. These changes were characterized by a down-regulation of the energetic metabolism and the induction of virulence-related factors like biosynthesis of LPS, production of hemolysins and lysozyme, membrane transport, heme uptake, or production of siderophores. However, the expression pattern of virulence factors under iron limitation showed interesting differences at warm and cold temperatures. Chemotaxis, motility, as well as the T6SS1 genes are expressed at higher levels at 25°C than at 15°C. By contrast, hemolysin RTX pore-forming toxin, T6SS2, and the genes associated with exopolysaccharides synthesis were preferentially expressed at 15°C. Notably, at this temperature, the siderophore piscibactin system was strongly up-regulated. In contrast, at 25°C, piscibactin genes were down-regulated and the vanchrobactin siderophore system seems to supply all the necessary iron to the cell. The results showed that V. anguillarum adjusts the expression of virulence factors responding to two environmental signals, iron levels and temperature. Thus, the relative relevance of each virulence factor for each fish species could vary depending on the water temperature. The results give clues about the physiological adaptations that allow V. anguillarum to cause infections in different fishes and could be relevant for vaccine development against fish vibriosis.

GlaR (YugA)-a novel RpiR-family transcription activator of the Leloir pathway of galactose utilization in Lactococcus lactis IL1403

Bacteria can utilize diverse sugars as carbon and energy source, but the regulatory mechanisms directing the choice of the preferred substrate are often poorly understood. Here, we analyzed the role of the YugA protein (now designated GlaR—Galactose–lactose operon Regulatory protein) of the RpiR family as a transcriptional activator of galactose (gal genes) and lactose (lac genes) utilization genes in Lactococcus lactis IL1403. In this bacterium, gal genes forming the Leloir operon are combined with lac genes in a single so‐called gal–lac operon. The first gene of this operon is the lacS gene encoding galactose permease. The glaR gene encoding GlaR lies directly upstream of the gal–lac gene cluster and is transcribed in the same direction. This genetic layout and the presence of glaR homologues in the closest neighborhood to the Leloir or gal–lac operons are highly conserved only among Lactococcus species. Deletion of glaR disabled galactose utilization and abrogated or decreased expression of the gal–lac genes. The GlaR‐dependent regulation of the gal–lac operon depends on its specific binding to a DNA region upstream of the lacS gene activating lacS expression and increasing the expression of the operon genes localized downstream. Notably, expression of lacS‐downstream genes, namely galMKTE, thgA and lacZ, is partially independent of the GlaR‐driven activation likely due to the presence of additional promoters. The glaR transcription itself is not subject to catabolite control protein A (CcpA) carbon catabolite repression (CRR) and is induced by galactose. Up to date, no similar mechanism has been reported in other lactic acid bacteria species. These results reveal a novel regulatory protein and shed new light on the regulation of carbohydrate catabolism in L. lactis IL1403, and by similarity, probably also in other lactococci.

Shewanella oneidensis MR-1 Utilizes both Sodium- and Proton-Pumping NADH Dehydrogenases during Aerobic Growth

Shewanella oneidensis MR-1 is a metal-reducing bacterium with the ability to utilize many different terminal electron acceptors, including oxygen and solid-metal oxides. Both metal oxide reduction and aerobic respiration have been studied extensively in this organism. However, electron transport chain processes upstream of the terminal oxidoreductases have been relatively understudied in this organism, especially electron transfer from NADH to respiratory quinones. Genome annotation indicates that S. oneidensis MR-1 encodes four NADH dehydrogenases, a proton-translocating dehydrogenase (Nuo), two sodium ion-translocating dehydrogenases (Nqr1 and Nqr2), and an “uncoupling” dehydrogenase (Ndh), but none of these complexes have been studied. Therefore, we conducted a study specifically focused on the effects of individual NADH dehydrogenase knockouts in S. oneidensis MR-1. We observed that two of the single-mutant strains, the ΔnuoN and ΔnqrF1 mutants, exhibited significant growth defects compared with the wild type. However, the defects were minor and only apparent under certain growth conditions. Further testing of the ΔnuoN ΔnqrF1 double-mutant strain yielded no growth in minimal medium under oxic conditions, indicating that Nuo and Nqr1 have overlapping functions, but at least one is necessary for aerobic growth. Coutilization of proton- and sodium ion-dependent energetics has important implications for the growth of this organism in environments with varied pH and salinity, including microbial electrochemical systems.

IMPORTANCE Bacteria utilize a wide variety of metabolic pathways that allow them to take advantage of different energy sources, and to do so with varied efficiency. The efficiency of a metabolic process determines the growth yield of an organism, or the amount of biomass it produces per amount of substrate consumed. This parameter has important implications in biotechnology and wastewater treatment, where low growth yields are often preferred to minimize the production of microbial biomass. In this study, we investigated respiratory pathways containing NADH dehydrogenases with varied efficiency (i.e., the number of ions translocated per NADH oxidized) in the metal-reducing bacterium Shewanella oneidensis MR-1. We observed that two different respiratory pathways are used concurrently, and at least one pathway must be functional for growth under oxic conditions.

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.

The effect of global transcriptional regulators on the anaerobic fermentative metabolism of Escherichia coli

Global transcription factors are known to regulate the anaerobic growth of Escherichia coli on glucose. These transcription factors help the organism to sense oxygen and accordingly regulate the synthesis of mixed acid producing enzymes. Five global transcription factors, namely ArcA, Fnr, IhfA-B, Crp and Fis, are known to play an important role in the growth phenotype of the organism in the transition from anaerobic to aerobic conditions. The effect of deletion of most of these global transcription factors on the growth phenotype has not been characterized under strict anaerobic fermentation conditions. In order to enumerate the role of global transcription factors in central carbon metabolism, experiments were performed using single deletion mutants of the above mentioned global transcription regulators. The mutants demonstrated lower growth rates, ranging from 3-75% lower growth as compared to the wild-type strain along with varying glucose uptake rates. Global transcription regulators help in lowering formate and acetate synthesis, thereby effectively channeling the carbon towards redox balance (through ethanol formation) and biomass synthesis. Flux analysis of mutant strains indicated that deletion of a single transcription factor alone does not play a significant role in the normalized flux distribution of the central carbon metabolism.

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.

One ligand, two regulators and three binding sites: How KDPG controls primary carbon metabolism in Pseudomonas

Effective regulation of primary carbon metabolism is critically important for bacteria to successfully adapt to different environments. We have identified an uncharacterised transcriptional regulator; RccR, that controls this process in response to carbon source availability. Disruption of rccR in the plant-associated microbe Pseudomonas fluorescens inhibits growth in defined media, and compromises its ability to colonise the wheat rhizosphere. Structurally, RccR is almost identical to the Entner-Doudoroff (ED) pathway regulator HexR, and both proteins are controlled by the same ED-intermediate; 2-keto-3-deoxy-6-phosphogluconate (KDPG). Despite these similarities, HexR and RccR control entirely different aspects of primary metabolism, with RccR regulating pyruvate metabolism (aceEF), the glyoxylate shunt (aceA, glcB, pntAA) and gluconeogenesis (pckA, gap). RccR displays complex and unusual regulatory behaviour; switching repression between the pyruvate metabolism and glyoxylate shunt/gluconeogenesis loci depending on the available carbon source. This regulatory complexity is enabled by two distinct pseudo-palindromic binding sites, differing only in the length of their linker regions, with KDPG binding increasing affinity for the 28 bp aceA binding site but decreasing affinity for the 15 bp aceE site. Thus, RccR is able to simultaneously suppress and activate gene expression in response to carbon source availability. Together, the RccR and HexR regulators enable the rapid coordination of multiple aspects of primary carbon metabolism, in response to levels of a single key intermediate.

Here we show how Pseudomonas controls multiple different primary carbon metabolism pathways by sensing levels of KDPG, an Entner Doudoroff (ED) pathway intermediate. KDPG binds to two highly similar transcription factors; the ED regulator HexR and the previously uncharacterised protein RccR. RccR inversely controls the glyoxylate shunt, gluconeogenesis and pyruvate metabolism, suppressing the first two pathways as pyruvate metabolism genes are expressed, and vice versa. This complex regulation is enabled by two distinct RccR-binding consensus sequences in the RccR regulon promoters. KDPG binding simultaneously increases RccR affinity for the glyoxylate shunt and gluconeogenesis promoters, and releases repression of pyruvate metabolism. This elegant two-regulator circuit allows Pseudomonas to rapidly respond to carbon source availability by sensing a single key intermediate, KDPG.

Transport and Catabolism of Carbohydrates by

We identified the genes encoding the proteins for the transport of glucose and maltose in <i>Neisseria meningitidis</i> strain 2C4-3. A mutant deleted for <i>NMV_1892</i><i>(glcP)</i> no longer grew on glucose and deletion of <i>NMV_0424</i><i>(malY)</i> prevented the utilization of maltose. We also purified and characterized glucokinase and α-phosphoglucomutase, which catalyze early catabolic steps of the two carbohydrates. <i>N. meningitidis</i> catabolizes the two carbohydrates either via the Entner-Doudoroff (ED) pathway or the pentose phosphate pathway, thereby forming glyceraldehyde-3-P and either pyruvate or fructose-6-P, respectively. We purified and characterized several key enzymes of the two pathways. The genes required for the transformation of glucose into gluconate-6-P and its further catabolism via the ED pathway are organized in two adjacent operons. <i>N. meningitidis</i> also contains genes encoding proteins which exhibit similarity to the gluconate transporter <i>(NMV_2230)</i> and gluconate kinase <i>(NMV_2231)</i> of Enterobacteriaceae and Firmicutes. However, gluconate might not be the real substrate of <i>NMV_2230</i> because <i>N. meningitidi</i>s was not able to grow on gluconate as the sole carbon source. Surprisingly, deletion of <i>NMV_2230</i> stimulated growth in minimal medium in the presence and absence of glucose and drastically slowed the clearance of <i>N. meningitidis</i> cells from transgenic mice after intraperitoneal challenge.

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.

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.

HexR Controls Glucose-Responsive Genes and Central Carbon Metabolism in Neisseria meningitidis

Neisseria meningitidis, an exclusively human pathogen and the leading cause of bacterial meningitis, must adapt to different host niches during human infection. N. meningitidis can utilize a restricted range of carbon sources, including lactate, glucose, and pyruvate, whose concentrations vary in host niches. Microarray analysis of N. meningitidis grown in a chemically defined medium in the presence or absence of glucose allowed us to identify genes regulated by carbon source availability. Most such genes are implicated in energy metabolism and transport, and some are implicated in virulence. In particular, genes involved in glucose catabolism were upregulated, whereas genes involved in the tricarboxylic acid cycle were downregulated. Several genes encoding surface-exposed proteins, including the MafA adhesins and Neisseria surface protein A, were upregulated in the presence of glucose. Our microarray analysis led to the identification of a glucose-responsive hexR-like transcriptional regulator that controls genes of the central carbon metabolism of N. meningitidis in response to glucose. We characterized the HexR regulon and showed that the hexR gene is accountable for some of the glucose-responsive regulation; in vitro assays with the purified protein showed that HexR binds to the promoters of the central metabolic operons of the bacterium. Based on DNA sequence alignment of the target sites, we propose a 17-bp pseudopalindromic consensus HexR binding motif. Furthermore, N. meningitidis strains lacking hexR expression were deficient in establishing successful bacteremia in an infant rat model of infection, indicating the importance of this regulator for the survival of this pathogen in vivo.

IMPORTANCE

Neisseria meningitidis grows on a limited range of nutrients during infection. We analyzed the gene expression of N. meningitidis in response to glucose, the main energy source available in human blood, and we found that glucose regulates many genes implicated in energy metabolism and nutrient transport, as well as some implicated in virulence. We identified and characterized a transcriptional regulator (HexR) that controls metabolic genes of N. meningitidis in response to glucose. We generated a mutant lacking HexR and found that the mutant was impaired in causing systemic infection in animal models. Since N. meningitidis lacks known bacterial regulators of energy metabolism, our findings suggest that HexR plays a major role in its biology by regulating metabolism in response to environmental signals.

Novel Functions and Regulation of Cryptic Cellobiose Operons in Escherichia coli

Presence of cellobiose as a sole carbon source induces mutations in the chb and asc operons of Escherichia coli and allows it to grow on cellobiose. We previously engineered these two operons with synthetic constitutive promoters and achieved efficient cellobiose metabolism through adaptive evolution. In this study, we characterized two mutations observed in the efficient cellobiose metabolizing strain: duplication of RBS of ascB gene, (β-glucosidase of asc operon) and nonsense mutation in yebK, (an uncharacterized transcription factor). Mutations in yebK play a dominant role by modulating the length of lag phase, relative to the growth rate of the strain when transferred from a rich medium to minimal cellobiose medium. Mutations in ascB, on the other hand, are specific for cellobiose and help in enhancing the specific growth rate. Taken together, our results show that ascB of the asc operon is controlled by an internal putative promoter in addition to the native cryptic promoter, and the transcription factor yebK helps to remodel the host physiology for cellobiose metabolism. While previous studies characterized the stress-induced mutations that allowed growth on cellobiose, here, we characterize the adaptation-induced mutations that help in enhancing cellobiose metabolic ability. This study will shed new light on the regulatory changes and factors that are needed for the functional coupling of the host physiology to the activated cryptic cellobiose metabolism.

Global transcriptome analysis reveals small RNAs affecting Neisseria meningitidis bacteremia

Most bacterial small RNAs (sRNAs) are post-transcriptional regulators involved in adaptive responses, controlling gene expression by modulating translation or stability of their target mRNAs often in concert with the RNA chaperone Hfq. Neisseria meningitides, the leading cause of bacterial meningitis, is able to adapt to different host niches during human infection. However, only a few sRNAs and their functions have been fully described to date. Recently, transcriptional expression profiling of N. meningitides in human blood ex vivo revealed 91 differentially expressed putative sRNAs. Here we expanded this analysis by performing a global transcriptome study after exposure of N. meningitides to physiologically relevant stress signals (e.g. heat shock, oxidative stress, iron and carbon source limitation). and we identified putative sRNAs that were differentially expressed in vitro. A set of 98 putative sRNAs was obtained by analyzing transcriptome data and 8 new sRNAs were validated, both by Northern blot and by primer extension techniques. Deletion of selected sRNAs caused attenuation of N. meningitides infection in the in vivo infant rat model, leading to the identification of the first sRNAs influencing meningococcal bacteremia. Further analysis indicated that one of the sRNAs affecting bacteremia responded to carbon source availability through repression by a GntR-like transcriptional regulator. Both the sRNA and the GntR-like regulator are implicated in the control of gene expression from a common network involved in energy metabolism.

Data-driven integration of genome-scale regulatory and metabolic network models

Microbes are diverse and extremely versatile organisms that play vital roles in all ecological niches. Understanding and harnessing microbial systems will be key to the sustainability of our planet. One approach to improving our knowledge of microbial processes is through data-driven and mechanism-informed computational modeling. Individual models of biological networks (such as metabolism, transcription, and signaling) have played pivotal roles in driving microbial research through the years. These networks, however, are highly interconnected and function in concert-a fact that has led to the development of a variety of approaches aimed at simulating the integrated functions of two or more network types. Though the task of integrating these different models is fraught with new challenges, the large amounts of high-throughput data sets being generated, and algorithms being developed, means that the time is at hand for concerted efforts to build integrated regulatory-metabolic networks in a data-driven fashion. In this perspective, we review current approaches for constructing integrated regulatory-metabolic models and outline new strategies for future development of these network models for any microbial system.

An integrated approach to reconstructing genome-scale transcriptional regulatory networks

Transcriptional regulatory networks (TRNs) program cells to dynamically alter their gene expression in response to changing internal or environmental conditions. In this study, we develop a novel workflow for generating large-scale TRN models that integrates comparative genomics data, global gene expression analyses, and intrinsic properties of transcription factors (TFs). An assessment of this workflow using benchmark datasets for the well-studied γ-proteobacterium Escherichia coli showed that it outperforms expression-based inference approaches, having a significantly larger area under the precision-recall curve. Further analysis indicated that this integrated workflow captures different aspects of the E. coli TRN than expression-based approaches, potentially making them highly complementary. We leveraged this new workflow and observations to build a large-scale TRN model for the α-Proteobacterium Rhodobacter sphaeroides that comprises 120 gene clusters, 1211 genes (including 93 TFs), 1858 predicted protein-DNA interactions and 76 DNA binding motifs. We found that ~67% of the predicted gene clusters in this TRN are enriched for functions ranging from photosynthesis or central carbon metabolism to environmental stress responses. We also found that members of many of the predicted gene clusters were consistent with prior knowledge in R. sphaeroides and/or other bacteria. Experimental validation of predictions from this R. sphaeroides TRN model showed that high precision and recall was also obtained for TFs involved in photosynthesis (PpsR), carbon metabolism (RSP_0489) and iron homeostasis (RSP_3341). In addition, this integrative approach enabled generation of TRNs with increased information content relative to R. sphaeroides TRN models built via other approaches. We also show how this approach can be used to simultaneously produce TRN models for each related organism used in the comparative genomics analysis. Our results highlight the advantages of integrating comparative genomics of closely related organisms with gene expression data to assemble large-scale TRN models with high-quality predictions.

The ever growing amount of genomic data enables the assembly of large-scale network models that can provide important new insights into living systems. However, assembly and validation of such large-scale models can be challenging, since we often lack sufficient information to make accurate predictions. This work describes a new approach for constructing large-scale transcriptional regulatory networks of individual cells. We show that the reconstructed network captures a significantly larger fraction of cellular regulatory processes than networks generated by other existing approaches. We predict this approach, with appropriate refinements, will allow reconstruction of large-scale transcriptional network models for a variety of other organisms. As we work towards modeling the function of cells or complex ecosystems, individually reconstructed network models of signaling, information transfer and metabolism, can be integrated to provide high information predictions and insights not otherwise obtainable.

Expression of extra-cellular levansucrase in Pseudomonas syringae is controlled by the in planta fitness-promoting metabolic repressor HexR

BACKGROUND

Pseudomonas syringae pv. glycinea PG4180 causes bacterial blight on soybean plants and enters the leaf tissue through stomata or open wounds, where it encounters a sucrose-rich milieu. Sucrose is utilized by invading bacteria via the secreted enzyme, levansucrase (Lsc), liberating glucose and forming the polyfructan levan. P. syringae PG4180 possesses two functional lsc alleles transcribed at virulence-promoting low temperatures.

RESULTS

We hypothesized that transcription of lsc is controlled by the hexose metabolism repressor, HexR, since potential HexR binding sites were identified upstream of both lsc genes. A hexR mutant of PG4180 was significantly growth-impaired when incubated with sucrose or glucose as sole carbon source, but exhibited wild type growth when arabinose was provided. Analyses of lsc expression resulted in higher transcript and protein levels in the hexR mutant as compared to the wild type. The hexR mutant's ability to multiply in planta was reduced. HexR did not seem to impact hrp gene expression as evidenced by the hexR mutant's unaltered hypersensitive response in tobacco and its unmodified protein secretion pattern as compared to the wild type under hrp-inducing conditions.

CONCLUSIONS

Our data suggested a co-regulation of genes involved in extra-cellular sugar acquisition with those involved in intra-cellular energy-providing metabolic pathways in P. syringae.

CceR and AkgR regulate central carbon and energy metabolism in alphaproteobacteria

Many pathways of carbon and energy metabolism are conserved across the phylogeny, but the networks that regulate their expression or activity often vary considerably among organisms. In this work, we show that two previously uncharacterized transcription factors (TFs) are direct regulators of genes encoding enzymes of central carbon and energy metabolism in the alphaproteobacterium Rhodobacter sphaeroides. The LacI family member CceR (RSP_1663) directly represses genes encoding enzymes in the Entner-Doudoroff pathway, while activating those encoding the F1F0 ATPase and enzymes of the tricarboxylic acid (TCA) cycle and gluconeogenesis, providing a direct transcriptional network connection between carbon and energy metabolism. We identified bases that are important for CceR DNA binding and showed that DNA binding by this TF is inhibited by 6-phosphogluconate. We also showed that the GntR family TF AkgR (RSP_0981) directly activates genes encoding several TCA cycle enzymes, and we identified conditions where its activity is increased. The properties of single and double ΔCceR and ΔAkgR mutants illustrate that these 2 TFs cooperatively regulate carbon and energy metabolism. Comparative genomic analysis indicates that CceR and AkgR orthologs are found in other alphaproteobacteria, where they are predicted to have a conserved function in regulating central carbon metabolism. Our characterization of CceR and AkgR has provided important new insight into the networks that control central carbon and energy metabolism in alphaproteobacteria that can be exploited to modify or engineer new traits in these widespread and versatile bacteria.

IMPORTANCE

To extract and conserve energy from nutrients, cells coordinate a set of metabolic pathways into integrated networks. Many pathways that conserve energy or interconvert metabolites are conserved across cells, but the networks regulating these processes are often highly variable. In this study, we characterize two previously unknown transcriptional regulators of carbon and energy metabolism that are conserved in alphaproteobacteria, a group of abundant, environmentally and biotechnologically important organisms. We identify the genes they regulate, the DNA sequences they recognize, the metabolite that controls the activity of one of the regulators, and conditions where they are required for growth. We provide important new insight into conserved cellular networks that can also be used to improve a variety of hosts for converting feedstock into valuable products.

Insights into Vibrio parahaemolyticus CHN25 response to artificial gastric fluid stress by transcriptomic analysis

Vibrio parahaemolyticus is the causative agent of food-borne gastroenteritis disease. Once consumed, human acid gastric fluid is perhaps one of the most important environmental stresses imposed on the bacterium. Herein, for the first time, we investigated Vibrio parahaemolyticus CHN25 response to artificial gastric fluid (AGF) stress by transcriptomic analysis. The bacterium at logarithmic growth phase (LGP) displayed lower survival rates than that at stationary growth phase (SGP) under a sub-lethal acid condition (pH 4.9). Transcriptome data revealed that 11.6% of the expressed genes in Vibrio parahaemolyticus CHN25 was up-regulated in LGP cells after exposed to AGF (pH 4.9) for 30 min, including those involved in sugar transport, nitrogen metabolism, energy production and protein biosynthesis, whereas 14.0% of the genes was down-regulated, such as ATP-binding cassette (ABC) transporter and flagellar biosynthesis genes. In contrast, the AGF stress only elicited 3.4% of the genes from SGP cells, the majority of which were attenuated in expression. Moreover, the number of expressed regulator genes was also substantially reduced in SGP cells. Comparison of transcriptome profiles further revealed forty-one growth-phase independent genes in the AGF stress, however, half of which displayed distinct expression features between the two growth phases. Vibrio parahaemolyticus seemed to have evolved a number of molecular strategies for coping with the acid stress. The data here will facilitate future studies for environmental stresses and pathogenicity of the leading seafood-borne pathogen worldwide.

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.

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.

Genomic reconstruction of the transcriptional regulatory network in Bacillus subtilis

The adaptation of microorganisms to their environment is controlled by complex transcriptional regulatory networks (TRNs), which are still only partially understood even for model species. Genome scale annotation of regulatory features of genes and TRN reconstruction are challenging tasks of microbial genomics. We used the knowledge-driven comparative-genomics approach implemented in the RegPredict Web server to infer TRN in the model Gram-positive bacterium Bacillus subtilis and 10 related Bacillales species. For transcription factor (TF) regulons, we combined the available information from the DBTBS database and the literature with bioinformatics tools, allowing inference of TF binding sites (TFBSs), comparative analysis of the genomic context of predicted TFBSs, functional assignment of target genes, and effector prediction. For RNA regulons, we used known RNA regulatory motifs collected in the Rfam database to scan genomes and analyze the genomic context of new RNA sites. The inferred TRN in B. subtilis comprises regulons for 129 TFs and 24 regulatory RNA families. First, we analyzed 66 TF regulons with previously known TFBSs in B. subtilis and projected them to other Bacillales genomes, resulting in refinement of TFBS motifs and identification of novel regulon members. Second, we inferred motifs and described regulons for 28 experimentally studied TFs with previously unknown TFBSs. Third, we discovered novel motifs and reconstructed regulons for 36 previously uncharacterized TFs. The inferred collection of regulons is available in the RegPrecise database (http://regprecise.lbl.gov/) and can be used in genetic experiments, metabolic modeling, and evolutionary analysis.

Somewhat in control--the role of transcription in regulating microbial metabolic fluxes

The most common way for microbes to control their metabolism is by controlling enzyme levels through transcriptional regulation. Yet recent studies have shown that in many cases, perturbations to the transcriptional regulatory network do not result in altered metabolic phenotypes on the level of the flux distribution. We suggest that this may be a consequence of cells protecting their metabolism against stochastic fluctuations in expression as well as enabling a fast response for those fluxes that may need to be changed quickly. Furthermore, it is impossible for a regulatory program to guarantee optimal expression levels in all conditions. Several studies have found examples of demonstrably suboptimal regulation of gene expression, and improvements to the regulatory network have been investigated in laboratory evolution experiments.

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.

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.

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.