Untargeted metabolomics suffers from incomplete raw data processing

Introduction

Untargeted metabolomics is a powerful tool for biological discoveries. To analyze the complex raw data, significant advances in computational approaches have been made, yet it is not clear how exhaustive and reliable the data analysis results are.

Objectives

Assessment of the quality of raw data processing in untargeted metabolomics.

Methods

Five published untargeted metabolomics studies, were reanalyzed.

Results

Omissions of at least 50 relevant compounds from the original results as well as examples of representative mistakes were reported for each study.

Conclusion

Incomplete raw data processing shows unexplored potential of current and legacy data.

     

Urinary antihypertensive drug metabolite screening using molecular networking coupled to high-resolution mass spectrometry fragmentation

Introduction

Mass spectrometry is the current technique of choice in studying drug metabolism. High-resolution mass spectrometry in combination with MS/MS gas-phase experiments has the potential to contribute to rapid advances in this field. However, the data emerging from such fragmentation spectral files pose challenges to downstream analysis, given their complexity and size.

Objectives

This study aims to detect and visualize antihypertensive drug metabolites in untargeted metabolomics experiments based on the spectral similarity of their fragmentation spectra. Furthermore, spectral clusters of endogenous metabolites were also examined.

Methods

Here we apply a molecular networking approach to seek drugs and their metabolites, in fragmentation spectra from urine derived from a cohort of 26 patients on antihypertensive therapy. The mass spectrometry data was collected on a Thermo Q-Exactive coupled to pHILIC chromatography using data dependent analysis (DDA) MS/MS gas-phase experiments.

Results

In total, 165 separate drug metabolites were found and structurally annotated (17 by spectral matching and 122 by classification based on a clustered fragmentation pattern). The clusters could be traced to 13 drugs including the known antihypertensives verapamil, losartan and amlodipine. The molecular networking approach also generated clusters of endogenous metabolites, including carnitine derivatives, and conjugates containing glutamine, glutamate and trigonelline.

Conclusions

The approach offers unprecedented capability in the untargeted identification of drugs and their metabolites at the population level and has great potential to contribute to understanding stratified responses to drugs where differences in drug metabolism may determine treatment outcome.

     

A metabolomics guided exploration of marine natural product chemical space

Introduction

Natural products from culture collections have enormous impact in advancing discovery programs for metabolites of biotechnological importance. These discovery efforts rely on the metabolomic characterization of strain collections.

Objective

Many emerging approaches compare metabolomic profiles of such collections, but few enable the analysis and prioritization of thousands of samples from diverse organisms while delivering chemistry specific read outs.

Method

In this work we utilize untargeted LC–MS/MS based metabolomics together with molecular networking to inventory the chemistries associated with 1000 marine microorganisms.

Result

This approach annotated 76 molecular families (a spectral match rate of 28 %), including clinically and biotechnologically important molecules such as valinomycin, actinomycin D, and desferrioxamine E. Targeting a molecular family produced primarily by one microorganism led to the isolation and structure elucidation of two new molecules designated maridric acids A and B.

Conclusion

Molecular networking guided exploration of large culture collections allows for rapid dereplication of know molecules and can highlight producers of uniques metabolites. These methods, together with large culture collections and growing databases, allow for data driven strain prioritization with a focus on novel chemistries.

     

An untargeted metabolomics method for archived newborn dried blood spots in epidemiologic studies

Introduction

For pediatric diseases like childhood leukemia, a short latency period points to in-utero exposures as potentially important risk factors. Untargeted metabolomics of small molecules in archived newborn dried blood spots (DBS) offers an avenue for discovering early-life exposures that contribute to disease risks.

Objectives

The purpose of this study was to develop a quantitative method for untargeted analysis of archived newborn DBS for use in an epidemiological study (California Childhood Leukemia Study, CCLS).

Methods

Using experimental DBS from the blood of an adult volunteer, we optimized extraction of small molecules and integrated measurement of potassium as a proxy for blood hematocrit. We then applied this extraction method to 4.7-mm punches from 106 control DBS samples from the CCLS. Sample extracts were analyzed with liquid chromatography—high resolution mass spectrometry (LC-HRMS) and an untargeted workflow was used to screen for metabolites that discriminate population characteristics such as sex, ethnicity, and birth weight.

Results

Thousands of small molecules were measured in extracts of archived DBS. Normalizing for potassium levels removed variability related to varying hematocrit across DBS punches. Of the roughly 1000 prevalent small molecules that were tested, multivariate linear regression detected significant associations with ethnicity (three metabolites) and birth weight (15 metabolites) after adjusting for multiple testing.

Conclusions

This untargeted workflow can be used for analysis of small molecules in archived DBS to discover novel biomarkers, to provide insights into the initiation and progression of diseases, and to provide guidance for disease prevention.

     

Large-scale untargeted LC-MS metabolomics data correction using between-batch feature alignment and cluster-based within-batch signal intensity drift correction

Introduction

Liquid chromatography-mass spectrometry (LC-MS) is a commonly used technique in untargeted metabolomics owing to broad coverage of metabolites, high sensitivity and simple sample preparation. However, data generated from multiple batches are affected by measurement errors inherent to alterations in signal intensity, drift in mass accuracy and retention times between samples both within and between batches. These measurement errors reduce repeatability and reproducibility and may thus decrease the power to detect biological responses and obscure interpretation.

Objective

Our aim was to develop procedures to address and correct for within- and between-batch variability in processing multiple-batch untargeted LC-MS metabolomics data to increase their quality.

Methods

Algorithms were developed for: (i) alignment and merging of features that are systematically misaligned between batches, through aggregating feature presence/missingness on batch level and combining similar features orthogonally present between batches; and (ii) within-batch drift correction using a cluster-based approach that allows multiple drift patterns within batch. Furthermore, a heuristic criterion was developed for the feature-wise choice of reference-based or population-based between-batch normalisation.

Results

In authentic data, between-batch alignment resulted in picking 15 % more features and deconvoluting 15 % of features previously erroneously aligned. Within-batch correction provided a decrease in median quality control feature coefficient of variation from 20.5 to 15.1 %. Algorithms are open source and available as an R package (‘batchCorr’).

Conclusions

The developed procedures provide unbiased measures of improved data quality, with implications for improved data analysis. Although developed for LC-MS based metabolomics, these methods are generic and can be applied to other data suffering from similar limitations.

     

Optimization of fecal sample preparation for untargeted LC-HRMS based metabolomics

Introduction

Collecting feces is easy. It offers direct outcome to endogenous and microbial metabolites.

Objectives

In a context of lack of consensus about fecal sample preparation, especially in animal species, we developed a robust protocol allowing untargeted LC-HRMS fingerprinting.

Methods

The conditions of extraction (quantity, preparation, solvents, dilutions) were investigated in bovine feces.

Results

A rapid and simple protocol involving feces extraction with methanol (1/3, M/V) followed by centrifugation and a step filtration (10 kDa) was developed.

Conclusion

The workflow generated repeatable and informative fingerprints for robust metabolome characterization.

     

ASICS: an automatic method for identification and quantification of metabolites in complex 1D <Superscript>1</Superscript>H NMR spectra

Introduction

Experiments in metabolomics rely on the identification and quantification of metabolites in complex biological mixtures. This remains one of the major challenges in NMR/mass spectrometry analysis of metabolic profiles. These features are mandatory to make metabolomics asserting a general approach to test a priori formulated hypotheses on the basis of exhaustive metabolome characterization rather than an exploratory tool dealing with unknown metabolic features.

Objectives

In this article we propose a method, named ASICS, based on a strong statistical theory that handles automatically the metabolites identification and quantification in proton NMR spectra.

Methods

A statistical linear model is built to explain a complex spectrum using a library containing pure metabolite spectra. This model can handle local or global chemical shift variations due to experimental conditions using a warping function. A statistical lasso-type estimator identifies and quantifies the metabolites in the complex spectrum. This estimator shows good statistical properties and handles peak overlapping issues.

Results

The performances of the method were investigated on known mixtures (such as synthetic urine) and on plasma datasets from duck and human. Results show noteworthy performances, outperforming current existing methods.

Conclusion

ASICS is a completely automated procedure to identify and quantify metabolites in ¹H NMR spectra of biological mixtures. It will enable empowering NMR-based metabolomics by quickly and accurately helping experts to obtain metabolic profiles.

     

Untargeted metabolomics approach using UPLC-ESI-QTOF-MS to explore the metabolome of fresh-cut iceberg lettuce

Introduction

The shelf-life of fresh-cut lettuce after storage is limited by several factors that affect its quality and lead to consumer rejection. Different metabolic events occur after cutting as an abiotic stress response.

Objectives

This study aims to explore the metabolome of iceberg lettuce and to understand the changes related to storage time and genetics applying an untargeted metabolomics approach.

Methods

Two cultivars with different browning susceptibility, fast-browning (FB) and slow-browning (SB), were analyzed by UPLC-ESI-QTOF-MS just after cutting (d0) and after five days of storage (d5). Extraction, metabolic profiling, and data-pretreatment procedures were optimized to obtain a robust and reliable data set.

Results

Preliminary principal component analysis and hierarchical cluster analysis of the full dataset [around 8551 extracted, aligned and filtered molecular features (MFs)] showed a clear separation between the different samples (FB-d0, FB-d5, SB-d0, and SB-d5), highlighting a clear storage time-dependent effect. After statistical analysis applying Student’s t test, 536 MFs were detected as significantly different between d0 and d5 of storage in FB and 633 in SB. Some of them (221) were common to both cultivars. Out of these significant compounds, 22 were tentatively identified by matching their molecular formulae with those previously reported in the literature. Five families of metabolites were detected: amino acids, phenolic compounds, sesquiterpene lactones, fatty acids, and lysophospholipids. All compounds showed a clear trend to decrease at d5 except phenolic compounds that increased after storage.

Conclusion

The untargeted metabolomics analysis is a powerful tool for characterizing the changes on lettuce metabolome associated with cultivar and especially with storage time. Some families of compounds affected by storage time were reported to be closely related to quality loss.

     

Untargeted and stable isotope-assisted metabolomic analysis of MDA-MB-231 cells under hypoxia

Introduction

Hypoxia commonly occurs in cancers and is highly related with the occurrence, development and metastasis of cancer. Treatment of triple negative breast cancer remains challenge. Knowledge about the metabolic status of triple negative breast cancer cell lines in hypoxia is valuable for the understanding of molecular mechanisms of this tumor subtype to develop effective therapeutics.

Objectives

Comprehensively characterize the metabolic profiles of triple negative breast cancer cell line MDA-MB-231 in normoxia and hypoxia and the pathways involved in metabolic changes in hypoxia.

Methods

Differences in metabolic profiles affected pathways of MDA-MB-231 cells in normoxia and hypoxia were characterized using GC–MS based untargeted and stable isotope assisted metabolomic techniques.

Results

Thirty-three metabolites were significantly changed in hypoxia and nine pathways were involved. Hypoxia increased glycolysis, inhibited TCA cycle, pentose phosphate pathway and pyruvate carboxylation, while increased glutaminolysis in MDA-MB-231 cells.

Conclusion

The current results provide metabolic differences of MDA-MB-231 cells in normoxia and hypoxia conditions as well as the involved metabolic pathways, demonstrating the power of combined use of untargeted and stable isotope-assisted metabolomic methods in comprehensive metabolomic analysis.

     

Bayesian estimation of the number of protonation sites for urinary metabolites from NMR spectroscopic data

Introduction

To aid the development of better algorithms for \(^1\)H NMR data analysis, such as alignment or peak-fitting, it is important to characterise and model chemical shift changes caused by variation in pH. The number of protonation sites, a key parameter in the theoretical relationship between pH and chemical shift, is traditionally estimated from the molecular structure, which is often unknown in untargeted metabolomics applications.

Objective

We aim to use observed NMR chemical shift titration data to estimate the number of protonation sites for a range of urinary metabolites.

Methods

A pool of urine from healthy subjects was titrated in the range pH 2–12, standard \(^1\)H NMR spectra were acquired and positions of 51 peaks (corresponding to 32 identified metabolites) were recorded. A theoretical model of chemical shift was fit to the data using a Bayesian statistical framework, using model selection procedures in a Markov Chain Monte Carlo algorithm to estimate the number of protonation sites for each molecule.

Results

The estimated number of protonation sites was found to be correct for 41 out of 51 peaks. In some cases, the number of sites was incorrectly estimated, due to very close pKa values or a limited amount of data in the required pH range.

Conclusions

Given appropriate data, it is possible to estimate the number of protonation sites for many metabolites typically observed in \(^1\)H NMR metabolomics without knowledge of the molecular structure. This approach may be a valuable resource for the development of future automated metabolite alignment, annotation and peak fitting algorithms.

     

KniMet: a pipeline for the processing of chromatography–mass spectrometry metabolomics data

Introduction

Data processing is one of the biggest problems in metabolomics, given the high number of samples analyzed and the need of multiple software packages for each step of the processing workflow.

Objectives

Merge in the same platform the steps required for metabolomics data processing.

Methods

KniMet is a workflow for the processing of mass spectrometry-metabolomics data based on the KNIME Analytics platform.

Results

The approach includes key steps to follow in metabolomics data processing: feature filtering, missing value imputation, normalization, batch correction and annotation.

Conclusion

KniMet provides the user with a local, modular and customizable workflow for the processing of both GC–MS and LC–MS open profiling data.

     

A causal network analysis in an observational study identifies metabolomics pathways influencing plasma triglyceride levels

Introduction

Plasma triglyceride levels are a risk factor for coronary heart disease. Triglyceride metabolism is well characterized, but challenges remain to identify novel paths to lower levels. A metabolomics analysis may help identify such novel pathways and, therefore, provide hints about new drug targets.

Objectives

In an observational study, causal relationships in the metabolomics level of granularity are taken into account to distinguish metabolites and pathways having a direct effect on plasma triglyceride levels from those which are only associated with or have indirect effect on triglyceride.

Method

The analysis began by leveraging near-complete information from the genome level of granularity using the GDAG algorithm to identify a robust causal network over 122 metabolites in an upper level of granularity. Knowing the metabolomics causal relationships, we enter the triglyceride variable in the model to identify metabolites with direct effect on plasma triglyceride levels. We carried out the same analysis on triglycerides measured over five different visits spanning 24 years.

Result

Nine metabolites out of 122 metabolites under consideration influenced directly plasma triglyceride levels. Given these nine metabolites, the rest of metabolites in the study do not have a significant effect on triglyceride levels at significance level alpha = 0.001. Therefore, for the further analysis and interpretations about triglyceride levels, the focus should be on these nine metabolites out of 122 metabolites in the study. The metabolites with the strongest effects at the baseline visit were arachidonate and carnitine, followed by 9-hydroxy-octadecadenoic acid and palmitoylglycerophosphoinositol. The influence of arachidonate on triglyceride levels remained significant even at the fourth visit, which was 10 years after the baseline visit.

Conclusion

These results demonstrate the utility of integrating multi-omics data in a granularity framework to identify novel candidate pathways to lower risk factor levels.

     

Metabolic profiling of <Emphasis Type="Italic">Daphnia magna</Emphasis> exposed to environmental stressors by GC–MS and chemometric tools

Introduction

Climate change is a major concern for the scientific community, demanding novel information about the effects of environmental stressors on living organisms. Metabolic profiling is required for achieving the most extensive possible range of compounds and their concentration changes on stressed conditions.

Objectives

Individuals of the crustacean species Daphnia magna were exposed to three different abiotic factors linked to global climate change: high salinity, high temperature levels and hypoxia. Advanced chemometric tools were used to characterize the metabolites affected by the exposure.

Method

An exploratory analysis of gas chromatography-mass spectrometry (GC–MS) data was performed to discriminate between control and exposed daphnid samples. Due to the complexity of these GC–MS data sets, a comprehensive untargeted analysis of the full scan data was performed using multivariate curve resolution-alternating least squares (MCR-ALS) method. This approach enabled to resolve most of the metabolite signals from interference peaks caused by derivatization reactions. Metabolites with significant changes in their peak areas were tentatively identified and the involved metabolic pathways explored.

Results

D. magna metabolic biomarkers are proposed for the considered physical factors. Metabolites related with energy metabolic pathways including some amino acids, carbohydrates, organic acids and nucleosides were identified as potential biomarkers of the investigated treatments.

Conclusions

The proposed untargeted GC–MS metabolomics strategy and multivariate data analysis tools were useful to investigate D. magna metabolome under environmental stressed conditions.

     

Investigation of altered urinary metabolomic profiles of invasive ductal carcinoma of breast using targeted and untargeted approaches

Introduction

Invasive ductal carcinoma (IDC) is a type of breast cancer, usually detected in advanced stages due to its asymptomatic nature which ultimately leads to low survival rate. Identification of urinary metabolic adaptations induced by IDC to understand the disease pathophysiology and monitor therapy response would be a helpful approach in clinical settings. Moreover, its non-invasive and cost effective strategy better suited to minimize apprehension among high risk population.

Objective

This study aims toward investigating the urinary metabolic alterations of IDC by targeted (LC-MRM/MS) and untargeted (GC–MS) approaches for the better understanding of the disease pathophysiology and monitoring therapy response.

Methods

Urinary metabolic alterations of IDC subjects (63) and control subjects (63) were explored by targeted (LC-MRM/MS) and untargeted (GC–MS) approaches. IDC specific urinary metabolomics signature was extracted by applying both univariate and multivariate statistical tools.

Results

Statistical analysis identified 39 urinary metabolites with the highest contribution to metabolomic alterations specific to IDC. Out of which, 19 metabolites were identified from targeted LC-MRM/MS analysis, while 20 were identified from the untargeted GC–MS analysis. Receiver operator characteristic (ROC) curve analysis evidenced 6 most discriminatory metabolites from each type of approach that could differentiate between IDC subjects and controls with higher sensitivity and specificity. Furthermore, metabolic pathway analysis depicted several dysregulated pathways in IDC including sugar, amino acid, nucleotide metabolism, TCA cycle etc.

Conclusions

Overall, this study provides valuable inputs regarding altered urinary metabolites which improved our knowledge on urinary metabolomic alterations induced by IDC. Moreover, this study identified several dysregulated metabolic pathways which offer further insight into the disease pathophysiology.

     

A review of metabolism-associated biomarkers in lung cancer diagnosis and treatment

Introduction

Lung cancer continues to be the leading cause of cancer-related mortality worldwide. Early detection has proven essential to extend survival. Genomic and proteomic advances have provided impetus to the effort dedicated to detect and diagnose the disease at an earlier stage. Recently, the study of metabolites associated with tumor formation and progression has inaugurated the era of cancer metabolomics to aid in this effort.

Objectives

This review summarizes recent work regarding novel metabolites with the potential to serve as biomarkers for early lung tumor detection, evaluation of disease progression, and prediction of patient outcomes.

Method

We compare the metabolite profiling of cancer patients with that of healthy individuals, and the metabolites identified in tissue and biofluid samples and their usefulness as lung cancer biomarkers. We discuss metabolite alterations in tumor versus paired non-tumor lung tissues, as well as metabolite alterations in different stages of lung cancers and their usefulness as indicators of disease progression and overall survival. We evaluate metabolite dysregulation in different types of lung cancers, and those associated with lung cancer versus other lung diseases. We also examine metabolite differences between lung cancer patients and smokers/risk-factor individuals.

Result

Although an extensive list of metabolites has been evaluated to distinguish between these cases, refinement of methods is further required for adequate patient diagnosis and treatment.

Conclusion

We conclude that with technological advancement, metabolomics may be able to replace more invasive and costly diagnostic procedures while also providing the means to more effectively tailor treatment to patient-specific tumors.

     

Urinary Metabolomic Signatures in Obstructive Sleep Apnea through Targeted Metabolomic Analysis: A Pilot Study

Introduction

Obstructive sleep apnea (OSA) is very common sleep problem, and it is associated with serious morbidities such as cardiovascular diseases and metabolic diseases. Overnight polysomnography (PSG) is the gold standard test for OSA, but it is expensive and requires specific facilities and equipment. Thus, novel screening methods are needed for effective diagnosis and follow-up in OSA.

Objectives

The aims of the study were to investigate the urinary metabolic signatures and identify potential urine markers for OSA using a mass spectrometry (MS)-based assay for targeted metabolomics.

Methods

Urine samples were collected from 48 male subjects who visited a sleep clinic for suspicious OSA. All underwent overnight in-laboratory polysomnography. The Biocrates AbsoluteIDQ p180 kit was used for targeted metabolomics.

Results

Among the 86 metabolites quantified, three acylcarnitines, one biogenic amine, two glycerophospholipids, and two sphingomyelins were differently expressed in OSA patients [apnea-hypopnea index (AHI) ≥5] compared with control groups (AHI <5 and/or simple snoring with no other sleep disorders). Additional partial correlation and multivariate logistic regression analysis revealed that long-chain acylcarnitine C14:1, symmetric dimethylarginine, and sphingomyelin C18:1 might be potential biomarkers for OSA. Receiver operating characteristic analysis showed favorable predictive properties of these metabolites. Furthermore, a combination of the metabolites exceeding cutoff values yielded further improved sensitivity or specificity.

Conclusions

MS-based targeted metabolomics identified specific classes of urinary metabolites that were up-regulated in OSA patients. Further assessments in large populations are required to clarify the screening values of these metabolite markers.

     

From correlation to causation: analysis of metabolomics data using systems biology approaches

Introduction

Metabolomics is a well-established tool in systems biology, especially in the top–down approach. Metabolomics experiments often results in discovery studies that provide intriguing biological hypotheses but rarely offer mechanistic explanation of such findings. In this light, the interpretation of metabolomics data can be boosted by deploying systems biology approaches.

Objectives

This review aims to provide an overview of systems biology approaches that are relevant to metabolomics and to discuss some successful applications of these methods.

Methods

We review the most recent applications of systems biology tools in the field of metabolomics, such as network inference and analysis, metabolic modelling and pathways analysis.

Results

We offer an ample overview of systems biology tools that can be applied to address metabolomics problems. The characteristics and application results of these tools are discussed also in a comparative manner.

Conclusions

Systems biology-enhanced analysis of metabolomics data can provide insights into the molecular mechanisms originating the observed metabolic profiles and enhance the scientific impact of metabolomics studies.

     

Metabolite secretion in microorganisms: the theory of metabolic overflow put to the test

Introduction

Microbial cells secrete many metabolites during growth, including important intermediates of the central carbon metabolism. This has not been taken into account by researchers when modeling microbial metabolism for metabolic engineering and systems biology studies.

Materials and Methods

The uptake of metabolites by microorganisms is well studied, but our knowledge of how and why they secrete different intracellular compounds is poor. The secretion of metabolites by microbial cells has traditionally been regarded as a consequence of intracellular metabolic overflow.

Conclusions

Here, we provide evidence based on time-series metabolomics data that microbial cells eliminate some metabolites in response to environmental cues, independent of metabolic overflow. Moreover, we review the different mechanisms of metabolite secretion and explore how this knowledge can benefit metabolic modeling and engineering.

     

Characterizing alternative feeds for rainbow trout (<Emphasis Type="Italic">O. mykiss</Emphasis>) by <Superscript>1</Superscript>H NMR metabolomics

Introduction

Fish feed formulations are constantly evolving to improve the quality of diets for farmed fish and to ensure the sustainability of the aquaculture sector. Nowadays, insect, microalgae and yeast are feedstuff candidates for new feeds. However, the characterization of aquafeed is still based on proximate and targeted analyses which may not be sufficient to assess feed quality.

Objectives

Our aim was to highlight the soluble compounds that specifically differ between selected plant-based feeds complemented with alternative feedstuffs and discuss their origin and potential for fish nutrition.

Methods

A growth trial was carried out to evaluate growth performances and feed conversion ratios of fish fed plant-based, commercial, insect, spirulina and yeast feeds. ¹H NMR metabolomics profiling of each feed was performed using a CPMG sequence on polar extracts. Spectra were processed, and data were analyzed using multivariate and univariate analyses to compare alternative feeds to a plant-based feed.

Results

Fish fed insect or yeast feed showed the best growth performances associated with the lowest feed conversion ratios compared to plant-based feed. Soluble compound ¹H NMR profiles of insect and spirulina alternative feeds differed significantly from the plant-based one that clustered with yeast feed. In insect and spirulina feeds, specific differences compared to plant-based feed concerned glycerol and 3-hydroxybutyrate, respectively.

Conclusion

This strategy based on compositional differences between plant-based and alternative feeds can be useful for detecting compounds unsuspected until now that could impact fish metabolism.