Search

Advances in omics technologies have allowed spatially resolved molecular profiling of single cells, providing a window not only into the diversity and distribution of cell types within a tissue, but also into the effects of interactions between cells in shaping the transcriptional landscape. Cells send chemical and mechanical signals which are received by other cells, where they can subsequently initiate context-specific gene regulatory responses. These interactions and their responses shape the individual molecular phenotype of a cell in a given microenvironment. RNAs or proteins measured in individual cells, together with the cells' spatial distribution, provide invaluable information about these mechanisms and the regulation of genes beyond processes occurring independently in each individual cell. "SpaCeNet" is a method designed to elucidate both the intracellular molecular networks (how molecular variables affect each other within the cell) and the intercellular molecular networks (how cells affect molecular variables in their neighbors). This is achieved by estimating conditional independence (CI) relations between captured variables within individual cells and by disentangling these from CI relations between variables of different cells., (© 2024 Schrod et al.; Published by Cold Spring Harbor Laboratory Press.)

Gene regulatory networks (GRNs) are effective tools for inferring complex interactions between molecules that regulate biological processes and hence can provide insights into drivers of biological systems. Inferring coexpression networks is a critical element of GRN inference, as the correlation between expression patterns may indicate that genes are coregulated by common factors. However, methods that estimate coexpression networks generally derive an aggregate network representing the mean regulatory properties of the population and so fail to fully capture population heterogeneity. Bayesian optimized networks obtained by assimilating omic data (BONOBO) is a scalable Bayesian model for deriving individual sample-specific coexpression matrices that recognizes variations in molecular interactions across individuals. For each sample, BONOBO assumes a Gaussian distribution on the log-transformed centered gene expression and a conjugate prior distribution on the sample-specific coexpression matrix constructed from all other samples in the data. Combining the sample-specific gene coexpression with the prior distribution, BONOBO yields a closed-form solution for the posterior distribution of the sample-specific coexpression matrices, thus allowing the analysis of large data sets. We demonstrate BONOBO's utility in several contexts, including analyzing gene regulation in yeast transcription factor knockout studies, the prognostic significance of miRNA-mRNA interaction in human breast cancer subtypes, and sex differences in gene regulation within human thyroid tissue. We find that BONOBO outperforms other methods that have been used for sample-specific coexpression network inference and provides insight into individual differences in the drivers of biological processes., (© 2024 Saha et al.; Published by Cold Spring Harbor Laboratory Press.)

There is increasing recognition that the sex chromosomes, X and Y, play an important role in health and disease that goes beyond the determination of biological sex. Loss of the Y chromosome (LOY) in blood, which occurs naturally in aging men, has been found to be a driver of cardiac fibrosis and heart failure mortality. LOY also occurs in most solid tumors in males and is often associated with worse survival, suggesting that LOY may give tumor cells a growth or survival advantage. We analyzed LOY in lung adenocarcinoma (LUAD) using both bulk and single-cell expression data and found evidence suggesting that LOY affects the tumor immune environment by altering cancer/testis antigen expression and consequently facilitating tumor immune evasion. Analyzing immunotherapy data, we show that LOY and changes in expression of particular cancer/testis antigens are associated with response to pembrolizumab treatment and outcome, providing a new and powerful biomarker for predicting immunotherapy response in LUAD tumors in males.

Background: Posttraumatic stress disorder (PTSD) is characterized by severe distress and associated with cardiometabolic diseases. Studies in military and clinical populations suggest dysregulated metabolomic processes may be a key mechanism. Prior work identified and validated a metabolite-based distress score (MDS) linked with depression and anxiety and subsequent cardiometabolic diseases. Here, we assessed whether PTSD shares metabolic alterations with depression and anxiety and also if additional metabolites are related to PTSD., Methods: We leveraged plasma metabolomics data from three subsamples nested within the Nurses' Health Study II, including 2835 women with 2950 blood samples collected across three timepoints (1996-2014) and 339 known metabolites consistently assayed by mass spectrometrybased techniques. Trauma and PTSD exposures were assessed in 2008 and characterized as follows: lifetime trauma without PTSD, lifetime PTSD in remission, and persistent PTSD symptoms. Associations between the exposures and the MDS or individual metabolites were estimated within each subsample adjusting for potential confounders and combined in random-effects meta-analyses., Results: Persistent PTSD symptoms were associated with higher levels of the previously developed MDS for depression and anxiety. Out of 339 metabolites, we identified nine metabolites (primarily elevated glycerophospholipids) associated with persistent symptoms (false discovery rate<0.05). No metabolite associations were found with the other PTSD-related exposures., Conclusions: As the first large-scale, population-based metabolomics analysis of PTSD, our study highlighted shared and distinct metabolic differences linked to PTSD versus depression or anxiety. We identified novel metabolite markers associated with PTSD symptom persistence, suggesting further connections with metabolic dysregulation that may have downstream consequences for health., Competing Interests: Conflict of interest: None declared.

Background: Lung adenocarcinoma (LUAD) has been observed to have significant sex differences in incidence, prognosis, and response to therapy. However, the molecular mechanisms responsible for these disparities have not been investigated extensively., Methods: Sample-specific gene regulatory network methods were used to analyze RNA sequencing data from non-cancerous human lung samples from The Genotype Tissue Expression Project (GTEx) and lung adenocarcinoma primary tumor samples from The Cancer Genome Atlas (TCGA); results were validated on independent data., Results: We found that genes associated with key biological pathways including cell proliferation, immune response and drug metabolism are differentially regulated between males and females in both healthy lung tissue and tumor, and that these regulatory differences are further perturbed by tobacco smoking. We also discovered significant sex bias in transcription factor targeting patterns of clinically actionable oncogenes and tumor suppressor genes, including AKT2 and KRAS. Using differentially regulated genes between healthy and tumor samples in conjunction with a drug repurposing tool, we identified several small-molecule drugs that might have sex-biased efficacy as cancer therapeutics and further validated this observation using an independent cell line database., Conclusions: These findings underscore the importance of including sex as a biological variable and considering gene regulatory processes in developing strategies for disease prevention and management., (© 2024. The Author(s).)

Compared to men, women often develop COPD at an earlier age with worse respiratory symptoms despite lower smoking exposure. However, most preventive, and therapeutic strategies ignore biological sex differences in COPD. Our goal was to better understand sex-specific gene regulatory processes in lung tissue and the molecular basis for sex differences in COPD onset and severity. We analyzed lung tissue gene expression and DNA methylation data from 747 individuals in the Lung Tissue Research Consortium (LTRC), and 85 individuals in an independent dataset. We identified sex differences in COPD-associated gene regulation using gene regulatory networks. We used linear regression to test for sex-biased associations of methylation with lung function, emphysema, smoking, and age. Analyzing gene regulatory networks in the control group, we identified that genes involved in the extracellular matrix (ECM) have higher transcriptional factor targeting in females than in males. However, this pattern is reversed in COPD, with males showing stronger regulatory targeting of ECM-related genes than females. Smoking exposure, age, lung function, and emphysema were all associated with sex-specific differential methylation of ECM-related genes. We identified sex-based gene regulatory patterns of ECM-related genes associated with lung function and emphysema. Multiple factors including epigenetics, smoking, aging, and cell heterogeneity influence sex-specific gene regulation in COPD. Our findings underscore the importance of considering sex as a key factor in disease susceptibility and severity.

Aging is the primary risk factor for many individual cancer types, including lung adenocarcinoma (LUAD). To understand how aging-related alterations in the regulation of key cellular processes might affect LUAD risk and survival outcomes, we built individual (person)-specific gene regulatory networks integrating gene expression, transcription factor protein-protein interaction, and sequence motif data, using PANDA/LIONESS algorithms, for both non-cancerous lung tissue samples from the Genotype Tissue Expression (GTEx) project and LUAD samples from The Cancer Genome Atlas (TCGA). In GTEx, we found that pathways involved in cell proliferation and immune response are increasingly targeted by regulatory transcription factors with age; these aging-associated alterations are accelerated by tobacco smoking and resemble oncogenic shifts in the regulatory landscape observed in LUAD and suggests that dysregulation of aging pathways might be associated with an increased risk of LUAD. Comparing normal adjacent samples from individuals with LUAD with healthy lung tissue samples from those without LUAD, we found that aging-associated genes show greater aging-biased targeting patterns in younger individuals with LUAD compared to their healthy counterparts of similar age, a pattern suggestive of age acceleration. This implies that an accelerated aging process may be responsible for tumor incidence in younger individuals. Using drug repurposing tool CLUEreg, we found small molecule drugs with potential geroprotective effects that may alter the accelerating aging profiles we found. We also observed that, in contrast to chronological age, a network-informed aging signature was associated with survival and response to chemotherapy in LUAD., Competing Interests: Declaration of interests The authors declare no competing interests.

Context: Psychological distress has been linked to diabetes risk. Few population-based, epidemiologic studies have investigated the potential molecular mechanisms (eg, metabolic dysregulation) underlying this association., Objective: To evaluate the association between a metabolomic signature for psychological distress and diabetes risk., Methods: We conducted a nested case-control study of plasma metabolomics and diabetes risk in the Nurses' Health Study, including 728 women (mean age: 55.2 years) with incident diabetes and 728 matched controls. Blood samples were collected between 1989 and 1990 and incident diabetes was diagnosed between 1992 and 2008. Based on our prior work, we calculated a weighted plasma metabolite-based distress score (MDS) comprised of 19 metabolites. We used conditional logistic regression accounting for matching factors and other diabetes risk factors to estimate odds ratios (OR) and 95% confidence intervals (CI) for diabetes risk according to MDS., Results: After adjusting for sociodemographic factors, family history of diabetes, and health behaviors, the OR (95% CI) for diabetes risk across quintiles of the MDS was 1.00 (reference) for Q1, 1.16 (0.77, 1.73) for Q2, 1.30 (0.88, 1.91) for Q3, 1.99 (1.36, 2.92) for Q4, and 2.47 (1.66, 3.67) for Q5. Each SD increase in MDS was associated with 36% higher diabetes risk (95% CI: 1.21, 1.54; P-trend <.0001). This association was moderately attenuated after additional adjustment for body mass index (comparable OR: 1.17; 95% CI: 1.02, 1.35; P-trend = .02). The MDS explained 17.6% of the association between self-reported psychological distress (defined as presence of depression or anxiety symptoms) and diabetes risk (P = .04)., Conclusion: MDS was significantly associated with diabetes risk in women. These results suggest that differences in multiple lipid and amino acid metabolites may underlie the observed association between psychological distress and diabetes risk., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Endocrine Society. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: Chronic psychological distress is associated with increased risk of cardiovascular disease (CVD) and investigators have posited inflammatory factors may be centrally involved in these relationships. However, mechanistic evidence and molecular underpinnings of these processes remain unclear, and data are particularly sparse among women. This study examined if a metabolite profile linked with distress was associated with increased CVD risk and inflammation-related risk factors., Methods: A plasma metabolite-based distress score (MDS) of twenty chronic psychological distress-related metabolites was developed in cross-sectional, 1:1 matched case-control data comprised of 558 women from the Nurses' Health Study (NHS; 279 women with distress, 279 controls). This MDS was then evaluated in two other cohorts: the Women's Health Initiative Observational Cohort (WHI-OS) and the Prevención con Dieta Mediterránea (PREDIMED) trial. We tested the MDS's association with risk of future CVD in each sample and with levels of C-reactive protein (CRP) in the WHI-OS. The WHI-OS subsample included 944 postmenopausal women (472 CHD cases; mean time to event = 5.8 years); the PREDIMED subsample included 980 men and women (224 CVD cases, mean time to event = 3.1 years)., Results: In the WHI-OS, a 1-SD increase in the plasma MDS was associated with a 20% increased incident CHD risk (odds ratio [OR] = 1.20, 95% CI: 1.04 - 1.38), adjusting for known CVD risk factors excluding total and HDL cholesterol. This association was attenuated after including total and HDL cholesterol. CRP mediated an average 12.9% (95% CI: 4.9% - 28%, p < 10 -15 ) of the total effect of MDS on CHD risk when adjusting for matching factors. This effect was attenuated after adjusting for known CVD risk factors. Of the metabolites in the MDS, tryptophan and threonine were inversely associated with incident CHD risk in univariate models. In PREDIMED, each one SD increase in the MDS was associated with an OR of 1.19 (95% CI: 1.00 - 1.41) for incident CVD risk, after adjusting all risk factors. Similar associations were observed in men and women. Four metabolites in the MDS were associated with incident CVD risk in PREDIMED in univariate models. Biliverdin and C36:5 phosphatidylcholine (PC) plasmalogen had inverse associations; C16:0 ceramide and C18:0 lysophosphatidylethanolamine(LPE) each had positive associations with CVD risk., Conclusions: Our study points to molecular alterations that may underlie the association between chronic distress and subsequent risk of cardiovascular disease in adults., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Metabolomic epidemiology is the high-throughput study of the relationship between metabolites and health-related traits. This emerging and rapidly growing field has improved our understanding of disease aetiology and contributed to advances in precision medicine. As the field continues to develop, metabolomic epidemiology could lead to the discovery of diagnostic biomarkers predictive of disease risk, aiding in earlier disease detection and better prognosis. In this Review, we discuss key advances facilitated by the field of metabolomic epidemiology for a range of conditions, including cardiometabolic diseases, cancer, Alzheimer's disease and COVID-19, with a focus on potential clinical utility. Core principles in metabolomic epidemiology, including study design, causal inference methods and multi-omic integration, are briefly discussed. Future directions required for clinical translation of metabolomic epidemiology findings are summarized, emphasizing public health implications. Further work is needed to establish which metabolites reproducibly improve clinical risk prediction in diverse populations and are causally related to disease progression., (© 2023. Springer Nature Limited.)

Metabolomics holds great promise for uncovering insights around biological processes impacting disease in human epidemiological studies. Metabolites can be measured across biological samples, including plasma, serum, saliva, urine, stool, and whole organs and tissues, offering a means to characterize metabolic processes relevant to disease etiology and traits of interest. Metabolomic epidemiology studies face unique challenges, such as identifying metabolites from targeted and untargeted assays, defining standards for quality control, harmonizing results across platforms that often capture different metabolites, and developing statistical methods for high-dimensional and correlated metabolomic data. In this review, we introduce metabolomic epidemiology to the broader scientific community, discuss opportunities and challenges presented by these studies, and highlight emerging innovations that hold promise to uncover new biological insights., Competing Interests: Declaration of interests None declared by authors., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Gaussian graphical models (GGMs) are a popular form of network model in which nodes represent features in multivariate normal data and edges reflect conditional dependencies between these features. GGM estimation is an active area of research. Currently available tools for GGM estimation require investigators to make several choices regarding algorithms, scoring criteria, and tuning parameters. An estimated GGM may be highly sensitive to these choices, and the accuracy of each method can vary based on structural characteristics of the network such as topology, degree distribution, and density. Because these characteristics are a priori unknown, it is not straightforward to establish universal guidelines for choosing a GGM estimation method. We address this problem by introducing SpiderLearner, an ensemble method that constructs a consensus network from multiple estimated GGMs. Given a set of candidate methods, SpiderLearner estimates the optimal convex combination of results from each method using a likelihood-based loss function. K $$ K $$ -fold cross-validation is applied in this process, reducing the risk of overfitting. In simulations, SpiderLearner performs better than or comparably to the best candidate methods according to a variety of metrics, including relative Frobenius norm and out-of-sample likelihood. We apply SpiderLearner to publicly available ovarian cancer gene expression data including 2013 participants from 13 diverse studies, demonstrating our tool's potential to identify biomarkers of complex disease. SpiderLearner is implemented as flexible, extensible, open-source code in the R package ensembleGGM at https://github.com/katehoffshutta/ensembleGGM., (© 2023 John Wiley & Sons Ltd.)

Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods., (© 2023. The Author(s).)

The increasing quantity of multi-omic data, such as methylomic and transcriptomic profiles collected on the same specimen or even on the same cell, provides a unique opportunity to explore the complex interactions that define cell phenotype and govern cellular responses to perturbations. We propose a network approach based on Gaussian Graphical Models (GGMs) that facilitates the joint analysis of paired omics data. This method, called DRAGON (Determining Regulatory Associations using Graphical models on multi-Omic Networks), calibrates its parameters to achieve an optimal trade-off between the network's complexity and estimation accuracy, while explicitly accounting for the characteristics of each of the assessed omics 'layers.' In simulation studies, we show that DRAGON adapts to edge density and feature size differences between omics layers, improving model inference and edge recovery compared to state-of-the-art methods. We further demonstrate in an analysis of joint transcriptome - methylome data from TCGA breast cancer specimens that DRAGON can identify key molecular mechanisms such as gene regulation via promoter methylation. In particular, we identify Transcription Factor AP-2 Beta (TFAP2B) as a potential multi-omic biomarker for basal-type breast cancer. DRAGON is available as open-source code in Python through the Network Zoo package (netZooPy v0.8; netzoo.github.io)., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Psychological distress can be conceptualized as an umbrella term encompassing symptoms of depression, anxiety, posttraumatic stress disorder (PTSD), or stress more generally. A systematic review of metabolomic markers associated with distress has the potential to reveal underlying molecular mechanisms linking distress to adverse health outcomes. The current systematic review extends prior reviews of clinical depressive disorders by synthesizing 39 existing studies that examined metabolomic markers for PTSD, anxiety disorders, and subclinical psychological distress in biological specimens. Most studies were based on small sets of pre-selected candidate metabolites, with few metabolites overlapping between studies. Vast heterogeneity was observed in study design and inconsistent patterns of association emerged between distress and metabolites. To gain a more robust understanding of distress and its metabolomic signatures, future research should include 1) large, population-based samples and longitudinal assessments, 2) replication and validation in diverse populations, 3) and agnostic metabolomic strategies profiling hundreds of targeted and nontargeted metabolites. Addressing these research priorities will improve the scope and reproducibility of future metabolomic studies of psychological distress., Competing Interests: Conflict of interest None., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Gaussian graphical models (GGMs) provide a framework for modeling conditional dependencies in multivariate data. In this tutorial, we provide an overview of GGM theory and a demonstration of various GGM tools in R. The mathematical foundations of GGMs are introduced with the goal of enabling the researcher to draw practical conclusions by interpreting model results. Background literature is presented, emphasizing methods recently developed for high-dimensional applications such as genomics, proteomics, or metabolomics. The application of these methods is illustrated using a publicly available dataset of gene expression profiles from 578 participants with ovarian cancer in The Cancer Genome Atlas. Stand-alone code for the demonstration is available as an RMarkdown file at https://github.com/katehoffshutta/ggmTutorial., (© 2022 John Wiley & Sons Ltd.)

Background and Objectives: Women have higher lifetime risk of stroke than men, and metabolic factors seem more strongly associated with stroke for women than men. However, few studies in either men or women have evaluated metabolomic profiles and incident stroke., Methods: We applied liquid chromatography-tandem mass spectrometry to measure 519 plasma metabolites in a discovery set of women in the Nurses' Health Study (NHS; 454 incident ischemic stroke cases, 454 controls) with validation in 2 independent, prospective cohorts: Prevención con Dieta Mediterránea (PREDIMED; 118 stroke cases, 791 controls) and Nurses' Health Study 2 (NHS2; 49 ischemic stroke cases, 49 controls). We applied logistic regression models with stroke as the outcome to adjust for multiple risk factors; the false discovery rate was controlled through the q value method., Results: Twenty-three metabolites were significantly associated with incident stroke in NHS after adjustment for traditional risk factors ( q < 0.05). Of these, 14 metabolites were available in PREDIMED and 3 were significantly associated with incident stroke: methionine sulfoxide, N6-acetyllysine, and sucrose ( q < 0.05). In NHS2, one of the 23 metabolites (glucuronate) was significantly associated with incident stroke ( q < 0.05). For all 4 metabolites, higher levels were associated with increased risk. These 4 metabolites were used to create a stroke metabolite score (SMS) in the NHS and tested in PREDIMED. Per unit of standard deviation of SMS, the odds ratio for incident stroke was 4.12 (95% confidence interval [CI] 2.26-7.51) in PREDIMED, after adjustment for risk factors. In PREDIMED, the area under the receiver operating characteristic curve (AUC) for the model including SMS and traditional risk factors was 0.70 (95% CI 0.75-0.79) vs the AUC for the model including the traditional risk factors only of 0.65 (95% CI 0.70-0.75), corresponding to a 5% improvement in risk prediction with SMS ( p < 0.005)., Discussion: Metabolites associated with stroke included 2 amino acids, a carboxylic acid, and sucrose. A composite SMS including these metabolites was associated with ischemic stroke and showed improvement in risk prediction beyond traditional risk factors., Classification of Evidence: This study provides Class II evidence that a SMS accurately predicts incident ischemic stroke risk., (© 2021 American Academy of Neurology.)

Background: Higher circulating carotenoids are associated with lower breast cancer risk. The underlying biology remains under-explored., Methods: We profiled 293 prediagnostic plasma metabolites in a nested case-control study ( n = 887 cases) within the Nurses' Health Studies. Associations between circulating carotenoids and metabolites were identified using linear-mixed models (FDR ≤ 0.05), and we further selected metabolites most predictive of carotenoids with LASSO. Metabolic signatures for carotenoids were calculated as weighted sums of LASSO selected metabolites. We further evaluated the metabolic signatures in relation to breast cancer risk using conditional logistic-regression., Results: We identified 48 to 110 metabolites associated with plasma levels of α-carotene, β-carotene, β-cryptoxanthin, estimated-vitamin-A-potential, lutein/zeaxanthin, and lycopene, which included primarily positively associated metabolites implicated in immune regulation (tryptophan), redox balance (plasmalogens, glutamine), epigenetic regulations (acetylated-/methylated-metabolites), and primarily inversely associated metabolites involved in β-oxidation (carnitines; FDR ≤ 0.05). The metabolomic signatures derived for β-carotene (Q4 vs. Q1 relative risk RR = 0.74, P trend = 0.02), and estimated-vitamin-A-potential (Q4 vs. Q1 RR = 0.74, P trend = 0.02)-measured ≥10 years before diagnosis-were associated with lower breast cancer risk. Modest attenuations of RR for measured levels of β-carotene and estimated-vitamin-A-potential were seen when we adjusted for their corresponding metabolic signatures., Conclusions: Metabolites involved in immune regulation, redox balance, membrane signaling, and β-oxidation were associated with plasma carotenoids. Although some metabolites may reflect shared common food sources or compartmental colocalization with carotenoids, others may signal the underlying pathways of carotenoids-associated lowered breast cancer risk., Impact: Consumption of carotenoid-rich diet is associated with a wide-range of metabolic changes which may help to reduce breast cancer risk., (©2021 American Association for Cancer Research.)

Several forms of chronic distress including anxiety and depression are associated with adverse cardiometabolic outcomes. Metabolic alterations may underlie these associations. Whether these forms of distress are associated with metabolic alterations even after accounting for comorbid conditions and other factors remains unclear. Using an agnostic approach, this study examines a broad range of metabolites in relation to chronic distress among women. For this cross-sectional study of chronic distress and 577 plasma metabolites, data are from different substudies within the Women's Health Initiative (WHI) and Nurses' Health Studies (NHSI, NHSII). Chronic distress was characterized by depressive symptoms and other depression indicators in the WHI and NHSII substudies, and by combined indicators of anxiety and depressive symptoms in the NHSI substudy. We used a two-phase discovery-validation framework, with WHI (N = 1317) and NHSII (N = 218) substudies in the discovery phase (identifying metabolites associated with distress) and NHSI (N = 558) substudy in the validation phase. A differential network analysis provided a systems-level assessment of metabolomic alterations under chronic distress. Analyses adjusted for potential confounders and mediators (demographics, comorbidities, medications, lifestyle factors). In the discovery phase, 46 metabolites were significantly associated with depression measures. In validation, six of these metabolites demonstrated significant associations with chronic distress after adjustment for potential confounders. Among women with high distress, we found lower gamma-aminobutyric acid (GABA), threonine, biliverdin, and serotonin and higher C16:0 ceramide and 3-methylxanthine. Our findings suggest chronic distress is associated with metabolomic alterations and provide specific targets for future study of biological pathways in chronic diseases., (Copyright © 2021 Elsevier Ltd. All rights reserved.)