Easy-Amanida: An R Shiny application for the meta-analysis of aggregate results in clinical metabolomics using Amanida and Webchem.

Meta-analysis is a useful tool in clinical research, as it combines the results of multiple clinical studies to improve precision when answering a particular scientific question. While there has been a substantial increase in publications using meta-analysis in various clinical research topics, the number of published meta-analyses in metabolomics is significantly lower compared to other omics disciplines. Metabolomics is the study of small chemical compounds in living organisms, which provides important insights into an organism's phenotype. However, the wide variety of compounds and the different experimental methods used in metabolomics make it challenging to perform a thorough meta-analysis. Additionally, there is a lack of consensus on reporting statistical estimates, and the high number of compound naming synonyms further complicates the process. Easy-Amanida is a new tool that combines two R packages, "amanida" and "webchem", to enable meta-analysis of aggregate statistical data, like p-value and fold-change, while ensuring the compounds naming harmonization. The Easy-Amanida app is implemented in Shiny, an R package add-on for interactive web apps, and provides a workflow to optimize the naming combination. This article describes all the steps to perform the meta-analysis using Easy-Amanida, including an illustrative example for interpreting the results. The use of aggregate statistics metrics extends the use of Easy-Amanida beyond the metabolomics field.

A Metabolites Merging Strategy (MMS): Harmonization to Enable Studies' Intercomparison.

Metabolomics encounters challenges in cross-study comparisons due to diverse metabolite nomenclature and reporting practices. To bridge this gap, we introduce the Metabolites Merging Strategy (MMS), offering a systematic framework to harmonize multiple metabolite datasets for enhanced interstudy comparability. MMS has three steps. Step 1: Translation and merging of the different datasets by employing InChIKeys for data integration, encompassing the translation of metabolite names (if needed). Followed by Step 2: Attributes' retrieval from the InChIkey, including descriptors of name (title name from PubChem and RefMet name from Metabolomics Workbench), and chemical properties (molecular weight and molecular formula), both systematic (InChI, InChIKey, SMILES) and non-systematic identifiers (PubChem, CheBI, HMDB, KEGG, LipidMaps, DrugBank, Bin ID and CAS number), and their ontology. Finally, a meticulous three-step curation process is used to rectify disparities for conjugated base/acid compounds (optional step), missing attributes, and synonym checking (duplicated information). The MMS procedure is exemplified through a case study of urinary asthma metabolites, where MMS facilitated the identification of significant pathways hidden when no dataset merging strategy was followed. This study highlights the need for standardized and unified metabolite datasets to enhance the reproducibility and comparability of metabolomics studies.

The untargeted urine volatilome for biomedical applications: methodology and volatilome database.

Chemically diverse in compounds, urine can give us an insight into metabolic breakdown products from foods, drinks, drugs, environmental contaminants, endogenous waste metabolites, and bacterial by-products. Hundreds of them are volatile compounds; however, their composition has never been provided in detail, nor has the methodology used for urine volatilome untargeted analysis. Here, we summarize key elements for the untargeted analysis of urine volatilome from a comprehensive compilation of literature, including the latest reports published. Current achievements and limitations on each process step are discussed and compared. 34 studies were found retrieving all information from the urine treatment to the final results obtained. In this report, we provide the first specific urine volatilome database, consisting of 841 compounds from 80 different chemical classes.

Urine NMR Metabolomics for Precision Oncology in Colorectal Cancer.

Metabolomics is a fundamental approach to discovering novel biomarkers and their potential use for precision medicine. When applied for population screening, NMR-based metabolomics can become a powerful clinical tool in precision oncology. Urine tests can be more widely accepted due to their intrinsic non-invasiveness. Our review provides the first exhaustive evaluation of NMR metabolomics for the determination of colorectal cancer (CRC) in urine. A specific search in PubMed, Web of Science, and Scopus was performed, and 10 studies met the required criteria. There were no restrictions on the query for study type, leading to not only colorectal cancer samples versus control comparisons, but also prospective studies of surgical effects. With this review, all compounds in the included studies were merged into a database. In doing so, we identified up to 100 compounds in urine samples, and 11 were found in at least three articles. Results were analyzed in three groups: case (CRC and adenomas)/control, pre-/post-surgery, and combining both groups. When combining the case-control and the pre-/post-surgery groups, up to twelve compounds were found to be relevant. Seven down-regulated metabolites in CRC were identified, creatinine, 4-hydroxybenzoic acid, acetone, carnitine, d-glucose, hippuric acid, l-lysine, l-threonine, and pyruvic acid, and three up-regulated compounds in CRC were identified, acetic acid, phenylacetylglutamine, and urea. The pathways and enrichment analysis returned only two pathways significantly expressed: the pyruvate metabolism and the glycolysis/gluconeogenesis pathway. In both cases, only the pyruvic acid (down-regulated in urine of CRC patients, with cancer cell proliferation effect in the tissue) and acetic acid (up-regulated in urine of CRC patients, with chemoprotective effect) were present.

Amanida: an R package for meta-analysis of metabolomics non-integral data.

SUMMARY: The combination, analysis and evaluation of different studies which try to answer or solve the same scientific question, also known as a meta-analysis, plays a crucial role in answering relevant clinical relevant questions. Unfortunately, metabolomics studies rarely disclose all the statistical information needed to perform a meta-analysis. Here, we present a meta-analysis approach using only the most reported statistical parameters in this field: P-value and fold-change. The P-values are combined via Fisher's method and fold-changes by averaging, both weighted by the study size (n). The amanida package includes several visualization options: a volcano plot for quantitative results, a vote plot for total regulation behaviours (up/down regulations) for each compound, and a explore plot of the vote-counting results with the number of times a compound is found upregulated or downregulated. In this way, it is very easy to detect discrepancies between studies at a first glance. AVAILABILITY AND IMPLEMENTATION: Amanida code and documentation are at CRAN and https://github.com/mariallr/amanida. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comprehensive Volatilome and Metabolome Signatures of Colorectal Cancer in Urine: A Systematic Review and Meta-Analysis.

To increase compliance with colorectal cancer screening programs and to reduce the recommended screening age, cheaper and easy non-invasiveness alternatives to the fecal immunochemical test should be provided. Following the PRISMA procedure of studies that evaluated the metabolome and volatilome signatures of colorectal cancer in human urine samples, an exhaustive search in PubMed, Web of Science, and Scopus found 28 studies that met the required criteria. There were no restrictions on the query for the type of study, leading to not only colorectal cancer samples versus control comparison but also polyps versus control and prospective studies of surgical effects, CRC staging and comparisons of CRC with other cancers. With this systematic review, we identified up to 244 compounds in urine samples (3 shared compounds between the volatilome and metabolome), and 10 of them were relevant in more than three articles. In the meta-analysis, nine studies met the criteria for inclusion, and the results combining the case-control and the pre-/post-surgery groups, eleven compounds were found to be relevant. Four upregulated metabolites were identified, 3-hydroxybutyric acid, L-dopa, L-histidinol, and N1, N12-diacetylspermine and seven downregulated compounds were identified, pyruvic acid, hydroquinone, tartaric acid, and hippuric acid as metabolites and butyraldehyde, ether, and 1,1,6-trimethyl-1,2-dihydronaphthalene as volatiles.

rMSIproc: an R package for mass spectrometry imaging data processing.

SUMMARY: Mass spectrometry imaging (MSI) can reveal biochemical information directly from a tissue section. MSI generates a large quantity of complex spectral data which is still challenging to translate into relevant biochemical information. Here, we present rMSIproc, an open-source R package that implements a full data processing workflow for MSI experiments performed using TOF or FT-based mass spectrometers. The package provides a novel strategy for spectral alignment and recalibration, which allows to process multiple datasets simultaneously. This enables to perform a confident statistical analysis with multiple datasets from one or several experiments. rMSIproc is designed to work with files larger than the computer memory capacity and the algorithms are implemented using a multi-threading strategy. rMSIproc is a powerful tool able to take full advantage of modern computer systems to completely develop the whole MSI potential. AVAILABILITY AND IMPLEMENTATION: rMSIproc is freely available at https://github.com/prafols/rMSIproc. CONTACT: pere.rafols@urv.cat. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Assessing the potential of sputtered gold nanolayers in mass spectrometry imaging for metabolomics applications.

Mass spectrometry imaging (MSI) is a molecular imaging technique that maps the distribution of molecules in biological tissues with high spatial resolution. The most widely used MSI modality is matrix-assisted laser desorption/ionization (MALDI), mainly due to the large variety of analyte classes amenable for MALDI analysis. However, the organic matrices used in classical MALDI may impact the quality of the molecular images due to limited lateral resolution and strong background noise in the low mass range, hindering its use in metabolomics. Here we present a matrix-free laser desorption/ionization (LDI) technique based on the deposition of gold nanolayers on tissue sections by means of sputter-coating. This gold coating method is quick, fully automated, reproducible, and allows growing highly controlled gold nanolayers, necessary for high quality and high resolution MS image acquisition. The performance of the developed method has been tested through the acquisition of MS images of brain tissues. The obtained spectra showed a high number of MS peaks in the low mass region (m/z below 1000 Da) with few background peaks, demonstrating the ability of the sputtered gold nanolayers of promoting the desorption/ionization of a wide range of metabolites. These results, together with the reliable MS spectrum calibration using gold peaks, make the developed method a valuable alternative for MSI applications.

Novel automated workflow for spectral alignment and mass calibration in MS imaging using a sputtered Ag nanolayer.

Mass spectrometry imaging (MSI) is a technique that can map analyte spatial distribution directly onto a tissue section. This enables the spatial correlation of molecular entities with a tissue morphology to be investigated. Analyte annotation in MSI is intrinsically linked to the mass accuracy of the data. Mass accuracy and analyte identification are determined by such factors as the experimental set up and the data processing workflow. We present an MSI data processing workflow that uses a label-free approach to compensate for mass shifts. The algorithms developed were designed to perform efficiently even for datasets much larger than computer's memory. Herein, we present the application of the developed processing workflow to a dataset with more than 13.000 pixels and ∼50.000 mass channels. We assessed the overall mass accuracy in the range m/z 400 to 1200 using silver and gold sputtered nanolayers. With our novel processing workflow we were able to obtain mass errors as low as 5 ppm using a TOF instrument.

Signal preprocessing, multivariate analysis and software tools for MA(LDI)-TOF mass spectrometry imaging for biological applications.

Mass spectrometry imaging (MSI) is a label-free analytical technique capable of molecularly characterizing biological samples, including tissues and cell lines. The constant development of analytical instrumentation and strategies over the previous decade makes MSI a key tool in clinical research. Nevertheless, most MSI studies are limited to targeted analysis or the mere visualization of a few molecular species (proteins, peptides, metabolites, or lipids) in a region of interest without fully exploiting the possibilities inherent in the MSI technique, such as tissue classification and segmentation or the identification of relevant biomarkers from an untargeted approach. MSI data processing is challenging due to several factors. The large volume of mass spectra involved in a MSI experiment makes choosing the correct computational strategies critical. Furthermore, pixel to pixel variation inherent in the technique makes choosing the correct preprocessing steps critical. The primary aim of this review was to provide an overview of the data-processing steps and tools that can be applied to an MSI experiment, from preprocessing the raw data to the more advanced strategies for image visualization and segmentation. This review is particularly aimed at researchers performing MSI experiments and who are interested in incorporating new data-processing features, improving their computational strategy, and/or desire access to data-processing tools currently available. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 37:281-306, 2018.

rMSI: an R package for MS imaging data handling and visualization.

SUMMARY: R platform provides some packages that are useful to process mass spectrometry imaging (MSI) data; however, none of them provide an easy to use graphical user interface (GUI). Here, we introduce rMSI, an R package for MSI data analysis focused on providing an efficient way to manage MSI data together with a GUI integrated in R environment. MS data is loaded in rMSI custom format optimized to minimize the memory footprint yet maintaining a fast spectra access. The rMSI GUI is designed for simple and effective data exploration and visualization. Moreover, rMSI is designed to be integrated in the R environment through a library of functions that can be used to share MS data across others R packages. The release of rMSI for R environment establishes a novel and flexible platform for MSI data analysis, completely free and open-source. AVAILABILITY AND IMPLEMENTATION: The code, the documentation, a tutorial and example data are available open-source at: github.com/prafols/rMSI. CONTACT: jesus.brezmes@urv.cat. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Avoiding hard chromatographic segmentation: A moving window approach for the automated resolution of gas chromatography-mass spectrometry-based metabolomics signals by multivariate methods.

Gas chromatography-mass spectrometry (GC-MS) produces large and complex datasets characterized by co-eluted compounds and at trace levels, and with a distinct compound ion-redundancy as a result of the high fragmentation by the electron impact ionization. Compounds in GC-MS can be resolved by taking advantage of the multivariate nature of GC-MS data by applying multivariate resolution methods. However, multivariate methods have to be applied in small regions of the chromatogram, and therefore chromatograms are segmented prior to the application of the algorithms. The automation of this segmentation process is a challenging task as it implies separating between informative data and noise from the chromatogram. This study demonstrates the capabilities of independent component analysis-orthogonal signal deconvolution (ICA-OSD) and multivariate curve resolution-alternating least squares (MCR-ALS) with an overlapping moving window implementation to avoid the typical hard chromatographic segmentation. Also, after being resolved, compounds are aligned across samples by an automated alignment algorithm. We evaluated the proposed methods through a quantitative analysis of GC-qTOF MS data from 25 serum samples. The quantitative performance of both moving window ICA-OSD and MCR-ALS-based implementations was compared with the quantification of 33 compounds by the XCMS package. Results shown that most of the R2 coefficients of determination exhibited a high correlation (R2>0.90) in both ICA-OSD and MCR-ALS moving window-based approaches.

eRah: A Computational Tool Integrating Spectral Deconvolution and Alignment with Quantification and Identification of Metabolites in GC/MS-Based Metabolomics.

Gas chromatography coupled to mass spectrometry (GC/MS) has been a long-standing approach used for identifying small molecules due to the highly reproducible ionization process of electron impact ionization (EI). However, the use of GC-EI MS in untargeted metabolomics produces large and complex data sets characterized by coeluting compounds and extensive fragmentation of molecular ions caused by the hard electron ionization. In order to identify and extract quantitative information on metabolites across multiple biological samples, integrated computational workflows for data processing are needed. Here we introduce eRah, a free computational tool written in the open language R composed of five core functions: (i) noise filtering and baseline removal of GC/MS chromatograms, (ii) an innovative compound deconvolution process using multivariate analysis techniques based on compound match by local covariance (CMLC) and orthogonal signal deconvolution (OSD), (iii) alignment of mass spectra across samples, (iv) missing compound recovery, and (v) identification of metabolites by spectral library matching using publicly available mass spectra. eRah outputs a table with compound names, matching scores and the integrated area of compounds for each sample. The automated capabilities of eRah are demonstrated by the analysis of GC-time-of-flight (TOF) MS data from plasma samples of adolescents with hyperinsulinaemic androgen excess and healthy controls. The quantitative results of eRah are compared to centWave, the peak-picking algorithm implemented in the widely used XCMS package, MetAlign, and ChromaTOF software. Significantly dysregulated metabolites are further validated using pure standards and targeted analysis by GC-triple quadrupole (QqQ) MS, LC-QqQ, and NMR. eRah is freely available at http://CRAN.R-project.org/package=erah .

Automated resolution of chromatographic signals by independent component analysis-orthogonal signal deconvolution in comprehensive gas chromatography/mass spectrometry-based metabolomics.

Comprehensive gas chromatography-mass spectrometry (GC×GC-MS) provides a different perspective in metabolomics profiling of samples. However, algorithms for GC×GC-MS data processing are needed in order to automatically process the data and extract the purest information about the compounds appearing in complex biological samples. This study shows the capability of independent component analysis-orthogonal signal deconvolution (ICA-OSD), an algorithm based on blind source separation and distributed in an R package called osd, to extract the spectra of the compounds appearing in GC×GC-MS chromatograms in an automated manner. We studied the performance of ICA-OSD by the quantification of 38 metabolites through a set of 20 Jurkat cell samples analyzed by GC×GC-MS. The quantification by ICA-OSD was compared with a supervised quantification by selective ions, and most of the R(2) coefficients of determination were in good agreement (R(2)>0.90) while up to 24 cases exhibited an excellent linear relation (R(2)>0.95). We concluded that ICA-OSD can be used to resolve co-eluted compounds in GC×GC-MS.

Performance Comparison of Fuzzy ARTMAP and LDA in Qualitative Classification of Iranian Rosa damascena Essential Oils by an Electronic Nose.

Quality control of essential oils is an important topic in industrial processing of medicinal and aromatic plants. In this paper, the performance of Fuzzy Adaptive Resonant Theory Map (ARTMAP) and linear discriminant analysis (LDA) algorithms are compared in the specific task of quality classification of Rosa damascene essential oil samples (one of the most famous and valuable essential oils in the world) using an electronic nose (EN) system based on seven metal oxide semiconductor (MOS) sensors. First, with the aid of a GC-MS analysis, samples of Rosa damascene essential oils were classified into three different categories (low, middle, and high quality, classes C1, C2, and C3, respectively) based on the total percent of the most crucial qualitative compounds. An ad-hoc electronic nose (EN) system was implemented to sense the samples and acquire signals. Forty-nine features were extracted from the EN sensor matrix (seven parameters to describe each sensor curve response). The extracted features were ordered in relevance by the intra/inter variance criterion (Vr), also known as the Fisher discriminant. A leave-one-out cross validation technique was implemented for estimating the classification accuracy reached by both algorithms. Success rates were calculated using 10, 20, 30, and the entire selected features from the response of the sensor array. The results revealed a maximum classification accuracy of 99% when applying the Fuzzy ARTMAP algorithm and 82% for LDA, using the first 10 features in both cases. Further classification results explained that sub-optimal performance is likely to occur when all the response features are applied. It was found that an electronic nose system employing a Fuzzy ARTMAP classifier could become an accurate, easy, and inexpensive alternative tool for qualitative control in the production of Rosa damascene essential oil.

Compound identification in gas chromatography/mass spectrometry-based metabolomics by blind source separation.

Metabolomics GC-MS samples involve high complexity data that must be effectively resolved to produce chemically meaningful results. Multivariate curve resolution-alternating least squares (MCR-ALS) is the most frequently reported technique for that purpose. More recently, independent component analysis (ICA) has been reported as an alternative to MCR. Those algorithms attempt to infer a model describing the observed data and, therefore, the least squares regression used in MCR assumes that the data is a linear combination of that model. However, due to the high complexity of real data, the construction of a model to describe optimally the observed data is a critical step and these algorithms should prevent the influence from outlier data. This study proves independent component regression (ICR) as an alternative for GC-MS compound identification. Both ICR and MCR though require least squares regression to correctly resolve the mixtures. In this paper, a novel orthogonal signal deconvolution (OSD) approach is introduced, which uses principal component analysis to determine the compound spectra. The study includes a compound identification comparison between the results by ICA-OSD, MCR-OSD, ICR and MCR-ALS using pure standards and human serum samples. Results shows that ICR may be used as an alternative to multivariate curve methods, as ICR efficiency is comparable to MCR-ALS. Also, the study demonstrates that the proposed OSD approach achieves greater spectral resolution accuracy than the traditional least squares approach when compounds elute under undue interference of biological matrices.

Dolphin: a tool for automatic targeted metabolite profiling using 1D and 2D (1)H-NMR data.

One of the main challenges in nuclear magnetic resonance (NMR) metabolomics is to obtain valuable metabolic information from large datasets of raw NMR spectra in a high throughput, automatic, and reproducible way. To date, established software packages used to match and quantify metabolites in NMR spectra remain mostly manually operated, leading to low resolution results and subject to inconsistencies not attributable to the NMR technique itself. Here, we introduce a new software package, called Dolphin, able to automatically quantify a set of target metabolites in multiple sample measurements using an approach based on 1D and 2D NMR techniques to overcome the inherent limitations of 1D (1)H-NMR spectra in metabolomics. Dolphin takes advantage of the 2D J-resolved NMR spectroscopy signal dispersion to avoid inconsistencies in signal position detection, enhancing the reliability and confidence in metabolite matching. Furthermore, in order to improve accuracy in quantification, Dolphin uses 2D NMR spectra to obtain additional information on all neighboring signals surrounding the target metabolite. We have compared the targeted profiling results of Dolphin, recorded from standard biological mixtures, with those of two well established approaches in NMR metabolomics. Overall, Dolphin produced more accurate results with the added advantage of being a fully automated and high throughput processing package.

Metabolomic assessment of the effect of dietary cholesterol in the progressive development of fatty liver disease.

Nonalcoholic fatty liver disease is considered to be the hepatic manifestation of metabolic syndrome and is usually related to high-fat, high-cholesterol diets. With the rationale that the identification and quantification of metabolites in different metabolic pathways may facilitate the discovery of clinically accessible biomarkers, we report the use of (1)H NMR metabolomics for quantitative profiling of liver extracts from LDLr(-/-) mice, a well-documented mouse model of fatty liver disease. A total of 55 metabolites were identified, and multivariate analyses in a diet- and time-comparative strategy were performed. Dietary cholesterol increased the hepatic concentrations of cholesterol, triglycerides, and oleic acid but also decreased the [PUFA/MUFA] ratio as well as the relative amount of long-chain polyunsaturated fatty acids in the liver. This was also accompanied by variations of the hepatic concentration of taurine, glutathione, methionine, and carnitine. Heat-map correlation analyses demonstrated that hepatic inflammation and development of steatosis correlated with cholesterol and triglyceride NMR derived signals, respectively. We conclude that dietary cholesterol is a causal factor in the development of both liver steatosis and hepatic inflammation.

Metabolic phenotyping of genetically modified mice: An NMR metabonomic approach.

Monocyte chemoattractant protein-1 (MCP-1) plays a relevant role in macrophage migration but recent findings suggest an additional role in lipid and glucose metabolism. We report the use of (1)H NMR spectroscopy as a useful complementary method to assess the metabolic function of this gene in a comparative strategy. This metabonomic analysis was rapid, simple, quantitative and reproducible, and revealed a suggestive relationship between the expression of the MCP-1 gene and hepatic glucose and taurine concentrations. This approach should be considered in genetically modified mice when a metabolic alteration is suspected, or in routine assessment of metabolic phenotype.