Timothy M. D. Ebbels

Metabolic profiling, metabolomic and metabonomic procedures for NMR spectroscopy of urine, plasma, serum and tissue extracts

Metabolic profiling, metabolomic and metabonomic studies mainly involve the multicomponent analysis of biological fluids, tissue and cell extracts using NMR spectroscopy and/or mass spectrometry (MS). We summarize the main NMR spectroscopic applications in modern metabolic research, and provide detailed protocols for biofluid (urine, serum/plasma) and tissue sample collection and preparation, including the extraction of polar and lipophilic metabolites from tissues. 1H NMR spectroscopic techniques such as standard 1D spectroscopy, relaxation-edited, diffusion-edited and 2D J-resolved pulse sequences are widely used at the analysis stage to monitor different groups of metabolites and are described here. They are often followed by more detailed statistical analysis or additional 2D NMR analysis for biomarker discovery. The standard acquisition time per sample is 4–5 min for a simple 1D spectrum, and both preparation and analysis can be automated to allow application to high-throughput screening for clinical diagnostic and toxicological studies, as well as molecular phenotyping and functional genomics.

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.

Contemporary issues in toxicology - The role of metabonomics in toxicology and its evaluation by the COMET project

The role that metabonomics has in the evaluation of xenobiotic toxicity studies is presented here together with a brief summary of published studies. To provide a comprehensive assessment of this approach, the Consortium for Metabonomic Toxicology (COMET) has been formed between six pharmaceutical companies and Imperial College of Science, Technology and Medicine (IC), London, UK. The objective of this group is to define methodologies and to apply metabonomic data generated using (1)H NMR spectroscopy of urine and blood serum for preclinical toxicological screening of candidate drugs. This is being achieved by generating databases of results for a wide range of model toxins which serve as the raw material for computer-based expert systems for toxicity prediction. The project progress on the generation of comprehensive metabonomic databases and multivariate statistical models for prediction of toxicity, initially for liver and kidney toxicity in the rat and mouse, is reported. Additionally, both the analytical and biological variation which might arise through the use of metabonomics has been evaluated. An evaluation of intersite NMR analytical reproducibility has revealed a high degree of robustness. Second, a detailed comparison has been made of the ability of the six companies to provide consistent urine and serum samples using a study of the toxicity of hydrazine at two doses in the male rat, this study showing a high degree of consistency between samples from the various companies in terms of spectral patterns and biochemical composition. Differences between samples from the various companies were small compared to the biochemical effects of the toxin. A metabonomic model has been constructed for urine from control rats, enabling identification of outlier samples and the metabolic reasons for the deviation. Building on this success, and with the completion of studies on approximately 80 model toxins, first expert systems for prediction of liver and kidney toxicity have been generated.

Proposed minimum reporting standards for data analysis in metabolomics

The goal of this group is to define the reporting requirements associated with the statistical analysis (including univariate, multivariate, informatics, machine learning etc.) of metabolite data with respect to other measured/collected experimental data (often called meta-data). These definitions will embrace as many aspects of a complete metabolomics study as possible at this time. In chronological order this will include: Experimental Design, both in terms of sample collection/matching, and data acquisition scheduling of samples through whichever spectroscopic technology used; Deconvolution (if required); Pre-processing, for example, data cleaning, outlier detection, row/column scaling, or other transformations; Definition and parameterization of subsequent visualizations and Statistical/Machine learning Methods applied to the dataset; If required, a clear definition of the Model Validation Scheme used (including how data are split into training/validation/test sets); Formal indication on whether the data analysis has been Independently Tested (either by experimental reproduction, or blind hold out test set). Finally, data interpretation and the visual representations and hypotheses obtained from the data analyses.

High-resolution magic-angle-spinning NMR spectroscopy for metabolic profiling of intact tissues

Metabolic profiling, metabolomic and metabonomic studies require robust study protocols for any large-scale comparisons and evaluations. Detailed methods for solution-state NMR spectroscopy have been summarized in an earlier protocol. This protocol details the analysis of intact tissue samples by means of high-resolution magic-angle-spinning (HR-MAS) NMR spectroscopy and we provide a detailed description of sample collection, preparation and analysis. Described here are 1H NMR spectroscopic techniques such as the standard one-dimensional, relaxation-edited, diffusion-edited and two-dimensional J-resolved pulse experiments, as well as one-dimensional 31P NMR spectroscopy. These are used to monitor different groups of metabolites, e.g., sugars, amino acids and osmolytes as well as larger molecules such as lipids, non-invasively. Through the use of NMR-based diffusion coefficient and relaxation times measurements, information on molecular compartmentation and mobility can be gleaned. The NMR methods are often combined with statistical analysis for further metabonomics analysis and biomarker identification. The standard acquisition time per sample is 8–10 min for a simple one-dimensional 1H NMR spectrum, giving access to metabolite information while retaining tissue integrity and hence allowing direct comparison with histopathology and MRI/MRS findings or the evaluation together with biofluid metabolic-profiling data.

Recursive Segment-Wise Peak Alignment of Biological 1H 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.

Improved analysis of multivariate data by variable stability scaling: application to NMR-based metabolic profiling

Abstract Variable scaling alters the covariance structure of data, affecting the outcome of multivariate analysis and calibration. Here we present a new method, variable stability (VAST) scaling, which weights each variable according to a metric of its stability. The beneficial effect of VAST scaling is demonstrated for a data set of 1 H NMR spectra of urine acquired as part of a metabonomic study into the effects of unilateral nephrectomy in an animal model. The application of VAST scaling improved the class distinction and predictive power of partial least squares discriminant analysis (PLS-DA) models. The effects of other data scaling and pre-processing methods, such as orthogonal signal correction (OSC), were also tested. VAST scaling produced the most robust models in terms of class prediction, outperforming OSC in this aspect. As a result the subtle, but consistent, metabolic perturbation caused by unilateral nephrectomy could be accurately characterised despite the presence of much greater biological differences caused by normal physiological variation. VAST scaling presents itself as an interpretable, robust and easily implemented data treatment for the enhancement of multivariate data analysis.

Optimization and Evaluation of Metabolite Extraction Protocols for Untargeted Metabolic Profiling of Liver Samples by UPLC-MS

A series of six protocols were evaluated for UPLC-MS based untargeted metabolic profiling of liver extracts in terms of reproducibility and number of metabolite features obtained. These protocols, designed to extract both polar and nonpolar metabolites, were based on (i) a two stage extraction approach or (ii) a simultaneous extraction in a biphasic mixture, employing different volumes and combinations of extraction and resuspension solvents. A multivariate statistical strategy was developed to allow comparison of the multidimensional variation between the methods. The optimal protocol for profiling both polar and nonpolar metabolites was found to be an aqueous extraction with methanol/water followed by an organic extraction with dichloromethane/methanol, with resuspension of the dried extracts in methanol/water before UPLC-MS analysis. This protocol resulted in a median CV of feature intensities among experimental replicates of <20% for aqueous extracts and <30% for organic extracts. These data demonstrate the robustness of the proposed protocol for extracting metabolites from liver samples and make it well suited for untargeted liver profiling in studies exploring xenobiotic hepatotoxicity and clinical investigations of liver disease. The generic nature of this protocol facilitates its application to other tissues, for example, brain or lung, enhancing its utility in clinical and toxicological studies.

Metabolome-wide association study identifies multiple biomarkers that discriminate north and south Chinese populations at differing risks of cardiovascular disease: INTERMAP study.

Rates of heart disease and stroke vary markedly between north and south China. A (1)H NMR spectroscopy-based metabolome-wide association approach was used to identify urinary metabolites that discriminate between southern and northern Chinese population samples, to investigate population biomarkers that might relate to the difference in cardiovascular disease risk. NMR spectra were acquired from two 24-h urine specimens per person for 523 northern and 244 southern Chinese participants in the INTERMAP Study of macro/micronutrients and blood pressure. Discriminating metabolites were identified using orthogonal partial least squares discriminant analysis and assessed for statistical significance with conservative family wise error rate < 0.01 to minimize false positive findings. Urinary metabolites significantly (P < 1.2 × 10(-16) to 2.9 × 10(-69)) higher in northern than southern Chinese populations included dimethylglycine, alanine, lactate, branched-chain amino acids (isoleucine, leucine, valine), N-acetyls of glycoprotein fragments (including uromodulin), N-acetyl neuraminic acid, pentanoic/heptanoic acid, and methylguanidine; metabolites significantly (P < 1.1 × 10(-12) to 2 × 10(-127)) higher in the south were gut microbial cometabolites (hippurate, 4-cresyl sulfate, phenylacetylglutamine, 2-hydroxyisobutyrate), succinate, creatine, scyllo-inositol, prolinebetaine, and trans-aconitate. These findings indicate the importance of environmental influences (e.g., diet), endogenous metabolism, and mammalian-gut microbial cometabolism, which may help explain north-south China differences in cardiovascular disease risk.

The evolution of partial least squares models and related chemometric approaches in metabonomics and metabolic phenotyping

Metabonomics is a key element in systems biology, and with current analytical methods, generates vast amounts of quantitative or qualitative metabolic data. Understanding of the global function of the living organism can be achieved by integration of ‘omics’ approaches including metabonomics, genomics, transcriptomics and proteomics, increasing the complexity of the full data sets. Multivariate statistical approaches are well suited to extract the characterizing metabolic information associated with each level of dynamic process. In this review, we discuss techniques that have evolved from principal component analysis and partial least squares (PLS) methods with a focus on improved interpretation and modeling with respect to biomarker recovery and data visualization in the context of metabonomic applications. Visualization is of paramount importance to investigate complex metabolic signatures, the power and potential of which is illustrated with key papers. Recent improvements based on the removal of orthogonal variation are discussed in terms of interpretation enhancement, and are supported by relevant applications. Flexibility of PLS methods in general and of O‐PLS in particular allows implementation of derivative methods such as O2‐PLS, O‐PLS‐variance components, nonlinear methods, and batch modeling to improve analysis of complex data sets, which facilitates extraction of information related to subtle biological processes. These approaches can be used to address issues present in complex multi‐factorial data sets. Thus, we highlight the key advantages and limitations of the different latent variable applications for top‐down systems biology and assess the differences between the methods available. Copyright