Focus: a robust workflow for one-dimensional NMR spectral analysis

SPA-STOCSY: an automated tool for identifying annotated and non-annotated metabolites in high-throughput NMR spectra

MOTIVATION

Nuclear magnetic resonance spectroscopy (NMR) is widely used to analyze metabolites in biological samples, but the analysis requires specific expertise, it is time-consuming, and can be inaccurate. Here, we present a powerful automate tool, SPatial clustering Algorithm-Statistical TOtal Correlation SpectroscopY (SPA-STOCSY), which overcomes challenges faced when analyzing NMR data and identifies metabolites in a sample with high accuracy.

RESULTS

As a data-driven method, SPA-STOCSY estimates all parameters from the input dataset. It first investigates the covariance pattern among datapoints and then calculates the optimal threshold with which to cluster datapoints belonging to the same structural unit, i.e. the metabolite. Generated clusters are then automatically linked to a metabolite library to identify candidates. To assess SPA-STOCSY's efficiency and accuracy, we applied it to synthesized spectra and spectra acquired on Drosophila melanogaster tissue and human embryonic stem cells. In the synthesized spectra, SPA outperformed Statistical Recoupling of Variables (SRV), an existing method for clustering spectral peaks, by capturing a higher percentage of the signal regions and the close-to-zero noise regions. In the biological data, SPA-STOCSY performed comparably to the operator-based Chenomx analysis while avoiding operator bias, and it required <7 min of total computation time. Overall, SPA-STOCSY is a fast, accurate, and unbiased tool for untargeted analysis of metabolites in the NMR spectra. It may thus accelerate the use of NMR for scientific discoveries, medical diagnostics, and patient-specific decision making.

AVAILABILITY AND IMPLEMENTATION

The codes of SPA-STOCSY are available at https://github.com/LiuzLab/SPA-STOCSY.

SPA-STOCSY: An Automated Tool for Identification of Annotated and Non-Annotated Metabolites in High-Throughput NMR Spectra

Nuclear Magnetic Resonance (NMR) spectroscopy is widely used to analyze metabolites in biological samples, but the analysis can be cumbersome and inaccurate. Here, we present a powerful automated tool, SPA-STOCSY (Spatial Clustering Algorithm - Statistical Total Correlation Spectroscopy), which overcomes the challenges by identifying metabolites in each sample with high accuracy. As a data-driven method, SPA-STOCSY estimates all parameters from the input dataset, first investigating the covariance pattern and then calculating the optimal threshold with which to cluster data points belonging to the same structural unit, i.e. metabolite. The generated clusters are then automatically linked to a compound library to identify candidates. To assess SPA-STOCSY’s efficiency and accuracy, we applied it to synthesized and real NMR data obtained from Drosophila melanogaster brains and human embryonic stem cells. In the synthesized spectra, SPA outperforms Statistical Recoupling of Variables, an existing method for clustering spectral peaks, by capturing a higher percentage of the signal regions and the close-to-zero noise regions. In the real spectra, SPA-STOCSY performs comparably to operator-based Chenomx analysis but avoids operator bias and performs the analyses in less than seven minutes of total computation time. Overall, SPA-STOCSY is a fast, accurate, and unbiased tool for untargeted analysis of metabolites in the NMR spectra. As such, it might accelerate the utilization of NMR for scientific discoveries, medical diagnostics, and patient-specific decision making.

Statistical mediation of the relationships between chronological age and lipoproteins by nonessential amino acids in healthy men

Aging is a major risk factor for metabolic impairment that may lead to age-related diseases such as cardiovascular disease. Different mechanisms that may explain the interplay between aging and lipoproteins, and between aging and low-molecular-weight metabolites (LMWMs), in the metabolic dysregulation associated with age-related diseases have been described separately. Here, we statistically evaluated the possible mediation effects of LMWMs on the relationships between chronological age and lipoprotein concentrations in healthy men ranging from 19 to 75 years of age. Relative and absolute concentrations of LMWMs and lipoproteins, respectively, were assessed by nuclear magnetic resonance (NMR) spectroscopy. Multivariate linear regression and mediation analysis were conducted to explore the associations between age, lipoproteins and LMWMs. The statistical significance of the identified mediation effects was evaluated using the bootstrapping technique, and the identified mediation effects were validated on a publicly available dataset. Chronological age was statistically associated with five lipoprotein classes and subclasses. The mediation analysis showed that serine mediated 24.1% (95% CI: 22.9 – 24.7) of the effect of age on LDL-P, and glutamate mediated 17.9% (95% CI: 17.6 – 18.5) of the effect of age on large LDL-P. In the publicly available data, glutamate mediated the relationship between age and an NMR-derived surrogate of cholesterol. Our results suggest that the age-related increase in LDL particles may be mediated by a decrease in the nonessential amino acid glutamate. Future studies may contribute to a better understanding of the potential biological role of glutamate and LDL particles in aging mechanisms and age-related diseases.

Detection of highly overlapping peaks via adaptive apodization

Accurate peak detection is an essential component of many NMR tasks such as peak alignment, compound identification, and global spectral deconvolution. However, current peak detection approaches are generally limited by their ability to deal with spectral overlap, which has a deleterious effect on downstream data processing. In this work, we present the use of an adaptive apodization strategy that improves the detection of highly overlapping peaks. Sensitivity enhancement is used to identify general regions of interest and resolution enhancement is used to separate overlapping peaks, with parameters for both calculated directly from the data. Further limits on peak width help to reduce false positives. The method proposed in this work has been implemented in an open-source R package called rnmrfind that is available for download on GitHub (https://github.com/ssokolen/rnmrfind). A set of default parameters have been chosen to provide effective peak detection while keeping false positives to a minimum; however, application-specific tuning is possible through the modification of minimum peak width at half height (in Hz) and noise cutoff threshold (as a factor of estimated standard deviation). Comparison to existing packages rNMR and speaq on a series of ¹H NMR spectra demonstrates improved peak resolution with little to no apparent drawbacks.

The exposome paradigm to predict environmental health in terms of systemic homeostasis and resource balance based on NMR data science

The environment, from microbial ecosystems to recycled resources, fluctuates dynamically due to many physical, chemical and biological factors, the profile of which reflects changes in overall state, such as environmental illness caused by a collapse of homeostasis. To evaluate and predict environmental health in terms of systemic homeostasis and resource balance, a comprehensive understanding of these factors requires an approach based on the "exposome paradigm", namely the totality of exposure to all substances. Furthermore, in considering sustainable development to meet global population growth, it is important to gain an understanding of both the circulation of biological resources and waste recycling in human society. From this perspective, natural environment, agriculture, aquaculture, wastewater treatment in industry, biomass degradation and biodegradable materials design are at the forefront of current research. In this respect, nuclear magnetic resonance (NMR) offers tremendous advantages in the analysis of samples of molecular complexity, such as crude bio-extracts, intact cells and tissues, fibres, foods, feeds, fertilizers and environmental samples. Here we outline examples to promote an understanding of recent applications of solution-state, solid-state, time-domain NMR and magnetic resonance imaging (MRI) to the complex evaluation of organisms, materials and the environment. We also describe useful databases and informatics tools, as well as machine learning techniques for NMR analysis, demonstrating that NMR data science can be used to evaluate the exposome in both the natural environment and human society towards a sustainable future.

Analysis of NMR Metabolomics Data

In this chapter, we summarize data preprocessing and data analysis strategies used for analysis of NMR data for metabolomics studies. Metabolomics consists of the analysis of the low molecular weight compounds in cells, tissues, or biological fluids, and has been used to reveal biomarkers for early disease detection and diagnosis, to monitor interventions, and to provide information on pathway perturbations to inform mechanisms and identifying targets. Metabolic profiling (also termed metabotyping) involves the analysis of hundreds to thousands of molecules using mainly state-of-the-art mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy technologies. While NMR is less sensitive than mass spectrometry, NMR does provide a wealth of complex and information rich metabolite data. NMR data together with the use of conventional statistics, modeling methods, and bioinformatics tools reveals biomarker and mechanistic information. A typical NMR spectrum, with up to 64k data points, of a complex biological fluid or an extract of cells and tissues consists of thousands of sharp signals that are mainly derived from small molecules. In addition, a number of advanced NMR spectroscopic methods are available for extracting information on high molecular weight compounds such as lipids or lipoproteins. There are numerous data preprocessing, data reduction, and analysis methods developed and evolving in the field of NMR metabolomics. Our goal is to provide an extensive summary of NMR data preprocessing and analysis strategies by providing examples and open source and commercially available analysis software and bioinformatics tools.

Automated Analysis of Large-Scale NMR Data Generates Metabolomic Signatures and Links Them to Candidate Metabolites

Identification of metabolites in large-scale ¹H NMR data from human biofluids remains challenging due to the complexity of the spectra and their sensitivity to pH and ionic concentrations. In this work, we tested the capacity of three analysis tools to extract metabolite signatures from 968 NMR profiles of human urine samples. Specifically, we studied sets of covarying features derived from principal component analysis (PCA), the iterative signature algorithm (ISA), and averaged correlation profiles (ACP), a new method we devised inspired by the STOCSY approach. We used our previously developed metabomatching method to match the sets generated by these algorithms to NMR spectra of individual metabolites available in public databases. On the basis of the number and quality of the matches, we concluded that ISA and ACP can robustly identify ten and nine metabolites, respectively, half of which were shared, while PCA did not produce any signatures with robust matches.

Similarity and differential NMR spectroscopy in metabolomics: application to the analysis of vegetable oils with 1H and 13C NMR

INTRODUCTION

In NMR based metabolomics there is a need for tools to easily compare spectra and to extract the maximum of information from the data.

OBJECTIVES

The calculation of similarity and performing differential NMR spectroscopy provides important additional information for classification and validation in metabolomics experiments.

METHODS

From 13 different vegetable oils samples were analysed by ¹H and ¹³C NMR. The similarity between spectra was calculated and differential NMR spectroscopy was used to discover marker compounds.

RESULTS

The similarity between the individual spectra was calculated for the spectra of all samples. The similarity was used to verify and improve the alignment. For vegetable oils which showed a high similarity, e.g. chia seed oil and linseed oil, differential NMR spectroscopy was used to discover marker compounds.

CONCLUSIONS

The calculation of similarity is an important tool to reveal variability between samples and spectra and can be used to verify data sets and improve alignment or binning procedures. With differential spectroscopy marker compounds are easily discovered. The methods can be seen as an important addition to the routine procedures of metabolomics experiments.

Metabomatching: Using genetic association to identify metabolites in proton NMR spectroscopy

A metabolome-wide genome-wide association study (mGWAS) aims to discover the effects of genetic variants on metabolome phenotypes. Most mGWASes use as phenotypes concentrations of limited sets of metabolites that can be identified and quantified from spectral information. In contrast, in an untargeted mGWAS both identification and quantification are forgone and, instead, all measured metabolome features are tested for association with genetic variants. While the untargeted approach does not discard data that may have eluded identification, the interpretation of associated features remains a challenge. To address this issue, we developed metabomatching to identify the metabolites underlying significant associations observed in untargeted mGWASes on proton NMR metabolome data. Metabomatching capitalizes on genetic spiking, the concept that because metabolome features associated with a genetic variant tend to correspond to the peaks of the NMR spectrum of the underlying metabolite, genetic association can allow for identification. Applied to the untargeted mGWASes in the SHIP and CoLaus cohorts and using 180 reference NMR spectra of the urine metabolome database, metabomatching successfully identified the underlying metabolite in 14 of 19, and 8 of 9 associations, respectively. The accuracy and efficiency of our method make it a strong contender for facilitating or complementing metabolomics analyses in large cohorts, where the availability of genetic, or other data, enables our approach, but targeted quantification is limited.

New Strategies and Challenges in Lung Proteomics and Metabolomics. An Official American Thoracic Society Workshop Report

This document presents the proceedings from the workshop entitled, "New Strategies and Challenges in Lung Proteomics and Metabolomics" held February 4th-5th, 2016, in Denver, Colorado. It was sponsored by the National Heart Lung Blood Institute, the American Thoracic Society, the Colorado Biological Mass Spectrometry Society, and National Jewish Health. The goal of this workshop was to convene, for the first time, relevant experts in lung proteomics and metabolomics to discuss and overcome specific challenges in these fields that are unique to the lung. The main objectives of this workshop were to identify, review, and/or understand: (1) emerging technologies in metabolomics and proteomics as applied to the study of the lung; (2) the unique composition and challenges of lung-specific biological specimens for metabolomic and proteomic analysis; (3) the diverse informatics approaches and databases unique to metabolomics and proteomics, with special emphasis on the lung; (4) integrative platforms across genetic and genomic databases that can be applied to lung-related metabolomic and proteomic studies; and (5) the clinical applications of proteomics and metabolomics. The major findings and conclusions of this workshop are summarized at the end of the report, and outline the progress and challenges that face these rapidly advancing fields.

NMR window of molecular complexity showing homeostasis in superorganisms

NMR offers tremendous advantages in the analyses of molecular complexity. The “big-data” are produced during the acquisition of fingerprints that must be stored and shared for posterior analysis and verifications.

Navigating freely-available software tools for metabolomics analysis

INTRODUCTION

The field of metabolomics has expanded greatly over the past two decades, both as an experimental science with applications in many areas, as well as in regards to data standards and bioinformatics software tools. The diversity of experimental designs and instrumental technologies used for metabolomics has led to the need for distinct data analysis methods and the development of many software tools.

OBJECTIVES

To compile a comprehensive list of the most widely used freely available software and tools that are used primarily in metabolomics.

METHODS

The most widely used tools were selected for inclusion in the review by either ≥ 50 citations on Web of Science (as of 08/09/16) or the use of the tool being reported in the recent Metabolomics Society survey. Tools were then categorised by the type of instrumental data (i.e. LC-MS, GC-MS or NMR) and the functionality (i.e. pre- and post-processing, statistical analysis, workflow and other functions) they are designed for.

RESULTS

A comprehensive list of the most used tools was compiled. Each tool is discussed within the context of its application domain and in relation to comparable tools of the same domain. An extended list including additional tools is available at https://github.com/RASpicer/MetabolomicsTools which is classified and searchable via a simple controlled vocabulary.

CONCLUSION

This review presents the most widely used tools for metabolomics analysis, categorised based on their main functionality. As future work, we suggest a direct comparison of tools' abilities to perform specific data analysis tasks e.g. peak picking.

Urine metabolome profiling of immune-mediated inflammatory diseases

BACKGROUND

Immune-mediated inflammatory diseases (IMIDs) are a group of complex and prevalent diseases where disease diagnostic and activity monitoring is highly challenging. The determination of the metabolite profiles of biological samples is becoming a powerful approach to identify new biomarkers of clinical utility. In order to identify new metabolite biomarkers of diagnosis and disease activity, we have performed the first large-scale profiling of the urine metabolome of the six most prevalent IMIDs: rheumatoid arthritis, psoriatic arthritis, psoriasis, systemic lupus erythematosus, Crohn's disease, and ulcerative colitis.

METHODS

Using nuclear magnetic resonance, we analyzed the urine metabolome in a discovery cohort of 1210 patients and 100 controls. Within each IMID, two patient subgroups were recruited representing extreme disease activity (very high vs. very low). Metabolite association analysis with disease diagnosis and disease activity was performed using multivariate linear regression in order to control for the effects of clinical, epidemiological, or technical variability. After multiple test correction, the most significant metabolite biomarkers were validated in an independent cohort of 1200 patients and 200 controls.

RESULTS

In the discovery cohort, we identified 28 significant associations between urine metabolite levels and disease diagnosis and three significant metabolite associations with disease activity (P FDR < 0.05). Using the validation cohort, we validated 26 of the diagnostic associations and all three metabolite associations with disease activity (P FDR < 0.05). Combining all diagnostic biomarkers using multivariate classifiers we obtained a good disease prediction accuracy in all IMIDs and particularly high in inflammatory bowel diseases. Several of the associated metabolites were found to be commonly altered in multiple IMIDs, some of which can be considered as hub biomarkers. The analysis of the metabolic reactions connecting the IMID-associated metabolites showed an over-representation of citric acid cycle, phenylalanine, and glycine-serine metabolism pathways.

CONCLUSIONS

This study shows that urine is a source of biomarkers of clinical utility in IMIDs. We have found that IMIDs show similar metabolic changes, particularly between clinically similar diseases and we have found, for the first time, the presence of hub metabolites. These findings represent an important step in the development of more efficient and less invasive diagnostic and disease monitoring methods in IMIDs.

Genome-Wide Association Study with Targeted and Non-targeted NMR Metabolomics Identifies 15 Novel Loci of Urinary Human Metabolic Individuality

Genome-wide association studies with metabolic traits (mGWAS) uncovered many genetic variants that influence human metabolism. These genetically influenced metabotypes (GIMs) contribute to our metabolic individuality, our capacity to respond to environmental challenges, and our susceptibility to specific diseases. While metabolic homeostasis in blood is a well investigated topic in large mGWAS with over 150 known loci, metabolic detoxification through urinary excretion has only been addressed by few small mGWAS with only 11 associated loci so far. Here we report the largest mGWAS to date, combining targeted and non-targeted ¹H NMR analysis of urine samples from 3,861 participants of the SHIP-0 cohort and 1,691 subjects of the KORA F4 cohort. We identified and replicated 22 loci with significant associations with urinary traits, 15 of which are new (HIBCH, CPS1, AGXT, XYLB, TKT, ETNPPL, SLC6A19, DMGDH, SLC36A2, GLDC, SLC6A13, ACSM3, SLC5A11, PNMT, SLC13A3). Two-thirds of the urinary loci also have a metabolite association in blood. For all but one of the 6 loci where significant associations target the same metabolite in blood and urine, the genetic effects have the same direction in both fluids. In contrast, for the SLC5A11 locus, we found increased levels of myo-inositol in urine whereas mGWAS in blood reported decreased levels for the same genetic variant. This might indicate less effective re-absorption of myo-inositol in the kidneys of carriers. In summary, our study more than doubles the number of known loci that influence urinary phenotypes. It thus allows novel insights into the relationship between blood homeostasis and its regulation through excretion. The newly discovered loci also include variants previously linked to chronic kidney disease (CPS1, SLC6A13), pulmonary hypertension (CPS1), and ischemic stroke (XYLB). By establishing connections from gene to disease via metabolic traits our results provide novel hypotheses about molecular mechanisms involved in the etiology of diseases.

Human metabolism is influenced by genetic and environmental factors defining a person’s metabolic individuality. This individuality is linked to personal differences in the ability to react on metabolic challenges and in the susceptibility to specific diseases. By investigating how common variants in genetic regions (loci) affect individual blood metabolite levels, the substantial contribution of genetic inheritance to metabolic individuality has been demonstrated previously. Meanwhile, more than 150 loci influencing metabolic homeostasis in blood are known. Here we shift the focus to genetic variants that modulate urinary metabolite excretion, for which only 11 loci were reported so far. In the largest genetic study on urinary metabolites to date, we identified 15 additional loci. Most of the 26 loci also affect blood metabolite levels. This shows that the metabolic individuality seen in blood is also reflected in urine, which is expected when urine is regarded as “diluted blood”. Nonetheless, we also found loci that appear to primarily influence metabolite excretion. For instance, we identified genetic variants near a gene of a transporter that change the capability for renal re-absorption of the transporter’s substrate. Thus, our findings could help to elucidate molecular mechanisms influencing kidney function and the body’s detoxification capabilities.

mQTL.NMR: an integrated suite for genetic mapping of quantitative variations of (1)H NMR-based metabolic profiles

High-throughput (1)H nuclear magnetic resonance (NMR) is an increasingly popular robust approach for qualitative and quantitative metabolic profiling, which can be used in conjunction with genomic techniques to discover novel genetic associations through metabotype quantitative trait locus (mQTL) mapping. There is therefore a crucial necessity to develop specialized tools for an accurate detection and unbiased interpretability of the genetically determined metabolic signals. Here we introduce and implement a combined chemoinformatic approach for objective and systematic analysis of untargeted (1)H NMR-based metabolic profiles in quantitative genetic contexts. The R/Bioconductor mQTL.NMR package was designed to (i) perform a series of preprocessing steps restoring spectral dependency in collinear NMR data sets to reduce the multiple testing burden, (ii) carry out robust and accurate mQTL mapping in human cohorts as well as in rodent models, (iii) statistically enhance structural assignment of genetically determined metabolites, and (iv) illustrate results with a series of visualization tools. Built-in flexibility and implementation in the powerful R/Bioconductor framework allow key preprocessing steps such as peak alignment, normalization, or dimensionality reduction to be tailored to specific problems. The mQTL.NMR package is freely available with its source code through the Comprehensive R/Bioconductor repository and its own website ( http://www.ican-institute.org/tools/ ). It represents a significant advance to facilitate untargeted metabolomic data processing and quantitative analysis and their genetic mapping.

Analytical methods in untargeted metabolomics: state of the art in 2015

Metabolomics comprises the methods and techniques that are used to measure the small molecule composition of biofluids and tissues, and is actually one of the most rapidly evolving research fields. The determination of the metabolomic profile - the metabolome - has multiple applications in many biological sciences, including the developing of new diagnostic tools in medicine. Recent technological advances in nuclear magnetic resonance and mass spectrometry are significantly improving our capacity to obtain more data from each biological sample. Consequently, there is a need for fast and accurate statistical and bioinformatic tools that can deal with the complexity and volume of the data generated in metabolomic studies. In this review, we provide an update of the most commonly used analytical methods in metabolomics, starting from raw data processing and ending with pathway analysis and biomarker identification. Finally, the integration of metabolomic profiles with molecular data from other high-throughput biotechnologies is also reviewed.

MVAPACK: a complete data handling package for NMR metabolomics

Data handling in the field of NMR metabolomics has historically been reliant on either in-house mathematical routines or long chains of expensive commercial software. Thus, while the relatively simple biochemical protocols of metabolomics maintain a low barrier to entry, new practitioners of metabolomics experiments are forced to either purchase expensive software packages or craft their own data handling solutions from scratch. This inevitably complicates the standardization and communication of data handling protocols in the field. We report a newly developed open-source platform for complete NMR metabolomics data handling, MVAPACK, and describe its application on an example metabolic fingerprinting data set.

Quantitative Quantum Mechanical Spectral Analysis (qQMSA) of (1)H NMR spectra of complex mixtures and biofluids

The quantitative interpretation of (1)H NMR spectra of mixtures like the biofluids is a demanding task due to spectral complexity and overlap. Complications may arise also from water suppression, T2-editing, protein interactions, relaxation differences of the species, experimental artifacts and, furthermore, the spectra may contain unknown components and macromolecular background which cannot be easily separated from baseline. In this work, tools and strategies for quantitative Quantum Mechanical Spectral Analysis (qQMSA) of (1)H NMR spectra from complex mixtures were developed and systematically assessed. In the present approach, the signals of well-defined, stoichiometric components are described by a QM model, while the background is described by a multiterm baseline function and the unknown signals using optimizable and adjustable lines, regular multiplets or any spectral structures which can be composed from spectral lines. Any prior knowledge available from the spectrum can also be added to the model. Fitting strategies for weak and strongly overlapping spectral systems were developed and assessed using two basic model systems, the metabolite mixtures without and with macromolecular (serum) background. The analyses show that if the spectra are measured in high-throughput manner, the consistent absolute quantification demands some calibration to compensate the different response factors of the protons and compounds. On the other hand, the results show that also the T2-edited spectra can be measured so that they obey well the QM rules. In general, qQMSA exploits and interprets the spectral information in maximal way taking full advantage from the QM properties of the spectra and, at the same time, offers chemical confidence which means that individual components can be identified with high confidence on the basis of their accurate spectral parameters.