Karan Uppal

The Gene Expression Barcode: leveraging public data repositories to begin cataloging the human and murine transcriptomes

Various databases have harnessed the wealth of publicly available microarray data to address biological questions ranging from across-tissue differential expression to homologous gene expression. Despite their practical value, these databases rely on relative measures of expression and are unable to address the most fundamental question—which genes are expressed in a given cell type. The Gene Expression Barcode is the first database to provide reliable absolute measures of expression for most annotated genes for 131 human and 89 mouse tissue types, including diseased tissue. This is made possible by a novel algorithm that leverages information from the GEO and ArrayExpress public repositories to build statistical models that permit converting data from a single microarray into expressed/unexpressed calls for each gene. For selected platforms, users may upload data and obtain results in a matter of seconds. The raw data, curated annotation, and code used to create our resource are also available at http://rafalab.jhsph.edu/barcode.

xMSanalyzer: automated pipeline for improved feature detection and downstream analysis of large-scale, non-targeted metabolomics data

BackgroundDetection of low abundance metabolites is important for de novo mapping of metabolic pathways related to diet, microbiome or environmental exposures. Multiple algorithms are available to extract m/z features from liquid chromatography-mass spectral data in a conservative manner, which tends to preclude detection of low abundance chemicals and chemicals found in small subsets of samples. The present study provides software to enhance such algorithms for feature detection, quality assessment, and annotation.ResultsxMSanalyzer is a set of utilities for automated processing of metabolomics data. The utilites can be classified into four main modules to: 1) improve feature detection for replicate analyses by systematic re-extraction with multiple parameter settings and data merger to optimize the balance between sensitivity and reliability, 2) evaluate sample quality and feature consistency, 3) detect feature overlap between datasets, and 4) characterize high-resolution m/z matches to small molecule metabolites and biological pathways using multiple chemical databases. The package was tested with plasma samples and shown to more than double the number of features extracted while improving quantitative reliability of detection. MS/MS analysis of a random subset of peaks that were exclusively detected using xMSanalyzer confirmed that the optimization scheme improves detection of real metabolites.ConclusionsxMSanalyzer is a package of utilities for data extraction, quality control assessment, detection of overlapping and unique metabolites in multiple datasets, and batch annotation of metabolites. The program was designed to integrate with existing packages such as apLCMS and XCMS, but the framework can also be used to enhance data extraction for other LC/MS data software.

Plasma Metabolomics in Human Pulmonary Tuberculosis Disease: A Pilot Study

We aimed to characterize metabolites during tuberculosis (TB) disease and identify new pathophysiologic pathways involved in infection as well as biomarkers of TB onset, progression and resolution. Such data may inform development of new anti-tuberculosis drugs. Plasma samples from adults with newly diagnosed pulmonary TB disease and their matched, asymptomatic, sputum culture-negative household contacts were analyzed using liquid chromatography high-resolution mass spectrometry (LC-MS) to identify metabolites. Statistical and bioinformatics methods were used to select accurate mass/charge (m/z) ions that were significantly different between the two groups at a false discovery rate (FDR) of q<0.05. Two-way hierarchical cluster analysis (HCA) was used to identify clusters of ions contributing to separation of cases and controls, and metabolomics databases were used to match these ions to known metabolites. Identity of specific D-series resolvins, glutamate and Mycobacterium tuberculosis (Mtb)-derived trehalose-6-mycolate was confirmed using LC-MS/MS analysis. Over 23,000 metabolites were detected in untargeted metabolomic analysis and 61 metabolites were significantly different between the two groups. HCA revealed 8 metabolite clusters containing metabolites largely upregulated in patients with TB disease, including anti-TB drugs, glutamate, choline derivatives, Mycobacterium tuberculosis-derived cell wall glycolipids (trehalose-6-mycolate and phosphatidylinositol) and pro-resolving lipid mediators of inflammation, known to stimulate resolution, efferocytosis and microbial killing. The resolvins were confirmed to be RvD1, aspirin-triggered RvD1, and RvD2. This study shows that high-resolution metabolomic analysis can differentiate patients with active TB disease from their asymptomatic household contacts. Specific metabolites upregulated in the plasma of patients with active TB disease, including Mtb-derived glycolipids and resolvins, have potential as biomarkers and may reveal pathways involved in TB disease pathogenesis and resolution.

Metabolomics of ADSOL (AS-1) Red Blood Cell Storage

Population-based investigations suggest that red blood cells (RBCs) are therapeutically effective when collected, processed, and stored for up to 42 days under validated conditions before transfusion. However, some retrospective clinical studies have shown worse patient outcomes when transfused RBCs have been stored for the longest times. Furthermore, studies of RBC persistence in the circulation after transfusion have suggested that considerable donor-to-donor variability exists and may affect transfusion efficacy. To understand the limitations of current blood storage technologies and to develop approaches to improve RBC storage and transfusion efficacy, we investigated the global metabolic alterations that occur when RBCs are stored in AS-1 (AS1-RBC). Leukoreduced AS1-RBC units prepared from 9 volunteer research donors (12 total donated units) were serially sampled for metabolomics analysis over 42 days of refrigerated storage. Samples were tested by gas chromatography/mass spectrometry and liquid chromatography/tandem mass spectrometry, and specific biochemical compounds were identified by comparison to a library of purified standards. Over 3 experiments, 185 to 264 defined metabolites were quantified in stored RBC samples. Kinetic changes in these biochemicals confirmed known alterations in glycolysis and other pathways previously identified in RBCs stored in saline, adenine, glucose and mannitol solution (SAGM-RBC). Furthermore, we identified additional alterations not previously seen in SAGM-RBCs (eg, stable pentose phosphate pathway flux, progressive decreases in oxidized glutathione), and we delineated changes occurring in other metabolic pathways not previously studied (eg, S-adenosyl methionine cycle). These data are presented in the context of a detailed comparison with previous studies of SAGM-RBCs from human donors and murine AS1-RBCs. Global metabolic profiling of AS1-RBCs revealed a number of biochemical alterations in stored blood that may affect RBC viability during storage as well as therapeutic effectiveness of stored RBCs in transfusion recipients. These results provide future opportunities to more clearly pinpoint the metabolic defects during RBC storage, to identify biomarkers for donor screening and prerelease RBC testing, and to develop improved RBC storage solutions and methodologies.

Metabolome-Wide Association Study of Neovascular Age-Related Macular Degeneration

Purpose To determine if plasma metabolic profiles can detect differences between patients with neovascular age-related macular degeneration (NVAMD) and similarly-aged controls. Methods Metabolomic analysis using liquid chromatography with Fourier-transform mass spectrometry (LC-FTMS) was performed on plasma samples from 26 NVAMD patients and 19 controls. Data were collected from mass/charge ratio (m/z) 85 to 850 on a Thermo LTQ-FT mass spectrometer, and metabolic features were extracted using an adaptive processing software package. Both non-transformed and log2 transformed data were corrected using Benjamini and Hochberg False Discovery Rate (FDR) to account for multiple testing. Orthogonal Partial Least Squares-Discriminant Analysis was performed to determine metabolic features that distinguished NVAMD patients from controls. Individual m/z features were matched to the Kyoto Encyclopedia of Genes and Genomes database and the Metlin metabolomics database, and metabolic pathways associated with NVAMD were identified using MetScape. Results Of the 1680 total m/z features detected by LC-FTMS, 94 unique m/z features were significantly different between NVAMD patients and controls using FDR (q = 0.05). A comparison of these features to those found with log2 transformed data (n = 132, q = 0.2) revealed 40 features in common, reaffirming the involvement of certain metabolites. Such metabolites included di- and tripeptides, covalently modified amino acids, bile acids, and vitamin D-related metabolites. Correlation analysis revealed associations among certain significant features, and pathway analysis demonstrated broader changes in tyrosine metabolism, sulfur amino acid metabolism, and amino acids related to urea metabolism. Conclusions These data suggest that metabolomic analysis can identify a panel of individual metabolites that differ between NVAMD cases and controls. Pathway analysis can assess the involvement of certain metabolic pathways, such as tyrosine and urea metabolism, and can provide further insight into the pathophysiology of AMD.

Serum Metabolomics of Slow vs. Rapid Motor Progression Parkinson’s Disease: a Pilot Study

Progression of Parkinson’s disease (PD) is highly variable, indicating that differences between slow and rapid progression forms could provide valuable information for improved early detection and management. Unfortunately, this represents a complex problem due to the heterogeneous nature of humans in regards to demographic characteristics, genetics, diet, environmental exposures and health behaviors. In this pilot study, we employed high resolution mass spectrometry-based metabolic profiling to investigate the metabolic signatures of slow versus rapidly progressing PD present in human serum. Archival serum samples from PD patients obtained within 3 years of disease onset were analyzed via dual chromatography-high resolution mass spectrometry, with data extraction by xMSanalyzer and used to predict rapid or slow motor progression of these patients during follow-up. Statistical analyses, such as false discovery rate analysis and partial least squares discriminant analysis, yielded a list of statistically significant metabolic features and further investigation revealed potential biomarkers. In particular, N8-acetyl spermidine was found to be significantly elevated in the rapid progressors compared to both control subjects and slow progressors. Our exploratory data indicate that a fast motor progression disease phenotype can be distinguished early in disease using high resolution mass spectrometry-based metabolic profiling and that altered polyamine metabolism may be a predictive marker of rapidly progressing PD.

Reference Standardization for Mass Spectrometry and High-resolution Metabolomics Applications to Exposome Research

The exposome is the cumulative measure of environmental influences and associated biological responses throughout the lifespan, including exposures from the environment, diet, behavior, and endogenous processes. A major challenge for exposome research lies in the development of robust and affordable analytic procedures to measure the broad range of exposures and associated biologic impacts occurring over a lifetime. Biomonitoring is an established approach to evaluate internal body burden of environmental exposures, but use of biomonitoring for exposome research is often limited by the high costs associated with quantification of individual chemicals. High-resolution metabolomics (HRM) uses ultra-high resolution mass spectrometry with minimal sample preparation to support high-throughput relative quantification of thousands of environmental, dietary, and microbial chemicals. HRM also measures metabolites in most endogenous metabolic pathways, thereby providing simultaneous measurement of biologic responses to environmental exposures. The present research examined quantification strategies to enhance the usefulness of HRM data for cumulative exposome research. The results provide a simple reference standardization protocol in which individual chemical concentrations in unknown samples are estimated by comparison to a concurrently analyzed, pooled reference sample with known chemical concentrations. The approach was tested using blinded analyses of amino acids in human samples and was found to be comparable to independent laboratory results based on surrogate standardization or internal standardization. Quantification was reproducible over a 13-month period and extrapolated to thousands of chemicals. The results show that reference standardization protocol provides an effective strategy that will enhance data collection for cumulative exposome research. In principle, the approach can be extended to other types of mass spectrometry and other analytical methods.

Correlation of the lung microbiota with metabolic profiles in bronchoalveolar lavage fluid in HIV infection

BackgroundWhile 16S ribosomal RNA (rRNA) sequencing has been used to characterize the lung’s bacterial microbiota in human immunodeficiency virus (HIV)-infected individuals, taxonomic studies provide limited information on bacterial function and impact on the host. Metabolic profiles can provide functional information on host-microbe interactions in the lungs. We investigated the relationship between the respiratory microbiota and metabolic profiles in the bronchoalveolar lavage fluid of HIV-infected and HIV-uninfected outpatients.ResultsTargeted sequencing of the 16S rRNA gene was used to analyze the bacterial community structure and liquid chromatography-high-resolution mass spectrometry was used to detect features in bronchoalveolar lavage fluid. Global integration of all metabolic features with microbial species was done using sparse partial least squares regression. Thirty-nine HIV-infected subjects and 20 HIV-uninfected controls without acute respiratory symptoms were enrolled. Twelve mass-to-charge ratio (m/z) features from C18 analysis were significantly different between HIV-infected individuals and controls (false discovery rate (FDR) = 0.2); another 79 features were identified by network analysis. Further metabolite analysis demonstrated that four features were significantly overrepresented in the bronchoalveolar lavage (BAL) fluid of HIV-infected individuals compared to HIV-uninfected, including cystine, two complex carbohydrates, and 3,5-dibromo-l-tyrosine. There were 231 m/z features significantly associated with peripheral blood CD4 cell counts identified using sparse partial least squares regression (sPLS) at a variable importance on projection (VIP) threshold of 2. Twenty-five percent of these 91 m/z features were associated with various microbial species. Bacteria from families Caulobacteraceae, Staphylococcaceae, Nocardioidaceae, and genus Streptococcus were associated with the greatest number of features. Glycerophospholipid and lineolate pathways correlated with these bacteria.ConclusionsIn bronchoalveolar lavage fluid, specific metabolic profiles correlated with bacterial organisms known to play a role in the pathogenesis of pneumonia in HIV-infected individuals. These findings suggest that microbial communities and their interactions with the host may have functional metabolic impact in the lung.

Novel metabolic markers for the risk of diabetes development in American Indians

OBJECTIVE To identify novel metabolic markers for diabetes development in American Indians. RESEARCH DESIGN AND METHODS Using an untargeted high-resolution liquid chromatography–mass spectrometry, we conducted metabolomics analysis of study participants who developed incident diabetes (n = 133) and those who did not (n = 298) from 2,117 normoglycemic American Indians followed for an average of 5.5 years in the Strong Heart Family Study. Relative abundances of metabolites were quantified in baseline fasting plasma of all 431 participants. Prospective association of each metabolite with risk of developing type 2 diabetes (T2D) was examined using logistic regression adjusting for established diabetes risk factors. RESULTS Seven metabolites (five known and two unknown) significantly predict the risk of T2D. Notably, one metabolite matching 2-hydroxybiphenyl was significantly associated with an increased risk of diabetes, whereas four metabolites matching PC (22:6/20:4), (3S)-7-hydroxy-2′,3′,4′,5′,8-pentamethoxyisoflavan, or tetrapeptides were significantly associated with decreased risk of diabetes. A multimarker score comprising all seven metabolites significantly improved risk prediction beyond established diabetes risk factors including BMI, fasting glucose, and insulin resistance. CONCLUSIONS The findings suggest that these newly detected metabolites may represent novel prognostic markers of T2D in American Indians, a group suffering from a disproportionately high rate of T2D.

MetabNet: An R Package for Metabolic Association Analysis of High-Resolution Metabolomics Data.

Liquid-chromatography high-resolution mass spectrometry provides capability to measure >40,000 ions derived from metabolites in biologic samples. This presents challenges to confirm identities of known chemicals and delineate potential metabolic pathway associations of unidentified chemicals. We provide an R package for metabolic network analysis, MetabNet, to perform targeted metabolome-wide association study of specific metabolites to facilitate detection of their related metabolic pathways and network structures.