CORAL: A framework for rigorous self-validated data modeling and integrative, reproducible data analysis.

BACKGROUND: Many organizations face challenges in managing and analyzing data, especially when relevant datasets arise from multiple sources and methods. Analyzing heterogeneous datasets and additional derived data requires rigorous tracking of their interrelationships and provenance. This task has long been a Grand Challenge of data science and has more recently been formalized in the FAIR principles: that all data objects be Findable, Accessible, Interoperable, and Reusable, both for machines and for people. Adherence to these principles is necessary for proper stewardship of information, for testing regulatory compliance, for measuring the efficiency of processes, and for facilitating reuse of data-analytical frameworks. FINDINGS: We present the Contextual Ontology-based Repository Analysis Library (CORAL), a platform that greatly facilitates adherence to all 4 of the FAIR principles, including the especially difficult challenge of making heterogeneous datasets Interoperable and Reusable across all parts of a large, long-lasting organization. To achieve this, CORAL's data model requires that data generators extensively document the context for all data, and our tools maintain that context throughout the entire analysis pipeline. CORAL also features a web interface for data generators to upload and explore data, as well as a Jupyter notebook interface for data analysts, both backed by a common API. CONCLUSIONS: CORAL enables organizations to build FAIR data types on the fly as they are needed, avoiding the expense of bespoke data modeling. CORAL provides a uniquely powerful platform to enable integrative cross-dataset analyses, generating deeper insights than are possible using traditional analysis tools.

Binary Metabolic Phenotypes and Phenotype Diversity Metrics for the Functional Characterization of Microbial Communities.

The profiling of 16S rRNA revolutionized the exploration of microbiomes, allowing to describe community composition by enumerating relevant taxa and their abundances. However, taxonomic profiles alone lack interpretability in terms of bacterial metabolism, and their translation into functional characteristics of microbiomes is a challenging task. This bottom-up approach minimally requires a reference collection of major metabolic traits deduced from the complete genomes of individual organisms, an accurate method of projecting these traits from a reference collection to the analyzed amplicon sequence variants (ASVs), and, ultimately, an approach to a microbiome-wide aggregation of predicted individual traits into physiologically relevant cumulative metrics to characterize and compare multiple microbiome samples. In this study, we extended a previously introduced computational approach for the functional profiling of complex microbial communities, which is based on the concept of binary metabolic phenotypes encoding the presence ("1") or absence ("0") of various measurable physiological properties in individual organisms that are termed phenotype carriers or non-carriers, respectively. Derived from complete genomes via metabolic reconstruction, binary phenotypes provide a foundation for the prediction of functional traits for each ASV identified in a microbiome sample. Here, we introduced three distinct mapping schemes for a microbiome-wide phenotype prediction and assessed their accuracy on the 16S datasets of mock bacterial communities representing human gut microbiome (HGM) as well as on two large HGM datasets, the American Gut Project and the UK twins study. The 16S sequence-based scheme yielded a more accurate phenotype predictions, while the taxonomy-based schemes demonstrated a reasonable performance to warrant their application for other types of input data (e.g., from shotgun metagenomics or qPCR). In addition to the abundance-weighted Community Phenotype Indices (CPIs) reflecting the fractional representation of various phenotype carriers in microbiome samples, we employ metrics capturing the diversity of phenotype carriers, Phenotype Alpha Diversity (PAD) and Phenotype Beta Diversity (PBD). In combination with CPI, PAD allows to classify the robustness of metabolic phenotypes by their anticipated stability in the face of potential environmental perturbations. PBD provides a promising approach for detecting the metabolic features potentially contributing to disease-associated metabolic traits as illustrated by a comparative analysis of HGM samples from healthy and Crohn's disease cohorts.

Characterization of a Metal-Resistant Bacillus Strain With a High Molybdate Affinity ModA From Contaminated Sediments at the Oak Ridge Reservation.

A nitrate- and metal-contaminated site at the Oak Ridge Reservation (ORR) was previously shown to contain the metal molybdenum (Mo) at picomolar concentrations. This potentially limits microbial nitrate reduction, as Mo is required by the enzyme nitrate reductase, which catalyzes the first step of nitrate removal. Enrichment for anaerobic nitrate-reducing microbes from contaminated sediment at the ORR yielded Bacillus strain EB106-08-02-XG196. This bacterium grows in the presence of multiple metals (Cd, Ni, Cu, Co, Mn, and U) but also exhibits better growth compared to control strains, including Pseudomonas fluorescens N2E2 isolated from a pristine ORR environment under low molybdate concentrations (<1 nM). Molybdate is taken up by the molybdate binding protein, ModA, of the molybdate ATP-binding cassette transporter. ModA of XG196 is phylogenetically distinct from those of other characterized ModA proteins. The genes encoding ModA from XG196, P. fluorescens N2E2 and Escherichia coli K12 were expressed in E. coli and the recombinant proteins were purified. Isothermal titration calorimetry analysis showed that XG196 ModA has a higher affinity for molybdate than other ModA proteins with a molybdate binding constant (K D ) of 2.2 nM, about one order of magnitude lower than those of P. fluorescens N2E2 (27.0 nM) and E. coli K12 (25.0 nM). XG196 ModA also showed a fivefold higher affinity for molybdate than for tungstate (11 nM), whereas the ModA proteins from P. fluorescens N2E2 [K D (Mo) 27.0 nM, K D (W) 26.7 nM] and E. coli K12[(K D (Mo) 25.0 nM, K D (W) 23.8 nM] had similar affinities for the two oxyanions. We propose that high molybdate affinity coupled with resistance to multiple metals gives strain XG196 a competitive advantage in Mo-limited environments contaminated with high concentrations of metals and nitrate, as found at ORR.

Draft Genome Sequence of Bacillus sp. Strain EB106-08-02-XG196, Isolated from High-Nitrate-Contaminated Sediment.

Bacillus sp. strain EB106-08-02-XG196 was isolated from a high-nitrate- and heavy metal-contaminated site at the Oak Ridge Reservation in Tennessee. We report the draft genome sequence of this strain to provide insights into the genomic basis for surviving in this unique environment.

High-throughput mapping of the phage resistance landscape in E. coli.

Bacteriophages (phages) are critical players in the dynamics and function of microbial communities and drive processes as diverse as global biogeochemical cycles and human health. Phages tend to be predators finely tuned to attack specific hosts, even down to the strain level, which in turn defend themselves using an array of mechanisms. However, to date, efforts to rapidly and comprehensively identify bacterial host factors important in phage infection and resistance have yet to be fully realized. Here, we globally map the host genetic determinants involved in resistance to 14 phylogenetically diverse double-stranded DNA phages using two model Escherichia coli strains (K-12 and BL21) with known sequence divergence to demonstrate strain-specific differences. Using genome-wide loss-of-function and gain-of-function genetic technologies, we are able to confirm previously described phage receptors as well as uncover a number of previously unknown host factors that confer resistance to one or more of these phages. We uncover differences in resistance factors that strongly align with the susceptibility of K-12 and BL21 to specific phage. We also identify both phage-specific mechanisms, such as the unexpected role of cyclic-di-GMP in host sensitivity to phage N4, and more generic defenses, such as the overproduction of colanic acid capsular polysaccharide that defends against a wide array of phages. Our results indicate that host responses to phages can occur via diverse cellular mechanisms. Our systematic and high-throughput genetic workflow to characterize phage-host interaction determinants can be extended to diverse bacteria to generate datasets that allow predictive models of how phage-mediated selection will shape bacterial phenotype and evolution. The results of this study and future efforts to map the phage resistance landscape will lead to new insights into the coevolution of hosts and their phage, which can ultimately be used to design better phage therapeutic treatments and tools for precision microbiome engineering.

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.

Micronutrient Requirements and Sharing Capabilities of the Human Gut Microbiome.

The human gut microbiome harbors a diverse array of metabolic pathways contributing to its development and homeostasis via a complex web of diet-dependent metabolic interactions within the microbial community and host. Genomics-based reconstruction and predictive modeling of these interactions would provide a framework for diagnostics and treatment of dysbiosis-related syndromes via rational selection of therapeutic prebiotics and dietary nutrients. Of particular interest are micronutrients, such as B-group vitamins, precursors of indispensable metabolic cofactors, that are produced de novo by some gut bacteria (prototrophs) but must be provided exogenously in the diet for many other bacterial species (auxotrophs) as well as for the mammalian host. Cross-feeding of B vitamins between prototrophic and auxotrophic species is expected to strongly contribute to the homeostasis of microbial communities in the distal gut given the efficient absorption of dietary vitamins in the upper gastrointestinal tract. To confidently estimate the balance of microbiome micronutrient biosynthetic capabilities and requirements using available genomic data, we have performed a subsystems-based reconstruction of biogenesis, salvage and uptake for eight B vitamins (B1, B2, B3, B5, B6, B7, B9, and B12) and queuosine (essential factor in tRNA modification) over a reference set of 2,228 bacterial genomes representing 690 cultured species of the human gastrointestinal microbiota. This allowed us to classify the studied organisms with respect to their pathway variants and infer their prototrophic vs. auxotrophic phenotypes. In addition to canonical vitamin pathways, several conserved partial pathways were identified pointing to alternative routes of syntrophic metabolism and expanding a microbial vitamin "menu" by several pathway intermediates (vitamers) such as thiazole, quinolinate, dethiobiotin, pantoate. A cross-species comparison was applied to assess the extent of conservation of vitamin phenotypes at distinct taxonomic levels (from strains to families). The obtained reference collection combined with 16S rRNA gene-based phylogenetic profiles was used to deduce phenotype profiles of the human gut microbiota across in two large cohorts. This analysis provided the first estimate of B-vitamin requirements, production and sharing capabilities in the human gut microbiome establishing predictive phenotype profiling as a new approach to classification of microbiome samples. Future expansion of our reference genomic collection of metabolic phenotypes will allow further improvement in coverage and accuracy of predictive phenotype profiling of the human microbiome.

LurR is a regulator of the central lactate oxidation pathway in sulfate-reducing Desulfovibrio species.

The central carbon/lactate utilization pathway in the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough, is encoded by the highly conserved operon DVU3025-3033. Our earlier in vitro genome-wide study had suggested a network of four two-component system regulators that target this large operon; however, how these four regulators control this operon was not known. Here, we probe the regulation of the lactate utilization operon with mutant strains and DNA-protein binding assays. We show that the LurR response regulator is required for optimal growth and complete lactate utilization, and that it activates the DVU3025-3033 lactate oxidation operon as well as DVU2451, a lactate permease gene, in the presence of lactate. We show by electrophoretic mobility shift assays that LurR binds to three sites in the upstream region of DVU3025, the first gene of the operon. NrfR, a response regulator that is activated under nitrite stress, and LurR share similar binding site motifs and bind the same sites upstream of DVU3025. The DVU3025 promoter also has a binding site motif (Pho box) that is bound by PhoB, a two-component response regulator activated under phosphate limitation. The lactate utilization operon, the regulator LurR, and LurR binding sites are conserved across the order Desulfovibrionales whereas possible modulation of the lactate utilization genes by additional regulators such as NrfR and PhoB appears to be limited to D. vulgaris.

Dual-barcoded shotgun expression library sequencing for high-throughput characterization of functional traits in bacteria.

A major challenge in genomics is the knowledge gap between sequence and its encoded function. Gain-of-function methods based on gene overexpression are attractive avenues for phenotype-based functional screens, but are not easily applied in high-throughput across many experimental conditions. Here, we present Dual Barcoded Shotgun Expression Library Sequencing (Dub-seq), a method that uses random DNA barcodes to greatly increase experimental throughput. As a demonstration of this approach, we construct a Dub-seq library with Escherichia coli genomic DNA, performed 155 genome-wide fitness assays in 52 experimental conditions, and identified overexpression phenotypes for 813 genes. We show that Dub-seq data is reproducible, accurately recapitulates known biology, and identifies hundreds of novel gain-of-function phenotypes for E. coli genes, a subset of which we verified with assays of individual strains. Dub-seq provides complementary information to loss-of-function approaches and will facilitate rapid and systematic functional characterization of microbial genomes.

Author Correction: Jungle Express is a versatile repressor system for tight transcriptional control.

In the original version of this Article, an incorrect URL was provided in the Data Availability Statement regarding the deposition of plasmids listed in Supplementary Table 4. The correct URL is https://public-registry.jbei.org/folders/378 . This error has been corrected in both the PDF and HTML versions of the Article.

Jungle Express is a versatile repressor system for tight transcriptional control.

Tightly regulated promoters are essential for numerous biological applications, where strong inducibility, portability, and scalability are desirable. Current systems are often incompatible with large-scale fermentations due to high inducer costs and strict media requirements. Here, we describe the bottom-up engineering of 'Jungle Express', an expression system that enables efficient gene regulation in diverse proteobacteria. This system is guided by EilR, a multidrug-binding repressor with high affinity to its optimized operator and cationic dyes that act as powerful inducers at negligible costs. In E. coli, the engineered promoters exhibit minimal basal transcription and are inducible over four orders of magnitude by 1 µM crystal violet, reaching expression levels exceeding those of the strongest current bacterial systems. Further, we provide molecular insights into specific interactions of EilR with its operator and with two inducers. The versatility of Jungle Express opens the way for tightly controlled and efficient gene expression that is not restricted to host organism, substrate, or scale.

Multiple signaling systems target a core set of transition metal homeostasis genes using similar binding motifs.

Bacterial response to metals can require complex regulation. We report an overlapping regulation for copper and zinc resistance genes in the denitrifying bacterium, Pseudomonas stutzeri RCH2, by three two-component regulatory proteins CopR1, CopR2 and CzcR. We conducted genome-wide evaluations to identify gene targets of two paralogous regulators, CopR1 and CopR2, annotated for copper signaling, and compared the results with the gene targets for CzcR, implicated in zinc signaling. We discovered that the CopRs and CzcR have largely common targets, and crossregulate a core set of P. stutzeri copper and zinc responsive genes. We established that this crossregulation is enabled by a conserved binding motif in the upstream regulatory regions of the target genes. The crossregulation is physiologically relevant as these regulators synergistically and antagonistically target multicopper oxidases, metal efflux and sequestration systems. CopR1 and CopR2 upregulate two cop operons encoding copper tolerance genes, while all three regulators downregulate a putative copper chaperone, Psest_1595. CzcR also upregulated the oprD gene and the CzcIABC Zn2+ efflux system, while CopR1 and CopR2 downregulated these genes. Our study suggests that crossregulation of copper and zinc homeostasis can be advantageous, and in P. stutzeri this is enabled by shared binding motifs for multiple response regulators.

Draft Genome Sequences of Two Janthinobacteriumlividum Strains, Isolated from Pristine Groundwater Collected from the Oak Ridge Field Research Center.

We present here the draft genome sequences of two Janthinobacterium lividum strains, GW456P and GW458P, isolated from groundwater samples collected from a background site at the Oak Ridge Field Research Center. Production of a purple pigment by these two strains was observed when grown on diluted (1/10) LB agar plates.

Complete genome sequence of Pseudomonas stutzeri strain RCH2 isolated from a Hexavalent Chromium [Cr(VI)] contaminated site.

Hexavalent Chromium [Cr(VI)] is a widespread contaminant found in soil, sediment, and ground water in several DOE sites, including Hanford 100 H area. In order to stimulate microbially mediated reduction of Cr(VI) at this site, a poly-lactate hydrogen release compound was injected into the chromium contaminated aquifer. Targeted enrichment of dominant nitrate-reducing bacteria post injection resulted in the isolation of Pseudomonas stutzeri strain RCH2. P. stutzeri strain RCH2 was isolated using acetate as the electron donor and is a complete denitrifier. Experiments with anaerobic washed cell suspension of strain RCH2 revealed it could reduce Cr(VI) and Fe(III). The genome of strain RCH2 was sequenced using a combination of Illumina and 454 sequencing technologies and contained a circular chromosome of 4.6 Mb and three plasmids. Global genome comparisons of strain RCH2 with six other fully sequenced P. stutzeri strains revealed most genomic regions are conserved, however strain RCH2 has an additional 244 genes, some of which are involved in chemotaxis, Flp pilus biogenesis and pyruvate/2-oxogluturate complex formation.

Novel Metal Cation Resistance Systems from Mutant Fitness Analysis of Denitrifying Pseudomonas stutzeri.

UNLABELLED: Metal ion transport systems have been studied extensively, but the specificity of a given transporter is often unclear from amino acid sequence data alone. In this study, predicted Cu(2+) and Zn(2+) resistance systems in Pseudomonas stutzeri strain RCH2 are compared with those experimentally implicated in Cu(2+) and Zn(2+) resistance, as determined by using a DNA-barcoded transposon mutant library. Mutant fitness data obtained under denitrifying conditions are combined with regulon predictions to yield a much more comprehensive picture of Cu(2+) and Zn(2+) resistance in strain RCH2. The results not only considerably expand what is known about well-established metal ion exporters (CzcCBA, CzcD, and CusCBA) and their accessory proteins (CzcI and CusF), they also reveal that isolates with mutations in some predicted Cu(2+) resistance systems do not show decreased fitness relative to the wild type when exposed to Cu(2+) In addition, new genes are identified that have no known connection to Zn(2+) (corB, corC, Psest_3226, Psest_3322, and Psest_0618) or Cu(2+) resistance (Mrp antiporter subunit gene, Psest_2850, and Psest_0584) but are crucial for resistance to these metal cations. Growth of individual deletion mutants lacking corB, corC, Psest_3226, or Psest_3322 confirmed the observed Zn-dependent phenotypes. Notably, to our knowledge, this is the first time a bacterial homolog of TMEM165, a human gene responsible for a congenital glycosylation disorder, has been deleted and the resulting strain characterized. Finally, the fitness values indicate Cu(2+)- and Zn(2+)-based inhibition of nitrite reductase and interference with molybdenum cofactor biosynthesis for nitrate reductase. These results extend the current understanding of Cu(2+) and Zn(2+) efflux and resistance and their effects on denitrifying metabolism. IMPORTANCE: In this study, genome-wide mutant fitness data in P. stutzeri RCH2 combined with regulon predictions identify several proteins of unknown function that are involved in resisting zinc and copper toxicity. For zinc, these include a member of the UPF0016 protein family that was previously implicated in Ca(2+)/H(+) antiport and a human congenital glycosylation disorder, CorB and CorC, which were previously linked to Mg(2+) transport, and Psest_3322 and Psest_0618, two proteins with no characterized homologs. Experiments using mutants lacking Psest_3226, Psest_3322, corB, corC, or czcI verified their proposed functions, which will enable future studies of these little-characterized zinc resistance determinants. Likewise, Psest_2850, annotated as an ion antiporter subunit, and the conserved hypothetical protein Psest_0584 are implicated in copper resistance. Physiological connections between previous studies and phenotypes presented here are discussed. Functional and mechanistic understanding of transport proteins improves the understanding of systems in which members of the same protein family, including those in humans, can have different functions.

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.

σ54-dependent regulome in Desulfovibrio vulgaris Hildenborough.

BACKGROUND: The σ(54) subunit controls a unique class of promoters in bacteria. Such promoters, without exception, require enhancer binding proteins (EBPs) for transcription initiation. Desulfovibrio vulgaris Hildenborough, a model bacterium for sulfate reduction studies, has a high number of EBPs, more than most sequenced bacteria. The cellular processes regulated by many of these EBPs remain unknown. RESULTS: To characterize the σ(54)-dependent regulome of D. vulgaris Hildenborough, we identified EBP binding motifs and regulated genes by a combination of computational and experimental techniques. These predictions were supported by our reconstruction of σ(54)-dependent promoters by comparative genomics. We reassessed and refined the results of earlier studies on regulation in D. vulgaris Hildenborough and consolidated them with our new findings. It allowed us to reconstruct the σ(54) regulome in D. vulgaris Hildenborough. This regulome includes 36 regulons that consist of 201 coding genes and 4 non-coding RNAs, and is involved in nitrogen, carbon and energy metabolism, regulation, transmembrane transport and various extracellular functions. To the best of our knowledge, this is the first report of direct regulation of alanine dehydrogenase, pyruvate metabolism genes and type III secretion system by σ(54)-dependent regulators. CONCLUSIONS: The σ(54)-dependent regulome is an important component of transcriptional regulatory network in D. vulgaris Hildenborough and related free-living Deltaproteobacteria. Our study provides a representative collection of σ(54)-dependent regulons that can be used for regulation prediction in Deltaproteobacteria and other taxa.

Regulation of Nitrite Stress Response in Desulfovibrio vulgaris Hildenborough, a Model Sulfate-Reducing Bacterium.

UNLABELLED: Sulfate-reducing bacteria (SRB) are sensitive to low concentrations of nitrite, and nitrite has been used to control SRB-related biofouling in oil fields. Desulfovibrio vulgaris Hildenborough, a model SRB, carries a cytochrome c-type nitrite reductase (nrfHA) that confers resistance to low concentrations of nitrite. The regulation of this nitrite reductase has not been directly examined to date. In this study, we show that DVU0621 (NrfR), a sigma54-dependent two-component system response regulator, is the positive regulator for this operon. NrfR activates the expression of the nrfHA operon in response to nitrite stress. We also show that nrfR is needed for fitness at low cell densities in the presence of nitrite because inactivation of nrfR affects the rate of nitrite reduction. We also predict and validate the binding sites for NrfR upstream of the nrfHA operon using purified NrfR in gel shift assays. We discuss possible roles for NrfR in regulating nitrate reductase genes in nitrate-utilizing Desulfovibrio spp. IMPORTANCE: The NrfA nitrite reductase is prevalent across several bacterial phyla and required for dissimilatory nitrite reduction. However, regulation of the nrfA gene has been studied in only a few nitrate-utilizing bacteria. Here, we show that in D. vulgaris, a bacterium that does not respire nitrate, the expression of nrfHA is induced by NrfR upon nitrite stress. This is the first report of regulation of nrfA by a sigma54-dependent two-component system. Our study increases our knowledge of nitrite stress responses and possibly of the regulation of nitrate reduction in SRB.

Comparative genomics of transcriptional regulation of methionine metabolism in Proteobacteria.

Methionine metabolism and uptake genes in Proteobacteria are controlled by a variety of RNA and DNA regulatory systems. We have applied comparative genomics to reconstruct regulons for three known transcription factors, MetJ, MetR, and SahR, and three known riboswitch motifs, SAH, SAM-SAH, and SAM_alpha, in ∼ 200 genomes from 22 taxonomic groups of Proteobacteria. We also identified two novel regulons: a SahR-like transcription factor SamR controlling various methionine biosynthesis genes in the Xanthomonadales group, and a potential RNA regulatory element with terminator-antiterminator mechanism controlling the metX or metZ genes in beta-proteobacteria. For each analyzed regulator we identified the core, taxon-specific and genome-specific regulon members. By analyzing the distribution of these regulators in bacterial genomes and by comparing their regulon contents we elucidated possible evolutionary scenarios for the regulation of the methionine metabolism genes in Proteobacteria.

Chlorate reduction in Shewanella algae†ACDC is a recently acquired metabolism characterized by gene loss, suboptimal regulation and oxidative stress.

Previous work on respiratory chlorate reduction has biochemically identified the terminal reductase ClrABC and the chlorite detoxifying enzyme Cld. In Shewanella algae ACDC, genes encoding these enzymes reside on composite transposons whose core we refer to as the chlorate reduction composite transposon interior (CRI). To better understand this metabolism in ACDC, we used RNA-seq and proteomics to predict carbon and electron flow during chlorate reduction and posit that formate is an important electron carrier with lactate as the electron donor, but that NADH predominates on acetate. Chlorate-specific transcription of electron transport chain components or the CRI was not observed, but clr and cld transcription was attenuated by oxygen. The major chlorate-specific response related to oxidative stress and was indicative of reactive chlorine species production. A genetic system based on rpsL-streptomycin counter selection was developed to further dissect the metabolism, but ACDC readily lost the CRI via homologous recombination of the composite transposon's flanking insertion sequences. An engineered strain containing a single chromosomal CRI did not grow on chlorate, but overexpression of cld and its neighbouring cytochrome c restored growth. We postulate that the recently acquired CRI underwent copy-number expansion to circumvent insufficient expression of key genes in the pathway.