- 1Medical Scientist Training Program, University of Alabama at Birmingham, Birmingham, AL, United States
- 2Department of Epidemiology, University of Alabama at Birmingham, Birmingham, AL, United States
- 3Department of Neurology, Massachusetts General Hospital, Boston, MA, United States
- 4Department of Biostatistics, University of Alabama at Birmingham, Birmingham, AL, United States
- 5Department of Epidemiology, University of Texas Health Science Center, Houston, TX, United States
- 6Division of Neurology, Department of Medicine and Division of Academic Affairs, Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand
- 7HudsonAlpha Institute for Biotechnology, Huntsville, AL, United States
- 8Office of the Provost, University of South Carolina, Columbia, SC, United States
Introduction: Metabolic syndrome (MetS) increases the risk of cardiovascular disease and death. Previous ‘-omics’ studies have identified dysregulated serum metabolites and aberrant DNA methylation in the setting of MetS. However, the relationship between the metabolome and epigenome have not been elucidated. In this study, we identified serum metabolites associated with MetS and DNA methylation, and we conducted bidirectional Mendelian randomization (MR) to assess causal relationships between metabolites and methylation.
Methods: We leveraged metabolomic and genomic data from a national United States cohort of older adults (REGARDS), as well as metabolomic, epigenomic, and genomic data from a family-based study of hypertension (HyperGEN). We conducted metabolite profiling for MetS in REGARDS using weighted logistic regression models and validated them in HyperGEN. Validated metabolites were selected for methylation studies which fit linear mixed models between metabolites and six CpG sites previously linked to MetS. Statistically significant metabolite-CpG pairs were selected for two-sample, bidirectional MR.
Results: Forward MR indicated that glucose and serine metabolites were causal on CpG methylation near CPT1A [B(SE): −0.003 (0.002), p = 0.028 and B(SE): 0.029 (0.011), p = 0.030, respectively] and that serine metabolites were causal on ABCG1 [B(SE): −0.008(0.003), p = 0.006] and SREBF1 [B(SE): −0.009(0.004), p = 0.018] methylation, which suggested a protective effect of serine. Reverse MR showed a bidirectional relationship between cg06500161 (ABCG1) and serine [B(SE): −1.534 (0.668), p = 0.023].
Discussion: The metabolome may contribute to the relationship between MetS and epigenetic modifications.
Introduction
Metabolic syndrome (MetS), defined by a cluster of cardiometabolic risk factors (hypertension, elevated blood glucose, abnormal lipid levels, and abdominal obesity), is described as “more than the sum of its parts” (Reilly and Rader, 2003; Huang, 2009). In the last 4 decades, the prevalence of MetS has risen by 35% overall, with consistent increases across all sociodemographic groups; further, more than one-third of all United States adults has MetS (Moore et al., 2017). MetS increases the risk of cardiovascular diseases more than the risk of these traits individually (Fahed et al., 2022). Previous metabolomics studies have discovered dysregulated serum metabolite concentrations in the setting of MetS, but mechanisms of how the metabolome confers MetS risk are unclear (Roberts et al., 2014; Azab et al., 2021; Carioca et al., 2021; Comte et al., 2021; Jialal et al., 2021; Wu et al., 2021; Zhang et al., 2021; Lind et al., 2022). Similarly, epigenome-wide association studies (EWAS) have identified, with consistent validation, alterations in DNA methylation of cytosine-phosphate-guanine (CpG) sites near genes that contribute to lipid metabolism, inflammation pathways, and hormone signaling in MetS. Importantly, the epigenome is sensitive to the environment (which may include the metabolome). Few studies have integrated the metabolome with MetS-related epigenomic sites, which could provide additional information on how the metabolome may alter individual risk for MetS.
Longitudinal data on ‘-omics’ are generally unavailable to understand causal relationships, particularly between -omic measurements (e.g., epigenomics and metabolomics). When only cross-sectional data are available, Mendelian randomization (MR) analyses provide an alternative. Importantly, past MR analyses have suggested that cardiometabolic traits are causal on methylation changes (Akinyemiju et al., 2018; Hidalgo et al., 2021; Jones et al., 2021). Also, metabolite reactions, such as the methionine cycle, may contribute to DNA methylation as methyl donors via DNA methyltransferase reactions (Zhang, 2018; Maddocks et al., 2016; Y et al., 2016). Yet, pathways of how MetS-related traits drive methylation variations have not been fully elucidated (Zaghlool et al., 2018). Methylation patterns associated with MetS that have been observed via EWAS could be downstream of dysregulated metabolite levels. Therefore, understanding the relationship between the metabolome and epigenome may help identify potential mechanisms of DNA methylation in cardiometabolic dysfunction. The objective of this study was to utilize genomic and metabolomic data from the Reasons for Geographic and Racial Differences in Stroke (REGARDS) Study, as well as genomic, epigenomic, and metabolomic data from the Hypertension Genetic Epidemiology Network (HyperGEN) to, first, conduct a multi-omics analysis of MetS and second, establish the directionality of the observed associations through bidirectional MR analyses.
Materials and methods
Study populations
We leveraged data from a subset of participants from the REGARDS Study. The REGARDS objectives and design have previously been described (Howard et al., 2005). Briefly, 30,239 community-dwelling, non-Hispanic Black and White adults aged 45 and older were recruited between 2003 and 2007 across the contiguous United States with oversampling from the Southeast region. Metabolomics were conducted in a subset of individuals who were part of a case-cohort study for incident strokes (N = 2,165) (Ament et al., 2022). A subset of these participants also had available genotype data (N = 1,865). Metabolite profiling, genotyping, and quality control metrics are detailed in the Supplementary Methods.
We also utilized available genomic (N = 1,898), epigenomic (N = 614), and metabolomic data (N = 300) from the Hypertension Genetic Epidemiology Network (HyperGEN). Recruitment methods for HyperGEN have been previously described (Williams et al., 2000). Briefly, HyperGEN—coordinated under the National Heart, Lung, and Blood Institute (NHLBI) Family Blood Pressure Program (FBPP)—was enriched for hypertensive adults and their siblings and/or adult offspring to evaluate genetic contributors to hypertension and target organ damage. Later, African American (AA) participants’ samples underwent deep coverage, whole genome sequencing (WGS) through the NLHBI Trans-Omics in Precision Medicine (TOPMed) Consortium. In separate ancillary studies, AA participant samples were selected for epigenotyping and metabolite measurements for studies of left ventricular hypertrophy (under an extreme phenotype design) and chronic kidney disease, respectively. Complete descriptions of data collection in HyperGEN are available in the Supplementary Methods. The different datasets used for this analysis are summarized in Supplementary Figure S1.
MetS phenotyping
MetS was defined as having at least 3 of the following at the baseline visit according the United States National Cholesterol Education Program Adult Treatment Panel III (NCEP-ATP III) 2005 criteria: hypertension (SBP ≥ 130 mmHg or DBP ≥ 85 or self-reported use of anti-hypertensive medications); hypertriglyceridemia (triglyceride ≥ 150 mg/dL or self-reported use of lipid-lowering medication); hypoalphalipoproteinemia (HDL cholesterol <40 mg/dL for males or <50 mg/dL for females or self-reported use of lipid-lowering medication); elevated blood sugar (fasting blood glucose ≥ 100 mg/dL) or self-reported use of glucose-lowering medication or insulin); and abdominal obesity (waist circumference >40 cm for males or >35 cm for females (Huang, 2009).
Statistical analyses
MetS-metabolite association and validation
All analyses are summarized in Figure 1. In the REGARDS case-cohort (N = 2,039), weighted logistic regression models were fit to test the association between prevalent MetS (outcome) and a panel of 162 serum metabolites (exposure). Models were adjusted for age, gender, race, cigarette smoking, alcohol use, and stroke case status. In sensitivity analyses, we also adjusted models for fasting status and restricted the analysis to the randomly selected controls (N = 872). Metabolites with association p-values less than the Bonferroni-corrected significance threshold of p = 0.05/162 = 3.09e-04 were selected for validation in HyperGEN participants with metabolomics data (N = 300). We also selected 12 metabolites that were upstream of DNA methylation pathways irrespective of statistical significance in REGARDS for validation. HyperGEN metabolite models were similarly adjusted for age, gender, smoking status, alcohol use, and recruitment center. Race was not included as a covariate, as all HyperGEN metabolomics participants were AA. Validated MetS metabolites (p <0.05) were selected for further analyses among HyperGEN participants who also had both metabolomic and methylation data (N = 134); none of these individuals were related. Logistic regression modeling was completed in SAS version 9.4.
Methylation studies
We evaluated the relationship between HyperGEN-validated metabolites and six CpG sites that have previously been linked to MetS (Table 1): cg06500161, annotated to ATP-Binding Cassette Subfamily G Member 1 (ABCG1); cg00574958, annotated to Carnitine Palmitoyltransferase 1A (CPT1A); cg02650017, annotated to Phosphoethanolamine/Phosphocholine Phosphatase 1 (PHOSPHO1); cg11024682, annotated to Sterol Regulatory Element-Binding Protein 1 (SREBF1); cg18181703, annotated to Suppressor Of Cytokine Signaling 3 (SOCS3); and cg19693031, annotated to Thioredoxin-Interacting Protein (TXNIP) (Hidalgo et al., 2021). We fit linear mixed models for the CpG beta score (outcome) and metabolite (predictor), adjusted for age, gender, recruitment center, left ventricular mass index, principal components of ancestry (PCs), Houseman-estimated cell counts, and family relatedness (random effect). For the linear mixed models, we used the lme4 package in R (version 4.2.0). CpG-metabolite pairs with association p <0.05 were eligible for MR studies.
Two-sample, bidirectional Mendelian randomization
We conducted bidirectional MR using the TwoSampleMR package (version 0.5.6) in R (version 4.2.0) for statistically significant CpG-metabolite associations (Hemani et al., 2018). We hypothesized that metabolite levels were causal on CpG methylation (Figure 2). This MR approach uses single nucleotide polymorphisms (SNPs) associated with the trait of interest as the instrumental variables (IVs). Under the assumptions of the MR, if the SNPs that are strongly associated with the exposure are also associated with the outcome, then there is evidence to suggest a causal relationship. To obtain the exposure instruments (metabolite quantitative trait loci SNPs or mQTL SNPs), we conducted genome-wide association studies (GWAS) in PLINK 2.0 among participants with both metabolite and genetic data in REGARDS for selected metabolites, adjusted for age, sex, and the first ten PCs. We similarly conducted GWAS for selected CpG sites and corresponding WGS data in HyperGEN (N = 557) to obtain the outcome instruments (methylation quantitative loci SNPs or meQTL SNPs). HyperGEN meQTL models were adjusted for left ventricular mass index. Because the HyperGEN sample was all African American, we restricted the REGARDS GWAS to African American participants (N = 846) for SNP IV selection.
There were few SNPs in REGARDS with association p <5.00e-08, which corresponds to the traditional Bonferroni-corrected threshold for GWAS; thus, we considered a lower p-value threshold in order to increase the available IVs for forward MR analyses. Summary statistics for SNPs with association p <1.00e-06 in REGARDS metabolite GWAS (exposure IVs) and their corresponding estimates in HyperGEN (outcome IVs) served as the inputs for the MR in TwoSampleMR. We then harmonized allele orientation and effect estimates across the two cohorts. For CpG-metabolite pairs in which there was evidence to suggest a causal relationship in the forward direction (metabolite causal on the CpG, MR p <0.05), we investigated reverse causation as well. In reverse MR, the exposure IVs were SNPs that were strongly associated with the respective CpG in HyperGEN with p <1.00e-06 (methylation quantitative trait loci SNPs or meQTL SNPs), and the corresponding SNP associations from the metabolite GWAS were the outcome IVs.
We obtained MR estimates via multiple methods: Egger, weighted median, inverse variance-weighted (IVW), simple, and weighted. Generally, the IVW method provides the greatest statistical power, but it assumes balanced pleiotropy. Thus, when there was significant pleiotropy (Egger regression intercept p <0.05), we assessed the Egger estimate; otherwise, we considered the IVW estimate. We also conducted multiple sensitivity analyses for each CpG-metabolite pair in both forward and reverse MR (Burgess et al., 2017). These included tests for heterogeneity, single SNP MR (which computes MR estimates using each SNP as an individual IV instead of a combined effect of all the SNPs selected as IVs), and leave-one-out analyses—which computes MR estimates multiple times while leaving out each SNP from the set of exposure IVs to determine whether an individual SNP is the primary driver of the observed effect. We also conducted the Steiger test for causal direction, which determined whether the proportion of the variance in the exposure explained by the exposure IVs was significantly greater than the proportion of the variance in the outcome explained by those same IVs. In all analyses, two-sided p-values are assessed.
Results
Baseline characteristics of the REGARDS and HyperGEN sub-study cohorts are presented in Table 2. In discovery models, 28 metabolites were significantly associated with prevalent MetS in REGARDS after correcting for multiple comparisons (Supplementary Figure S2). Of these, two metabolites (cystathionine and S-adenosylhomocysteine (SAH)) were known to be upstream of DNA methylation pathways. Of the 10 remaining methylation-related metabolites, the following were marginally associated with MetS (p <0.05): cystine, S-adenosylmethionine (SAM), sarcosine, serine. Unlike other methylation metabolites, higher serine levels were associated with a lower odds of MetS in REGARDS.
TABLE 2. Baseline characteristics of REGARDS and HyperGEN participants with metabolite measurements.
Moreover, associations were robust when we restricted the analysis to the controls in the stroke case-cohort and when we accounted for fasting status (Supplementary Table S1). Of the 38 total metabolites that were either statistically significant and/or are involved in DNA methylation pathways, 31 metabolites were available for analysis in HyperGEN. In validation models, D-gluconic acid (DGA), glucose, isoleucine, and leucine were strongly associated with MetS (p <1.00e-05). Additionally, 2-aminoadipic acid (2-AAA), alanine, glutamate, serine, and taurine were modestly associated with MetS (p <0.05). Associations for MetS metabolites in REGARDS and HyperGEN are presented in Supplementary Table S2.
In HyperGEN, 44.7% of participants in the metabolomics cohort also had available methylation data. The prevalence of MetS in this subset (74.6%) was concordant with the full metabolomics cohort (71.7%). In the methylation analysis, we assessed the relationships between the 9 validated MetS metabolites and the selected CpG sites. In the fully adjusted models, DGA, glucose, glutamate, isoleucine, serine, and taurine were associated with methylation of at least one CpG site. However, 2-AAA, alanine, and leucine, were not associated with any CpG site. Furthermore, cg18181703 (SOCS3) was not associated with any metabolite.
CpG-metabolite pairs that were statistically significant in linear mixed modeling were eligible for MR analysis (Table 3). All CpG-metabolite associations are presented in Supplementary Table S3. Of the 13 statistically significant CpG-metabolite associations, only glucose (CPT1A) and serine (ABCG1, CPT1A, SREBF1) metabolites showed causal effects on methylation in the forward MR (Table 4). Glucose SNPs were inversely associated methylation of cg00574958 (CPT1A). And serine-associated SNPs were inversely associated with methylation of CpG sites at ABCG1 and SREBF1, though positively associated with CpG sites at CPT1A and TXNIP.
Metabolite effects on methylation were robust in sensitivity analyses for the selection of SNPs as IVs, including leave-one-out analyses and single SNP MR (Figure 3). Moreover, tests for heterogeneity were not statistically significant for any CpG-metabolite pair, indicating that the SNPs that were used as IVs did not violate MR assumptions. And the Steiger tests confirmed the causal directionality of metabolite influence on methylation.
FIGURE 3. MR sensitivity analyses for serine and glucose metabolites. Single SNP MR (left) and leave-one-out MR (right) for (A) serine-ABCG1 and (B) glucose-CPT1A.
For metabolites with significant causal effects on methylation in the forward MR, reverse MR showed that CpG methylation was not causal on metabolites, apart from cg06500161 (ABCG1) on serine (B [SE]: −1.534 [0.668], p = 0.023). Results were consistent in leave-one-out and single SNP MR estimations. Because there were more than 700 IVs for ABCG1 in reverse MR, we considered a p-value threshold of p <5.00e-08 to reduce the number of SNPs that were potentially in linkage disequilibrium (LD) with causal IVs. Complete reverse MR results are presented in Supplementary Table S4.
Discussion
In this study, we leveraged individual-level metabolomic, epigenomic, and genomic data to, first, identify metabolites associated with MetS and, second, evaluate the relationship between these metabolites and CpG sites that have been linked to MetS in prior studies. We determined that serine and glucose may causally affect DNA methylation and that there may be a bidirectional relationship between serine and methylation of the ABCG1 gene region. Notably, serine is also upstream of DNA methyltransferase reactions, unlike the metabolites (e.g., DGA) that were strongly associated with MetS and methylation but were null in MR analyses.
Previous metabolome-wide association studies (MWAS) for MetS have identified metabolites of statistical significance (Roberts et al., 2020; Azab et al., 2021; Carioca et al., 2021; Sun et al., 2021; Lind et al., 2022). However, replication of these metabolites across cohorts is lacking, and integration of additional ‘-omics’ data has been recommended to further explicate the pathophysiology of these variations in metabolite levels (Monnerie et al., 2020; Comte et al., 2021). For example, an integrative proteomics-metabolomics analysis identified metabolic pathways that were perturbed in the setting of MetS, including serine metabolism (Chen et al., 2022). Another study that combined metabolomic, proteomic, and transcriptomic data showed that mitochondrial respiratory chain dysfunction alters serine biosynthesis pathways (Bao et al., 2016). Although mechanisms are not fully understood, multi-omics approaches may illuminate biochemical pathways that contribute to MetS.
In this analysis, we detected causal relationships between serine and methylation of multiple MetS CpG sites. The inverse associations with CpG sites at ABCG1 and SREBF1, as well as positive associations cg00574958 at CPT1A, suggest a protective relationship (e.g., increased CPT1A methylation is associated with less risk for MetS). Serine contributes to epigenetic modifications, as its metabolism supports the methionine cycle and subsequent methylation of DNA, RNA, and histones via methyltransferase reactions (Kadayifci et al., 2018; Gadecka and Bielak-Zmijewska, 2019; Lin et al., 2022; Pardo et al., 2022). The effect of serine metabolism on epigenetic variations has primarily been weighed in the context of cancer, but this metabolite could similarly affect cardiometabolic health via these same methylation pathways. For example, animal studies suggest that dietary supplementation with serine and taurine may prevent MetS (Ji et al., 2016; Rafiee et al., 2022; Wang et al., 2022). And a recent study in REGARDS found that serine levels were positively associated with a plant-based dietary pattern (Bhave et al., 2022). While mechanistic studies linking serine to DNA methylation and cardiometabolic dysfunction are lacking, current evidence provides a theoretical framework for the observed associations in this study.
EWAS for MetS and related traits have shown, with consistent validation across diverse populations, associations with the CpG sites we selected for this analysis (Akinyemiju et al., 2018; Hidalgo et al., 2021). Although this study was not able to link methylation to gene expression, previous studies have demonstrated this, as well as identified causal relationships with cardiometabolic traits in Mendelian randomization analyses (Cameron et al., 2023). Previously, Hidalgo et al. demonstrated that the selected CpG sites could be combined to form a methylation risk score for MetS in the HyperGEN (African ancestry) and GOLDN (European ancestry) cohorts (Hidalgo et al., 2021). Further, across multiple studies, cg00574958 (CPT1A) has demonstrated inverse associations with triglycerides and very low-density lipoprotein (VLDL) cholesterol (Irvin et al., 2014; Lai et al., 2016; Braun et al., 2017; Zheng et al., 2022). Similarly, inverse relationships between cg19693031 (TXNIP) and diabetes have been observed (Nuotio et al., 2020; Miller et al., 2023). These sites are located at or near CpG islands and other regulatory regions (Table 1). Thus, variations in their methylation may be linked to altered expression of genes that contribute to lipid metabolism, hyperglycemia, and inflammation (Jones et al., 2021; Lu et al., 2022). While studies are limited, MR for methylation and cardiometabolic traits suggest that altered methylation is the consequence, not the cause, of metabolic dysfunction (Dekkers et al., 2016; Sayols-Baixeras et al., 2018; Zaghlool et al., 2018). It is possible that, in the setting of MetS, the catabolism of certain metabolites (e.g., serine) is altered and may induce methylation changes; and as such, DNA methylation functions as a biomarker of an existing metabolic disturbance. At the same time, the bidirectional influence of methylation on metabolic traits cannot be fully ruled out. For example, a nested case-control study found that methylation of cg06500161 at ABCG1 was associated with incident diabetes in a population of adults in rural China (Qie et al., 2021). In our intermediate analyses of methylation and metabolites, we observed significant associations in the expected direction of effect based on prior studies, as well as a significant causal relationship of ABCG1 methylation on serine. Thus, there may be biological plausibility for bidirectional relationships between the environment, metabolome, methylation, and MetS.
While this study was neither the first MWAS of MetS nor the MR of metabolites and methylation, it is one of the few to validate MetS-associated metabolites in an external cohort. And, to our knowledge, it is the first to evaluate bidirectional causality. In a multi-omics analysis, Zaghlool et al. detected a causal, inverse relationship of VLDL-A on cg00574958 at CPT1A (Zaghlool et al., 2018). In our findings, glucose was also inversely associated with this CpG site, whereas there was a positive causal association of serine. Furthermore, ABCG1 methylation was causal on serine levels in reverse MR, whereas Zaghlool et al. were unable to investigate reverse causation. Overall, a lack of individual-level metabolite, methylation, and SNP data within a single cohort has limited investigators to one-sample and/or unidirectional MR approaches (Zaghlool et al., 2018; Wu et al., 2021). The former can lead to instrument bias, and in the latter, one cannot assess bidirectionality (Hartwig et al., 2016). We conducted MR with data from separate studies—thereby reducing the likelihood of instrument bias—and assessed reverse causation, unlike prior MR for DNA methylation and metabolites. Still, the number of participants in HyperGEN with available methylation and metabolite measurements were limited (N<150). This lack of power, coupled with the conservative Bonferroni-corrected significance threshold, may have led to false negatives. However, we also selected a group of metabolites with biological plausibility (i.e., potential methyl donors via DNA methyltransferase reactions) for causal testing even if they did not meet the corrected p-value threshold. Other limitations include between-study comparability, as REGARDS is a study of older adults (>45 years at baseline), and HyperGEN was enriched for hypertension. However, the robustness of our findings despite these cohort differences may suggest good generalizability.
In conclusion, we 1) conducted metabolite profiling for MetS; 2) evaluated the relationship between MetS metabolites and methylation of known MetS CpG sites; and 3) determined the causal direction of these relationships. Not only did we apply metabolite profiling in one of the largest and most diverse studies of MetS to date, but we also validated our findings in an external cohort for a total sample size of more than 2,000 individuals. The results showed that metabolites were causal on methylation, and one CpG-metabolite pair (cg00650161 (ABCG1)-serine) may also exhibit bidirectional influence. Our top findings suggested protective effects of serine in MetS, as well as causal relationships between serine and methylation of CpG sites linked to MetS. These findings were robust to multiple sensitivity analyses.
Longitudinal and mechanistic studies are needed to confirm these findings and further elucidate biological pathways. Further validation in more populations is also needed, as these relationships may vary in the context of different environments. Expanding our understanding of these relationships could potentially lead to the identification of precise biomarkers of MetS risk and cardiometabolic health overall.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: The metabolomics datasets used and/or analyzed during the current study are not publicly available but are available from the study upon reasonable request. However, genetic data and methylation data from HyperGEN can be accessed upon request via the TOPMed Consortium. Raw genotype data from REGARDS can be accessed via dbGaP (accession phs002719.v1.p1). Requests to access these datasets should be directed to REGARDS Study, cmVnYXJkc2FkbWluQHVhYi5lZHU=; TOPMed Consortium, https://dbgap.ncbi.nlm.nih.gov/aa/wga.cgi?page=login.
Ethics statement
The studies involving humans were approved by the University of Alabama at Birmingham Institutional Review Board. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
Conceptualization: AJ, ZA, WK, and MI. Methodology: AJ, ZA, AP, HT, and VS. Data collection/curation: ZA, NC, VS, DKA, and DMA. Data analysis: AJ and AP. Writing—original draft preparation: AJ and MI. Writing—review and editing: ZA, AP, NC, NK, HT, VS, WK, and MI. Supervision: WK and MI. Funding: DA, WK, and MI. All authors contributed to the article and approved the submitted version.
Funding
This research project is supported by the cooperative agreement U01 NS041588 co-funded by the National Institute of Neurological Disorders and Stroke (NINDS) and the National Institute on Aging (NIA), National Institutes of Health, Department of Health and Human Service. Additional funding was provided by the NIH/NINDS R01NS099209 and NIH/NHLBI R35HL155466. HyperGEN was also supported by the NIH/NHLBI R01HL055673. AJ is also supported by the National Institute on Diabetes and Digestive and Kidney Diseases (NIH/NIDDK) through F31DK128990.
Acknowledgments
The authors thank the other investigators, the staff, and the participants of the REGARDS and HyperGEN studies for their valuable contributions. A full list of participating REGARDS investigators and institutions can be found at: https://www.uab.edu/soph/regardsstudy/.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Author disclaimer
The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH and NINDS or the NIA. Representatives from NIH and NINDS did not have any role in the design and conduct of the study, the collection, management, analysis, and interpretation of the data, or the preparation or approval of the manuscript.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1184661/full#supplementary-material
References
Akinyemiju, T., Do, A. N., Patki, A., Aslibekyan, S., Zhi, D., Hidalgo, B., et al. (2018). Epigenome-wide association study of metabolic syndrome in African-American adults. Clin. Epigenetics 10, 49. doi:10.1186/s13148-018-0483-2
Ament, Z., Patki, A., Chaudhary, N., Bhave, V. M., Garcia Guarniz, A. L., Gao, Y., et al. (2022). Nucleosides associated with incident ischemic stroke in the REGARDS and JHS cohorts. Neurology 98 (21), e2097–e2107. doi:10.1212/WNL.0000000000200262
Azab, S. M., de Souza, R. J., Lamri, A., Shanmuganathan, M., Kroezen, Z., Schulze, K. M., et al. (2021). Metabolite profiles and the risk of metabolic syndrome in early childhood: A case-control study. BMC Med. 19 (1), 292. doi:10.1186/s12916-021-02162-7
Bao, X. R., Ong, S. E., Goldberger, O., Peng, J., Sharma, R., Thompson, D. A., et al. (2016). Mitochondrial dysfunction remodels one-carbon metabolism in human cells. Elife 5, e10575. doi:10.7554/eLife.10575
Bhave, V. M., Ament, Z., Patki, A., Gao, Y., Kijpaisalratana, N., Guo, B., et al. (2022). Plasma metabolites link dietary patterns to stroke risk. Ann. Neurol. 93, 500–510. doi:10.1002/ana.26552
Braun, K. V. E., Dhana, K., de Vries, P. S., Voortman, T., van Meurs Jbj, , Uitterlinden, A. G., et al. (2017). Epigenome-wide association study (EWAS) on lipids: the rotterdam study. Clin. Epigenetics 9, 15. doi:10.1186/s13148-016-0304-4
Burgess, S., Bowden, J., Fall, T., Ingelsson, E., and Thompson, S. G. (2017). Sensitivity analyses for robust causal inference from mendelian randomization analyses with multiple genetic variants. Epidemiology 28 (1), 30–42. doi:10.1097/EDE.0000000000000559
Cameron, V. A., Jones, G. T., Horwood, L. J., Pilbrow, A. P., Martin, J., Frampton, C., et al. (2023). DNA methylation patterns at birth predict health outcomes in young adults born very low birthweight. Clin. Epigenetics 15 (1), 47. doi:10.1186/s13148-023-01463-3
Carioca, A. A. F., Steluti, J., Carvalho, A. M., Silva, A. M., Silva, I. D. C. G., Fisberg, R. M., et al. (2021). Plasma metabolomics are associated with metabolic syndrome: A targeted approach. Nutrition 83, 111082. doi:10.1016/j.nut.2020.111082
Chen, M., Yang, Z., Gan, H., Wang, Y., Li, C., and Gao, Y. (2022). Investigation into potential mechanisms of metabolic syndrome by integrative analysis of metabolomics and proteomics. PLoS One 17 (7), e0270593. doi:10.1371/journal.pone.0270593
Comte, B., Monnerie, S., Brandolini-Bunlon, M., Canlet, C., Castelli, F., Chu-Van, E., et al. (2021). Multiplatform metabolomics for an integrative exploration of metabolic syndrome in older men. EBioMedicine 69, 103440. doi:10.1016/j.ebiom.2021.103440
Dekkers, K. F., van Iterson, M., Slieker, R. C., Moed, M. H., Bonder, M. J., van Galen, M., et al. (2016). Blood lipids influence DNA methylation in circulating cells. Genome Biol. 17 (1), 138. doi:10.1186/s13059-016-1000-6
Fahed, G., Aoun, L., Bou Zerdan, M., Allam, S., Bouferraa, Y., Assi, H. I., et al. (2022). Metabolic syndrome: updates on pathophysiology and management in 2021. Int. J. Mol. Sci. 23 (2), 786. doi:10.3390/ijms23020786
Gadecka, A., and Bielak-Zmijewska, A. (2019). Slowing down ageing: the role of nutrients and microbiota in modulation of the epigenome. Nutrients 11 (6), 1251. doi:10.3390/nu11061251
Hartwig, F. P., Davies, N. M., Hemani, G., and Davey Smith, G. (2016). Two-sample mendelian randomization: avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int. J. Epidemiol. 45 (6), 1717–1726. doi:10.1093/ije/dyx028
Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base platform supports systematic causal inference across the human phenome. Elife 7, e34408. doi:10.7554/eLife.34408
Hidalgo, B. A., Minniefield, B., Patki, A., Tanner, R., Bagheri, M., Tiwari, H. K., et al. (2021). A 6-CpG validated methylation risk score model for metabolic syndrome: the HyperGEN and GOLDN studies. PLoS One 16 (11), e0259836. doi:10.1371/journal.pone.0259836
Howard, V. J., Cushman, M., Pulley, L., Gomez, C. R., Go, R. C., Prineas, R. J., et al. (2005). The reasons for geographic and racial differences in stroke study: objectives and design. Neuroepidemiology 25 (3), 135–143. doi:10.1159/000086678
Huang, P. L. (2009). A comprehensive definition for metabolic syndrome. Dis. Model Mech. 2 (5-6), 231–237. doi:10.1242/dmm.001180
Irvin, M. R., Zhi, D., Joehanes, R., Mendelson, M., Aslibekyan, S., Claas, S. A., et al. (2014). Epigenome-wide association study of fasting blood lipids in the Genetics of Lipid-lowering Drugs and Diet Network study. Circulation 130 (7), 565–572. doi:10.1161/CIRCULATIONAHA.114.009158
Ji, Y., Wu, Z., Dai, Z., Sun, K., Wang, J., and Wu, G. (2016). Nutritional epigenetics with a focus on amino acids: implications for the development and treatment of metabolic syndrome. J. Nutr. Biochem. 27, 1–8. doi:10.1016/j.jnutbio.2015.08.003
Jialal, I., Patel, A., Devaraj, S., and Adams-Huet, B. (2021). Metabolites that activate the inflammasome in nascent metabolic syndrome. J. Diabetes Complicat. 35 (3), 107836. doi:10.1016/j.jdiacomp.2020.107836
Jones, A. C., Irvin, M. R., Claas, S. A., and Arnett, D. K. (2021). Lipid phenotypes and DNA methylation: A review of the literature. Curr. Atheroscler. Rep. 23 (11), 71. doi:10.1007/s11883-021-00965-w
Kadayifci, F. Z., Zheng, S., and Pan, Y. X. (2018). Molecular mechanisms underlying the link between diet and DNA methylation. Int. J. Mol. Sci. 19 (12), 4055. doi:10.3390/ijms19124055
Lai, C. Q., Wojczynski, M. K., Parnell, L. D., Hidalgo, B. A., Irvin, M. R., Aslibekyan, S., et al. (2016). Epigenome-wide association study of triglyceride postprandial responses to a high-fat dietary challenge. J. Lipid Res. 57 (12), 2200–2207. doi:10.1194/jlr.M069948
Lin, J. G., Kourtis, S., Ghose, R., Pardo Lorente, N., Kubicek, S., and Sdelci, S. (2022). Metabolic modulation of transcription: the role of one-carbon metabolism. Cell Chem. Biol. 29, 1664–1679. doi:10.1016/j.chembiol.2022.11.009
Lind, L., Salihovic, S., Sundström, J., Elmståhl, S., Hammar, U., Dekkers, K., et al. (2022). Metabolic profiling of obesity with and without the metabolic syndrome - a multi-sample evaluation. J. Clin. Endocrinol. Metab. 107, 1337–1345. doi:10.1210/clinem/dgab922
Lu, B., Chen, J., Xu, G., Grayson, T. B., Jing, G., Jo, S., et al. (2022). Alpha cell thioredoxin-interacting Protein deletion improves diabetes-associated hyperglycemia and hyperglucagonemia. Endocrinology 163, bqac133. doi:10.1210/endocr/bqac133
Maddocks, O. D., Labuschagne, C. F., Adams, P. D., and Vousden, K. H. (2016). Serine metabolism supports the methionine cycle and DNA/RNA methylation through de novo ATP synthesis in cancer cells. Mol. Cell 61 (2), 210–221. doi:10.1016/j.molcel.2015.12.014
Miller, R. G., Mychaleckyj, J. C., Onengut-Gumuscu, S., Orchard, T. J., and Costacou, T. (2023). TXNIP DNA methylation is associated with glycemic control over 28 years in type 1 diabetes: findings from the Pittsburgh epidemiology of diabetes complications (EDC) study. BMJ Open Diabetes Res. Care 11 (1), e003068. doi:10.1136/bmjdrc-2022-003068
Monnerie, S., Comte, B., Ziegler, D., Morais, J. A., Pujos-Guillot, E., and Gaudreau, P. (2020). Metabolomic and lipidomic signatures of metabolic syndrome and its physiological components in adults: A systematic review. Sci. Rep. 10 (1), 669. doi:10.1038/s41598-019-56909-7
Moore, J. X., Chaudhary, N., and Akinyemiju, T. (2017). Metabolic syndrome prevalence by race/ethnicity and sex in the United States, national health and nutrition examination survey, 1988-2012. Prev. Chronic Dis. 14, E24. doi:10.5888/pcd14.160287
Nuotio, M. L., Pervjakova, N., Joensuu, A., Karhunen, V., Hiekkalinna, T., Milani, L., et al. (2020). An epigenome-wide association study of metabolic syndrome and its components. Sci. Rep. 10 (1), 20567. doi:10.1038/s41598-020-77506-z
Pardo, J. C., Ruiz de Porras, V., Gil, J., Font, A., Puig-Domingo, M., and Jordà, M. (2022). Lipid metabolism and epigenetics crosstalk in prostate cancer. Nutrients 14 (4), 851. doi:10.3390/nu14040851
Qie, R., Chen, Q., Wang, T., Chen, X., Wang, J., Cheng, R., et al. (2021). Association of ABCG1 gene methylation and its dynamic change status with incident type 2 diabetes mellitus: the rural Chinese cohort study. J. Hum. Genet. 66 (4), 347–357. doi:10.1038/s10038-020-00848-z
Rafiee, Z., García-Serrano, A. M., and Duarte, J. M. N. (2022). Taurine supplementation as a neuroprotective strategy upon brain dysfunction in metabolic syndrome and diabetes. Nutrients 14 (6), 1292. doi:10.3390/nu14061292
Reilly, M. P., and Rader, D. J. (2003). The metabolic syndrome: more than the sum of its parts? Circulation 108 (13), 1546–1551. doi:10.1161/01.CIR.0000088846.10655.E0
Roberts, J. A., Varma, V. R., Huang, C. W., An, Y., Oommen, A., Tanaka, T., et al. (2020). Blood metabolite signature of metabolic syndrome implicates alterations in amino acid metabolism: findings from the Baltimore longitudinal study of aging (BLSA) and the tsuruoka metabolomics cohort study (TMCS). Int. J. Mol. Sci. 21 (4), 1249. doi:10.3390/ijms21041249
Roberts, L. D., Koulman, A., and Griffin, J. L. (2014). Towards metabolic biomarkers of insulin resistance and type 2 diabetes: progress from the metabolome. Lancet Diabetes Endocrinol. 2 (1), 65–75. doi:10.1016/S2213-8587(13)70143-8
Sayols-Baixeras, S., Tiwari, H. K., and Aslibekyan, S. W. (2018). Disentangling associations between DNA methylation and blood lipids: A mendelian randomization approach. BMC Proc. 12 (9), 23. doi:10.1186/s12919-018-0119-8
Sun, S., He, D., Luo, C., Lin, X., Wu, J., Yin, X., et al. (2021). Metabolic syndrome and its components are associated with altered amino acid profile in Chinese han population. Front. Endocrinol. (Lausanne) 12, 795044. doi:10.3389/fendo.2021.795044
Wang, S., Jung, S., and Ko, K. S. (2022). Effects of amino acids supplementation on lipid and glucose metabolism in HepG2 cells. Nutrients 14 (15), 3050. doi:10.3390/nu14153050
Williams, R. R., Rao, D. C., Ellison, R. C., Arnett, D. K., Heiss, G., Oberman, A., et al. (2000). NHLBI family blood pressure program: methodology and recruitment in the HyperGEN network. Hypertension genetic epidemiology network. Ann. Epidemiol. 10 (6), 389–400. doi:10.1016/s1047-2797(00)00063-6
Wu, Q., Li, J., Sun, X., He, D., Cheng, Z., Zhang, X., et al. (2021). Multi-stage metabolomics and genetic analyses identified metabolite biomarkers of metabolic syndrome and their genetic determinants. EBioMedicine 74, 103707. doi:10.1016/j.ebiom.2021.103707
Yang, M., and Vousden, K. H. (2016). Serine and one-carbon metabolism in cancer. Nat. Rev. Cancer 16 (10), 650–662. doi:10.1038/nrc.2016.81
Zaghlool, S. B., Mook-Kanamori, D. O., Kader, S., Stephan, N., Halama, A., Engelke, R., et al. (2018). Deep molecular phenotypes link complex disorders and physiological insult to CpG methylation. Hum. Mol. Genet. 27 (6), 1106–1121. doi:10.1093/hmg/ddy006
Zhang, N. (2018). Role of methionine on epigenetic modification of DNA methylation and gene expression in animals. Anim. Nutr. 4 (1), 11–16. doi:10.1016/j.aninu.2017.08.009
Zhang, Y., Zhang, H., Rong, S., Bian, C., Yang, Y., and Pan, H. (2021). NMR spectroscopy based metabolomics confirms the aggravation of metabolic disorder in metabolic syndrome combined with hyperuricemia. Nutr. Metab. Cardiovasc Dis. 31 (8), 2449–2457. doi:10.1016/j.numecd.2021.05.015
Keywords: Mendelian randomization, metabolic syndrome, metabolomics, DNA methylation, multiomics
Citation: Jones AC, Ament Z, Patki A, Chaudhary NS, Srinivasasainagendra V, Kijpaisalratana N, Absher DM, Tiwari HK, Arnett DK, Kimberly WT and Irvin MR (2023) Metabolite profiles and DNA methylation in metabolic syndrome: a two-sample, bidirectional Mendelian randomization. Front. Genet. 14:1184661. doi: 10.3389/fgene.2023.1184661
Received: 12 March 2023; Accepted: 07 September 2023;
Published: 15 September 2023.
Edited by:
Michael E. Symonds, University of Nottingham, United KingdomReviewed by:
Swapan Kumar Das, Wake Forest University, United StatesMartin Broberg, University of Helsinki, Finland
Copyright © 2023 Jones, Ament, Patki, Chaudhary, Srinivasasainagendra, Kijpaisalratana, Absher, Tiwari, Arnett, Kimberly and Irvin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alana C. Jones, YWNqb25lc0B1YWIuZWR1