Denoising Search doubles the number of metabolite and exposome annotations in human plasma using an Orbitrap Astral mass spectrometer.

Chemical exposures may impact human metabolism and contribute to the etiology of neurodegenerative disorders like Alzheimer's Disease (AD). Identifying these small metabolites involves matching experimental spectra to reference spectra in databases. However, environmental chemicals or physiologically active metabolites are usually present at low concentrations in human specimens. The presence of noise ions can significantly degrade spectral quality, leading to false negatives and reduced identification rates. In response to this challenge, the Spectral Denoising algorithm removes both chemical and electronic noise. Spectral Denoising outperformed alternative methods in benchmarking studies on 240 tested metabolites. It improved high confident compound identifications at an average 35-fold lower concentrations than previously achievable. Spectral Denoising proved highly robust against varying levels of both chemical and electronic noise even with >150-fold higher intensity of noise ions than true fragment ions. For human plasma samples of AD patients that were analyzed on the Orbitrap Astral mass spectrometer, Denoising Search detected 2.3-fold more annotated compounds compared to the Exploris 240 Orbitrap instrument, including drug metabolites, household and industrial chemicals, and pesticides. This combination of advanced instrumentation with a superior denoising algorithm opens the door for precision medicine in exposome research.

An exploratory study on lipidomic profiles in a cohort of individuals with posttraumatic stress disorder.

Posttraumatic stress disorder (PTSD) can develop after trauma exposure. Some studies report that women develop PTSD at twice the rate of men, despite greater trauma exposure in men. Lipids and their metabolites (lipidome) regulate a myriad of key biological processes and pathways such as membrane integrity, oxidative stress, and neuroinflammation in the brain by maintaining neuronal connectivity and homeostasis. In this study, we analyzed the lipidome of 40 adults with PTSD and 40 trauma-exposed non-PTSD individuals (n = 20/sex/condition; 19-39 years old). Plasma samples were analyzed for lipidomics using Quadrupole Time-of-Flight (QToF) mass spectrometry. Additionally, ~ 90 measures were collected, on sleep, and mental and physical health indices. Poorer sleep quality was associated with greater PTSD severity in both sexes. The lipidomics analysis identified a total of 348 quantifiable known lipid metabolites and 1951 lipid metabolites that are yet unknown; known metabolites were part of 13 lipid subclasses. After adjusting for BMI and sleep quality, in women with PTSD, only one lipid subclass, phosphatidylethanolamine (PE) was altered, whereas, in men with PTSD, 9 out of 13 subclasses were altered compared to non-PTSD women and men, respectively. Severe PTSD was associated with 22% and 5% of altered lipid metabolites in men and women, respectively. Of the changed metabolites, only 0.5% measures (2 PEs and cholesterol) were common between women and men with PTSD. Several sphingomyelins, PEs, ceramides, and triglycerides were increased in men with severe PTSD. The correlations between triglycerides and ceramide metabolites with cholesterol metabolites and systolic blood pressure were dependent upon sex and PTSD status. Alterations in triglycerides and ceramides are linked with cardiac health and metabolic function in humans. Thus, disturbed sleep and higher body mass may have contributed to changes in the lipidome found in PTSD.

Resolving multi-image spatial lipidomic responses to inhaled toxicants by machine learning.

Regional responses to inhaled toxicants are essential to understand the pathogenesis of lung disease under exposure to air pollution. We evaluated the effect of combined allergen sensitization and ozone exposure on eliciting spatial differences in lipid distribution in the mouse lung that may contribute to ozone-induced exacerbations in asthma. Lung lobes from male and female BALB/c mice were cryosectioned and acquired by high resolution mass spectrometry imaging (MSI). Processed MSI peak annotations were validated by LC-MS/MS data from scraped tissue slides and microdissected lung tissue. Images were normalized and segmented into clusters. Interestingly, segmented clusters overlapped with stained serial tissue sections, enabling statistical analysis across biological replicates for morphologically relevant lung regions. Spatially distinct lipids had higher overall degree of unsaturated fatty acids in distal lung regions compared to proximal regions. Furthermore, the airway and alveolar epithelium exhibited significantly decreased sphingolipid and glycerophospholipid abundance in females, but not in males. We demonstrate the potential role of lipid saturation in healthy lung function and highlight sex differences in regional lung lipid distribution following ozone exposure. Our study provides a framework for future MSI experiments capable of relative quantification across biological replicates and expansion to multiple sample types, including human tissue.

High-Resolution Mass Spectrometry for Human Exposomics: Expanding Chemical Space Coverage.

In the modern "omics" era, measurement of the human exposome is a critical missing link between genetic drivers and disease outcomes. High-resolution mass spectrometry (HRMS), routinely used in proteomics and metabolomics, has emerged as a leading technology to broadly profile chemical exposure agents and related biomolecules for accurate mass measurement, high sensitivity, rapid data acquisition, and increased resolution of chemical space. Non-targeted approaches are increasingly accessible, supporting a shift from conventional hypothesis-driven, quantitation-centric targeted analyses toward data-driven, hypothesis-generating chemical exposome-wide profiling. However, HRMS-based exposomics encounters unique challenges. New analytical and computational infrastructures are needed to expand the analysis coverage through streamlined, scalable, and harmonized workflows and data pipelines that permit longitudinal chemical exposome tracking, retrospective validation, and multi-omics integration for meaningful health-oriented inferences. In this article, we survey the literature on state-of-the-art HRMS-based technologies, review current analytical workflows and informatic pipelines, and provide an up-to-date reference on exposomic approaches for chemists, toxicologists, epidemiologists, care providers, and stakeholders in health sciences and medicine. We propose efforts to benchmark fit-for-purpose platforms for expanding coverage of chemical space, including gas/liquid chromatography-HRMS (GC-HRMS and LC-HRMS), and discuss opportunities, challenges, and strategies to advance the burgeoning field of the exposome.

Longitudinal changes in DNA methylation during the onset of islet autoimmunity differentiate between reversion versus progression of islet autoimmunity.

Background: Type 1 diabetes (T1D) is preceded by a heterogenous pre-clinical phase, islet autoimmunity (IA). We aimed to identify pre vs. post-IA seroconversion (SV) changes in DNAm that differed across three IA progression phenotypes, those who lose autoantibodies (reverters), progress to clinical T1D (progressors), or maintain autoantibody levels (maintainers). Methods: This epigenome-wide association study (EWAS) included longitudinal DNAm measurements in blood (Illumina 450K and EPIC) from participants in Diabetes Autoimmunity Study in the Young (DAISY) who developed IA, one or more islet autoantibodies on at least two consecutive visits. We compared reverters - individuals who sero-reverted, negative for all autoantibodies on at least two consecutive visits and did not develop T1D (n=41); maintainers - continued to test positive for autoantibodies but did not develop T1D (n=60); progressors - developed clinical T1D (n=42). DNAm data were measured before (pre-SV visit) and after IA (post-SV visit). Linear mixed models were used to test for differences in pre- vs post-SV changes in DNAm across the three groups. Linear mixed models were also used to test for group differences in average DNAm. Cell proportions, age, and sex were adjusted for in all models. Median follow-up across all participants was 15.5 yrs. (interquartile range (IQR): 10.8-18.7). Results: The median age at the pre-SV visit was 2.2 yrs. (IQR: 0.8-5.3) in progressors, compared to 6.0 yrs. (IQR: 1.3-8.4) in reverters, and 5.7 yrs. (IQR: 1.4-9.7) in maintainers. Median time between the visits was similar in reverters 1.4 yrs. (IQR: 1-1.9), maintainers 1.3 yrs. (IQR: 1.0-2.0), and progressors 1.8 yrs. (IQR: 1.0-2.0). Changes in DNAm, pre- vs post-SV, differed across the groups at one site (cg16066195) and 11 regions. Average DNAm (mean of pre- and post-SV) differed across 22 regions. Conclusion: Differentially changing DNAm regions were located in genomic areas related to beta cell function, immune cell differentiation, and immune cell function.

Metabolic Responses to an Acute Glucose Challenge: The Differential Effects of Eight Weeks of Almond vs. Cracker Consumption in Young Adults.

This study investigated the dynamic responses to an acute glucose challenge following chronic almond versus cracker consumption for 8 weeks (clinicaltrials.gov ID: NCT03084003). Seventy-three young adults (age: 18-19 years, BMI: 18-41 kg/m2) participated in an 8-week randomized, controlled, parallel-arm intervention and were randomly assigned to consume either almonds (2 oz/d, n=38) or an isocaloric control snack of graham crackers (325 kcal/d, n=35) daily for 8 weeks. Twenty participants from each group underwent a 2-hour oral glucose tolerance test (oGTT) at the end of the 8-week intervention. Metabolite abundances in the oGTT serum samples were quantified using untargeted metabolomics, and targeted analyses for free PUFAs, total fatty acids, oxylipins, and endocannabinoids. Multivariate, univariate, and chemical enrichment analyses were conducted to identify significant metabolic shifts. Findings exhibit a biphasic lipid response distinguished by higher levels of unsaturated triglycerides in the earlier periods of the oGTT followed by lower levels in the latter period in the almond versus cracker group (p-value<0.05, chemical enrichment analyses). Almond (vs. cracker) consumption was also associated with higher AUC120 min of aminomalonate, and oxylipins (p-value<0.05), but lower AUC120 min of L-cystine, N-acetylmannosamine, and isoheptadecanoic acid (p-value<0.05). Additionally, the Matsuda Index in the almond group correlated with AUC120 min of CE 22:6 (r=-0.46; p-value<0.05) and 12,13 DiHOME (r=0.45; p-value<0.05). Almond consumption for 8 weeks leads to dynamic, differential shifts in response to an acute glucose challenge, marked by alterations in lipid and amino acid mediators involved in metabolic and physiological pathways.

Diagnosing Parkinson's disease and monitoring its progression: Biomarkers from combined GC-TOF MS and LC-MS/MS untargeted metabolomics.

Parkinson's disease (PD) is a prevalent neurodegenerative disorder with a poorly understood etiology. An accurate diagnosis of idiopathic PD remains challenging as misdiagnosis is common in routine clinical practice. Moreover, current therapeutics focus on symptomatic management rather than curing or slowing down disease progression. Therefore, identification of potential PD biomarkers and providing a better understanding of the underlying disease pathophysiology are urgent. Herein, hydrophilic interaction liquid chromatography-mass spectrometry (LC-MS/MS) and gas chromatography-mass spectrometry (GC-TOF MS) based metabolomics approaches were used to profile the serum metabolome of 50 patients with different stages of idiopathic PD (early, mid and advanced) and 45 age-matched controls. Levels of 57 metabolites including cysteine-S-sulfate and N-acetyl tryptophan were significantly higher in patients with PD compared to controls, with lower amounts of additional 51 metabolites including vanillic acid, and N-acetylaspartic acid. Xanthines, including caffeine and its downstream metabolites, were lowered in patients with PD relative to controls indicating a potential role caffeine and its metabolites against neuronal damage. Seven metabolites, namely cysteine-S-sulfate, 1-methylxanthine, vanillic acid, N-acetylaspartic acid, 3-N-acetyl tryptophan, 5-methoxytryptophol, and 13-HODE yielded a ROC curve with a high classification accuracy (AUC 0.977). Comparison between different PD stages showed that cysteine-S-sulfate levels were significantly increasing with the advancement of PD stages while LPI 20:4 was significantly decreasing with disease progression. Our findings provide new biomarker candidates to assist in the diagnosis of PD and monitor its progression. Unusual metabolites like cysteine-S-sulfate might point to therapeutic targets that could enhance the development of novel PD treatments, such as NMDA antagonists.

Data processing solutions to render metabolomics more quantitative: case studies in food and clinical metabolomics using Metabox 2.0.

In classic semiquantitative metabolomics, metabolite intensities are affected by biological factors and other unwanted variations. A systematic evaluation of the data processing methods is crucial to identify adequate processing procedures for a given experimental setup. Current comparative studies are mostly focused on peak area data but not on absolute concentrations. In this study, we evaluated data processing methods to produce outputs that were most similar to the corresponding absolute quantified data. We examined the data distribution characteristics, fold difference patterns between 2 metabolites, and sample variance. We used 2 metabolomic datasets from a retail milk study and a lupus nephritis cohort as test cases. When studying the impact of data normalization, transformation, scaling, and combinations of these methods, we found that the cross-contribution compensating multiple standard normalization (ccmn) method, followed by square root data transformation, was most appropriate for a well-controlled study such as the milk study dataset. Regarding the lupus nephritis cohort study, only ccmn normalization could slightly improve the data quality of the noisy cohort. Since the assessment accounted for the resemblance between processed data and the corresponding absolute quantified data, our results denote a helpful guideline for processing metabolomic datasets within a similar context (food and clinical metabolomics). Finally, we introduce Metabox 2.0, which enables thorough analysis of metabolomic data, including data processing, biomarker analysis, integrative analysis, and data interpretation. It was successfully used to process and analyze the data in this study. An online web version is available at http://metsysbio.com/metabox.

Heated tobacco product IQOS induces unique metabolic signatures in human bronchial epithelial cells.

Metabolic signatures are lacking for heated tobacco products, making it crucial to identify new biosignatures of lung damage. This will enable the establishment of product-specific guidelines and an understanding of associated toxicity. https://bit.ly/3TkhBox.

Hydrophilic Interaction Liquid Chromatography-Hydrogen/Deuterium Exchange-Mass Spectrometry (HILIC-HDX-MS) for Untargeted Metabolomics.

Liquid chromatography with mass spectrometry (LC-MS)-based metabolomics detects thousands of molecular features (retention time-m/z pairs) in biological samples per analysis, yet the metabolite annotation rate remains low, with 90% of signals classified as unknowns. To enhance the metabolite annotation rates, researchers employ tandem mass spectral libraries and challenging in silico fragmentation software. Hydrogen/deuterium exchange mass spectrometry (HDX-MS) may offer an additional layer of structural information in untargeted metabolomics, especially for identifying specific unidentified metabolites that are revealed to be statistically significant. Here, we investigate the potential of hydrophilic interaction liquid chromatography (HILIC)-HDX-MS in untargeted metabolomics. Specifically, we evaluate the effectiveness of two approaches using hypothetical targets: the post-column addition of deuterium oxide (D2O) and the on-column HILIC-HDX-MS method. To illustrate the practical application of HILIC-HDX-MS, we apply this methodology using the in silico fragmentation software MS-FINDER to an unknown compound detected in various biological samples, including plasma, serum, tissues, and feces during HILIC-MS profiling, subsequently identified as N1-acetylspermidine.

The lipidome of posttraumatic stress disorder.

Posttraumatic stress disorder (PTSD) can develop after trauma exposure. Some studies report that women develop PTSD at twice the rate of men, despite greater trauma exposure in men. Lipids and their metabolites (lipidome) regulate a myriad of key biological processes and pathways such as membrane integrity, oxidative stress, and neuroinflammation in the brain by maintaining neuronal connectivity and homeostasis. In this study, we analyzed the lipidome of 40 individuals with PTSD and 40 trauma-exposed non-PTSD individuals. Plasma samples were analyzed for lipidomics using Quadrupole Time-of-Flight (QToF) mass spectrometry. Additionally, ~ 90 measures were collected, on sleep, mental and physical health indices. Sleep quality worsened as PTSD severity increased in both sexes. The lipidomics analysis identified a total of 348 quantifiable known lipid metabolites and 1951 lipid metabolites that are yet unknown; known metabolites were part of 13 classes of lipids. After adjusting for sleep quality, in women with PTSD, only one lipid subclass, phosphatidylethanolamine (PE) was altered, whereas, in men with PTSD, 9 out of 13 subclasses were altered compared to non-PTSD women and men, respectively. Severe PTSD was associated with 22% and 5% of altered lipid metabolites in men and women, respectively. Of the changed metabolites, only 0.5% measures (2 PEs and cholesterol) were common between women and men with PTSD. Several sphingomyelins, PEs, ceramides, and triglycerides were increased in men with severe PTSD. The triglycerides and ceramide metabolites that were most highly increased were correlated with cholesterol metabolites and systolic blood pressure in men but not always in women with PTSD. Alterations in triglycerides and ceramides are linked with cardiac health and metabolic function in humans. Thus, disturbed sleep and higher weight may have contributed to changes in the lipidome found in PTSD.

Unbiased serum metabolomic analysis in cats with naturally occurring chronic enteropathies before and after medical intervention.

Chronic enteropathies (CE) are common disorders in cats and the differentiation between the two main underlying diseases, inflammatory bowel disease (IBD) and low-grade intestinal T-cell lymphoma (LGITL), can be challenging. Characterization of the serum metabolome could provide further information on alterations of disease-associated metabolic pathways and may identify diagnostic or therapeutic targets. Unbiased metabolomics analysis of serum from 28 cats with CE (14 cats with IBD, 14 cats with LGITL) and 14 healthy controls identified 1,007 named metabolites, of which 129 were significantly different in cats with CE compared to healthy controls at baseline. Random Forest analysis revealed a predictive accuracy of 90% for differentiating controls from cats with chronic enteropathy. Metabolic pathways found to be significantly altered included phospholipids, amino acids, thiamine, and tryptophan metabolism. Several metabolites were found to be significantly different between cats with IBD versus LGITL, including several sphingolipids, phosphatidylcholine 40:7, uridine, pinitol, 3,4-dihydroxybenzoic acid, and glucuronic acid. However, random forest analysis revealed a poor group predictive accuracy of 60% for the differentiation of IBD from LGITL. Of 129 compounds found to be significantly different between healthy cats and cats with CE at baseline, 58 remained different following treatment.

A terminal metabolite of niacin promotes vascular inflammation and contributes to cardiovascular disease risk.

Despite intensive preventive cardiovascular disease (CVD) efforts, substantial residual CVD risk remains even for individuals receiving all guideline-recommended interventions. Niacin is an essential micronutrient fortified in food staples, but its role in CVD is not well understood. In this study, untargeted metabolomics analysis of fasting plasma from stable cardiac patients in a prospective discovery cohort (n = 1,162 total, n = 422 females) suggested that niacin metabolism was associated with incident major adverse cardiovascular events (MACE). Serum levels of the terminal metabolites of excess niacin, N1-methyl-2-pyridone-5-carboxamide (2PY) and N1-methyl-4-pyridone-3-carboxamide (4PY), were associated with increased 3-year MACE risk in two validation cohorts (US n = 2,331 total, n = 774 females; European n = 832 total, n = 249 females) (adjusted hazard ratio (HR) (95% confidence interval) for 2PY: 1.64 (1.10-2.42) and 2.02 (1.29-3.18), respectively; for 4PY: 1.89 (1.26-2.84) and 1.99 (1.26-3.14), respectively). Phenome-wide association analysis of the genetic variant rs10496731, which was significantly associated with both 2PY and 4PY levels, revealed an association of this variant with levels of soluble vascular adhesion molecule 1 (sVCAM-1). Further meta-analysis confirmed association of rs10496731 with sVCAM-1 (n = 106,000 total, n = 53,075 females, P = 3.6 × 10-18). Moreover, sVCAM-1 levels were significantly correlated with both 2PY and 4PY in a validation cohort (n = 974 total, n = 333 females) (2PY: rho = 0.13, P = 7.7 × 10-5; 4PY: rho = 0.18, P = 1.1 × 10-8). Lastly, treatment with physiological levels of 4PY, but not its structural isomer 2PY, induced expression of VCAM-1 and leukocyte adherence to vascular endothelium in mice. Collectively, these results indicate that the terminal breakdown products of excess niacin, 2PY and 4PY, are both associated with residual CVD risk. They also suggest an inflammation-dependent mechanism underlying the clinical association between 4PY and MACE.

metabCombiner 2.0: Disparate Multi-Dataset Feature Alignment for LC-MS Metabolomics.

Liquid chromatography-high-resolution mass spectrometry (LC-HRMS), as applied to untargeted metabolomics, enables the simultaneous detection of thousands of small molecules, generating complex datasets. Alignment is a crucial step in data processing pipelines, whereby LC-MS features derived from common ions are assembled into a unified matrix amenable to further analysis. Variability in the analytical factors that influence liquid chromatography separations complicates data alignment. This is prominent when aligning data acquired in different laboratories, generated using non-identical instruments, or between batches from large-scale studies. Previously, we developed metabCombiner for aligning disparately acquired LC-MS metabolomics datasets. Here, we report significant upgrades to metabCombiner that enable the stepwise alignment of multiple untargeted LC-MS metabolomics datasets, facilitating inter-laboratory reproducibility studies. To accomplish this, a "primary" feature list is used as a template for matching compounds in "target" feature lists. We demonstrate this workflow by aligning four lipidomics datasets from core laboratories generated using each institution's in-house LC-MS instrumentation and methods. We also introduce batchCombine, an application of the metabCombiner framework for aligning experiments composed of multiple batches. metabCombiner is available as an R package on Github and Bioconductor, along with a new online version implemented as an R Shiny App.

Metabolomic Biomarkers of Dietary Approaches to Stop Hypertension (DASH) Dietary Patterns in Pregnant Women.

Objective: the aim of this study was to identify plasma metabolomic markers of Dietary Approaches to Stop Hypertension (DASH) dietary patterns in pregnant women. Methods: This study included 186 women who had both dietary intake and metabolome measured from a nested case-control study within the NICHD Fetal Growth Studies-Singletons cohort (FGS). Dietary intakes were ascertained at 8-13 gestational weeks (GW) using the Food Frequency Questionnaire (FFQ) and DASH scores were calculated based on eight food and nutrient components. Fasting plasma samples were collected at 15-26 GW and untargeted metabolomic profiling was performed. Multivariable linear regression models were used to examine the association of individual metabolites with the DASH score. Least absolute shrinkage and selection operator (LASSO) regression was used to select a panel of metabolites jointly associated with the DASH score. Results: Of the total 460 known metabolites, 92 were individually associated with DASH score in linear regressions, 25 were selected as a panel by LASSO regressions, and 18 were identified by both methods. Among the top 18 metabolites, there were 11 lipids and lipid-like molecules (i.e., TG (49:1), TG (52:2), PC (31:0), PC (35:3), PC (36:4) C, PC (36:5) B, PC (38:4) B, PC (42:6), SM (d32:0), gamma-tocopherol, and dodecanoic acid), 5 organic acids and derivatives (i.e., asparagine, beta-alanine, glycine, taurine, and hydroxycarbamate), 1 organic oxygen compound (i.e., xylitol), and 1 organoheterocyclic compound (i.e., maleimide). Conclusions: our study identified plasma metabolomic markers for DASH dietary patterns in pregnant women, with most of being lipids and lipid-like molecules.

Longitudinal Lipidomic Signature of Coronary Heart Disease in American Indian People.

BACKGROUND: Dyslipidemia is an independent risk factor for coronary heart disease (CHD). Standard lipid panel cannot capture the complexity of the blood lipidome (ie, all molecular lipids in the blood). To date, very few large-scale epidemiological studies have assessed the full spectrum of the blood lipidome on risk of CHD, especially in a longitudinal setting. METHODS AND RESULTS: Using an untargeted liquid chromatography-mass spectrometry, we repeatedly measured 1542 lipid species from 1835 unique American Indian participants who attended 2 clinical visits (≈5.5 years apart) and followed up to 17.8 years in the Strong Heart Family Study (SHFS). We first identified baseline lipid species associated with risk of CHD, followed by replication in a European population. The model adjusted for age, sex, body mass index, smoking, hypertension, diabetes, low-density lipoprotein cholesterol, estimated glomerular filtration rate, education, and physical activity at baseline. We then examined the longitudinal association between changes in lipid species and changes in cardiovascular risk factors during follow-up. Multiple testing was controlled by the false discovery rate. We found that baseline levels of multiple lipid species (eg, phosphatidylcholines, phosphatidylethanolamines, and ceramides) were associated with the risk of CHD and improved the prediction accuracy over conventional risk factors in American Indian people. Some identified lipids in American Indian people were replicated in European people. Longitudinal changes in multiple lipid species (eg, acylcarnitines, phosphatidylcholines, and triacylglycerols) were associated with changes in cardiovascular risk factors. CONCLUSIONS: Baseline plasma lipids and their longitudinal changes over time are associated with risk of CHD. These findings provide novel insights into the role of dyslipidemia in CHD.

IDSL.GOA: gene ontology analysis for interpreting metabolomic datasets.

Biological interpretation of metabolomic datasets often ends at a pathway analysis step to find the over-represented metabolic pathways in the list of statistically significant metabolites. However, definitions of biochemical pathways and metabolite coverage vary among different curated databases, leading to missed interpretations. For the lists of genes, transcripts and proteins, Gene Ontology (GO) terms over-presentation analysis has become a standardized approach for biological interpretation. But, GO analysis has not been achieved for metabolomic datasets. We present a new knowledgebase (KB) and the online tool, Gene Ontology Analysis by the Integrated Data Science Laboratory for Metabolomics and Exposomics (IDSL.GOA) to conduct GO over-representation analysis for a metabolite list. The IDSL.GOA KB covers 2393 metabolic GO terms and associated 3144 genes, 1,492 EC annotations, and 2621 metabolites. IDSL.GOA analysis of a case study of older versus young female brain cortex metabolome highlighted 82 GO terms being significantly overrepresented (FDR < 0.05). We showed how IDSL.GOA identified key and relevant GO metabolic processes that were not yet covered in other pathway databases. Overall, we suggest that interpretation of metabolite lists should not be limited to only pathway maps and can also leverage GO terms as well. IDSL.GOA provides a useful tool for this purpose, allowing for a more comprehensive and accurate analysis of metabolite pathway data. IDSL.GOA tool can be accessed at https://goa.idsl.me/ .

IDSL.GOA: Gene Ontology Analysis for Interpreting Metabolomic datasets.

Biological interpretation of metabolomic datasets often ends at a pathway analysis step to find the over-represented metabolic pathways in the list of statistically significant metabolites. However, definitions of biochemical pathways and metabolite coverage vary among different curated databases, leading to missed interpretations. For the lists of genes, transcripts and proteins, Gene Ontology (GO) terms over-presentation analysis has become a standardized approach for biological interpretation. But, GO analysis has not been achieved for metabolomic datasets. We present a new knowledgebase (KB) and the online tool, Gene Ontology Analysis by the Integrated Data Science Laboratory for Metabolomics and Exposomics (IDSL.GOA) to conduct GO over-representation analysis for a metabolite list. The IDSL.GOA KB covers 2,393 metabolic GO terms and associated 3,144 genes, 1,492 EC annotations, and 2,621 metabolites. IDSL.GOA analysis of a case study of older vs young female brain cortex metabolome highlighted 82 GO terms being significantly overrepresented (FDR <0.05). We showed how IDSL.GOA identified key and relevant GO metabolic processes that were not yet covered in other pathway databases. Overall, we suggest that interpretation of metabolite lists should not be limited to only pathway maps and can also leverage GO terms as well. IDSL.GOA provides a useful tool for this purpose, allowing for a more comprehensive and accurate analysis of metabolite pathway data. IDSL.GOA tool can be accessed at https://goa.idsl.me/.

Plasma lipidomic markers of diet quality are associated with incident coronary heart disease in American Indian adults: the Strong Heart Family Study.

BACKGROUND: Identifying lipidomic markers of diet quality is needed to inform the development of biomarkers of diet, and to understand the mechanisms driving the diet- coronary heart disease (CHD) association. OBJECTIVES: This study aimed to identify lipidomic markers of diet quality and examine whether these lipids are associated with incident CHD. METHODS: Using liquid chromatography-mass spectrometry, we measured 1542 lipid species from 1694 American Indian adults (aged 18-75 years, 62% female) in the Strong Heart Family Study. Participants were followed up for development of CHD through 2020. Information on the past year diet was collected using the Block Food Frequency Questionnaire, and diet quality was assessed using the Alternative Healthy Eating Index-2010 (AHEI). Mixed-effects linear regression was used to identify individual lipids cross-sectionally associated with AHEI. In prospective analysis, Cox frailty model was used to estimate the hazard ratio (HR) of each AHEI-related lipid for incident CHD. All models were adjusted for age, sex, center, education, body mass index, smoking, alcohol drinking, level of physical activity, energy intake, diabetes, hypertension, and use of lipid-lowering drugs. Multiple testing was controlled at a false discovery rate of <0.05. RESULTS: Among 1542 lipid species measured, 71 lipid species (23 known), including acylcarnitine, cholesterol esters, glycerophospholipids, sphingomyelins and triacylglycerols, were associated with AHEI. Most of the identified lipids were associated with consumption of ω-3 (n-3) fatty acids. In total, 147 participants developed CHD during a mean follow-up of 17.8 years. Among the diet-related lipids, 10 lipids [5 known: cholesterol ester (CE)(22:5)B, phosphatidylcholine (PC)(p-14:0/22:1)/PC(o-14:0/22:1), PC(p-38:3)/PC(o-38:4)B, phosphatidylethanolamine (PE)(p-18:0/20:4)/PE(o-18:0/20:4), and sphingomyelin (d36:2)A] were associated with incident CHD. On average, each standard deviation increase in the baseline level of these 5 lipids was associated with 17%-23% increased risk of CHD (from HR: 1.17; 95% CI: 1, 1.36; to HR: 1.23; 95% CI: 1.05, 1.43). CONCLUSIONS: In this study, lipidomic markers of diet quality in American Indian adults are found. Some diet-related lipids are associated with risk of CHD beyond established risk factors.