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.

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.

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.

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.

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.

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.

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.

Identification of pre-diagnostic lipid sets associated with liver cancer risk using untargeted lipidomics and chemical set analysis: A nested case-control study within the ATBC cohort

In pre-disposed individuals, a reprogramming of the hepatic lipid metabolism may support liver cancer initiation. We conducted a high-resolution mass spectrometry based untargeted lipidomics analysis of pre-diagnostic serum samples from a nested case-control study (219 liver cancer cases and 219 controls) within the Alpha-Tocopherol, Beta-Carotene Cancer Prevention (ATBC) Study. Out of 462 annotated lipids, 158 (34.2%) were associated with liver cancer risk in a conditional logistic regression analysis at a false discovery rate (FDR) <0.05. A chemical set enrichment analysis (ChemRICH) and co-regulatory set analysis suggested that 22/28 lipid classes and 47/83 correlation modules were significantly associated with liver cancer risk (FDR <0.05). Strong positive associations were observed for monounsaturated fatty acids (MUFA), triacylglycerols (TAGs) and phosphatidylcholines (PCs) having MUFA acyl chains. Negative associations were observed for sphingolipids (ceramides and sphingomyelins), lysophosphatidylcholines, cholesterol esters and polyunsaturated fatty acids (PUFA) containing TAGs and PCs. Stearoyl-CoA desaturase enzyme 1 (SCD1), a rate limiting enzyme in fatty acid metabolism and ceramidases seems to be critical in this reprogramming. In conclusion, our study reports pre-diagnostic lipid changes that provide novel insights into hepatic lipid metabolism reprogramming may contribute to a pro-cell growth and anti-apoptotic tissue environment and, in turn, support liver cancer initiation.

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/.

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.

HNF4α isoforms regulate the circadian balance between carbohydrate and lipid metabolism in the liver

Hepatocyte Nuclear Factor 4α (HNF4α), a master regulator of hepatocyte differentiation, is regulated by two promoters (P1 and P2) which drive the expression of different isoforms. P1-HNF4α is the major isoform in the adult liver while P2-HNF4α is thought to be expressed only in fetal liver and liver cancer. Here, we show that P2-HNF4α is indeed expressed in the normal adult liver at Zeitgeber time (ZT)9 and ZT21. Using exon swap mice that express only P2-HNF4α we show that this isoform orchestrates a distinct transcriptome and metabolome via unique chromatin and protein-protein interactions, including with different clock proteins at different times of the day leading to subtle differences in circadian gene regulation. Furthermore, deletion of the Clock gene alters the circadian oscillation of P2- (but not P1-)HNF4α RNA, revealing a complex feedback loop between the HNF4α isoforms and the hepatic clock. Finally, we demonstrate that while P1-HNF4α drives gluconeogenesis, P2-HNF4α drives ketogenesis and is required for elevated levels of ketone bodies in female mice. Taken together, we propose that the highly conserved two-promoter structure of the Hnf4a gene is an evolutionarily conserved mechanism to maintain the balance between gluconeogenesis and ketogenesis in the liver in a circadian fashion.

Longitudinal lipidomic profiles during pregnancy and associations with neonatal anthropometry: findings from a multiracial cohort

BACKGROUND

Maternal lipidomic profiling offers promise for characterizing lipid metabolites during pregnancy, but longitudinal data are limited. This study aimed to examine associations of longitudinal lipidomic profiles during pregnancy with multiple neonatal anthropometry using data from a multiracial cohort.

METHODS

We measured untargeted plasma lipidome profiles among 321 pregnant women from the NICHD Fetal Growth Study-Singletons using plasma samples collected longitudinally during four study visits at gestational weeks (GW) 10-14, 15-26, 23-31, and 33-39, respectively. We evaluated individual lipidomic metabolites at each study visit in association with neonatal anthropometry. We also evaluated the associations longitudinally by constructing lipid networks using weighted correlation network analysis and common networks using consensus network analysis across four visits using linear mixed-effects models with the adjustment of false discover rate.

FINDINGS

Multiple triglycerides (TG) were positively associated with birth weight (BW), BW Z-score, length and head circumference, while some cholesteryl ester (CE), phosphatidylcholine (PC), sphingomyelines (SM), phosphatidylethanolamines (PE), and lysophosphatidylcholines (LPC 20:3) families were inversely associated with BW, length, abdominal and head circumference at different GWs. Longitudinal trajectories of TG, PC, and glucosylcermides (GlcCer) were associated with BW, and CE (18:2) with BW z-score, length, and sum of skinfolds (SS), while some PC and PE were significantly associated with abdominal and head circumference. Modules of TG at GW 10-14 and 15-26 mainly were associated with BW. At GW 33-39, two networks of LPC (20:3) and of PC, TG, and CE, showed inverse associations with abdominal circumference. Distinct trajectories within two consensus modules with changes in TG, CE, PC, and LPC showed significant differences in BW and length.

INTERPRETATION

The results demonstrated that longitudinal changes of TGs during early- and mid-pregnancy and changes of PC, LPC, and CE during late-pregnancy were significantly associated with neonatal anthropometry.

FUNDING

National Institute of Child Health and Human Development intramural funding.

LibGen: Generating High Quality Spectral Libraries of Natural Products for EAD-, UVPD-, and HCD-High Resolution Mass Spectrometers

Compound annotation using spectral-matching algorithms is vital for (MS/MS)-based metabolomics research, but is hindered by the lack of high-quality reference MS/MS library spectra. Finding and removing errors from libraries, including noise ions, is mostly done manually. This process is both error-prone and time-consuming. To address these challenges, we have developed an automated library curation pipeline, LibGen, to universally build novel spectral libraries. This pipeline corrects mass errors, denoises spectra by subformula assignments, and performs quality control of the reference spectra by calculating explained intensity and spectral entropy. We employed LibGen to generate three high-quality libraries with chemical standards of 2241 natural products. To this end, we used an IQ-X orbital ion trap mass spectrometer to generate 1947 classic high-energy collision dissociation spectra (HCD) as well as 1093 ultraviolet-photodissociation (UVPD) mass spectra. The third library was generated by an electron-activated collision dissociation (EAD) 7600 ZenoTOF mass spectrometer yielding 3244 MS/MS spectra. The natural compounds covered 140 chemical classes from prenol lipids to benzypyrans with >97% of the compounds showing <0.2 Tanimoto-similarity, demonstrating a very high structural variance. Mass spectra showed much higher information content for both UVPD- and EAD-mass spectra compared to classic HCD spectra when using spectral entropy calculations. We validated the denoising algorithm by acquiring MS/MS spectra at high concentration and at 13-fold diluted chemical standards. At low concentrations, a higher proportion of spectra showed apparent fragment ions that could not be explained by subformula losses of the parent molecule. When more than 10% of the total intensity of MS/MS fragments was regarded as noise ions, spectra were considered as low quality and were not included in the libraries. As the overall process is fully automated, LibGen can be utilized by all researchers who create or curate mass spectral libraries. The libraries we created here are publicly available at MassBank.us.

Lipid metabolism drives allele-specific early-stage hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) results from pathogenic variants in sarcomeric protein genes, that increase myocyte energy demand and lead to cardiac hypertrophy. But it is unknown whether a common metabolic trait underlies the cardiac phenotype at early disease stage. This study characterized two HCM mouse models (R92W-TnT, R403Q-MyHC) that demonstrate differences in mitochondrial function at early disease stage. Using a combination of cardiac phenotyping, transcriptomics, mass spectrometry-based metabolomics and computational modeling, we discovered allele-specific differences in cardiac structure/function and metabolic changes. TnT-mutant hearts had impaired energy substrate metabolism and increased phospholipid remodeling compared to MyHC-mutants. TnT-mutants showed increased incorporation of saturated fatty acid residues into ceramides, cardiolipin, and increased lipid peroxidation, that could underlie allele-specific differences in mitochondrial function and cardiomyopathy.

Microbe-derived uremic solutes enhance thrombosis potential in the host

p-Cresol sulfate (pCS) and indoxyl sulfate (IS), gut microbiome-derived metabolites, are traditionally associated with cardiovascular disease (CVD) risks in the setting of impaired kidney function. While pharmacologic provision of pCS or IS can promote pro-thrombotic phenotypes, neither the microbial enzymes involved nor direct gut microbial production have been linked to CVD. Untargeted metabolomics was performed on a discovery cohort (n = 1,149) with relatively preserved kidney function, followed by stable isotope-dilution mass spectrometry quantification of pCS and IS in an independent validation cohort (n = 3,954). Genetic engineering of human commensals to produce p-cresol and indole gain-of-function and loss-of-function mutants, followed by colonization of germ-free mice, and studies on host thrombosis were performed. Systemic pCS and IS levels were independently associated with all-cause mortality. Both in vitro and within colonized germ-free mice p-cresol productions were recapitulated by collaboration of two organisms: a Bacteroides strain that converts tyrosine to 4-hydroxyphenylacetate, and a Clostridium strain that decarboxylates 4-hydroxyphenylacetate to p-cresol. We then engineered a single organism, Bacteroides thetaiotaomicron, to produce p-cresol, indole, or both metabolites. Colonizing germ-free mice with engineered strains, we show the gut microbial genes for p-cresol (hpdBCA) and indole (tryptophanase) are sufficient to confer a pro-thrombotic phenotype in vivo. Moreover, human fecal metagenomics analyses show that abundances of hpdBCA and tryptophanase are associated with CVD. These studies show that pCS and IS, two abundant microbiome-derived metabolites, play a broader potential role in CVD than was previously known. They also suggest that therapeutic targeting of gut microbial p-cresol- and indole-producing pathways represent rational targets for CVD.IMPORTANCEAlterations in gut microbial composition and function have been linked to numerous diseases. Identifying microbial pathways responsible for producing molecules that adversely impact the host is an important first step in the development of therapeutic interventions. Here, we first use large-scale clinical observations to link blood levels of defined microbial products to cardiovascular disease risks. Notably, the previously identified uremic toxins p-cresol sulfate and indoxyl sulfate were shown to predict 5-year mortality risks. After identifying the microbes and microbial enzymes involved in the generation of these uremic toxins, we used bioengineering technologies coupled with colonization of germ-free mice to show that the gut microbial genes that generate p-cresol and indole are sufficient to confer p-cresol sulfate and indoxyl sulfate formation, and a pro-thrombotic phenotype in vivo. The findings and tools developed serve as a critical step in both the study and targeting of these gut microbial pathways in vivo.

Lipidomics analysis reveals the effect of Sirex noctilio infestation on the lipid metabolism in Pinus radiata needles

The Sirex noctilio's climatic adaption and rapid proliferation have caused Pinus mortality worldwide. The infestation combines the early effect of female S. noctilio gland secretion and the spreading symbiotic fungus Amylostereum areolatum. 'Lipidomics' is the study of all non-water-soluble components of the metabolome. Most of these non-water-soluble compounds correspond to lipids which can provide information about a biological activity, an organelle, an organism, or a disease. Using HPLC-MS/MS based lipidomics, 122 lipids were identified in P. radiata needles during S. noctilio infestation. Phosphatidic acids, N-acylethanolamines, and phosphatidylinositol-ceramides accumulated in infested trees could suggest a high level of phospholipases activities. The phosphatidylcholines were the most down-regulated species during infection, which could also suggest that they may be used as a substrate for up-regulated lipids. The accumulation of very long-chain fatty acids and long-chain fatty acids during the infestation could imply the tree defense response to create a barrier in the drilled zone to avoid larvae development and fungus proliferation. Also, the growth arrest phase of the trees during the prolonged infestation suggests a resistance response, regulated by the accumulation of NAE, which potentially shifts the tree energy to respond to the infestation.

Awakening the natural capability of psicose production in Escherichia coli

Due to the rampant rise in obesity and diabetes, consumers are desperately seeking for ways to reduce their sugar intake, but to date there are no options that are both accessible and without sacrifice of palatability. One of the most promising new ingredients in the food system as a non-nutritive sugar substitute with near perfect palatability is D-psicose. D-psicose is currently produced using an in vitro enzymatic isomerization of D-fructose, resulting in low yield and purity, and therefore requiring substantial downstream processing to obtain a high purity product. This has made adoption of D-psicose into products limited and results in significantly higher per unit costs, reducing accessibility to those most in need. Here, we found that Escherichia coli natively possesses a thermodynamically favorable pathway to produce D-psicose from D-glucose through a series of phosphorylation-epimerization-dephosphorylation steps. To increase carbon flux towards D-psicose production, we introduced a series of genetic modifications to pathway enzymes, central carbon metabolism, and competing metabolic pathways. In an attempt to maximize both cellular viability and D-psicose production, we implemented methods for the dynamic regulation of key genes including clustered regularly interspaced short palindromic repeats inhibition (CRISPRi) and stationary-phase promoters. The engineered strains achieved complete consumption of D-glucose and production of D-psicose, at a titer of 15.3 g L-1, productivity of 2 g L-1 h-1, and yield of 62% under test tube conditions. These results demonstrate the viability of whole-cell catalysis as a sustainable alternative to in vitro enzymatic synthesis for the accessible production of D-psicose.

Longitudinal lipidomic signature of carotid atherosclerosis in American Indians: Findings from the Strong Heart Family Study

BACKGROUND AND AIMS

Dyslipidemia is an independent risk factor for atherosclerosis and atherosclerotic cardiovascular disease (ASCVD). To date, a comprehensive assessment of individual lipid species associated with atherosclerosis is lacking in large-scale epidemiological studies, especially in a longitudinal setting. We investigated the association of circulating lipid species and its longitudinal changes with carotid atherosclerosis.

METHODS

Using liquid chromatograph-mass spectrometry, we repeatedly measured 1542 lipid species in 3687 plasma samples from 1918 unique American Indians attending two visits (mean ∼5 years apart) in the Strong Heart Family Study. Carotid atherosclerotic plaques were assessed by ultrasonography at each visit. We identified lipids associated with prevalence or progression of carotid plaques, adjusting age, sex, BMI, smoking, hypertension, diabetes, and eGFR. Then we examined whether longitudinal changes in lipids were associated with changes in cardiovascular risk factors. Multiple testing was controlled at false discovery rate (FDR) < 0.05.

RESULTS

Higher levels of sphingomyelins, ether-phosphatidylcholines, and triacylglycerols were significantly associated with prevalence or progression of carotid plaques (odds ratios ranged from 1.15 to 1.34). Longitudinal changes in multiple lipid species (e.g., acylcarnitines, phosphatidylcholines, triacylglycerols) were associated with changes in cardiometabolic traits (e.g., BMI, blood pressure, fasting glucose, eGFR). Network analysis identified differential lipid networks associated with plaque progression.

CONCLUSIONS

Baseline and longitudinal changes in multiple lipid species were significantly associated with carotid atherosclerosis and its progression in American Indians. Some plaque-related lipid species were also associated with risk for CVD events.

Flash entropy search to query all mass spectral libraries in real time

Public repositories of metabolomics mass spectra encompass more than 1 billion entries. With open search, dot product or entropy similarity, comparisons of a single tandem mass spectrometry spectrum take more than 8 h. Flash entropy search speeds up calculations more than 10,000 times to query 1 billion spectra in less than 2 s, without loss in accuracy. It benefits from using multiple threads and GPU calculations. This algorithm can fully exploit large spectral libraries with little memory overhead for any mass spectrometry laboratory.

Serum-Based Lipid Panels for Diagnosis of Idiopathic Parkinson's Disease

Parkinson's disease (PD) is a highly prevalent neurodegenerative movement disorder with an unclear etiology and a lack of definite diagnostic tests and effective treatments. About 95% of PD cases are idiopathic, in which none of the well-known genes underlying familial parkinsonism are mutated. We used untargeted liquid chromatography-mass spectrometry (LC-MS/MS) to profile the serum lipidome of 50 patients with different stages of idiopathic PD (early, mid, or advanced) and 45 age-matched controls. When comparing the PD patients to the control subjects, 169 lipids were significantly altered in both a univariate analysis and a multivariate partial least-squares discriminant analysis (PLS-DA). Compared to the controls, the patients with PD had higher levels of unsaturated triacylglycerides (e.g., TG O-56:9 and TG 52:3), saturated lysophosphatidylcholines (LPC 17:0, 16:0, and 15:0), and hydroxyeicosatetraenoic acid (12-HETE), while lower levels of phosphatidylserines (e.g., PS 40:4 and PS 16:0_22:4), sphingomyelins (SM 42:1), and ceramides (e.g., Cer 40:0 and 42:0) were found between the PD patients and the controls. A panel of 10 significantly altered lipids (PS 40:0, Cer 40:0, Cer 42:0, LPC 17:0, LPC 15:0, PC 37:7, PE O-40:8, PC O-42:4, FA 23:0, and SM 42:1) resulted in a strong receiver operating characteristic curve with an AUC = 0.974. This panel may, therefore, be useful for diagnosing PD. In addition, lipid panels may prove useful for distinguishing among the progression stages of PD. Using one-way ANOVA, 155 lipid species were significantly altered among the PD stages. Parkinson's disease progressed from the early to advanced stages with decreasing levels of PC 31:1, PC 38:4, and LPE 22:5. Conversely, LPC-O 20:0, PC O-42:3, FA 19:0, and FA 22:2 showed an increase in their levels with disease progression. Overall, this study shows an intriguing number of robust changes in specific serum lipids that may become useful for diagnosing PD and its progression, once panels have been validated in larger clinical trials and prospective studies.

Predicting chronic postsurgical pain: current evidence and a novel program to develop predictive biomarker signatures

ABSTRACT

Chronic pain affects more than 50 million Americans. Treatments remain inadequate, in large part, because the pathophysiological mechanisms underlying the development of chronic pain remain poorly understood. Pain biomarkers could potentially identify and measure biological pathways and phenotypical expressions that are altered by pain, provide insight into biological treatment targets, and help identify at-risk patients who might benefit from early intervention. Biomarkers are used to diagnose, track, and treat other diseases, but no validated clinical biomarkers exist yet for chronic pain. To address this problem, the National Institutes of Health Common Fund launched the Acute to Chronic Pain Signatures (A2CPS) program to evaluate candidate biomarkers, develop them into biosignatures, and discover novel biomarkers for chronification of pain after surgery. This article discusses candidate biomarkers identified by A2CPS for evaluation, including genomic, proteomic, metabolomic, lipidomic, neuroimaging, psychophysical, psychological, and behavioral measures. Acute to Chronic Pain Signatures will provide the most comprehensive investigation of biomarkers for the transition to chronic postsurgical pain undertaken to date. Data and analytic resources generatedby A2CPS will be shared with the scientific community in hopes that other investigators will extract valuable insights beyond A2CPS's initial findings. This article will review the identified biomarkers and rationale for including them, the current state of the science on biomarkers of the transition from acute to chronic pain, gaps in the literature, and how A2CPS will address these gaps.

Adipocytes reprogram cancer cell metabolism by diverting glucose towards glycerol-3-phosphate thereby promoting metastasis

In the tumor microenvironment, adipocytes function as an alternate fuel source for cancer cells. However, whether adipocytes influence macromolecular biosynthesis in cancer cells is unknown. Here we systematically characterized the bidirectional interaction between primary human adipocytes and ovarian cancer (OvCa) cells using multi-platform metabolomics, imaging mass spectrometry, isotope tracing and gene expression analysis. We report that, in OvCa cells co-cultured with adipocytes and in metastatic tumors, a part of the glucose from glycolysis is utilized for the biosynthesis of glycerol-3-phosphate (G3P). Normoxic HIF1α protein regulates the altered flow of glucose-derived carbons in cancer cells, resulting in increased glycerophospholipids and triacylglycerol synthesis. The knockdown of HIF1α or G3P acyltransferase 3 (a regulatory enzyme of glycerophospholipid synthesis) reduced metastasis in xenograft models of OvCa. In summary, we show that, in an adipose-rich tumor microenvironment, cancer cells generate G3P as a precursor for critical membrane and signaling components, thereby promoting metastasis. Targeting biosynthetic processes specific to adipose-rich tumor microenvironments might be an effective strategy against metastasis.

Bayesian Statistics Improves Biological Interpretability of Metabolomics Data from Human Cohorts

Univariate analyses of metabolomics data currently follow a frequentist approach, using p-values to reject a null hypothesis. We here propose the use of Bayesian statistics to quantify evidence supporting different hypotheses and discriminate between the null hypothesis versus the lack of statistical power. We used metabolomics data from three independent human cohorts that studied the plasma signatures of subjects with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). The data are publicly available, covering 84-197 subjects in each study with 562-888 identified metabolites of which 777 were common between the two studies and 93 were compounds reported in all three studies. We show how Bayesian statistics incorporates results from one study as "prior information" into the next study, thereby improving the overall assessment of the likelihood of finding specific differences between plasma metabolite levels. Using classic statistics and Benjamini-Hochberg FDR-corrections, Study 1 detected 18 metabolic differences and Study 2 detected no differences. Using Bayesian statistics on the same data, we found a high likelihood that 97 compounds were altered in concentration in Study 2, after using the results of Study 1 as the prior distributions. These findings included lower levels of peroxisome-produced ether-lipids, higher levels of long-chain unsaturated triacylglycerides, and the presence of exposome compounds that are explained by the difference in diet and medication between healthy subjects and ME/CFS patients. Although Study 3 reported only 92 compounds in common with the other two studies, these major differences were confirmed. We also found that prostaglandin F2alpha, a lipid mediator of physiological relevance, was reduced in ME/CFS patients across all three studies. The use of Bayesian statistics led to biological conclusions from metabolomic data that were not found through frequentist approaches. We propose that Bayesian statistics is highly useful for studies with similar research designs if similar metabolomic assays are used.

Improving Quantitative Accuracy in Nontargeted Lipidomics by Evaluating Adduct Formation

For large-scale lipidomic analyses, accurate and reproducible quantification of endogenous lipids is crucial for comparing results within and across studies. Many lipids present in liquid chromatography-electrospray ionization-mass spectrometry form various adducts with buffer components. The mechanisms and conditions that dictate adduct formation are still poorly understood. In a positive mode, neutral lipids like mono-, di-, and triacylglycerides and cholesteryl esters typically generate [M + NH4]+ adduct ions, although [M + Na]+, [M + K]+, and other (more complex) species can also be significantly abundant in MS1 precursor ion spectra. Variations in the ratios of these adducts (within and between matrices) can lead to dramatic inaccuracies during quantification. Here, we examine 48 unique diacylglycerol (DAG) species across 2366 mouse samples for eight matrix-specific data sets of plasma, liver, kidney, brain, heart muscle, gastrocnemius muscle, gonadal, and inguinal fat. Typically, no single adduct ion species accounted for more than 60% of the total observed abundance across each data set. Even within a single matrix, DAGs showed a high variability of adduct ratios. The ratio of [M + NH4]+ adduct ions was increased for longer-chain DAGs and for polyunsaturated DAGs, at the expense of reduced ratios of [M + Na]+ adducts. When using three deuterated internal DAG standards, we found that absolute concentrations were estimated with up to 70% error when only one adduct ion was used instead of all adducts combined. Importantly, when combining [M + NH4]+ and [M + Na]+ adduct ions, quantification results were within 5% accuracy compared to all adduct ions combined. Additional variance can be caused by other factors, such as instrument conditions or matrix effects.

Exploring the Impact of Organic Solvent Quality and Unusual Adduct Formation during LC-MS-Based Lipidomic Profiling

Liquid chromatography-mass spectrometry (LC-MS) is the key technique for analyzing complex lipids in biological samples. Various LC-MS modes are used for lipid separation, including different stationary phases, mobile-phase solvents, and modifiers. Quality control in lipidomics analysis is crucial to ensuring the generated data's reliability, reproducibility, and accuracy. While several quality control measures are commonly discussed, the impact of organic solvent quality during LC-MS analysis is often overlooked. Additionally, the annotation of complex lipids remains prone to biases, leading to potential misidentifications and incomplete characterization of lipid species. In this study, we investigate how LC-MS-grade isopropanol from different vendors may influence the quality of the mobile phase used in LC-MS-based untargeted lipidomic profiling of biological samples. Furthermore, we report the occurrence of an unusual, yet highly abundant, ethylamine adduct [M+46.0651]+ that may form for specific lipid subclasses during LC-MS analysis in positive electrospray ionization mode when acetonitrile is part of the mobile phase, potentially leading to lipid misidentification. These findings emphasize the importance of considering solvent quality in LC-MS analysis and highlight challenges in lipid annotation.

Automatic Assignment of Molecular Ion Species to Elemental Formulas in Gas Chromatography/Methane Chemical Ionization Accurate Mass Spectrometry

Gas chromatography-mass spectrometry (GC-MS) usually employs hard electron ionization, leading to extensive fragmentations that are suitable to identify compounds based on library matches. However, such spectra are less useful to structurally characterize unknown compounds that are absent from libraries, due to the lack of readily recognizable molecular ion species. We tested methane chemical ionization on 369 trimethylsilylated (TMS) derivatized metabolites using a quadrupole time-of-flight detector (QTOF). We developed an algorithm to automatically detect molecular ion species and tested SIRIUS software on how accurate the determination of molecular formulas was. The automatic workflow correctly recognized 289 (84%) of all 345 detected derivatized standards. Specifically, strong [M - CH3]+ fragments were observed in 290 of 345 derivatized chemicals, which enabled the automatic recognition of molecular adduct patterns. Using Sirius software, correct elemental formulas were retrieved in 87% of cases within the top three hits. When investigating the cases for which the automatic pattern analysis failed, we found that several metabolites showed a previously unknown [M + TMS]+ adduct formed by rearrangement. Methane chemical ionization with GC-QTOF mass spectrometry is a suitable avenue to identify molecular formulas for abundant unknown peaks.

Sexual Dimorphism of the Mouse Plasma Metabolome Is Associated with Phenotypes of 30 Gene Knockout Lines

Although metabolic alterations are observed in many monogenic and complex genetic disorders, the impact of most mammalian genes on cellular metabolism remains unknown. Understanding the effect of mouse gene dysfunction on metabolism can inform the functions of their human orthologues. We investigated the effect of loss-of-function mutations in 30 unique gene knockout (KO) lines on plasma metabolites, including genes coding for structural proteins (11 of 30), metabolic pathway enzymes (12 of 30) and protein kinases (7 of 30). Steroids, bile acids, oxylipins, primary metabolites, biogenic amines and complex lipids were analyzed with dedicated mass spectrometry platforms, yielding 827 identified metabolites in male and female KO mice and wildtype (WT) controls. Twenty-two percent of 23,698 KO versus WT comparison tests showed significant genotype effects on plasma metabolites. Fifty-six percent of identified metabolites were significantly different between the sexes in WT mice. Many of these metabolites were also found to have sexually dimorphic changes in KO lines. We used plasma metabolites to complement phenotype information exemplified for Dhfr, Idh1, Mfap4, Nek2, Npc2, Phyh and Sra1. The association of plasma metabolites with IMPC phenotypes showed dramatic sexual dimorphism in wildtype mice. We demonstrate how to link metabolomics to genotypes and (disease) phenotypes. Sex must be considered as critical factor in the biological interpretation of gene functions.

Profile Characterization of Biogenic Amines in Glioblastoma Patients Undergoing Standard-of-Care Treatment

INTRODUCTION

Biogenic amines play important roles throughout cellular metabolism. This study explores a role of biogenic amines in glioblastoma pathogenesis. Here, we characterize the plasma levels of biogenic amines in glioblastoma patients undergoing standard-of-care treatment.

METHODS

We examined 138 plasma samples from 36 patients with isocitrate dehydrogenase (IDH) wild-type glioblastoma at multiple stages of treatment. Untargeted gas chromatography-time of flight mass spectrometry (GC-TOF MS) was used to measure metabolite levels. Machine learning approaches were then used to develop a predictive tool based on these datasets.

RESULTS

Surgery was associated with increased levels of 12 metabolites and decreased levels of 11 metabolites. Chemoradiation was associated with increased levels of three metabolites and decreased levels of three other metabolites. Ensemble learning models, specifically random forest (RF) and AdaBoost (AB), accurately classified treatment phases with high accuracy (RF: 0.81 ± 0.04, AB: 0.78 ± 0.05). The metabolites sorbitol and N-methylisoleucine were identified as important predictive features and confirmed via SHAP.

CONCLUSION

To our knowledge, this is the first study to describe plasma biogenic amine signatures throughout the treatment of patients with glioblastoma. A larger study is needed to confirm these results with hopes of developing a diagnostic algorithm.

Denoising Autoencoder Normalization for Large-Scale Untargeted Metabolomics by Gas Chromatography-Mass Spectrometry

Large-scale metabolomics assays are widely used in epidemiology for biomarker discovery and risk assessments. However, systematic errors introduced by instrumental signal drifting pose a big challenge in large-scale assays, especially for derivatization-based gas chromatography-mass spectrometry (GC-MS). Here, we compare the results of different normalization methods for a study with more than 4000 human plasma samples involved in a type 2 diabetes cohort study, in addition to 413 pooled quality control (QC) samples, 413 commercial pooled plasma samples, and a set of 25 stable isotope-labeled internal standards used for every sample. Data acquisition was conducted across 1.2 years, including seven column changes. In total, 413 pooled QC (training) and 413 BioIVT samples (validation) were used for normalization comparisons. Surprisingly, neither internal standards nor sum-based normalizations yielded median precision of less than 30% across all 563 metabolite annotations. While the machine-learning-based SERRF algorithm gave 19% median precision based on the pooled quality control samples, external cross-validation with BioIVT plasma pools yielded a median 34% relative standard deviation (RSD). We developed a new method: systematic error reduction by denoising autoencoder (SERDA). SERDA lowered the median standard deviations of the training QC samples down to 16% RSD, yielding an overall error of 19% RSD when applied to the independent BioIVT validation QC samples. This is the largest study on GC-MS metabolomics ever reported, demonstrating that technical errors can be normalized and handled effectively for this assay. SERDA was further validated on two additional large-scale GC-MS-based human plasma metabolomics studies, confirming the superior performance of SERDA over SERRF or sum normalizations.

SMetaS: A Sample Metadata Standardizer for Metabolomics

Metabolomics has advanced to an extent where it is desired to standardize and compare data across individual studies. While past work in standardization has focused on data acquisition, data processing, and data storage aspects, metabolomics databases are useless without ontology-based descriptions of biological samples and study designs. We introduce here a user-centric tool to automatically standardize sample metadata. Using such a tool in frontends for metabolomic databases will dramatically increase the FAIRness (Findability, Accessibility, Interoperability, and Reusability) of data, specifically for data reuse and for finding datasets that share comparable sets of metadata, e.g., study meta-analyses, cross-species analyses or large scale metabolomic atlases. SMetaS (Sample Metadata Standardizer) combines a classic database with an API and frontend and is provided in a containerized environment. The tool has two user-centric components. In the first component, the user designs a sample metadata matrix and fills the cells using natural language terminology. In the second component, the tool transforms the completed matrix by replacing freetext terms with terms from fixed vocabularies. This transformation process is designed to maximize simplicity and is guided by, among other strategies, synonym matching and typographical fixing in an n-grams/nearest neighbors model approach. The tool enables downstream analysis of submitted studies and samples via string equality for FAIR retrospective use.

Differences in the Stool Metabolome between Vegans and Omnivores: Analyzing the NIST Stool Reference Material

To gain confidence in results of omic-data acquisitions, methods must be benchmarked using validated quality control materials. We report data combining both untargeted and targeted metabolomics assays for the analysis of four new human fecal reference materials developed by the U.S. National Institute of Standards and Technologies (NIST) for metagenomics and metabolomics measurements. These reference grade test materials (RGTM) were established by NIST based on two different diets and two different samples treatments, as follows: firstly, homogenized fecal matter from subjects eating vegan diets, stored and submitted in either lyophilized (RGTM 10162) or aqueous form (RGTM 10171); secondly, homogenized fecal matter from subjects eating omnivore diets, stored and submitted in either lyophilized (RGTM 10172) or aqueous form (RGTM 10173). We used four untargeted metabolomics assays (lipidomics, primary metabolites, biogenic amines and polyphenols) and one targeted assay on bile acids. A total of 3563 compounds were annotated by mass spectrometry, including 353 compounds that were annotated in more than one assay. Almost half of all compounds were annotated using hydrophilic interaction chromatography/accurate mass spectrometry, followed by the lipidomics and the polyphenol assays. In total, 910 metabolites were found in at least 4-fold different levels in fecal matter from vegans versus omnivores, specifically for peptides, amino acids and lipids. In comparison, only 251 compounds showed 4-fold differences between lyophilized and aqueous fecal samples, including DG O-34:0 and methionine sulfoxide. A range of diet-specific metabolites were identified to be significantly different between vegans and omnivores, exemplified by citrinin and C17:0-acylcarnitine for omnivores, and curcumin and lenticin for vegans. Bioactive molecules like acyl alpha-hydroxy-fatty acids (AAHFA) were differentially regulated in vegan versus omnivore fecal materials, highlighting the importance of diet-specific reference materials for dietary biomarker studies.

Longitudinal Lipidomic Profile of Hypertension in American Indians: Findings From the Strong Heart Family Study

BACKGROUND

Dyslipidemia is an important risk factor for hypertension and cardiovascular disease. Standard lipid panel cannot reflect the complexity of blood lipidome. The associations of individual lipid species with hypertension remain to be determined in large-scale epidemiological studies, especially in a longitudinal setting.

METHODS

Using liquid chromatography-mass spectrometry, we repeatedly measured 1542 lipid species in 3699 fasting plasma samples at 2 visits (1905 at baseline, 1794 at follow-up, ~5.5 years apart) from 1905 unique American Indians in the Strong Heart Family Study. We first identified baseline lipids associated with prevalent and incident hypertension, followed by replication of top hits in Europeans. We then conducted repeated measurement analysis to examine the associations of changes in lipid species with changes in systolic blood pressure, diastolic blood pressure, and mean arterial pressure. Network analysis was performed to identify lipid networks associated with the risk of hypertension.

RESULTS

Baseline levels of multiple lipid species, for example, glycerophospholipids, cholesterol esters, sphingomyelins, glycerolipids, and fatty acids, were significantly associated with both prevalent and incident hypertension in American Indians. Some lipids were confirmed in Europeans. Longitudinal changes in multiple lipid species, for example, acylcarnitines, phosphatidylcholines, fatty acids, and triacylglycerols, were significantly associated with changes in blood pressure measurements. Network analysis identified distinct lipidomic patterns associated with the risk of hypertension.

CONCLUSIONS

Baseline plasma lipid species and their longitudinal changes are significantly associated with hypertension development in American Indians. Our findings shed light on the role of dyslipidemia in hypertension and may offer potential opportunities for risk stratification and early prediction of hypertension.

Longitudinal lipidomic signatures of all-cause and CVD mortality in American Indians: findings from the Strong Heart Study

Dyslipidemia is an independent and modifiable risk factor for aging and age-related disorders. Routine lipid panel cannot capture all individual lipid species in blood (i.e., blood lipidome). To date, a comprehensive assessment of the blood lipidome associated with mortality is lacking in large-scale community-dwelling individuals, especially in a longitudinal setting. Using liquid chromatograph-mass spectrometry, we repeatedly measured individual lipid species in 3,821 plasma samples collected at two visits (~ 5.5 years apart) from 1,930 unique American Indians in the Strong Heart Family Study. We first identified baseline lipids associated with risks for all-cause mortality and CVD mortality (mean follow-up period: 17.8 years) in American Indians, followed by replication of top hits in European Caucasians in the Malmö Diet and Cancer-Cardiovascular Cohort (n = 3,943, mean follow-up period: 23.7 years). The model adjusted age, sex, BMI, smoking, hypertension, diabetes, and LDL-c at baseline. We then examined the associations between changes in lipid species and risk of mortality. Multiple testing was controlled by false discovery rate (FDR). We found that baseline levels and longitudinal changes of multiple lipid species, e.g., cholesterol esters, glycerophospholipids, sphingomyelins, and triacylglycerols, were significantly associated with risks of all-cause or CVD mortality. Many lipids identified in American Indians could be replicated in European Caucasians. Network analysis identified differential lipid networks associated with risk of mortality. Our findings provide novel insight into the role of dyslipidemia in disease mortality and offer potential biomarkers for early prediction and risk reduction in American Indians and other ethnic groups.

The BinDiscover database: a biology-focused meta-analysis tool for 156,000 GC-TOF MS metabolome samples

Metabolomics by gas chromatography/mass spectrometry (GC/MS) provides a standardized and reliable platform for understanding small molecule biology. Since 2005, the West Coast Metabolomics Center at the University of California at Davis has collated GC/MS metabolomics data from over 156,000 samples and 2000 studies into the standardized BinBase database. We believe that the observations from these samples will provide meaningful insight to biologists and that our data treatment and webtool will provide insight to others who seek to standardize disparate metabolomics studies. We here developed an easy-to-use query interface, BinDiscover, to enable intuitive, rapid hypothesis generation for biologists based on these metabolomic samples. BinDiscover creates observation summaries and graphics across a broad range of species, organs, diseases, and compounds. Throughout the components of BinDiscover, we emphasize the use of ontologies to aggregate large groups of samples based on the proximity of their metadata within these ontologies. This adjacency allows for the simultaneous exploration of entire categories such as "rodents", "digestive tract", or "amino acids". The ontologies are particularly relevant for BinDiscover's ontologically grouped differential analysis, which, like other components of BinDiscover, creates clear graphs and summary statistics across compounds and biological metadata. We exemplify BinDiscover's extensive applicability in three showcases across biological domains.

Alternative Identification of Glycosides Using MS/MS Matching with an In Silico-Modified Aglycone Mass Spectra Library

Glycosylation of metabolites serves multiple purposes. Adding sugars makes metabolites more water soluble and improves their biodistribution, stability, and detoxification. In plants, the increase in melting points enables storing otherwise volatile compounds that are released by hydrolysis when needed. Classically, glycosylated metabolites were identified by mass spectrometry (MS/MS) using [M-sugar] neutral losses. Herein, we studied 71 pairs of glycosides with their respective aglycones, including hexose, pentose, and glucuronide moieties. Using liquid chromatography (LC) coupled to electrospray ionization high-resolution mass spectrometry, we detected the classic [M-sugar] product ions for only 68% of glycosides. Instead, we found that most aglycone MS/MS product ions were conserved in the MS/MS spectra of their corresponding glycosides, even when no [M-sugar] neutral losses were observed. We added pentose and hexose units to the precursor masses of an MS/MS library of 3057 aglycones to enable rapid identification of glycosylated natural products with standard MS/MS search algorithms. When searching unknown compounds in untargeted LC-MS/MS metabolomics data of chocolate and tea, we structurally annotated 108 novel glycosides in standard MS-DIAL data processing. We uploaded this new in silico-glycosylated product MS/MS library to GitHub to enable users to detect natural product glycosides without authentic chemical standards.

Dietary Docosahexaenoic Acid and Glucose Systemic Metabolic Changes in the Mouse

The endocannabinoid system (ECS) participates in regulating whole body energy balance. Overactivation of the ECS has been associated with the negative consequence of obesity and type 2 diabetes. Since activators of the ECS rely on lipid-derived ligands, an investigation was conducted to determine whether dietary PUFA could influence the ECS to affect glucose clearance by measuring metabolites of macronutrient metabolism. C57/blk6 mice were fed a control or DHA-enriched semi-purified diet for a period of 112 d. Plasma, skeletal muscle, and liver were collected after 56 d and 112 d of feeding the diets for metabolomics analysis. Key findings characterized a shift in glucose metabolism and greater catabolism of fatty acids in mice fed the DHA diet. Glucose use and promotion of fatty acids as substrate were found based on levels of metabolic pathway intermediates and altered metabolic changes related to pathway flux with DHA feeding. Greater levels of DHA-derived glycerol lipids were found subsequently leading to the decrease of arachidonate-derived endocannabinoids (eCB). Levels of 1- and 2-arachidonylglcerol eCB in muscle and liver were lower in the DHA diet group compared to controls. These findings demonstrate that DHA feeding in mice alters macronutrient metabolism and may restore ECS tone by lowering arachidonic acid derived eCB.

Plasma lipidomic profile of depressive symptoms: a longitudinal study in a large sample of community-dwelling American Indians in the strong heart study

Dyslipidemia has been associated with depression, but individual lipid species associated with depression remain largely unknown. The temporal relationship between lipid metabolism and the development of depression also remains to be determined. We studied 3721 fasting plasma samples from 1978 American Indians attending two exams (2001-2003, 2006-2009, mean ~5.5 years apart) in the Strong Heart Family Study. Plasma lipids were repeatedly measured by untargeted liquid chromatography-mass spectrometry (LC-MS). Depressive symptoms were assessed using the 20-item Center for Epidemiologic Studies for Depression (CES-D). Participants at risk for depression were defined as total CES-D score ≥16. Generalized estimating equation (GEE) was used to examine the associations of lipid species with incident or prevalent depression, adjusting for covariates. The associations between changes in lipids and changes in depressive symptoms were additionally adjusted for baseline lipids. We found that lower levels of sphingomyelins and glycerophospholipids and higher level of lysophospholipids were significantly associated with incident and/or prevalent depression. Changes in sphingomyelins, glycerophospholipids, acylcarnitines, fatty acids and triacylglycerols were associated with changes in depressive symptoms and other psychosomatic traits. We also identified differential lipid networks associated with risk of depression. The observed alterations in lipid metabolism may affect depression through increasing the activities of acid sphingomyelinase and phospholipase A2, disturbing neurotransmitters and membrane signaling, enhancing inflammation, oxidative stress, and lipid peroxidation, and/or affecting energy storage in lipid droplets or membrane formation. These findings illuminate the mechanisms through which dyslipidemia may contribute to depression and provide initial evidence for targeting lipid metabolism in developing preventive and therapeutic interventions for depression.

Metabolomic analysis of maternal mid-gestation plasma and cord blood in autism spectrum disorders

The discovery of prenatal and neonatal molecular biomarkers has the potential to yield insights into autism spectrum disorder (ASD) and facilitate early diagnosis. We characterized metabolomic profiles in ASD using plasma samples collected in the Norwegian Autism Birth Cohort from mothers at weeks 17-21 gestation (maternal mid-gestation, MMG, n = 408) and from children on the day of birth (cord blood, CB, n = 418). We analyzed associations using sex-stratified adjusted logistic regression models with Bayesian analyses. Chemical enrichment analyses (ChemRICH) were performed to determine altered chemical clusters. We also employed machine learning algorithms to assess the utility of metabolomics as ASD biomarkers. We identified ASD associations with a variety of chemical compounds including arachidonic acid, glutamate, and glutamine, and metabolite clusters including hydroxy eicospentaenoic acids, phosphatidylcholines, and ceramides in MMG and CB plasma that are consistent with inflammation, disruption of membrane integrity, and impaired neurotransmission and neurotoxicity. Girls with ASD have disruption of ether/non-ether phospholipid balance in the MMG plasma that is similar to that found in other neurodevelopmental disorders. ASD boys in the CB analyses had the highest number of dysregulated chemical clusters. Machine learning classifiers distinguished ASD cases from controls with area under the receiver operating characteristic (AUROC) values ranging from 0.710 to 0.853. Predictive performance was better in CB analyses than in MMG. These findings may provide new insights into the sex-specific differences in ASD and have implications for discovery of biomarkers that may enable early detection and intervention.

Zinc treatment reverses and anti-Zn-regulated miRs suppress esophageal carcinomas in vivo

Esophageal squamous cell carcinoma (ESCC) is a deadly disease with few prevention or treatment options. ESCC development in humans and rodents is associated with Zn deficiency (ZD), inflammation, and overexpression of oncogenic microRNAs: miR-31 and miR-21. In a ZD-promoted ESCC rat model with upregulation of these miRs, systemic antimiR-31 suppresses the miR-31-EGLN3/STK40-NF-κB-controlled inflammatory pathway and ESCC. In this model, systemic delivery of Zn-regulated antimiR-31, followed by antimiR-21, restored expression of tumor-suppressor proteins targeted by these specific miRs: STK40/EGLN3 (miR-31), PDCD4 (miR-21), suppressing inflammation, promoting apoptosis, and inhibiting ESCC development. Moreover, ESCC-bearing Zn-deficient (ZD) rats receiving Zn medication showed a 47% decrease in ESCC incidence vs. Zn-untreated controls. Zn treatment eliminated ESCCs by affecting a spectrum of biological processes that included downregulation of expression of the two miRs and miR-31-controlled inflammatory pathway, stimulation of miR-21-PDCD4 axis apoptosis, and reversal of the ESCC metabolome: with decrease in putrescine, increase in glucose, accompanied by downregulation of metabolite enzymes ODC and HK2. Thus, Zn treatment or miR-31/21 silencing are effective therapeutic strategies for ESCC in this rodent model and should be examined in the human counterpart exhibiting the same biological processes.

Human metabolome variation along the upper intestinal tract

Most processing of the human diet occurs in the small intestine. Metabolites in the small intestine originate from host secretions, plus the ingested exposome1 and microbial transformations. Here we probe the spatiotemporal variation of upper intestinal luminal contents during routine daily digestion in 15 healthy male and female participants. For this, we use a non-invasive, ingestible sampling device to collect and analyse 274 intestinal samples and 60 corresponding stool homogenates by combining five mass spectrometry assays2,3 and 16S rRNA sequencing. We identify 1,909 metabolites, including sulfonolipids and fatty acid esters of hydroxy fatty acids (FAHFA) lipids. We observe that stool and intestinal metabolomes differ dramatically. Food metabolites display trends in dietary biomarkers, unexpected increases in dicarboxylic acids along the intestinal tract and a positive association between luminal keto acids and fruit intake. Diet-derived and microbially linked metabolites account for the largest inter-individual differences. Notably, two individuals who had taken antibiotics within 6 months before sampling show large variation in levels of bioactive FAHFAs and sulfonolipids and other microbially related metabolites. From inter-individual variation, we identify Blautia species as a candidate to be involved in FAHFA metabolism. In conclusion, non-invasive, in vivo sampling of the human small intestine and ascending colon under physiological conditions reveals links between diet, host and microbial metabolism.

Profiling the human intestinal environment under physiological conditions

The spatiotemporal structure of the human microbiome1,2, proteome3 and metabolome4,5 reflects and determines regional intestinal physiology and may have implications for disease6. Yet, little is known about the distribution of microorganisms, their environment and their biochemical activity in the gut because of reliance on stool samples and limited access to only some regions of the gut using endoscopy in fasting or sedated individuals7. To address these deficiencies, we developed an ingestible device that collects samples from multiple regions of the human intestinal tract during normal digestion. Collection of 240 intestinal samples from 15 healthy individuals using the device and subsequent multi-omics analyses identified significant differences between bacteria, phages, host proteins and metabolites in the intestines versus stool. Certain microbial taxa were differentially enriched and prophage induction was more prevalent in the intestines than in stool. The host proteome and bile acid profiles varied along the intestines and were highly distinct from those of stool. Correlations between gradients in bile acid concentrations and microbial abundance predicted species that altered the bile acid pool through deconjugation. Furthermore, microbially conjugated bile acid concentrations exhibited amino acid-dependent trends that were not apparent in stool. Overall, non-invasive, longitudinal profiling of microorganisms, proteins and bile acids along the intestinal tract under physiological conditions can help elucidate the roles of the gut microbiome and metabolome in human physiology and disease.

Short-Term Stability of Serum and Liver Extracts for Untargeted Metabolomics and Lipidomics

Thermal reactions can significantly alter the metabolomic and lipidomic content of biofluids and tissues during storage. In this study, we investigated the stability of polar metabolites and complex lipids in dry human serum and mouse liver extracts over a three-day period under various temperature conditions. Specifically, we tested temperatures of -80 °C (freezer), -24 °C (freezer), -0.5 °C (polystyrene box with gel-based ice packs), +5 °C (refrigerator), +23 °C (laboratory, room temperature), and +30 °C (thermostat) to simulate the time between sample extraction and analysis, shipping dry extracts to different labs as an alternative to dry ice, and document the impact of higher temperatures on sample integrity. The extracts were analyzed using five fast liquid chromatography-mass spectrometry (LC-MS) methods to screen polar metabolites and complex lipids, and over 600 metabolites were annotated in serum and liver extracts. We found that storing dry extracts at -24 °C and partially at -0.5 °C provided comparable results to -80 °C (reference condition). However, increasing the storage temperatures led to significant changes in oxidized triacylglycerols, phospholipids, and fatty acids within three days. Polar metabolites were mainly affected at storage temperatures of +23 °C and +30 °C.

In utero metabolomic signatures of refined grain intake and risk of gestational diabetes: A metabolome-wide association study

BACKGROUND

Epidemiologic evidence has linked refined grain intake to a higher risk of gestational diabetes (GDM), but the biological underpinnings remain unclear.

OBJECTIVES

We aimed to identify and validate refined grain-related metabolomic biomarkers for GDM risk.

METHODS

In a metabolome-wide association study of 91 cases with GDM and 180 matched controls without GDM (discovery set) nested in the prospective Pregnancy Environment and Lifestyle Study (PETALS), refined grain intake during preconception and early pregnancy and serum untargeted metabolomics were assessed at gestational weeks 10-13. We identified refined grain-related metabolites using multivariable linear regression and examined their prospective associations with GDM risk using conditional logistic regression. We further examined the predictivity of refined grain-related metabolites selected by least absolute shrinkage and selection operator regression in the discovery set and validation set (a random PETALS subsample of 38 individuals with and 336 without GDM).

RESULTS

Among 821 annotated serum (87.4% fasting) metabolites, 42 were associated with refined grain intake, of which 17 (70.6% in glycerolipids, glycerophospholipids, and sphingolipids clusters) were associated with subsequent GDM risk (all false discovery rate-adjusted P values <0.05). Adding 7 of 17 metabolites to a conventional risk factor-based prediction model increased the C-statistic for GDM risk in the discovery set from 0.71 (95% CI: 0.64, 0.77) to 0.77 (95% CI: 0.71, 0.83) and in the validation set from 0.77 (95% CI: 0.69, 0.86) to 0.81 (95% CI: 0.74, 0.89), both with P-for-difference <0.05.

CONCLUSIONS

Clusters of glycerolipids, glycerophospholipids, and sphingolipids may be implicated in the association between refined grain intake and GDM risk, as demonstrated by the significant associations of these metabolites with both refined grains and GDM risk and the incremental predictive value of these metabolites for GDM risk beyond the conventional risk factors. These findings provide evidence on the potential biological underpinnings linking refined grain intake to the risk of GDM and help identify novel disease-related dietary biomarkers to inform diet-related preventive strategies for GDM.

Effect of pentobarbital as a euthanasia agent on equine in vitro embryo production

Postmortem and pre-euthanasia oocyte retrieval provides the last opportunity to preserve the genetic material in mares. Pentobarbital (PB) is the most common euthanasia agent; however, its effect on the developmental competence of oocytes has not been determined. Here, we evaluated the concentration of PB in equine follicular fluid (FF) and investigated its effect on the developmental competence of oocytes using a bovine IVF model to overcome the low availability of equine oocytes. The concentration of PB was measured by gas-chromatography/mass-spectrometry in FF collected from mare ovaries immediately after euthanasia (n = 10), 24 h post-euthanasia (n = 10), and from the ovaries collected by ovariectomy (negative control; n = 10). The serum concentration of PB was also evaluated as a positive control. PB was detected in all FF samples with an average concentration of 56.5 μg/ml. Next, bovine cumulus-oocyte complexes (COC) were held in holding media with PB for 6 h at 60 μg/ml (H60, n = 196), 164 μg/ml (H164, n = 215) or without PB (control; n = 212). After holding, the oocytes were matured and fertilized in vitro, followed by in vitro culture to the blastocyst stage. The cumulus expansion grade, cleavage rate, blastocyst rate, embryo kinetic rate and the blastocyst cell numbers were compared among the experimental groups of bovine COC. Higher rates of Grade 1 cumulus expansion were found in controls (54%, 32-76%; median, min-max) in comparison to H60 and H164 (24%,11-33% and 13%, 8-44%; P < 0.001). The cleavage rate was higher in the controls than in H164 (64% vs. 44%; P < 0.01). Blastocyst rates (blastocyst/cleaved oocytes) and total cell number were not different among the groups (control 29%, H60 25%, and H164 24%). In a preliminary study, equine oocytes (n = 28) were exposed to PB in vitro for 6 h followed by intracytoplasmic sperm injection (ICSI) and in vitro embryo production. Exposed oocytes showed a numerically lower maturation rate (43% Vs 52%; P > 0.05) in comparison to the laboratory-established rate during the same timepoints. Overall, we showed that PB reaches the FF immediately after euthanasia, exposing oocytes to this drug. This exposure affected cumulus expansion and cleavage rates in a bovine model, suggesting initial damage caused by PB that may not completely impede the formation of embryos, although lower overall embryo numbers might be obtained.

Cellular transformation promotes the incorporation of docosahexaenoic acid into the endolysosome-specific lipid bis(monoacylglycerol)phosphate in breast cancer

Bis(monoacylglycero)phosphates (BMPs), a class of lipids highly enriched within endolysosomal organelles, are key components of the lysosomal intraluminal vesicles responsible for activating sphingolipid catabolic enzymes. While BMPs are understudied relative to other phospholipids, recent reports associate BMP dysregulation with a variety of pathological states including neurodegenerative diseases and lysosomal storage disorders. Since the dramatic lysosomal remodeling characteristic of cellular transformation could impact BMP abundance and function, we employed untargeted lipidomics approaches to identify and quantify BMP species in several in vitro and in vivo models of breast cancer and comparative non-transformed cells and tissues. We observed lower BMP levels within transformed cells relative to normal cells, and consistent enrichment of docosahexaenoic acid (22:6) fatty acyl chain-containing BMP species in both human- and mouse-derived mammary tumorigenesis models. Our functional analysis points to a working model whereby 22:6 BMPs serve as reactive oxygen species scavengers in tumor cells, protecting lysosomes from oxidant-induced lysosomal membrane permeabilization. Our findings suggest that breast tumor cells might divert polyunsaturated fatty acids into BMP lipids as part of an adaptive response to protect their lysosomes from elevated reactive oxygen species levels, and raise the possibility that BMP-mediated lysosomal protection is a tumor-specific vulnerability that may be exploited therapeutically.

Genome-wide analysis of oxylipins and oxylipin profiles in a pediatric population

Background

Oxylipins are inflammatory biomarkers derived from omega-3 and-6 fatty acids implicated in inflammatory diseases but have not been studied in a genome-wide association study (GWAS). The aim of this study was to identify genetic loci associated with oxylipins and oxylipin profiles to identify biologic pathways and therapeutic targets for oxylipins.

Methods

We conducted a GWAS of plasma oxylipins in 316 participants in the Diabetes Autoimmunity Study in the Young (DAISY). DNA samples were genotyped using the TEDDY-T1D Exome array, and additional variants were imputed using the Trans-Omics for Precision Medicine (TOPMed) multi-ancestry reference panel. Principal components analysis of 36 plasma oxylipins was used to capture oxylipin profiles. PC1 represented linoleic acid (LA)- and alpha-linolenic acid (ALA)-related oxylipins, and PC2 represented arachidonic acid (ARA)-related oxylipins. Oxylipin PC1, PC2, and the top five loading oxylipins from each PC were used as outcomes in the GWAS (genome-wide significance: p < 5×10-8).

Results

The SNP rs143070873 was associated with (p < 5×10-8) the LA-related oxylipin 9-HODE, and rs6444933 (downstream of CLDN11) was associated with the LA-related oxylipin 13 S-HODE. A locus between MIR1302-7 and LOC100131146, rs10118380 and an intronic variant in TRPM3 were associated with the ARA-related oxylipin 11-HETE. These loci are involved in inflammatory signaling cascades and interact with PLA2, an initial step to oxylipin biosynthesis.

Conclusion

Genetic loci involved in inflammation and oxylipin metabolism are associated with oxylipin levels.