Jayashree Ray

Interactive XCMS Online: Simplifying Advanced Metabolomic Data Processing and Subsequent Statistical Analyses

XCMS Online (xcmsonline.scripps.edu) is a cloud-based informatic platform designed to process and visualize mass-spectrometry-based, untargeted metabolomic data. Initially, the platform was developed for two-group comparisons to match the independent, “control” versus “disease” experimental design. Here, we introduce an enhanced XCMS Online interface that enables users to perform dependent (paired) two-group comparisons, meta-analysis, and multigroup comparisons, with comprehensive statistical output and interactive visualization tools. Newly incorporated statistical tests cover a wide array of univariate analyses. Multigroup comparison allows for the identification of differentially expressed metabolite features across multiple classes of data while higher order meta-analysis facilitates the identification of shared metabolic patterns across multiple two-group comparisons. Given the complexity of these data sets, we have developed an interactive platform where users can monitor the statistical output of univariate (cloud plots) and multivariate (PCA plots) data analysis in real time by adjusting the threshold and range of various parameters. On the interactive cloud plot, metabolite features can be filtered out by their significance level (p-value), fold change, mass-to-charge ratio, retention time, and intensity. The variation pattern of each feature can be visualized on both extracted-ion chromatograms and box plots. The interactive principal component analysis includes scores, loadings, and scree plots that can be adjusted depending on scaling criteria. The utility of XCMS functionalities is demonstrated through the metabolomic analysis of bacterial stress response and the comparison of lymphoblastic leukemia cell lines.

The genetic basis of energy conservation in the sulfate-reducing bacterium Desulfovibrio alaskensis G20.

Sulfate-reducing bacteria play major roles in the global carbon and sulfur cycles, but it remains unclear how reducing sulfate yields energy. To determine the genetic basis of energy conservation, we measured the fitness of thousands of pooled mutants of Desulfovibrio alaskensis G20 during growth in 12 different combinations of electron donors and acceptors. We show that ion pumping by the ferredoxin:NADH oxidoreductase Rnf is required whenever substrate-level phosphorylation is not possible. The uncharacterized complex Hdr/flox-1 (Dde_1207:13) is sometimes important alongside Rnf and may perform an electron bifurcation to generate more reduced ferredoxin from NADH to allow further ion pumping. Similarly, during the oxidation of malate or fumarate, the electron-bifurcating transhydrogenase NfnAB-2 (Dde_1250:1) is important and may generate reduced ferredoxin to allow additional ion pumping by Rnf. During formate oxidation, the periplasmic [NiFeSe] hydrogenase HysAB is required, which suggests that hydrogen forms in the periplasm, diffuses to the cytoplasm, and is used to reduce ferredoxin, thus providing a substrate for Rnf. During hydrogen utilization, the transmembrane electron transport complex Tmc is important and may move electrons from the periplasm into the cytoplasmic sulfite reduction pathway. Finally, mutants of many other putative electron carriers have no clear phenotype, which suggests that they are not important under our growth conditions, although we cannot rule out genetic redundancy.

Functional Genomics with a Comprehensive Library of Transposon Mutants for the Sulfate-Reducing Bacterium Desulfovibrio alaskensis G20

ABSTRACT The genomes of sulfate-reducing bacteria remain poorly characterized, largely due to a paucity of experimental data and genetic tools. To meet this challenge, we generated an archived library of 15,477 mapped transposon insertion mutants in the sulfate-reducing bacterium Desulfovibrio alaskensis G20. To demonstrate the utility of the individual mutants, we profiled gene expression in mutants of six regulatory genes and used these data, together with 1,313 high-confidence transcription start sites identified by tiling microarrays and transcriptome sequencing (5′ RNA-Seq), to update the regulons of Fur and Rex and to confirm the predicted regulons of LysX, PhnF, PerR, and Dde_3000, a histidine kinase. In addition to enabling single mutant investigations, the D. alaskensis G20 transposon mutants also contain DNA bar codes, which enables the pooling and analysis of mutant fitness for thousands of strains simultaneously. Using two pools of mutants that represent insertions in 2,369 unique protein-coding genes, we demonstrate that the hypothetical gene Dde_3007 is required for methionine biosynthesis. Using comparative genomics, we propose that Dde_3007 performs a missing step in methionine biosynthesis by transferring a sulfur group to O-phosphohomoserine to form homocysteine. Additionally, we show that the entire choline utilization cluster is important for fitness in choline sulfate medium, which confirms that a functional microcompartment is required for choline oxidation. Finally, we demonstrate that Dde_3291, a MerR-like transcription factor, is a choline-dependent activator of the choline utilization cluster. Taken together, our data set and genetic resources provide a foundation for systems-level investigation of a poorly studied group of bacteria of environmental and industrial importance. IMPORTANCE Sulfate-reducing bacteria contribute to global nutrient cycles and are a nuisance for the petroleum industry. Despite their environmental and industrial significance, the genomes of sulfate-reducing bacteria remain poorly characterized. Here, we describe a genetic approach to fill gaps in our knowledge of sulfate-reducing bacteria. We generated a large collection of archived, transposon mutants in Desulfovibrio alaskensis G20 and used the phenotypes of these mutant strains to infer the function of genes involved in gene regulation, methionine biosynthesis, and choline utilization. Our findings and mutant resources will enable systematic investigations into gene function, energy generation, stress response, and metabolism for this important group of bacteria. Sulfate-reducing bacteria contribute to global nutrient cycles and are a nuisance for the petroleum industry. Despite their environmental and industrial significance, the genomes of sulfate-reducing bacteria remain poorly characterized. Here, we describe a genetic approach to fill gaps in our knowledge of sulfate-reducing bacteria. We generated a large collection of archived, transposon mutants in Desulfovibrio alaskensis G20 and used the phenotypes of these mutant strains to infer the function of genes involved in gene regulation, methionine biosynthesis, and choline utilization. Our findings and mutant resources will enable systematic investigations into gene function, energy generation, stress response, and metabolism for this important group of bacteria.

Deep Annotation of Protein Function across Diverse Bacteria from Mutant Phenotypes

Summary The function of nearly half of all protein-coding genes identified in bacterial genomes remains unknown. To systematically explore the functions of these proteins, we generated saturated transposon mutant libraries from 25 diverse bacteria and we assayed mutant phenotypes across hundreds of distinct conditions. From 3,903 genome-wide mutant fitness assays, we obtained 14.9 million gene phenotype measurements and we identified a mutant phenotype for 8,487 proteins with previously unknown functions. The majority of these hypothetical proteins (57%) had phenotypes that were either specific to a few conditions or were similar to that of another gene, thus enabling us to make informed predictions of protein function. For 1,914 of these hypothetical proteins, the functional associations are conserved across related proteins from different bacteria, which confirms that these associations are genuine. This comprehensive catalogue of experimentally-annotated protein functions also enables the targeted exploration of specific biological processes. For example, sensitivity to a DNA-damaging agent revealed 28 known families of DNA repair proteins and 11 putative novel families. Across all sequenced bacteria, 14% of proteins that lack detailed annotations have an ortholog with a functional association in our data set. Our study demonstrates the utility and scalability of high-throughput genetics for large-scale annotation of bacterial proteins and provides a vast compendium of experimentally-determined protein functions across diverse bacteria.

The energy-conserving electron transfer system used by Desulfovibrio alaskensis strain G20 during pyruvate fermentation involves reduction of endogenously formed fumarate and cytoplasmic and membrane-bound complexes, Hdr-Flox and Rnf.

The adaptation capability of Desulfovibrio to natural fluctuations in electron acceptor availability was evaluated by studying Desulfovibrio alaskensis strain G20 under varying respiratory, fermentative and methanogenic coculture conditions in chemostats. Transition from lactate to pyruvate in coculture resulted in a dramatic shift in the population structure and closer interspecies cell-to-cell interactions. Lower methane production rates in coculture than predicted from pyruvate input was attributed to redirection of electron flow to fumarate reduction. Without a methanogenic partner, accumulation of H₂and formate resulted in greater succinate production. Comparative transcript and gene fitness analysis in concert with physiological data of G20 wildtype and mutants demonstrated that pyruvate fermentation involves respiration of cytoplasmically formed fumarate using cytoplasmic and membrane-bound energy-conserving complexes, Rnf, Hdr-Flox-1 and Hmc. At the low H₂/formate levels maintained in coculture, Rnf likely functions as proton-pumping ferredoxin (Fd): type-I cytochrome c oxidoreductase, which transitions to a proton-pumping Fd(red):  nicotinamide adenine dinucleotide (NAD⁺) oxidoreductase at high H₂/formate levels during fermentation in monoculture. Hdr-Flox-1 is postulated to recycle Fd(red) via a flavin-based electron bifurcation involving NADH, Fdox and the thiol/disulphide-containing DsrC. In a menaquinone (MQ)-based electron confurcation reaction, the high-molecular-weight cytochrome-c₃complex, Hmc, is proposed to then couple DsrC(red) and periplasmic H₂/formate oxidation using the MQ pool to fuel a membrane-bound fumarate reductase.

Mutant phenotypes for thousands of bacterial genes of unknown function

One-third of all protein-coding genes from bacterial genomes cannot be annotated with a function. Here, to investigate the functions of these genes, we present genome-wide mutant fitness data from 32 diverse bacteria across dozens of growth conditions. We identified mutant phenotypes for 11,779 protein-coding genes that had not been annotated with a specific function. Many genes could be associated with a specific condition because the gene affected fitness only in that condition, or with another gene in the same bacterium because they had similar mutant phenotypes. Of the poorly annotated genes, 2,316 had associations that have high confidence because they are conserved in other bacteria. By combining these conserved associations with comparative genomics, we identified putative DNA repair proteins; in addition, we propose specific functions for poorly annotated enzymes and transporters and for uncharacterized protein families. Our study demonstrates the scalability of microbial genetics and its utility for improving gene annotations.A large-scale mutagenesis screen identifies mutant phenotypes for over 11,000 protein-coding genes in bacteria that had previously not been assigned a specific function.

Molybdenum Availability Is Key to Nitrate Removal in Contaminated Groundwater Environments

ABSTRACT The concentrations of molybdenum (Mo) and 25 other metals were measured in groundwater samples from 80 wells on the Oak Ridge Reservation (ORR) (Oak Ridge, TN), many of which are contaminated with nitrate, as well as uranium and various other metals. The concentrations of nitrate and uranium were in the ranges of 0.1 μM to 230 mM and <0.2 nM to 580 μM, respectively. Almost all metals examined had significantly greater median concentrations in a subset of wells that were highly contaminated with uranium (≥126 nM). They included cadmium, manganese, and cobalt, which were 1,300- to 2,700-fold higher. A notable exception, however, was Mo, which had a lower median concentration in the uranium-contaminated wells. This is significant, because Mo is essential in the dissimilatory nitrate reduction branch of the global nitrogen cycle. It is required at the catalytic site of nitrate reductase, the enzyme that reduces nitrate to nitrite. Moreover, more than 85% of the groundwater samples contained less than 10 nM Mo, whereas concentrations of 10 to 100 nM Mo were required for efficient growth by nitrate reduction for two Pseudomonas strains isolated from ORR wells and by a model denitrifier, Pseudomonas stutzeri RCH2. Higher concentrations of Mo tended to inhibit the growth of these strains due to the accumulation of toxic concentrations of nitrite, and this effect was exacerbated at high nitrate concentrations. The relevance of these results to a Mo-based nitrate removal strategy and the potential community-driving role that Mo plays in contaminated environments are discussed.

A metabolic pathway for catabolizing levulinic acid in bacteria.

Microorganisms can catabolize a wide range of organic compounds and therefore have the potential to perform many industrially relevant bioconversions. One barrier to realizing the potential of biorefining strategies lies in our incomplete knowledge of metabolic pathways, including those that can be used to assimilate naturally abundant or easily generated feedstocks. For instance, levulinic acid (LA) is a carbon source that is readily obtainable as a dehydration product of lignocellulosic biomass and can serve as the sole carbon source for some bacteria. Yet, the genetics and structure of LA catabolism have remained unknown. Here, we report the identification and characterization of a seven-gene operon that enables LA catabolism in Pseudomonas putida KT2440. When the pathway was reconstituted with purified proteins, we observed the formation of four acyl-CoA intermediates, including a unique 4-phosphovaleryl-CoA and the previously observed 3-hydroxyvaleryl-CoA product. Using adaptive evolution, we obtained a mutant of Escherichia coli LS5218 with functional deletions of fadE and atoC that was capable of robust growth on LA when it expressed the five enzymes from the P. putida operon. This discovery will enable more efficient use of biomass hydrolysates and metabolic engineering to develop bioconversions using LA as a feedstock.Levulinic acid (LA) is a value-added chemical easily obtained from biomass. The pathway enabling LA degradation in Pseudomonas putida requires five enzymes and can be engineered into Escherichia coli, which could enable further biotechnological applications.

Exometabolomics Assisted Design and Validation of Synthetic Obligate Mutualism

Synthetic microbial ecology has the potential to enhance the productivity and resiliency of biotechnology processes compared to approaches using single isolates. Engineering microbial consortia is challenging; however, one approach that has attracted significant attention is the creation of synthetic obligate mutualism using auxotrophic mutants that depend on each other for exchange or cross-feeding of metabolites. Here, we describe the integration of mutant library fitness profiling with mass spectrometry based exometabolomics as a method for constructing synthetic mutualism based on cross-feeding. Two industrially important species lacking known ecological interactions, Zymomonas mobilis and Escherichia coli, were selected as the test species. Amino acid exometabolites identified in the spent medium of Z. mobilis were used to select three corresponding E. coli auxotrophs (proA, pheA and IlvA), as potential E. coli counterparts for the coculture. A pooled mutant fitness assay with a Z. mobilis transposon mutant library was used to identify mutants with improved growth in the presence of E. coli. An auxotroph mutant in a gene (ZMO0748) with sequence similarity to cysteine synthase A (cysK), was selected as the Z. mobilis counterpart for the coculture. Exometabolomic analysis of spent E. coli medium identified glutathione related metabolites as potentially available for rescue of the Z. mobilis cysteine synthase mutant. Three sets of cocultures between the Z. mobilis auxotroph and each of the three E. coli auxotrophs were monitored by optical density for growth and analyzed by flow cytometry to confirm high cell counts for each species. Taken together, our methods provide a technological framework for creating synthetic mutualisms combining existing screening based methods and exometabolomics for both the selection of obligate mutualism partners and elucidation of metabolites involved in auxotroph rescue.

Complete genome sequence of Cupriavidus basilensis 4G11, isolated from the Oak Ridge Field Research Center site

ABSTRACT Cupriavidus basilensis 4G11 was isolated from groundwater at the Oak Ridge Field Research Center (FRC) site. Here, we report the complete genome sequence and annotation of Cupriavidus basilensis 4G11. The genome contains 8,421,483 bp, 7,661 predicted protein-coding genes, and a total GC content of 64.4%.