speaq 2.0: A complete workflow for high-throughput 1D NMR spectra processing and quantification

Towards objective measurements of habitual dietary intake patterns: comparing NMR metabolomics and food frequency questionnaire data in a population-based cohort

BACKGROUND

Low-quality, non-diverse diet is a main risk factor for premature death. Accurate measurement of habitual diet is challenging and there is a need for validated objective methods. Blood metabolite patterns reflect direct or enzymatically diet-induced metabolites. Here, we aimed to evaluate associations between blood metabolite patterns and a priori and data-driven food intake patterns.

METHODS

1, 895 participants in the Northern Sweden Health and Disease Study, a population-based prospective cohort study, were included. Fasting plasma samples were analyzed with 1H Nuclear Magnetic Resonance. Food intake data from a 64-item validated food frequency questionnaire were summarized into a priori Healthy Diet Score (HDS), relative Mediterranean Diet Score (rMDS) and a set of plant-based diet indices (PDI) as well as data driven clusters from latent class analyses (LCA). Orthogonal projections to latent structures (OPLS) were used to explore clustering patterns of metabolites and their relation to reported dietary intake patterns.

RESULTS

Age, sex, body mass index, education and year of study participation had significant influence on OPLS metabolite models. OPLS models for healthful PDI and LCA-clusters were not significant, whereas for HDS, rMDS, PDI and unhealthful PDI significant models were obtained (CV-ANOVA p < 0.001). Still, model statistics were weak and the ability of the models to correctly classify participants into highest and lowest quartiles of rMDS, PDI and unhealthful PDI was poor (50%/78%, 42%/75% and 59%/70%, respectively).

CONCLUSION

Associations between blood metabolite patterns and a priori as well as data-driven food intake patterns were poor. NMR metabolomics may not be sufficiently sensitive to small metabolites that distinguish between complex dietary intake patterns, like lipids.

Multifaceted Analysis of Cerebrospinal Fluid and Serum from Progressive Multiple Sclerosis Patients: Potential Role of Vitamin C and Metal Ion Imbalance in the Divergence of Primary Progressive Multiple Sclerosis and Secondary Progressive Multiple Sclerosis

The progressive forms of multiple sclerosis (MS) primary progressive MS (PPMS) and secondary progressive MS (SPMS) are clinically distinguished by the rate at which symptoms worsen. Little is however known about the pathological mechanisms underlying the differential rate of accumulation of pathological changes. In this study, ¹H NMR spectroscopy was used to measure low-molecular-weight metabolites in paired cerebrospinal fluid (CSF) and serum of PPMS, SPMS, and control patients, as well as to determine lipoproteins and glycoproteins in serum samples. Additionally, neurodegenerative and inflammatory markers, neurofilament light (NFL) and chitinase-3-like protein 1 (CHI3L1), and the concentration of seven metal elements, Mg, Mn, Cu, Fe, Pb, Zn, and Ca, were also determined in both CSF and serum. The results indicate that the pathological changes associated with progressive MS are mainly localized in the central nervous system (CNS). More so, PPMS and SPMS patients with comparable disability status are pathologically similar in relation to neurodegeneration, neuroinflammation, and some metabolites that distinguish them from controls. However, the rapid progression of PPMS from the onset may be driven by a combination of neurotoxicity induced by heavy metals coupled with diminished CNS antioxidative capacity associated with differential intrathecal ascorbate retention and imbalance of Mg and Cu.

GIPMA: Global Intensity-Guided Peak Matching and Alignment for 2D 1H-13C HSQC-Based Metabolomics

Two-dimensional (2D) ¹H-¹³C heteronuclear single quantum coherence (HSQC) has been increasingly applied to metabolomics studies because it can greatly improve the resolving capability compared with one-dimensional (1D) ¹H NMR. However, preprocessing methods such as peak matching and alignment tools for 2D NMR-based metabolomics have lagged behind similar methods for 1D ¹H NMR-based metabolomics. Correct matching and alignment of 2D NMR spectral features across multiple samples are particularly important for subsequent multivariate data analysis. Considering different intensity dynamic ranges of a variety of metabolites and the chemical shift variation across the spectra of multiple samples, here, we developed an efficient peak matching and alignment algorithm for 2D ¹H-¹³C HSQC-based metabolomics, called global intensity-guided peak matching and alignment (GIPMA). In GIPMA, peaks identified in all spectra are pooled together and sorted by intensity. Chemical shift of a stronger peak is regarded to be more accurate and reliable than that of a weaker peak. The strongest undesignated peak is chosen as the reference of a new cluster if it is not located within the chemical shift tolerance of any existing peak cluster (PC), or otherwise it is matched to an existing PC and the aligned chemical shift of the PC is updated as the intensity-weighted average of the chemical shifts of all peaks in the cluster. Setting an optimum chemical shift tolerance (Δδ_o) is critical for the peak matching and alignment across multiple samples. GIPMA dynamically searches for and intelligently selects the Δδ_o for peak matching to maximize the number of valid peak clusters (vPC), that is, spectral features, among multiple samples. By GIPMA, fully automatic peakwise matching and alignment do not require any spectrum as initial reference, while the chemical shift of each PC is updated as the intensity-weighted average of the chemical shifts of all peaks in the same PC, which is warranted to be statistically more accurate. Accurate chemical shifts for each representative spectral feature will facilitate subsequent peak assignment and are essential for correct metabolite identification and result interpretation. The proposed method was demonstrated successfully on the spectra of six model mixtures consisting of seven typical metabolites, yielding correct matching of all known spectral features. The performance of GIPMA was also demonstrated on 2D ¹H-¹³C HSQC spectra of 87 real extracts of 29 samples of five Dendrobium species. Hierarchical cluster analysis (HCA) and principal component analysis (PCA) of the 87 matched and aligned spectra by GIPMA generates correct classification of the 29 samples into five groups. In summary, the proposed algorithm of GIPMA provided a practical peak matching and alignment method to facilitate 2D NMR-based metabolomics studies.

Problems, principles and progress in computational annotation of NMR metabolomics data

BACKGROUND

Compound identification remains a critical bottleneck in the process of exploiting Nuclear Magnetic Resonance (NMR) metabolomics data, especially for ¹H 1-dimensional (¹H 1D) data. As databases of reference compound spectra have grown, workflows have evolved to rely heavily on their search functions to facilitate this process by generating lists of potential metabolites found in complex mixture data, facilitating annotation and identification. However, approaches for validating and communicating annotations are most often guided by expert knowledge, and therefore are highly variable despite repeated efforts to align practices and define community standards.

AIM OF REVIEW

This review is aimed at broadening the application of automated annotation tools by discussing the key ideas of spectral matching and beginning to describe a set of terms to classify this information, thus advancing standards for communicating annotation confidence. Additionally, we hope that this review will facilitate the growing collaboration between chemical data scientists, software developers and the NMR metabolomics community aiding development of long-term software solutions.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We begin with a brief discussion of the typical untargeted NMR identification workflow. We differentiate between annotation (hypothesis generation, filtering), and identification (hypothesis testing, verification), and note the utility of different NMR data features for annotation. We then touch on three parts of annotation: (1) generation of queries, (2) matching queries to reference data, and (3) scoring and confidence estimation of potential matches for verification. In doing so, we highlight existing approaches to automated and semi-automated annotation from the perspective of the structural information they utilize, as well as how this information can be represented computationally.

Dietary restriction and medical therapy drives PPARα-regulated improvements in early diabetic kidney disease in male rats

The attenuation of diabetic kidney disease (DKD) by metabolic surgery is enhanced by pharmacotherapy promoting renal fatty acid oxidation (FAO). Using the Zucker Diabetic Fatty and Zucker Diabetic Sprague Dawley rat models of DKD, we conducted studies to determine if these effects could be replicated with a non-invasive bariatric mimetic intervention. Metabolic control and renal injury were compared in rats undergoing a dietary restriction plus medical therapy protocol (DMT; fenofibrate, liraglutide, metformin, ramipril, and rosuvastatin) and ad libitum-fed controls. The global renal cortical transcriptome and urinary 1H-NMR metabolomic profiles were also compared. Kidney cell type-specific and medication-specific transcriptomic responses were explored through in silico deconvolution. Transcriptomic and metabolomic correlates of improvements in kidney structure were defined using a molecular morphometric approach. The DMT protocol led to ∼20% weight loss, normalized metabolic parameters and was associated with reductions in indices of glomerular and proximal tubular injury. The transcriptomic response to DMT was dominated by changes in fenofibrate- and peroxisome proliferator-activated receptor-α (PPARα)-governed peroxisomal and mitochondrial FAO transcripts localizing to the proximal tubule. DMT induced urinary excretion of PPARα-regulated metabolites involved in nicotinamide metabolism and reversed DKD-associated changes in the urinary excretion of tricarboxylic acid (TCA) cycle intermediates. FAO transcripts and urinary nicotinamide and TCA cycle metabolites were moderately to strongly correlated with improvements in glomerular and proximal tubular injury. Weight loss plus pharmacological PPARα agonism is a promising means of attenuating DKD.

Preanalytical Pitfalls in Untargeted Plasma Nuclear Magnetic Resonance Metabolomics of Endocrine Hypertension

Despite considerable morbidity and mortality, numerous cases of endocrine hypertension (EHT) forms, including primary aldosteronism (PA), pheochromocytoma and functional paraganglioma (PPGL), and Cushing's syndrome (CS), remain undetected. We aimed to establish signatures for the different forms of EHT, investigate potentially confounding effects and establish unbiased disease biomarkers. Plasma samples were obtained from 13 biobanks across seven countries and analyzed using untargeted NMR metabolomics. We compared unstratified samples of 106 PHT patients to 231 EHT patients, including 104 PA, 94 PPGL and 33 CS patients. Spectra were subjected to a multivariate statistical comparison of PHT to EHT forms and the associated signatures were obtained. Three approaches were applied to investigate and correct confounding effects. Though we found signatures that could separate PHT from EHT forms, there were also key similarities with the signatures of sample center of origin and sample age. The study design restricted the applicability of the corrections employed. With the samples that were available, no biomarkers for PHT vs. EHT could be identified. The complexity of the confounding effects, evidenced by their robustness to correction approaches, highlighted the need for a consensus on how to deal with variabilities probably attributed to preanalytical factors in retrospective, multicenter metabolomics studies.

The Long-Term Culture of Human Fibroblasts Reveals a Spectroscopic Signature of Senescence

Aging is a complex process which leads to progressive loss of fitness/capability/ability, increasing susceptibility to disease and, ultimately, death. Regardless of the organism, there are some features common to aging, namely, the loss of proteostasis and cell senescence. Mammalian cell lines have been used as models to study the aging process, in particular, cell senescence. Thus, the aim of this study was to characterize the senescence-associated metabolic profile of a long-term culture of human fibroblasts using Fourier Transform Infrared and Nuclear Magnetic Resonance spectroscopy. We sub-cultivated fibroblasts from a newborn donor from passage 4 to passage 17 and the results showed deep changes in the spectroscopic profile of cells over time. Late passage cells were characterized by a decrease in the length of fatty acid chains, triglycerides and cholesterol and an increase in lipid unsaturation. We also found an increase in the content of intermolecular β-sheets, possibly indicating an increase in protein aggregation levels in cells of later passages. Metabolic profiling by NMR showed increased levels of extracellular lactate, phosphocholine and glycine in cells at later passages. This study suggests that spectroscopy approaches can be successfully used to study changes concomitant with cell senescence and validate the use of human fibroblasts as a model to monitor the aging process.

Urinary Metabolomic Changes Accompanying Albuminuria Remission following Gastric Bypass Surgery for Type 2 Diabetic Kidney Disease

In the Microvascular Outcomes after Metabolic Surgery randomised clinical trial (MOMS RCT, NCT01821508), combined metabolic surgery (gastric bypass) plus medical therapy (CSM) was superior to medical therapy alone (MTA) as a means of achieving albuminuria remission at 2-year follow-up in patients with obesity and early diabetic kidney disease (DKD). In the present study, we assessed the urinary ¹H-NMR metabolome in a subgroup of patients from both arms of the MOMS RCT at baseline and 6-month follow-up. Whilst CSM and MTA both reduced the urinary excretion of sugars, CSM generated a distinctive urinary metabolomic profile characterised by increases in host–microbial co-metabolites (N-phenylacetylglycine, trimethylamine N-oxide, and 4-aminobutyrate (GABA)) and amino acids (arginine and glutamine). Furthermore, reductions in aromatic amino acids (phenylalanine and tyrosine), as well as branched-chain amino acids (BCAAs) and related catabolites (valine, leucine, 3-hydroxyisobutyrate, 3-hydroxyisovalerate, and 3-methyl-2-oxovalerate), were observed following CSM but not MTA. Improvements in BMI did not correlate with improvements in metabolic and renal indices following CSM. Conversely, urinary metabolites changed by CSM at 6 months were moderately to strongly correlated with improvements in blood pressure, glycaemia, triglycerides, and albuminuria up to 24 months following treatment initiation, highlighting the potential involvement of these shifts in the urinary metabolomic profile in the metabolic and renoprotective effects of CSM.

Medications Activating Tubular Fatty Acid Oxidation Enhance the Protective Effects of Roux-en-Y Gastric Bypass Surgery in a Rat Model of Early Diabetic Kidney Disease

Background

Roux-en-Y gastric bypass surgery (RYGB) improves biochemical and histological parameters of diabetic kidney disease (DKD). Targeted adjunct medical therapy may enhance renoprotection following RYGB.

Methods

The effects of RYGB and RYGB plus fenofibrate, metformin, ramipril, and rosuvastatin (RYGB-FMRR) on metabolic control and histological and ultrastructural indices of glomerular and proximal tubular injury were compared in the Zucker Diabetic Sprague Dawley (ZDSD) rat model of DKD. Renal cortical transcriptomic (RNA-sequencing) and urinary metabolomic (1H-NMR spectroscopy) responses were profiled and integrated. Transcripts were assigned to kidney cell types through in silico deconvolution in kidney single-nucleus RNA-sequencing and microdissected tubular epithelial cell proteomics datasets. Medication-specific transcriptomic responses following RYGB-FMRR were explored using a network pharmacology approach. Omic correlates of improvements in structural and ultrastructural indices of renal injury were defined using a molecular morphometric approach.

Results

RYGB-FMRR was superior to RYGB alone with respect to metabolic control, albuminuria, and histological and ultrastructural indices of glomerular injury. RYGB-FMRR reversed DKD-associated changes in mitochondrial morphology in the proximal tubule to a greater extent than RYGB. Attenuation of transcriptomic pathway level activation of pro-fibrotic responses was greater after RYGB-FMRR than RYGB. Fenofibrate was found to be the principal medication effector of gene expression changes following RYGB-FMRR, which led to the transcriptional induction of PPARα-regulated genes that are predominantly expressed in the proximal tubule and which regulate peroxisomal and mitochondrial fatty acid oxidation (FAO). After omics integration, expression of these FAO transcripts positively correlated with urinary levels of PPARα-regulated nicotinamide metabolites and negatively correlated with urinary tricarboxylic acid (TCA) cycle intermediates. Changes in FAO transcripts and nicotinamide and TCA cycle metabolites following RYGB-FMRR correlated strongly with improvements in glomerular and proximal tubular injury.

Conclusions

Integrative multi-omic analyses point to PPARα-stimulated FAO in the proximal tubule as a dominant effector of treatment response to combined surgical and medical therapy in experimental DKD. Synergism between RYGB and pharmacological stimulation of FAO represents a promising combinatorial approach to the treatment of DKD in the setting of obesity.

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.

Nuclear Magnetic Resonance Spectroscopy in Clinical Metabolomics and Personalized Medicine: Current Challenges and Perspectives

Personalized medicine is probably the most promising area being developed in modern medicine. This approach attempts to optimize the therapies and the patient care based on the individual patient characteristics. Its success highly depends on the way the characterization of the disease and its evolution, the patient’s classification, its follow-up and the treatment could be optimized. Thus, personalized medicine must combine innovative tools to measure, integrate and model data. Towards this goal, clinical metabolomics appears as ideally suited to obtain relevant information. Indeed, the metabolomics signature brings crucial insight to stratify patients according to their responses to a pathology and/or a treatment, to provide prognostic and diagnostic biomarkers, and to improve therapeutic outcomes. However, the translation of metabolomics from laboratory studies to clinical practice remains a subsequent challenge. Nuclear magnetic resonance spectroscopy (NMR) and mass spectrometry (MS) are the two key platforms for the measurement of the metabolome. NMR has several advantages and features that are essential in clinical metabolomics. Indeed, NMR spectroscopy is inherently very robust, reproducible, unbiased, quantitative, informative at the structural molecular level, requires little sample preparation and reduced data processing. NMR is also well adapted to the measurement of large cohorts, to multi-sites and to longitudinal studies. This review focus on the potential of NMR in the context of clinical metabolomics and personalized medicine. Starting with the current status of NMR-based metabolomics at the clinical level and highlighting its strengths, weaknesses and challenges, this article also explores how, far from the initial “opposition” or “competition”, NMR and MS have been integrated and have demonstrated a great complementarity, in terms of sample classification and biomarker identification. Finally, a perspective discussion provides insight into the current methodological developments that could significantly raise NMR as a more resolutive, sensitive and accessible tool for clinical applications and point-of-care diagnosis. Thanks to these advances, NMR has a strong potential to join the other analytical tools currently used in clinical settings.

Glycemic Control and Metabolic Adaptation in Response to High-Fat versus High-Carbohydrate Diets-Data from a Randomized Cross-Over Study in Healthy Subjects

Granular study of metabolic responses to alterations in the ratio of dietary macro-nutrients can enhance our understanding of how dietary modifications influence patients with impaired glycemic control. In order to study the effect of diets enriched in fat or carbohydrates, fifteen healthy, normal-weight volunteers received, in a cross-over design, and in a randomized unblinded order, two weeks of an iso-caloric high-fat diet (HFD: 60E% from fat) and a high-carbohydrate diet (HCD: 60E% from carbohydrates). A mixed meal test (MMT) was performed at the end of each dietary period to examine glucose clearance kinetics and insulin and incretin hormone levels, as well as plasma metabolomic profiles. The MMT induced almost identical glycemia and insulinemia following the HFD or HCD. GLP-1 levels were higher after the HFD vs. HCD, whereas GIP did not differ. The HFD, compared to the HCD, increased the levels of several metabolomic markers of risk for the development of insulin resistance, e.g., branched-chain amino acid (valine and leucine), creatine and α-hydroxybutyric acid levels. In normal-weight, healthy volunteers, two weeks of the HFD vs. HCD showed similar profiles of meal-induced glycemia and insulinemia. Despite this, the HFD showed a metabolomic pattern implying a risk for a metabolic shift towards impaired insulin sensitivity in the long run.

NMR Metabolite Profiles of the Bivalve Mollusc Mytilus galloprovincialis Before and After Immune Stimulation With Vibrio splendidus

The hemolymph metabolome of Mytilus galloprovincialis injected with live Vibrio splendidus bacteria was analyzed by ¹H-NMR spectrometry. Changes in spectral hemolymph profiles were already detected after mussel acclimation (3 days at 18 or 25 °C). A significant decrease of succinic acid was accompanied by an increase of most free amino acids, mytilitol, and, to a smaller degree, osmolytes. These metabolic changes are consistent with effective osmoregulation, and the restart of aerobic respiration after the functional anaerobiosis occurred during transport. The injection of Vibrio splendidus in mussels acclimated at 18°C caused a significant decrease of several amino acids, sugars, and unassigned chemical species, more pronounced at 24 than at 12 h postinjection. Correlation heatmaps indicated dynamic metabolic adjustments and the relevance of protein turnover in maintaining the homeostasis during the response to stressful stimuli. This study confirms NMR-based metabolomics as a feasible analytical approach complementary to other omics techniques in the investigation of the functional mussel responses to environmental challenges.

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.

Statistical analysis in metabolic phenotyping

Metabolic phenotyping is an important tool in translational biomedical research. The advanced analytical technologies commonly used for phenotyping, including mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy, generate complex data requiring tailored statistical analysis methods. Detailed protocols have been published for data acquisition by liquid NMR, solid-state NMR, ultra-performance liquid chromatography (LC-)MS and gas chromatography (GC-)MS on biofluids or tissues and their preprocessing. Here we propose an efficient protocol (guidelines and software) for statistical analysis of metabolic data generated by these methods. Code for all steps is provided, and no prior coding skill is necessary. We offer efficient solutions for the different steps required within the complete phenotyping data analytics workflow: scaling, normalization, outlier detection, multivariate analysis to explore and model study-related effects, selection of candidate biomarkers, validation, multiple testing correction and performance evaluation of statistical models. We also provide a statistical power calculation algorithm and safeguards to ensure robust and meaningful experimental designs that deliver reliable results. We exemplify the protocol with a two-group classification study and data from an epidemiological cohort; however, the protocol can be easily modified to cover a wider range of experimental designs or incorporate different modeling approaches. This protocol describes a minimal set of analyses needed to rigorously investigate typical datasets encountered in metabolic phenotyping.

Chronic Kidney Disease Cohort Studies: A Guide to Metabolome Analyses

Kidney diseases still pose one of the biggest challenges for global health, and their heterogeneity and often high comorbidity load seriously hinders the unraveling of their underlying pathomechanisms and the delivery of optimal patient care. Metabolomics, the quantitative study of small organic compounds, called metabolites, in a biological specimen, is gaining more and more importance in nephrology research. Conducting a metabolomics study in human kidney disease cohorts, however, requires thorough knowledge about the key workflow steps: study planning, sample collection, metabolomics data acquisition and preprocessing, statistical/bioinformatics data analysis, and results interpretation within a biomedical context. This review provides a guide for future metabolomics studies in human kidney disease cohorts. We will offer an overview of important a priori considerations for metabolomics cohort studies, available analytical as well as statistical/bioinformatics data analysis techniques, and subsequent interpretation of metabolic findings. We will further point out potential research questions for metabolomics studies in the context of kidney diseases and summarize the main results and data availability of important studies already conducted in this field.

Improved Automated Quantification Algorithm (AQuA) and Its Application to NMR-Based Metabolomics of EDTA-Containing Plasma

We have recently presented an Automated Quantification Algorithm (AQuA) and demonstrated its utility for rapid and accurate absolute metabolite quantification in ¹H NMR spectra in which positions and line widths of signals were predicted from a constant metabolite spectral library. The AQuA quantifies based on one preselected signal per metabolite and employs library spectra to model interferences from other metabolite signals. However, for some types of spectra, the interspectral deviations of signal positions and line widths can be pronounced; hence, interferences cannot be modeled using a constant spectral library. We here address this issue and present an improved AQuA that handles interspectral deviations. The improved AQuA monitors and characterizes the appearance of specific signals in each spectrum and automatically adjusts the spectral library to model interferences accordingly. The performance of the improved AQuA was tested on a large data set from plasma samples collected using ethylenediaminetetraacetic acid (EDTA) as an anticoagulant (n = 772). These spectra provided a suitable test system for the improved AQuA since EDTA signals (i) vary in intensity, position, and line width between spectra and (ii) interfere with many signals from plasma metabolites targeted for quantification (n = 54). Without the improvement, ca. 20 out of the 54 metabolites would have been overestimated. This included acetylcarnitine and ornithine, which are considered particularly difficult to quantify with ¹H NMR in EDTA-containing plasma. Furthermore, the improved AQuA performed rapidly (<10 s for all spectra). We believe that the improved AQuA provides a basis for automated quantification in other data sets where specific signals show interspectral deviations.

NMR spectroscopy in wine authentication: An official control perspective

Wine authentication is vital in identifying malpractice and fraud, and various physical and chemical analytical techniques have been employed for this purpose. Besides wet chemistry, these include chromatography, isotopic ratio mass spectrometry, optical spectroscopy, and nuclear magnetic resonance (NMR) spectroscopy, which have been applied in recent years in combination with chemometric approaches. For many years, ² H NMR spectroscopy was the method of choice and achieved official recognition in the detection of sugar addition to grape products. Recently, ¹ H NMR spectroscopy, a simpler and faster method (in terms of sample preparation), has gathered more and more attention in wine analysis, even if it still lacks official recognition. This technique makes targeted quantitative determination of wine ingredients and nontargeted detection of the metabolomic fingerprint of a wine sample possible. This review summarizes the possibilities and limitations of ¹ H NMR spectroscopy in analytical wine authentication, by reviewing its applications as reported in the literature. Examples of commercial and open-source solutions combining NMR spectroscopy and chemometrics are also examined herein, together with its opportunities of becoming an official method.

Joint Automatic Metabolite Identification and Quantification of a Set of 1H NMR Spectra

Metabolomics is a promising approach to characterize phenotypes or to identify biomarkers. It is also easily accessible through NMR, which can provide a comprehensive understanding of the metabolome of any living organisms. However, the analysis of ¹H NMR spectrum remains difficult, mainly due to the different problems encountered to perform automatic identification and quantification of metabolites in a reproducible way. In addition, methods that perform automatic identification and quantification of metabolites are often designed to process one given complex mixture spectrum at a time. Hence, when a set of complex mixture spectra coming from the same experiment has to be processed, the approach is simply repeated independently for every spectrum, despite their resemblance. Here, we present new methods that are the first to either align spectra or to identify and quantify metabolites by integrating information coming from several complex spectra of the same experiment. The performances of these new methods are then evaluated on both simulated and real datasets. The results show an improvement in the metabolite identification and in the accuracy of metabolite quantifications, especially when the concentration is low. This joint procedure is available in version 2.0 of ASICS package.

Affine Transformation of Negative Values for NMR Metabolomics Using the mrbin R Package

Data from untargeted metabolomics studies employing nuclear magnetic resonance (NMR) spectroscopy oftentimes contain negative values. These negative values hamper data processing and analysis algorithms and prevent the use of such data in multiomics integration settings. New methods to deal with such negative values are thus an urgent need in the metabolomics community. This study presents affine transformation of negative values (ATNV), a novel algorithm for replacement of negative values in NMR data sets. ATNV was implemented in the R package mrbin, which features interactive menus for user-friendly application and is available for free for various operating systems within the free R statistical programming language. The novel algorithms were tested on a set of human urinary NMR spectra and were able to successfully identify relevant metabolites.

AlpsNMR: an R package for signal processing of fully untargeted NMR-based metabolomics

SUMMARY

Nuclear magnetic resonance (NMR)-based metabolomics is widely used to obtain metabolic fingerprints of biological systems. While targeted workflows require previous knowledge of metabolites, prior to statistical analysis, untargeted approaches remain a challenge. Computational tools dealing with fully untargeted NMR-based metabolomics are still scarce or not user-friendly. Therefore, we developed AlpsNMR (Automated spectraL Processing System for NMR), an R package that provides automated and efficient signal processing for untargeted NMR metabolomics. AlpsNMR includes spectra loading, metadata handling, automated outlier detection, spectra alignment and peak-picking, integration and normalization. The resulting output can be used for further statistical analysis. AlpsNMR proved effective in detecting metabolite changes in a test case. The tool allows less experienced users to easily implement this workflow from spectra to a ready-to-use dataset in their routines.

AVAILABILITY AND IMPLEMENTATION

The AlpsNMR R package and tutorial is freely available to download from http://github.com/sipss/AlpsNMR under the MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exploring feature selection of St John's wort grown under different light spectra using 1 H-NMR spectroscopy

INTRODUCTION

Nuclear magnetic resonance (NMR) spectroscopy combined with multivariate statistical analysis can provide tools to help detect differences in plant chemistry when grown under varying conditions. Hypericum perforatum, or Saint John's wort, plants are a suitable model to explore methods of discrimination between early stage plants grown in different conditions.

OBJECTIVES

The purpose of this work was to develop a method for identifying differences in chemical profiles between young Hypericum perforatum plants grown under different lighting conditions.

MATERIAL AND METHODS

Cuttings were grown for 3 weeks under different light conditions. Plant extracts were prepared in MeOD-d₄ and analysed by ¹ H-NMR. A multivariate analysis method of the NMR data was developed in an effort to determine variations in chemical profiles.

RESULTS

The method identified specific metabolites as drivers of difference between the plants grown under different light conditions. STOCSY (statistical total correlation spectroscopy) and quantification of highlighted metabolites supported the findings of the multivariate analysis. Glutamine, sucrose and fructose were found to be chemical markers of light quality in this study.

CONCLUSION

NMR metabolomics using a medium field instrument could find differences in plant chemistry when grown in different conditions. This method could easily be extended to benchtop instruments and be used for crop monitoring and growth condition optimisation.

ASICS: an R package for a whole analysis workflow of 1D 1H NMR spectra

MOTIVATION

In metabolomics, the detection of new biomarkers from Nuclear Magnetic Resonance (NMR) spectra is a promising approach. However, this analysis remains difficult due to the lack of a whole workflow that handles spectra pre-processing, automatic identification and quantification of metabolites and statistical analyses, in a reproducible way.

RESULTS

We present ASICS, an R package that contains a complete workflow to analyse spectra from NMR experiments. It contains an automatic approach to identify and quantify metabolites in a complex mixture spectrum and uses the results of the quantification in untargeted and targeted statistical analyses. ASICS was shown to improve the precision of quantification in comparison to existing methods on two independent datasets. In addition, ASICS successfully recovered most metabolites that were found important to explain a two level condition describing the samples by a manual and expert analysis based on bucketing. It also found new relevant metabolites involved in metabolic pathways related to risk factors associated with the condition.

AVAILABILITY AND IMPLEMENTATION

ASICS is distributed as an R package, available on Bioconductor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A comparison of high-throughput plasma NMR protocols for comparative untargeted metabolomics

INTRODUCTION

When analyzing the human plasma metabolome with Nuclear Magnetic Resonance (NMR) spectroscopy, the Carr-Purcell-Meiboom-Gill (CPMG) experiment is commonly employed for large studies. However, this process can lead to compromised statistical analyses due to residual macromolecule signals. In addition, the utilization of Trimethylsilylpropanoic acid (TSP) as an internal standard often leads to quantification issues, and binning, as a spectral summarization step, can result in features not clearly assignable to metabolites.

OBJECTIVES

Our aim was to establish a new complete protocol for large plasma cohorts collected with the purpose of describing the comparative metabolic profile of groups of samples.

METHODS

We compared the conventional CPMG approach to a novel procedure that involves diffusion NMR, using the Longitudinal Eddy-Current Delay (LED) experiment, maleic acid (MA) as the quantification reference and peak picking for spectral reduction. This comparison was carried out using the ultrafiltration method as a gold standard in a simple sample classification experiment, with Partial Least Squares-Discriminant Analysis (PLS-DA) and the resulting metabolic signatures for multivariate data analysis. In addition, the quantification capabilities of the method were evaluated.

RESULTS

We found that the LED method applied was able to detect more metabolites than CPMG and suppress macromolecule signals more efficiently. The complete protocol was able to yield PLS-DA models with enhanced classification accuracy as well as a more reliable set of important features than the conventional CPMG approach. Assessment of the quantitative capabilities of the method resulted in good linearity, recovery and agreement with an established amino acid assay for the majority of the metabolites tested. Regarding repeatability, ~ 85% of all peaks had an adequately low coefficient of variation (< 30%) in replicate samples.

CONCLUSION

Overall, our comparison yielded a high-throughput untargeted plasma NMR protocol for optimized data acquisition and processing that is expected to be a valuable contribution in the field of metabolic biomarker discovery.

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.

The metaRbolomics Toolbox in Bioconductor and beyond

Metabolomics aims to measure and characterise the complex composition of metabolites in a biological system. Metabolomics studies involve sophisticated analytical techniques such as mass spectrometry and nuclear magnetic resonance spectroscopy, and generate large amounts of high-dimensional and complex experimental data. Open source processing and analysis tools are of major interest in light of innovative, open and reproducible science. The scientific community has developed a wide range of open source software, providing freely available advanced processing and analysis approaches. The programming and statistics environment R has emerged as one of the most popular environments to process and analyse Metabolomics datasets. A major benefit of such an environment is the possibility of connecting different tools into more complex workflows. Combining reusable data processing R scripts with the experimental data thus allows for open, reproducible research. This review provides an extensive overview of existing packages in R for different steps in a typical computational metabolomics workflow, including data processing, biostatistics, metabolite annotation and identification, and biochemical network and pathway analysis. Multifunctional workflows, possible user interfaces and integration into workflow management systems are also reviewed. In total, this review summarises more than two hundred metabolomics specific packages primarily available on CRAN, Bioconductor and GitHub.

Advances and challenges in development of precision psychiatry through clinical metabolomics on mood and psychotic disorders

Metabolomics is defined as the study of the global metabolite profile in a system under a given set of conditions. The objective of this review is to comprehensively assess the literature on metabolomics in mood disorders and schizophrenia and provide data for mental health researchers about the challenges and potentials of metabolomics. The majority of studies in metabolomics in Psychiatry uses peripheral blood or urine. The most widely used analytical techniques in metabolomics research are nuclear magnetic resonance (NMR) and mass spectrometry (MS). They are multiparametric and provide extensive structural and conformational information on multiple chemical classes. NMR is useful in untargeted analysis, which focuses on biosignatures or 'metabolic fingerprints' of illnesses. MS targeted metabolomics approach focuses on the identification and quantification of selected metabolites known to be involved in a particular metabolic pathway. The available studies of metabolomics in Schizophrenia, Bipolar Disorder and Major Depressive Disorder suggest a potential in investigating metabolic pathways involved in these diseases' pathophysiology and response to treatment, as well as its potential in biomarkers identification.

Tools and resources for metabolomics research community: A 2017-2018 update

The scale at which MS- and NMR-based platforms generate metabolomics datasets for both research, core, and clinical facilities to address challenges in the various sciences-ranging from biomedical to agricultural-is underappreciated. Thus, metabolomics efforts spanning microbe, environment, plant, animal, and human systems have led to continual and concomitant growth of in silico resources for analysis and interpretation of these datasets. These software tools, resources, and databases drive the field forward to help keep pace with the amount of data being generated and the sophisticated and diverse analytical platforms that are being used to generate these metabolomics datasets. To address challenges in data preprocessing, metabolite annotation, statistical interrogation, visualization, interpretation, and integration, the metabolomics and informatics research community comes up with hundreds of tools every year. The purpose of the present review is to provide a brief and useful summary of more than 95 metabolomics tools, software, and databases that were either developed or significantly improved during 2017-2018. We hope to see this review help readers, developers, and researchers to obtain informed access to these thorough lists of resources for further improvisation, implementation, and application in due course of time.