Improving peak detection in high-resolution LC/MS metabolomics data using preexisting knowledge and machine learning approach

TIHI Toolkit: A Peak Finder and Analyzer for Spectroscopic Data

Complex signal vectors, particularly spectra, are integral to many scientific domains. Interpreting these signals often involves decomposing them into contributions from independent components and subtraction or deconvolution of the channel and instrument noise. Despite the fundamental nature of this task, researchers frequently rely on costly commercial tools. To make such tools accessible to all, we present Tihi, interactive, open-source multiplatform software for interpolation, denoising, baseline correction, peak detection, and signal decomposition. Tihi provides a user-friendly graphical interface (GUI) that facilitates the analysis of spectroscopic data and more. It allows researchers to contribute to and freely distribute these tools, ensuring broad accessibility and fostering collaborative improvements. We present examples demonstrating the efficiency of the program using the spectra of different systems acquired by different spectroscopic techniques, including Raman (aspirin), IR (solid ammonia), XRD (anatase), and UV-vis (petal tip from the Puya alpestris flower). These examples showcase a variety of spectra that differ significantly, from signals with narrow profiles to signals with very broad profiles. This demonstrates the versatility of Tihi for peak identification in a wide range of spectroscopic techniques.

Prenatal and Early Life Exposure to Ambient Air Pollutants Is Associated with the Fecal Metabolome in the First Two Years of Life

Prenatal and early life air pollution exposure has been linked with several adverse health outcomes. However, the mechanisms underlying these relationships are not yet fully understood. Therefore, this study utilizes fecal metabolomics to determine if pre- and postnatal exposure to ambient air pollutants (i.e., PM10, PM2.5, and NO2) is associated with the fecal metabolome in the first 2 years of life in a Latino cohort from Southern California. The aims of this analysis were to estimate associations between (1) prenatal air pollution exposure with fecal metabolic features at 1-month of age, (2) prior month postnatal air pollution exposure with fecal metabolites from 1-month to 2 years of age, and (3) how postnatal air pollution exposure impacts the change over time of fecal metabolites in the first 2 years of life. Prenatal exposure to air pollutants was associated with several Level-1 metabolites, including those involved in vitamin B6 and tyrosine metabolism. Prior month air pollution exposure in the postnatal period was associated with Level-1 metabolites involved in histidine metabolism. Lastly, we found that pre- and postnatal ambient air pollution exposure was associated with changes in metabolic features involved in metabolic pathways including amino acid metabolism, histidine metabolism, and fatty acid metabolism.

Deep learning framework for peak detection at the intact level of therapeutic proteins

While automated peak detection functionalities are available in commercially accessible software, achieving optimal true positive rates frequently necessitates visual inspection and manual adjustments. In the initial phase of this study, hetero-variants (glycoforms) of a monoclonal antibody were distinguished using liquid chromatography-mass spectrometry, revealing discernible peaks at the intact level. To comprehensively identify each peak (hetero-variant) in the intact-level analysis, a deep learning approach utilizing convolutional neural networks (CNNs) was employed in the subsequent phase of the study. In the current case study, utilizing conventional software for peak identification, five peaks were detected using a 0.5 threshold, whereas seven peaks were identified using the CNN model. The model exhibited strong performance with a probability area under the curve (AUC) of 0.9949, surpassing that of partial least squares discriminant analysis (PLS-DA) (probability AUC of 0.8041), and locally weighted regression (LWR) (probability AUC of 0.6885) on the data acquired during experimentation in real-time. The AUC of the receiver operating characteristic curve also illustrated the superior performance of the CNN over PLS-DA and LWR.

Human Milk Composition Is Associated with Maternal Body Mass Index in a Cross-Sectional, Untargeted Metabolomics Analysis of Human Milk from Guatemalan Mothers

Background

Maternal overweight and obesity has been associated with poor lactation performance including delayed lactogenesis and reduced duration. However, the effect on human milk composition is less well understood.

Objectives

We evaluated the relationship of maternal BMI on the human milk metabolome among Guatemalan mothers.

Methods

We used data from 75 Guatemalan mothers who participated in the Household Air Pollution Intervention Network trial. Maternal BMI was measured between 9 and <20 weeks of gestation. Milk samples were collected at a single time point using aseptic collection from one breast at 6 mo postpartum and analyzed using high-resolution mass spectrometry. A cross-sectional untargeted high-resolution metabolomics analysis was performed by coupling hydrophilic interaction liquid chromatography (HILIC) and reverse phase C18 chromatography with mass spectrometry. Metabolic features associated with maternal BMI were determined by a metabolome-wide association study (MWAS), adjusting for baseline maternal age, education, and dietary diversity, and perturbations in metabolic pathways were identified by pathway enrichment analysis.

Results

The mean age of participants at baseline was 23.62 ± 3.81 y, and mean BMI was 24.27 ± 4.22 kg/m2. Of the total metabolic features detected by HILIC column (19,199 features) and by C18 column (11,594 features), BMI was associated with 1026 HILIC and 500 C18 features. Enriched pathways represented amino acid metabolism, galactose metabolism, and xenobiotic metabolic metabolism. However, no significant features were identified after adjusting for multiple comparisons using the Benjamini-Hochberg false discovery rate procedure (FDRBH < 0.2).

Conclusions

Findings from this untargeted MWAS indicate that maternal BMI is associated with metabolic perturbations of galactose metabolism, xenobiotic metabolism, and xenobiotic metabolism by cytochrome p450 and biosynthesis of amino acid pathways. Significant metabolic pathway alterations detected in human milk were associated with energy metabolism-related pathways including carbohydrate and amino acid metabolism.This trial was registered at clinicaltrials.gov as NCT02944682.

Metabolomic Effects of Folic Acid Supplementation in Adults: Evidence from the FACT Trial

BACKGROUND

Folic acid (FA) is the oxidized form of folate found in supplements and FA-fortified foods. Most FA is reduced by dihydrofolate reductase to 5-methyltetrahydrofolate (5mTHF); the latter is the form of folate naturally found in foods. Ingestion of FA increases the plasma levels of both 5mTHF and unmetabolized FA (UMFA). Limited information is available on the downstream metabolic effects of FA supplementation, including potential effects associated with UMFA.

OBJECTIVE

We aimed to assess the metabolic effects of FA-supplementation, and the associations of plasma 5mTHF and UMFA with the metabolome in FA-naïve Bangladeshi adults.

METHODS

Sixty participants were selected from the Folic Acid and Creatine Trial; half received 800 μg FA/day for 12 weeks and half placebo. Plasma metabolome profiles were measured by high-resolution mass spectrometry, including 170 identified metabolites and 26,541 metabolic features. Penalized regression methods were used to assess the associations of targeted metabolites with FA-supplementation, plasma 5mTHF, and plasma UMFA. Pathway analyses were conducted using Mummichog.

RESULTS

In penalized models of identified metabolites, FA-supplementation was associated with higher choline. Changes in 5mTHF concentrations were positively associated with metabolites involved in amino acid metabolism (5-hydroxyindoleacetic acid, acetylmethionine, creatinine, guanidinoacetate, hydroxyproline/n-acetylalanine) and 2 fatty acids (docosahexaenoic acid and linoleic acid). Changes in 5mTHF concentrations were negatively associated with acetylglutamate, acetyllysine, carnitine, propionyl carnitine, cinnamic acid, homogentisate, arachidonic acid, and nicotine. UMFA concentrations were associated with lower levels of arachidonic acid. Together, metabolites selected across all models were related to lipids, aromatic amino acid metabolism, and the urea cycle. Analyses of nontargeted metabolic features identified additional pathways associated with FA supplementation.

CONCLUSION

In addition to the recapitulation of several expected metabolic changes associated with 5mTHF, we observed additional metabolites/pathways associated with FA-supplementation and UMFA. Further studies are needed to confirm these associations and assess their potential implications for human health.

TRIAL REGISTRATION NUMBER

This trial was registered at https://clinicaltrials.gov as NCT01050556.

The potential role of early life feeding patterns in shaping the infant fecal metabolome: implications for neurodevelopmental outcomes

Infant fecal metabolomics can provide valuable insights into the associations of nutrition, dietary patterns, and health outcomes in early life. Breastmilk is typically classified as the best source of nutrition for nearly all infants. However, exclusive breastfeeding may not always be possible for all infants. This study aimed to characterize associations between levels of mixed breastfeeding and formula feeding, along with solid food consumption and the infant fecal metabolome at 1- and 6-months of age. As a secondary aim, we examined how feeding-associated metabolites may be associated with early life neurodevelopmental outcomes. Fecal samples were collected at 1- and 6-months, and metabolic features were assessed via untargeted liquid chromatography/high-resolution mass spectrometry. Feeding groups were defined at 1-month as 1) exclusively breastfed, 2) breastfed >50% of feedings, or 3) formula fed ≥50% of feedings. Six-month groups were defined as majority breastmilk (>50%) or majority formula fed (≥50%) complemented by solid foods. Neurodevelopmental outcomes were assessed using the Bayley Scales of Infant Development at 2 years. Changes in the infant fecal metabolome were associated with feeding patterns at 1- and 6-months. Feeding patterns were associated with the intensities of a total of 57 fecal metabolites at 1-month and 25 metabolites at 6-months, which were either associated with increased breastmilk or increased formula feeding. Most breastmilk-associated metabolites, which are involved in lipid metabolism and cellular processes like cell signaling, were associated with higher neurodevelopmental scores, while formula-associated metabolites were associated with lower neurodevelopmental scores. These findings offer preliminary evidence that feeding patterns are associated with altered infant fecal metabolomes, which may be associated with cognitive development later in life.

Fecal microbiota transplantation promotes reduction of antimicrobial resistance by strain replacement

Multidrug-resistant organism (MDRO) colonization is a fundamental challenge in antimicrobial resistance. Limited studies have shown that fecal microbiota transplantation (FMT) can reduce MDRO colonization, but its mechanisms are poorly understood. We conducted a randomized, controlled trial of FMT for MDRO decolonization in renal transplant recipients called PREMIX (NCT02922816). Eleven participants were enrolled and randomized 1:1 to FMT or an observation period followed by delayed FMT if stool cultures were MDRO positive at day 36. Participants who were MDRO positive after one FMT were treated with a second FMT. At last visit, eight of nine patients who completed all treatments were MDRO culture negative. FMT-treated participants had longer time to recurrent MDRO infection versus PREMIX-eligible controls who were not treated with FMT. Key taxa (Akkermansia muciniphila, Alistipes putredinis, Phocaeicola dorei, Phascolarctobacterium faecium, Alistipes species, Mesosutterella massiliensis, Barnesiella intestinihominis, and Faecalibacterium prausnitzii) from the single feces donor used in the study that engrafted in recipients and metabolites such as short-chain fatty acids and bile acids in FMT-responding participants uncovered leads for rational microbiome therapeutic and diagnostic development. Metagenomic analyses revealed a previously unobserved mechanism of MDRO eradication by conspecific strain competition in an FMT-treated subset. Susceptible Enterobacterales strains that replaced baseline extended-spectrum β-lactamase-producing strains were not detectable in donor microbiota manufactured as FMT doses but in one case were detectable in the recipient before FMT. These data suggest that FMT may provide a path to exploit strain competition to reduce MDRO colonization.

Metabolomic changes associated with chronic arsenic exposure in a Bangladeshi population

Chronic exposure to arsenic (As) remains a global public health concern and our understanding of the biological mechanisms underlying the adverse effects of As exposure remains incomplete. Here, we used a high-resolution metabolomics approach to examine how As affects metabolic pathways in humans. We selected 60 non-smoking adults from the Folic Acid and Creatine Trial (FACT). Inorganic (AsIII, AsV) and organic (monomethylarsonous acid [MMAs], dimethylarsinous Acid [DMAs]) As species were measured in blood and urine collected at baseline and at 12 weeks. Plasma metabolome profiles were measured using untargeted high-resolution mass spectrometry. Associations of blood and urinary As with 170 confirmed metabolites and >26,000 untargeted spectral features were modeled using a metabolome-wide association study (MWAS) approach. Models were adjusted for age, sex, visit, and BMI and corrected for false discovery rate (FDR). In the MWAS screening of confirmed metabolites, 17 were associated with ≥1 blood As species (FDR<0.05), including fatty acids, neurotransmitter metabolites, and amino acids. These results were consistent across blood As species and between blood and urine As. Untargeted MWAS identified 423 spectral features associated with ≥1 blood As species. Unlike the confirmed metabolites, untargeted model results were not consistent across As species, with AsV and DMAs showing distinct association patterns. Mummichog pathway analysis revealed 12 enriched metabolic pathways that overlapped with the 17 identified metabolites, including one carbon metabolism, tricarboxylic acid cycle, fatty acid metabolism, and purine metabolism. Exposure to As may affect numerous essential pathways that underlie the well-characterized associations of As with multiple chronic diseases.

Longitudinal profiles of the fecal metabolome during the first 2 years of life

During the first 2 years of life, the infant gut microbiome is rapidly developing, and gut bacteria may impact host health through the production of metabolites that can have systemic effects. Thus, the fecal metabolome represents a functional readout of gut bacteria. Despite the important role that fecal metabolites may play in infant health, the development of the infant fecal metabolome has not yet been thoroughly characterized using frequent, repeated sampling during the first 2 years of life. Here, we described the development of the fecal metabolome in a cohort of 101 Latino infants with data collected at 1-, 6-, 12-, 18-, and 24-months of age. We showed that the fecal metabolome is highly conserved across time and highly personalized, with metabolic profiles being largely driven by intra-individual variability. Finally, we also identified several novel metabolites and metabolic pathways that changed significantly with infant age, such as valerobetaine and amino acid metabolism, among others.

Use of high-resolution metabolomics to assess the biological perturbations associated with maternal exposure to Bisphenol A and Bisphenol F among pregnant African American women

BACKGROUND

Human and animal exposure to bisphenol A (BPA) has been associated with adverse developmental and reproductive effects. The molecular mechanisms by which BPA exposure exerts its effects are not well-understood, even less known about its analogues bisphenol F (BPF). To address these knowledge gaps, we conducted an untargeted metabolome-wide association study (MWAS) to identify metabolic perturbations associated with BPA/BPF exposures in a pregnant African American cohort.

METHODS

From a subset of study participants enrolled in the Atlanta African American Maternal-Child cohort, we collected both urine samples, for targeted exposure assessment of BPA (N = 230) and BPF (N = 48), and serum samples, for high-resolution metabolomics (HRM) profiling (N = 230), during early pregnancy (8-14 weeks' gestation). Using an established untargeted HRM workflow consisting of MWAS modeling, pathway enrichment analysis, and chemical annotation and confirmation, we investigated the potential metabolic pathways and features associated with BPA/BPF exposures.

RESULTS

The geometric mean creatinine-adjusted concentrations of urinary BPA and BPF were 0.85 ± 2.58 and 0.70 ± 4.71 µg/g creatinine, respectively. After false positive discovery rate correction at 20 % level, 264 and 733 unique metabolic features were significantly associated with urinary BPA and BPF concentrations, representing 10 and 12 metabolic pathways, respectively. Three metabolic pathways, including steroid hormones biosynthesis, lysine and lipoate metabolism, were significantly associated with both BPA and BPF exposure. Using chemical standards, we have confirmed the chemical identity of 16 metabolites significantly associated with BPA or BPF exposure.

CONCLUSIONS

Our findings support that exposure to BPA and BPF in pregnant women is associated with the perturbation of aromatic amino acid metabolism, xenobiotics metabolism, steroid biosynthesis, and other amino acid metabolism closely linked to stress responses, inflammation, neural development, reproduction, and weight regulation.

Clinical metabolomics for inborn errors of metabolism

Metabolism is a highly regulated process that provides nutrients to cells and essential building blocks for the synthesis of protein, DNA and other macromolecules. In healthy biological systems, metabolism maintains a steady state in which the concentrations of metabolites are relatively constant yet are subject to metabolic demands and environmental stimuli. Rare genetic disorders, such as inborn errors of metabolism (IEM), cause defects in regulatory enzymes or proteins leading to metabolic pathway disruption and metabolite accumulation or deficiency. Traditionally, the laboratory diagnosis of IEMs has been limited to analytical methods that target specific metabolites such as amino acids and acyl carnitines. This approach is effective as a screening method for the most common IEM disorders but lacks the comprehensive coverage of metabolites that is necessary to identify rare disorders that present with nonspecific clinical symptoms. Fortunately, advancements in technology and data analytics has introduced a new field of study called metabolomics which has allowed scientists to perform comprehensive metabolite profiling of biological systems to provide insight into mechanism of action and gene function. Since metabolomics seeks to measure all small molecule metabolites in a biological specimen, it provides an innovative approach to evaluating disease in patients with rare genetic disorders. In this review we provide insight into the appropriate application of metabolomics in clinical settings. We discuss the advantages and limitations of the method and provide details related to the technology, data analytics and statistical modeling required for metabolomic profiling of patients with IEMs.

Operationalizing the Exposome Using Passive Silicone Samplers

The exposome, which is defined as the cumulative effect of environmental exposures and corresponding biological responses, aims to provide a comprehensive measure for evaluating non-genetic causes of disease. Operationalization of the exposome for environmental health and precision medicine has been limited by the lack of a universal approach for characterizing complex exposures, particularly as they vary temporally and geographically. To overcome these challenges, passive sampling devices (PSDs) provide a key measurement strategy for deep exposome phenotyping, which aims to provide comprehensive chemical assessment using untargeted high-resolution mass spectrometry for exposome-wide association studies. To highlight the advantages of silicone PSDs, we review their use in population studies and evaluate the broad range of applications and chemical classes characterized using these samplers. We assess key aspects of incorporating PSDs within observational studies, including the need to preclean samplers prior to use to remove impurities that interfere with compound detection, analytical considerations, and cost. We close with strategies on how to incorporate measures of the external exposome using PSDs, and their advantages for reducing variability in exposure measures and providing a more thorough accounting of the exposome. Continued development and application of silicone PSDs will facilitate greater understanding of how environmental exposures drive disease risk, while providing a feasible strategy for incorporating untargeted, high-resolution characterization of the external exposome in human studies.

Per- and polyfluoroalkyl substance (PFAS) exposure, maternal metabolomic perturbation, and fetal growth in African American women: A meet-in-the-middle approach

BACKGROUND

Prenatal exposures to per- and polyfluoroalkyl substances (PFAS) have been linked to reduced fetal growth. However, the detailed molecular mechanisms remain largely unknown. This study aims to investigate biological pathways and intermediate biomarkers underlying the association between serum PFAS and fetal growth using high-resolution metabolomics in a cohort of pregnant African American women in the Atlanta area, Georgia.

METHODS

Serum perfluorohexane sulfonic acid (PFHxS), perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), and perfluorononanoic acid (PFNA) measurements and untargeted serum metabolomics profiling were conducted in 313 pregnant African American women at 8-14 weeks gestation. Multiple linear regression models were applied to assess the associations of PFAS with birth weight and small-for-gestational age (SGA) birth. A high-resolution metabolomics workflow including metabolome-wide association study, pathway enrichment analysis, and chemical annotation and confirmation with a meet-in-the-middle approach was performed to characterize the biological pathways and intermediate biomarkers of the PFAS-fetal growth relationship.

RESULTS

Each log2-unit increase in serum PFNA concentration was significantly associated with higher odds of SGA birth (OR = 1.32, 95% CI 1.07, 1.63); similar but borderline significant associations were found in PFOA (OR = 1.20, 95% CI 0.94, 1.49) with SGA. Among 25,516 metabolic features extracted from the serum samples, we successfully annotated and confirmed 10 overlapping metabolites associated with both PFAS and fetal growth endpoints, including glycine, taurine, uric acid, ferulic acid, 2-hexyl-3-phenyl-2-propenal, unsaturated fatty acid C18:1, androgenic hormone conjugate, parent bile acid, and bile acid-glycine conjugate. Also, we identified 21 overlapping metabolic pathways from pathway enrichment analyses. These overlapping metabolites and pathways were closely related to amino acid, lipid and fatty acid, bile acid, and androgenic hormone metabolism perturbations.

CONCLUSION

In this cohort of pregnant African American women, higher serum concentrations of PFOA and PFNA were associated with reduced fetal growth. Perturbations of biological pathways involved in amino acid, lipid and fatty acid, bile acid, and androgenic hormone metabolism were associated with PFAS exposures and reduced fetal growth, and uric acid was shown to be a potential intermediate biomarker. Our results provide opportunities for future studies to develop early detection and intervention for PFAS-induced fetal growth restriction.

Microbial metabolite delta-valerobetaine is a diet-dependent obesogen

Obesity and obesity-related metabolic disorders are linked to the intestinal microbiome. However, the causality of changes in the microbiome-host interaction affecting energy metabolism remains controversial. Here, we show the microbiome-derived metabolite δ-valerobetaine (VB) is a diet-dependent obesogen that is increased with phenotypic obesity and is correlated with visceral adipose tissue mass in humans. VB is absent in germ-free mice and their mitochondria but present in ex-germ-free conventionalized mice and their mitochondria. Mechanistic studies in vivo and in vitro show VB is produced by diverse bacterial species and inhibits mitochondrial fatty acid oxidation through decreasing cellular carnitine and mitochondrial long-chain acyl-coenzyme As. VB administration to germ-free and conventional mice increases visceral fat mass and exacerbates hepatic steatosis with a western diet but not control diet. Thus, VB provides a molecular target to understand and potentially manage microbiome-host symbiosis or dysbiosis in diet-dependent obesity.

Large scale enzyme based xenobiotic identification for exposomics

Advances in genomics have revealed many of the genetic underpinnings of human disease, but exposomics methods are currently inadequate to obtain a similar level of understanding of environmental contributions to human disease. Exposomics methods are limited by low abundance of xenobiotic metabolites and lack of authentic standards, which precludes identification using solely mass spectrometry-based criteria. Here, we develop and validate a method for enzymatic generation of xenobiotic metabolites for use with high-resolution mass spectrometry (HRMS) for chemical identification. Generated xenobiotic metabolites were used to confirm identities of respective metabolites in mice and human samples based upon accurate mass, retention time and co-occurrence with related xenobiotic metabolites. The results establish a generally applicable enzyme-based identification (EBI) for mass spectrometry identification of xenobiotic metabolites and could complement existing criteria for chemical identification.

Plasma metabolic phenotypes of HPV-associated vs smoking-associated head and neck cancer and patient survival

BACKGROUND

Metabolic differences between human papillomavirus (HPV)-associated head and neck squamous cell carcinoma (HNSCC) and smoking-associated HNSCC may partially explain differences in prognosis. The former relies on mitochondrial oxidative phosphorylation (OXPHOS) while the latter relies on glycolysis. These differences have not been studied in blood.

METHODS

We extracted metabolites using untargeted liquid chromatography high-resolution mass spectrometry from pretreatment plasma in a cohort of 55 HPV-associated and 82 smoking-associated HNSCC subjects. Metabolic pathway enrichment analysis of differentially expressed metabolites produced pathway-based signatures. Significant pathways (P<0.05) were reduced via principal components analysis and assessed with overall survival via Cox models. We classified each subject as glycolytic or OXPHOS phenotype and assessed it with survival.

RESULTS

Of 2,410 analyzed metabolites, 191 were differentially expressed. Relative to smoking-associated HNSCC, bile acid biosynthesis (P<0.0001) and octadecatrienoic acid beta-oxidation (P=0.01), were upregulated in HPV-associated HNSCC, while galactose metabolism (P=0.001) and vitamin B6 metabolism (P=0.01) were downregulated; the first two suggest an OXPHOS phenotype while the latter two suggest glycolytic. First principal components of bile acid biosynthesis (HR=0.52 per standard deviation, 95% CI:0.38-0.72, P<0.001) and octadecatrienoic acid beta-oxidation (HR=0.54 per sd, 95% CI:0.38-0.78, P<0.001) were significantly associated with overall survival independent of HPV and smoking. The glycolytic vs OXPHOS phenotype was also independently associated with survival (HR=3.17, 95% CI:1.07-9.35; P=0.04).

CONCLUSIONS

Plasma metabolites related to glycolysis and mitochondrial OXPHOS may be biomarkers of HNSCC patient prognosis independent of HPV or smoking. Future investigations should determine if they predict treatment efficacy.

IMPACT

Blood metabolomics may be a useful marker to aid HNSCC patient prognosis.

High sensitivity and specificity feature detection in liquid chromatography-mass spectrometry data: A deep learning framework

Feature detection is a crucial pre-processing step for high-resolution liquid chromatography-mass spectrometry (LC-MS) data analysis. Typical practices based on thresholds or rigid mathematical assumptions can cause ineffective performance in detecting low abundance and non-ideal distributed compounds. We herein introduce a novel feature detection method based on deep learning named SeA-M2Net that considers feature detection as an image-based object detection task. By fully employing raw data directly, and integrating all related factors (e.g., LC elution, charge state, and isotope distribution) with two-dimensional pseudo color images to calculate the probability of the presence of the compound, low abundance compounds can be well preserved and observed. More importantly, SeA-M2Net, with deep multilevel and multiscale structures focuses on compound pattern detection in a learned method instead of assuming a mathematical parametric model. All parameters in SeA-M2Net are learned from data in the training procedure, thus allowing for maximum flexibility of pattern distribution deformation. The algorithm is tested on several LC-MS datasets of multiple biological samples obtained from different instruments with varied experimental settings. We demonstrate the superiority of the new approach in handling complex compound patterns (e.g., low abundance, overlapping regions, LC shifts, and missing values). Our experiments indicate that SeA-M2Net outperforms widely used detection methods in terms of detection accuracy.

CPVA: a web-based metabolomic tool for chromatographic peak visualization and annotation

MOTIVATION

Liquid chromatography-mass spectrometry-based non-targeted metabolomics is routinely performed to qualitatively and quantitatively analyze a tremendous amount of metabolite signals in complex biological samples. However, false-positive peaks in the datasets are commonly detected as metabolite signals by using many popular software, resulting in non-reliable measurement.

RESULTS

To reduce false-positive calling, we developed an interactive web tool, termed CPVA, for visualization and accurate annotation of the detected peaks in non-targeted metabolomics data. We used a chromatogram-centric strategy to unfold the characteristics of chromatographic peaks through visualization of peak morphology metrics, with additional functions to annotate adducts, isotopes and contaminants. CPVA is a free, user-friendly tool to help users to identify peak background noises and contaminants, resulting in decrease of false-positive or redundant peak calling, thereby improving the data quality of non-targeted metabolomics studies.

AVAILABILITY AND IMPLEMENTATION

The CPVA is freely available at http://cpva.eastus.cloudapp.azure.com. Source code and installation instructions are available on GitHub: https://github.com/13479776/cpva.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Use of Integrative Interactomics for Improvement of Farm Animal Health and Welfare: An Example with Fescue Toxicosis

Rapid scientific advances are increasing our understanding of the way complex biological interactions integrate to maintain homeostatic balance and how seemingly small, localized perturbations can lead to systemic effects. The ‘omics movement, alongside increased throughput resulting from statistical and computational advances, has transformed our understanding of disease mechanisms and the multi-dimensional interaction between environmental stressors and host physiology through data integration into multi-dimensional analyses, i.e., integrative interactomics. This review focuses on the use of high-throughput technologies in farm animal research, including health- and toxicology-related papers. Although limited, we highlight recent animal agriculture-centered reports from the integrative multi-‘omics movement. We provide an example with fescue toxicosis, an economically costly disease affecting grazing livestock, and describe how integrative interactomics can be applied to a disease with a complex pathophysiology in the pursuit of novel treatment and management approaches. We outline how ‘omics techniques have been used thus far to understand fescue toxicosis pathophysiology, lay out a framework for the fescue toxicosis integrome, identify some challenges we foresee, and offer possible means for addressing these challenges. Finally, we briefly discuss how the example with fescue toxicosis could be used for other agriculturally important animal health and welfare problems.

Food Phenotyping: Recording and Processing of Non-Targeted Liquid Chromatography Mass Spectrometry Data for Verifying Food Authenticity

Experiments based on metabolomics represent powerful approaches to the experimental verification of the integrity of food. In particular, high-resolution non-targeted analyses, which are carried out by means of liquid chromatography-mass spectrometry systems (LC-MS), offer a variety of options. However, an enormous amount of data is recorded, which must be processed in a correspondingly complex manner. The evaluation of LC-MS based non-targeted data is not entirely trivial and a wide variety of strategies have been developed that can be used in this regard. In this paper, an overview of the mandatory steps regarding data acquisition is given first, followed by a presentation of the required preprocessing steps for data evaluation. Then some multivariate analysis methods are discussed, which have proven to be particularly suitable in this context in recent years. The publication closes with information on the identification of marker compounds.

GMSimpute: a generalized two-step Lasso approach to impute missing values in label-free mass spectrum analysis

Motivation

Missingness in label-free mass spectrometry is inherent to the technology. A computational approach to recover missing values in metabolomics and proteomics datasets is important. Most existing methods are designed under a particular assumption, either missing at random or under the detection limit. If the missing pattern deviates from the assumption, it may lead to biased results. Hence, we investigate the missing patterns in free mass spectrometry data and develop an omnibus approach GMSimpute, to allow effective imputation accommodating different missing patterns.

Results

Three proteomics datasets and one metabolomics dataset indicate missing values could be a mixture of abundance-dependent and abundance-independent missingness. We assess the performance of GMSimpute using simulated data (with a wide range of 80 missing patterns) and metabolomics data from the Cancer Genome Atlas breast cancer and clear cell renal cell carcinoma studies. Using Pearson correlation and normalized root mean square errors between the true and imputed abundance, we compare its performance to K-nearest neighbors’ type approaches, Random Forest, GSimp, a model-based method implemented in DanteR and minimum values. The results indicate GMSimpute provides higher accuracy in imputation and exhibits stable performance across different missing patterns. In addition, GMSimpute is able to identify the features in downstream differential expression analysis with high accuracy when applied to the Cancer Genome Atlas datasets.

Availability and implementation

GMSimpute is on CRAN: https://cran.r-project.org/web/packages/GMSimpute/index.html.

Supplementary information

Supplementary data are available at Bioinformatics online.

Reference Standardization for Quantification and Harmonization of Large-Scale Metabolomics

Reference standardization was developed to address quantification and harmonization challenges for high-resolution metabolomics (HRM) data collected across different studies or analytical methods. Reference standardization relies on the concurrent analysis of calibrated pooled reference samples at predefined intervals and enables a single-step batch correction and quantification for high-throughput metabolomics. Here, we provide quantitative measures of approximately 200 metabolites for each of three pooled reference materials (220 metabolites for Qstd3, 211 metabolites for CHEAR, 204 metabolites for NIST1950) and show application of this approach for quantification supports harmonization of metabolomics data collected from 3677 human samples in 17 separate studies analyzed by two complementary HRM methods over a 17-month period. The results establish reference standardization as a method suitable for harmonizing large-scale metabolomics data and extending capabilities to quantify large numbers of known and unidentified metabolites detected by high-resolution mass spectrometry methods.

Cardiovascular Risk and Resilience Among Black Adults: Rationale and Design of the MECA Study

Background

Cardiovascular disease incidence, prevalence, morbidity, and mortality have declined in the past several decades; however, disparities persist among subsets of the population. Notably, blacks have not experienced the same improvements on the whole as whites. Furthermore, frequent reports of relatively poorer health statistics among the black population have led to a broad assumption that black race reliably predicts relatively poorer health outcomes. However, substantial intraethnic and intraracial heterogeneity exists; moreover, individuals with similar risk factors and environmental exposures are often known to experience vastly different cardiovascular health outcomes. Thus, some individuals have good outcomes even in the presence of cardiovascular risk factors, a concept known as resilience.

Methods and Results

The MECA (Morehouse‐Emory Center for Health Equity) Study was designed to investigate the multilevel exposures that contribute to “resilience” in the face of risk for poor cardiovascular health among blacks in the greater Atlanta, GA, metropolitan area. We used census tract data to determine “at‐risk” and “resilient” neighborhoods with high or low prevalence of cardiovascular morbidity and mortality, based on cardiovascular death, hospitalization, and emergency department visits for blacks. More than 1400 individuals from these census tracts assented to demographic, health, and psychosocial questionnaires administered through telephone surveys. Afterwards, ≈500 individuals were recruited to enroll in a clinical study, where risk biomarkers, such as oxidative stress, and inflammatory markers, endothelial progenitor cells, metabolomic and microRNA profiles, and subclinical vascular dysfunction were measured. In addition, comprehensive behavioral questionnaires were collected and ideal cardiovascular health metrics were assessed using the American Heart Association's Life Simple 7 measure. Last, 150 individuals with low Life Simple 7 were recruited and randomized to a behavioral mobile health (eHealth) plus health coach or eHealth only intervention and followed up for improvement.

Conclusions

The MECA Study is investigating socioenvironmental and individual behavioral measures that promote resilience to cardiovascular disease in blacks by assessing biological, functional, and molecular mechanisms.

REGISTRATION

URL: https://www.clini​caltr​ials.gov. Unique identifier: NCT03308812.

An automated system for predicting detection limit and precision profile from a chromatogram

This paper presents a basic model of an automated system for predicting the detection limit and precision profile (plot of relative standard deviation (RSD) of measurements against concentration) in chromatography. The fundamental assumption is that the major source of response errors at low sample concentrations is background noise and at high concentrations, it is the volumes injected into an HPLC system by a sample injector. The noise is approximated by the mixed random processes of the first order autoregressive process AR(1) and white noise. The research procedures are: (1) the description of the standard deviation (SD) of measurements in terms of the parameters of the mixed random processes; (2) the algorithm for the parameter estimation of the mixed processes from actual background noise; (3) the mathematical distinction between noise and signal in a chromatogram. When compounds are chromatographically separated, each obtained signal is given the detection limit and precision profile on laboratory-made software. A file of a chromatogram is the only requirement for the theoretical prediction of measurement uncertainty and therefore the repeated measurements of real samples can be dispensed with. The theoretically predicted RSDs are verified by comparing them with the statistical RSDs obtained by repeated measurements. Signal shapes on noise are illustrated at the detection limit and quantitation limit, the signal-to-noise ratios of which are close to the widely adopted values, 3 and 10, respectively.

Ammonium Fluoride as Suitable Additive for HILIC-Based LC-HRMS Metabolomics

Hydrophilic Interaction Liquid Chromatography (HILIC) chromatography is widely applied in metabolomics as a complementary strategy to reverse phase chromatography. Nevertheless, it still faces several issues in terms of peak shape and compounds ionization, limiting the automatic de-convolution and data semi-quantification performed through dedicated software. A way to improve the chromatographic and ionization performance of a HILIC method is to modify the electrostatic interactions of the analytes with both mobile and stationary phases. In this study, using a ZIC-HILIC chromatographic phase, we evaluated the performance of ammonium fluoride (AF) as additive salt, comparing its performance to ammonium acetate (AA). Three comparative criteria were selected: (1) identification and peak quality of 34 standards following a metabolomics-specific evaluation approach, (2) an intraday repeatability test with real samples and (3) performing two real metabolomics fingerprints with the AF method to evaluate its inter-day repeatability. The AF method showed not only higher ionization efficiency and signal-to-noise ratio but also better repeatability and robustness than the AA approach. A tips and tricks section is then added, aiming at improving method replicability for further users. In conclusion, ammonium fluoride as additive salt presents several advantages and might be considered as a step forward in the application of robust HILIC methods in metabolomics.

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.

Metabolomics in the clinic: A review of the shared and unique features of untargeted metabolomics for clinical research and clinical testing

Metabolomics is the untargeted measurement of the metabolome, which is composed of the complement of small molecules detected in a biological sample. As such, metabolomic analysis produces a global biochemical phenotype. It is a technology that has been utilized in the research setting for over a decade. The metabolome is directly linked to and is influenced by genetics, epigenetics, environmental factors, and the microbiome-all of which affect health. Metabolomics can be applied to human clinical diagnostics and to other fields such as veterinary medicine, nutrition, exercise, physiology, agriculture/plant biochemistry, and toxicology. Applications of metabolomics in clinical testing are emerging, but several aspects of its use as a clinical test differ from applications focused on research or biomarker discovery and need to be considered for metabolomics clinical test data to have optimum impact, be meaningful, and be used responsibly. In this review, we deconstruct aspects and challenges of metabolomics for clinical testing by illustrating the significance of test design, accurate and precise data acquisition, quality control, data processing, n-of-1 comparison to a reference population, and biochemical pathway analysis. We describe how metabolomics technology is integral to defining individual biochemical phenotypes, elaborates on human health and disease, and fits within the precision medicine landscape. Finally, we conclude by outlining some future steps needed to bring metabolomics into the clinical space and to be recognized by the broader medical and regulatory fields.

Trace Phosphate Improves ZIC-pHILIC Peak Shape, Sensitivity, and Coverage for Untargeted Metabolomics

Existing hydrophilic interaction liquid chromatography (HILIC) methods, considered individually, each exhibit poor chromatographic performance for a substantial fraction of polar metabolites. In addition to limiting metabolome coverage, such deficiencies also complicate automated data processing. Here we show that some of these analytical challenges can be addressed for the ZIC-pHILIC, a zwitterionic stationary phase commonly used in metabolomics, with the addition of trace levels of phosphate. Specifically, micromolar phosphate extended metabolome coverage by hundreds of credentialed features, improved peak shapes, and reduced peak-detection errors during informatic processing. Although the addition of high levels of phosphate (millimolar) as a HILIC mobile phase buffer has been explored previously, such concentrations interfere with mass spectrometric (MS) detection. We show that using phosphate as a trace additive at micromolar concentrations improves analysis by electrospray MS, increasing signal for a diverse set of polar standards. Given the small amount of phosphate needed, comparable chromatographic improvements were also achieved by direct addition of phosphate to the sample during reconstitution. Our results suggest that defects in ZIC-pHILIC performance are predominantly driven by electrostatic interactions, which can be modulated by phosphate. These findings constitute both a methodological improvement for untargeted metabolomics and an advance in our understanding of the mechanisms limiting HILIC coverage.

Particulate metal exposures induce plasma metabolome changes in a commuter panel study

INTRODUCTION

Advances in liquid chromatography-mass spectrometry (LC-MS) have enabled high-resolution metabolomics (HRM) to emerge as a sensitive tool for measuring environmental exposures and corresponding biological response. Using measurements collected as part of a large, panel-based study of car commuters, the current analysis examines in-vehicle air pollution concentrations, targeted inflammatory biomarker levels, and metabolomic profiles to trace potential metabolic perturbations associated with on-road traffic exposures.

METHODS

A 60-person panel of adults participated in a crossover study, where each participant conducted a highway commute and randomized to either a side-street commute or clinic exposure session. In addition to in-vehicle exposure characterizations, participants contributed pre- and post-exposure dried blood spots for 2-hr changes in targeted proinflammatory and vascular injury biomarkers and 10-hr changes in the plasma metabolome. Samples were analyzed on a Thermo QExactive MS system in positive and negative electrospray ionization (ESI) mode. Data were processed and analyzed in R using apLCMS, xMSanalyzer, and limma. Features associated with environmental exposures or biological endpoints were identified with a linear mixed effects model and annotated through human metabolic pathway analysis in mummichog.

RESULTS

HRM detected 10-hr perturbations in 110 features associated with in-vehicle, particulate metal exposures (Al, Pb, and Fe) which reflect changes in arachidonic acid, leukotriene, and tryptophan metabolism. Two-hour changes in proinflammatory biomarkers hs-CRP, IL-6, IL-8, and IL-1β were also associated with 10-hr changes in the plasma metabolome, suggesting diverse amino acid, leukotriene, and antioxidant metabolism effects. A putatively identified metabolite, 20-OH-LTB4, decreased after in-vehicle exposure to particulate metals, suggesting a subclinical immune response.

CONCLUSIONS

Acute exposures to traffic-related air pollutants are associated with broad inflammatory response, including several traditional markers of inflammation.

High-Resolution Metabolomics Assessment of Military Personnel: Evaluating Analytical Strategies for Chemical Detection

OBJECTIVE

The aim of this study was to maximize detection of serum metabolites with high-resolution metabolomics (HRM).

METHODS

Department of Defense Serum Repository (DoDSR) samples were analyzed using ultrahigh resolution mass spectrometry with three complementary chromatographic phases and four ionization modes. Chemical coverage was evaluated by number of ions detected and accurate mass matches to a human metabolomics database.

RESULTS

Individual HRM platforms provided accurate mass matches for up to 58% of the KEGG metabolite database. Combining two analytical methods increased matches to 72% and included metabolites in most major human metabolic pathways and chemical classes. Detection and feature quality varied by analytical configuration.

CONCLUSIONS

Dual chromatography HRM with positive and negative electrospray ionization provides an effective generalized method for metabolic assessment of military personnel.

One Step Forward for Reducing False Positive and False Negative Compound Identifications from Mass Spectrometry Metabolomics Data: New Algorithms for Constructing Extracted Ion Chromatograms and Detecting Chromatographic Peaks

False positive and false negative peaks detected from extracted ion chromatograms (EIC) are an urgent problem with existing software packages that preprocess untargeted liquid or gas chromatography-mass spectrometry metabolomics data because they can translate downstream into spurious or missing compound identifications. We have developed new algorithms that carry out the sequential construction of EICs and detection of EIC peaks. We compare the new algorithms to two popular software packages XCMS and MZmine 2 and present evidence that these new algorithms detect significantly fewer false positives. Regarding the detection of compounds known to be present in the data, the new algorithms perform at least as well as XCMS and MZmine 2. Furthermore, we present evidence that mass tolerance in m/z should be favored rather than mass tolerance in ppm in the process of constructing EICs. The mass tolerance parameter plays a critical role in the EIC construction process and can have immense impact on the detection of EIC peaks.

Computational Metabolomics: A Framework for the Million Metabolome

"Sola dosis facit venenum." These words of Paracelsus, "the dose makes the poison", can lead to a cavalier attitude concerning potential toxicities of the vast array of low abundance environmental chemicals to which humans are exposed. Exposome research teaches that 80-85% of human disease is linked to environmental exposures. The human exposome is estimated to include >400,000 environmental chemicals, most of which are uncharacterized with regard to human health. In fact, mass spectrometry measures >200,000 m/z features (ions) in microliter volumes derived from human samples; most are unidentified. This crystallizes a grand challenge for chemical research in toxicology: to develop reliable and affordable analytical methods to understand health impacts of the extensive human chemical experience. To this end, there appears to be no choice but to abandon the limitations of measuring one chemical at a time. The present review looks at progress in computational metabolomics to provide probability-based annotation linking ions to known chemicals and serve as a foundation for unambiguous designation of unidentified ions for toxicologic study. We review methods to characterize ions in terms of accurate mass m/z, chromatographic retention time, correlation of adduct, isotopic and fragment forms, association with metabolic pathways and measurement of collision-induced dissociation products, collision cross section, and chirality. Such information can support a largely unambiguous system for documenting unidentified ions in environmental surveillance and human biomonitoring. Assembly of this data would provide a resource to characterize and understand health risks of the array of low-abundance chemicals to which humans are exposed.

High-resolution metabolomics of occupational exposure to trichloroethylene

Background: Occupational exposure to trichloroethylene (TCE) has been linked to adverse health outcomes including non-Hodgkin’s lymphoma and kidney and liver cancer; however, TCE’s mode of action for development of these diseases in humans is not well understood.

Methods: Non-targeted metabolomics analysis of plasma obtained from 80 TCE-exposed workers [full shift exposure range of 0.4 to 230 parts-per-million of air (ppm_a)] and 95 matched controls were completed by ultra-high resolution mass spectrometry. Biological response to TCE exposure was determined using a metabolome-wide association study (MWAS) framework, with metabolic changes and plasma TCE metabolites evaluated by dose-response and pathway enrichment. Biological perturbations were then linked to immunological, renal and exposure molecular markers measured in the same population.

Results: Metabolic features associated with TCE exposure included known TCE metabolites, unidentifiable chlorinated compounds and endogenous metabolites. Exposure resulted in a systemic response in endogenous metabolism, including disruption in purine catabolism and decreases in sulphur amino acid and bile acid biosynthesis pathways. Metabolite associations with TCE exposure included uric acid (β = 0.13, P-value = 3.6 × 10⁻⁵), glutamine (β = 0.08, P-value = 0.0013), cystine (β = 0.75, P-value = 0.0022), methylthioadenosine (β = −1.6, P-value = 0.0043), taurine (β = −2.4, P-value = 0.0011) and chenodeoxycholic acid (β = −1.3, P-value = 0.0039), which are consistent with known toxic effects of TCE, including immunosuppression, hepatotoxicity and nephrotoxicity. Correlation with additional exposure markers and physiological endpoints supported known disease associations.

Conclusions: High-resolution metabolomics correlates measured occupational exposure to internal dose and metabolic response, providing insight into molecular mechanisms of exposure-related disease aetiology.

General strategies to increase the repeatability in non-target screening by liquid chromatography-high resolution mass spectrometry

This article focuses on the data evaluation of non-target high-resolution LC-MS profiles of water samples. Taking into account multiple technical replicates, the difficulties in peak recognition and the related problems of false positive and false negative findings are systematically demonstrated. On the basis of a combinatorial approach, different models involving sophisticated workflows are evaluated, particularly with regard to the repeatability. In addition, the improvement resulting from data processing was systematically taken into consideration where the recovery of spiked standards emphasized that real peaks of interest were barely or not removed by the derived filter criteria. The comprehensive evaluation included different matrix types spiked with up to 263 analytical standards which were analyzed repeatedly leading to a total number of more than 250 injections that were incorporated in the assessment of different models of data processing. It was found that the analysis of multiple replicates is the key factor as, on the one hand, it provides the option of integrating valuable filters in order to minimize the false positive rate and, on the other hand, allows correcting partially false negative findings occurring during the peak recognition. The developed processing strategies including replicates clearly point to an enhanced data quality since both the repeatability as well as the peak recognition could be considerably improved. As proof of concept, four different matrix types, including a wastewater treatment plant (WWTP) effluent, were spiked with 130 isotopically labeled standards at different concentration levels. Despite the stringent filter criteria, at 100 ng L(-1) recovery rates of up to 93% were reached in the positive ionization mode. The proposed model, comprising three technical replicates, filters less than 5% and 2% of the standards recognized at 100 and 500 ng L(-1), respectively and thus indicates the general applicability of the presented strategies.

Local false discovery rate estimation using feature reliability in LC/MS metabolomics data

False discovery rate (FDR) control is an important tool of statistical inference in feature selection. In mass spectrometry-based metabolomics data, features can be measured at different levels of reliability and false features are often detected in untargeted metabolite profiling as chemical and/or bioinformatics noise. The traditional false discovery rate methods treat all features equally, which can cause substantial loss of statistical power to detect differentially expressed features. We propose a reliability index for mass spectrometry-based metabolomics data with repeated measurements, which is quantified using a composite measure. We then present a new method to estimate the local false discovery rate (lfdr) that incorporates feature reliability. In simulations, our proposed method achieved better balance between sensitivity and controlling false discovery, as compared to traditional lfdr estimation. We applied our method to a real metabolomics dataset and were able to detect more differentially expressed metabolites that were biologically meaningful.