A Unified Conceptual Framework for Metabolic Phenotyping in Diagnosis and Prognosis.

Understanding metabotype (multicomponent metabolic characteristics) variation can help to generate new diagnostic and prognostic biomarkers, as well as models, with potential to impact on patient management. We present a suite of conceptual approaches for the generation, analysis, and understanding of metabotypes from body fluids and tissues. We describe and exemplify four fundamental approaches to the generation and utilization of metabotype data via multiparametric measurement of (i) metabolite levels, (ii) metabolic trajectories, (iii) metabolic entropies, and (iv) metabolic networks and correlations in space and time. This conceptual framework can underpin metabotyping in the scenario of personalized medicine, with the aim of improving clinical outcomes for patients, but the framework will have value and utility in areas of metabolic profiling well beyond this exemplar.

ChemDistiller: an engine for metabolite annotation in mass spectrometry.

Motivation: High-resolution mass spectrometry permits simultaneous detection of thousands of different metabolites in biological samples; however, their automated annotation still presents a challenge due to the limited number of tailored computational solutions freely available to the scientific community. Results: Here, we introduce ChemDistiller, a customizable engine that combines automated large-scale annotation of metabolites using tandem MS data with a compiled database containing tens of millions of compounds with pre-calculated 'fingerprints' and fragmentation patterns. Our tests using publicly and commercially available tandem MS spectra for reference compounds show retrievals rates comparable to or exceeding the ones obtainable by the current state-of-the-art solutions in the field while offering higher throughput, scalability and processing speed. Availability and implementation: Source code freely available for download at https://bitbucket.org/iAnalytica/chemdistillerpython. Supplementary information: Supplementary data are available at Bioinformatics online.

Translational utility of a hierarchical classification strategy in biomolecular data analytics.

Hierarchical classification (HC) stratifies and classifies data from broad classes into more specific classes. Unlike commonly used data classification strategies, this enables the probabilistic prediction of unknown classes at different levels, minimizing the burden of incomplete databases. Despite these advantages, its translational application in biomedical sciences has been limited. We describe and demonstrate the implementation of a HC approach for "omics-driven" classification of 15 bacterial species at various taxonomic levels achieving 90-100% accuracy, and 9 cancer types into morphological types and 35 subtypes with 99% and 76% accuracy, respectively. Unknown bacterial species were probabilistically assigned with 100% accuracy to their respective genus or family using mass spectra (n = 284). Cancer types were predicted by mRNA data (n = 1960) for most subtypes with 95-100% accuracy. This has high relevance in clinical practice where complete datasets are difficult to compile with the continuous evolution of diseases and emergence of new strains, yet prediction of unknown classes, such as bacterial species, at upper hierarchy levels may be sufficient to initiate antimicrobial therapy. The algorithms presented here can be directly translated into clinical-use with any quantitative data, and have broad application potential, from unlabeled sample identification, to hierarchical feature selection, and discovery of new taxonomic variants.

Investigation of the Impact of Desorption Electrospray Ionization Sprayer Geometry on Its Performance in Imaging of Biological Tissue.

In this study, the impact of sprayer design and geometry on performance in desorption electrospray ionization mass spectrometry (DESI-MS) is assessed, as the sprayer is thought to be a major source of variability. Absolute intensity repeatability, spectral composition, and classification accuracy for biological tissues are considered. Marked differences in tissue analysis performance are seen between the commercially available and a lab-built sprayer. These are thought to be associated with the geometry of the solvent capillary and the resulting shape of the primary electrospray. Experiments with a sprayer with a fixed solvent capillary position show that capillary orientation has a crucial impact on tissue complex lipid signal and can lead to an almost complete loss of signal. Absolute intensity repeatability is compared for five lab-built sprayers using pork liver sections. Repeatability ranges from 1 to 224% for individual sprayers and peaks of different spectral abundance. Between sprayers, repeatability is 16%, 9%, 23%, and 34% for high, medium, low, and very low abundance peaks, respectively. To assess the impact of sprayer variability on tissue classification using multivariate statistical tools, nine human colorectal adenocarcinoma sections are analyzed with three lab-built sprayers, and classification accuracy for adenocarcinoma versus the surrounding stroma is assessed. It ranges from 80.7 to 94.5% between the three sprayers and is 86.5% overall. The presented results confirm that the sprayer setup needs to be closely controlled to obtain reliable data, and a new sprayer setup with a fixed solvent capillary geometry should be developed.

Mass Spectrometric Analysis of Exhaled Breath for the Identification of Volatile Organic Compound Biomarkers in Esophageal and Gastric Adenocarcinoma.

OBJECTIVE: The present study assessed whether exhaled breath analysis using Selected Ion Flow Tube Mass Spectrometry could distinguish esophageal and gastric adenocarcinoma from noncancer controls. BACKGROUND: The majority of patients with upper gastrointestinal cancer present with advanced disease, resulting in poor long-term survival rates. Novel methods are needed to diagnose potentially curable upper gastrointestinal malignancies. METHODS: A Profile-3 Selected Ion Flow Tube Mass Spectrometry instrument was used for analysis of volatile organic compounds (VOCs) within exhaled breath samples. All study participants had undergone upper gastrointestinal endoscopy on the day of breath sampling. Receiver operating characteristic analysis and a diagnostic risk prediction model were used to assess the discriminatory accuracy of the identified VOCs. RESULTS: Exhaled breath samples were analyzed from 81 patients with esophageal (N = 48) or gastric adenocarcinoma (N = 33) and 129 controls including Barrett's metaplasia (N = 16), benign upper gastrointestinal diseases (N = 62), or a normal upper gastrointestinal tract (N = 51). Twelve VOCs-pentanoic acid, hexanoic acid, phenol, methyl phenol, ethyl phenol, butanal, pentanal, hexanal, heptanal, octanal, nonanal, and decanal-were present at significantly higher concentrations (P < 0.05) in the cancer groups than in the noncancer controls. The area under the ROC curve using these significant VOCs to discriminate esophageal and gastric adenocarcinoma from those with normal upper gastrointestinal tracts was 0.97 and 0.98, respectively. The area under the ROC curve for the model and validation subsets of the diagnostic prediction model was 0.92 ±â€Š0.01 and 0.87 ±â€Š0.03, respectively. CONCLUSIONS: Distinct exhaled breath VOC profiles can distinguish patients with esophageal and gastric adenocarcinoma from noncancer controls.

Metabolic, immune, and gut microbial signals mount a systems response to Leishmania major infection.

Parasitic infections such as leishmaniasis induce a cascade of host physiological responses, including metabolic and immunological changes. Infection with Leishmania major parasites causes cutaneous leishmaniasis in humans, a neglected tropical disease that is difficult to manage. To understand the determinants of pathology, we studied L. major infection in two mouse models: the self-healing C57BL/6 strain and the nonhealing BALB/c strain. Metabolic profiling of urine, plasma, and feces via proton NMR spectroscopy was performed to discover parasite-specific imprints on global host metabolism. Plasma cytokine status and fecal microbiome were also characterized as additional metrics of the host response to infection. Results demonstrated differences in glucose and lipid metabolism, distinctive immunological phenotypes, and shifts in microbial composition between the two models. We present a novel approach to integrate such metrics using correlation network analyses, whereby self-healing mice demonstrated an orchestrated interaction between the biological measures shortly after infection. In contrast, the response observed in nonhealing mice was delayed and fragmented. Our study suggests that trans-system communication across host metabolism, the innate immune system, and gut microbiome is key for a successful host response to L. major and provides a new concept, potentially translatable to other diseases.

Characterization and identification of clinically relevant microorganisms using rapid evaporative ionization mass spectrometry.

Rapid evaporative ionization mass spectrometry (REIMS) was investigated for its suitability as a general identification system for bacteria and fungi. Strains of 28 clinically relevant bacterial species were analyzed in negative ion mode, and corresponding data was subjected to unsupervised and supervised multivariate statistical analyses. The created supervised model yielded correct cross-validation results of 95.9%, 97.8%, and 100% on species, genus, and Gram-stain level, respectively. These results were not affected by the resolution of the mass spectral data. Blind identification tests were performed for strains cultured on different culture media and analyzed using different instrumental platforms which led to 97.8-100% correct identification. Seven different Escherichia coli strains were subjected to different culture conditions and were distinguishable with 88% accuracy. In addition, the technique proved suitable to distinguish five pathogenic Candida species with 98.8% accuracy without any further modification to the experimental workflow. These results prove that REIMS is sufficiently specific to serve as a culture condition-independent tool for the identification and characterization of microorganisms.

Chemo-informatic strategy for imaging mass spectrometry-based hyperspectral profiling of lipid signatures in colorectal cancer.

Mass spectrometry imaging (MSI) provides the opportunity to investigate tumor biology from an entirely novel biochemical perspective and could lead to the identification of a new pool of cancer biomarkers. Effective clinical translation of histology-driven MSI in systems oncology requires precise colocalization of morphological and biochemical features as well as advanced methods for data treatment and interrogation. Currently proposed MSI workflows are subject to several limitations, including nonoptimized raw data preprocessing, imprecise image coregistration, and limited pattern recognition capabilities. Here we outline a comprehensive strategy for histology-driven MSI, using desorption electrospray ionization that covers (i) optimized data preprocessing for improved information recovery; (ii) precise image coregistration; and (iii) efficient extraction of tissue-specific molecular ion signatures for enhanced biochemical distinction of different tissue types. The proposed workflow has been used to investigate region-specific lipid signatures in colorectal cancer tissue. Unique lipid patterns were observed using this approach according to tissue type, and a tissue recognition system using multivariate molecular ion patterns allowed highly accurate (>98%) identification of pixels according to morphology (cancer, healthy mucosa, smooth muscle, and microvasculature). This strategy offers unique insights into tumor microenvironmental biochemistry and should facilitate compilation of a large-scale tissue morphology-specific MSI spectral database with which to pursue next-generation, fully automated histological approaches.

Analysis of intact bacteria using rapid evaporative ionisation mass spectrometry.

An identification system for microorganisms based on recently developed rapid evaporative ionisation mass spectrometry (REIMS) is presented. Nine bacterial species cultured on various growth media were correctly identified to family-, genus-, and species-level based on their different mass spectral fingerprints using a cross-validated maximum margin criterion model.

Optimized preprocessing of ultra-performance liquid chromatography/mass spectrometry urinary metabolic profiles for improved information recovery.

Ultra-performance liquid chromatography coupled to mass spectrometry (UPLC/MS) has been used increasingly for measuring changes of low molecular weight metabolites in biofluids/tissues in response to biological challenges such as drug toxicity and disease processes. Typically samples show high variability in concentration, and the derived metabolic profiles have a heteroscedastic noise structure characterized by increasing variance as a function of increased signal intensity. These sources of experimental and instrumental noise substantially complicate information recovery when statistical tools are used. We apply and compare several preprocessing procedures and introduce a statistical error model to account for these bioanalytical complexities. In particular, the use of total intensity, median fold change, locally weighted scatter plot smoothing, and quantile normalizations to reduce extraneous variance induced by sample dilution were compared. We demonstrate that the UPLC/MS peak intensities of urine samples should respond linearly to variable sample dilution across the intensity range. While all four studied normalization methods performed reasonably well in reducing dilution-induced variation of urine samples in the absence of biological variation, the median fold change normalization is least compromised by the biologically relevant changes in mixture components and is thus preferable. Additionally, the application of a subsequent log-based transformation was successful in stabilizing the variance with respect to peak intensity, confirming the predominant influence of multiplicative noise in peak intensities from UPLC/MS-derived metabolic profile data sets. We demonstrate that variance-stabilizing transformation and normalization are critical preprocessing steps that can benefit greatly metabolic information recovery from such data sets when widely applied chemometric methods are used.

A metabolic entropy approach for measurements of systemic metabolic disruptions in patho-physiological States.

Multicellular organisms maintain the stability of their internal environment using metabolic and physiological regulatory mechanisms that are disrupted during disease. The loss of homeostatic control results in a complex set of disordered states that may lead to metabolic network failure and irreversible system damage. We have applied a new statistical entropy-based approach to quantify temporal systemic disorder (divergence of metabolic responses) in experimental patho-physiological states, via NMR-spectroscopy generated metabolic profiles of urine. A recovery (R-) potential metric has also been developed to evaluate the relative extent to which defined metabolic processes are perturbed in the context of a global system in terms of multiple changes in concentrations of biofluid components accompanying the disrupted functional activity. This approach is sensitive to physiological as well as pathological interventions. We show that global disruptions of metabolic processes, lesion reversibility, and disorder in metabolic responses to a stressor can be visualized via metabolic entropy metrics, giving insights into biological robustness and thus providing a new tool for assessing deviation from homeostatic regulation.

Urinary metabolic phenotyping differentiates children with autism from their unaffected siblings and age-matched controls.

Autism is an early onset developmental disorder with a severe life-long impact on behavior and social functioning that has associated metabolic abnormalities. The urinary metabolic phenotypes of individuals (age range=3-9 years old) diagnosed with autism using the DSM-IV-TR criteria (n = 39; male = 35; female = 4), together with their nonautistic siblings (n = 28; male = 14; female = 14) and age-matched healthy volunteers (n = 34, male = 17; female = 17) have been characterized for the first time using (1)H NMR spectroscopy and pattern recognition methods. Novel findings associated with alterations in nicotinic acid metabolism within autistic individuals showing increased urinary excretion of N-methyl-2-pyridone-5-carboxamide, N-methyl nicotinic acid, and N-methyl nicotinamide indicate a perturbation in the tryptophan-nicotinic acid metabolic pathway. Multivariate statistical analysis indicated urinary patterns of the free amino acids, glutamate and taurine were significantly different between groups with the autistic children showing higher levels of urinary taurine and a lower level of urinary glutamate, indicating perturbation in sulfur and amino acid metabolism in these children. Additionally, metabolic phenotype (metabotype) differences were observed between autistic and control children, which were associated with perturbations in the relative patterns of urinary mammalian-microbial cometabolites including dimethylamine, hippurate, and phenyacetylglutamine. These biochemical changes are consistent with some of the known abnormalities of gut microbiota found in autistic individuals and the associated gastrointestinal dysfunction and may be of value in monitoring the success of therapeutic interventions.

Cluster analysis statistical spectroscopy using nuclear magnetic resonance generated metabolic data sets from perturbed biological systems.

We present a new approach for analysis, information recovery, and display of biological (1)H nuclear magnetic resonance (NMR) spectral data, cluster analysis statistical spectroscopy (CLASSY), which profiles qualitative and quantitative changes in biofluid metabolic composition by utilizing a novel local-global correlation clustering scheme to identify structurally related spectral peaks and arrange metabolites by similarity of temporal dynamic variation. Underlying spectral data sets are presented in a novel graphical format to represent high-dimensionality biochemical information conveying both statistical metabolite relationships and their responses to experimental perturbation simultaneously in a high-throughput and intuitive manner. The method is exemplified using multiple 600 MHz (1)H NMR spectra of rat (n = 40) urine samples collected over 160 h following the development of experimental pancreatitis induced by L-arginine (ARG) and a wider range of model toxins including acetaminophen, galactosamine, and 2-bromoethanamine. The CLASSY approach deconvolutes complex biofluid mixture spectra into quantitative fold-change metabolic trajectories and clusters metabolites by commonalities of coexpression patterns. We demonstrate that the developing pathological processes cause coordinated changes in the levels of many compounds which share similar pathway connectivities. Variability in individual responses to toxin exposure is also readily detected and visualized allowing the assessment of interanimal variability. As an untargeted, unsupervised approach, CLASSY provides significant advantages in biological information recovery in terms of increased throughput, interpretability, and robustness and has wide potential metabonomic/metabolomic applications in clinical, toxicological, and nutritional studies of biofluids as well as in studies of cellular biochemistry, microbial fermentation monitoring, and functional genomics.

Recursive segment-wise peak alignment of biological (1)h NMR spectra for improved metabolic biomarker recovery.

Chemical shift variation in small-molecule (1)H NMR signals of biofluids complicates biomarker information recovery in metabonomic studies when using multivariate statistical and pattern recognition tools. Current peak realignment methods are generally time-consuming or align major peaks at the expense of minor peak shift accuracy. We present a novel recursive segment-wise peak alignment (RSPA) method to reduce variability in peak positions across the multiple (1)H NMR spectra used in metabonomic studies. The method refines a segmentation of reference and test spectra in a top-down fashion, sequentially subdividing the initial larger segments, as required, to improve the local spectral alignment. We also describe a general procedure that allows robust comparison of realignment quality of various available methods for a range of peak intensities. The RSPA method is illustrated with respect to 140 (1)H NMR rat urine spectra from a caloric restriction study and is compared with several other widely used peak alignment methods. We demonstrate the superior performance of the RSPA alignment over a wide range of peaks and its capacity to enhance interpretability and robustness of multivariate statistical tools. The approach is widely applicable for NMR-based metabolic studies and is potentially suitable for many other types of data sets such as chromatographic profiles and MS data.

Human metabolic phenotype diversity and its association with diet and blood pressure.

Metabolic phenotypes are the products of interactions among a variety of factors-dietary, other lifestyle/environmental, gut microbial and genetic. We use a large-scale exploratory analytical approach to investigate metabolic phenotype variation across and within four human populations, based on 1H NMR spectroscopy. Metabolites discriminating across populations are then linked to data for individuals on blood pressure, a major risk factor for coronary heart disease and stroke (leading causes of mortality worldwide). We analyse spectra from two 24-hour urine specimens for each of 4,630 participants from the INTERMAP epidemiological study, involving 17 population samples aged 40-59 in China, Japan, UK and USA. We show that urinary metabolite excretion patterns for East Asian and western population samples, with contrasting diets, diet-related major risk factors, and coronary heart disease/stroke rates, are significantly differentiated (P < 10(-16)), as are Chinese/Japanese metabolic phenotypes, and subgroups with differences in dietary vegetable/animal protein and blood pressure. Among discriminatory metabolites, we quantify four and show association (P < 0.05 to P < 0.0001) of mean 24-hour urinary formate excretion with blood pressure in multiple regression analyses for individuals. Mean 24-hour urinary excretion of alanine (direct) and hippurate (inverse), reflecting diet and gut microbial activities, are also associated with blood pressure of individuals. Metabolic phenotyping applied to high-quality epidemiological data offers the potential to develop an area of aetiopathogenetic knowledge involving discovery of novel biomarkers related to cardiovascular disease risk.