Skip to main content
  • Research article
  • Open access
  • Published:

Circulating metabolites may illustrate relationship of alcohol consumption with cardiovascular disease

Abstract

Background

Metabolite signatures of long-term alcohol consumption are lacking. To better understand the molecular basis linking alcohol drinking and cardiovascular disease (CVD), we investigated circulating metabolites associated with long-term alcohol consumption and examined whether these metabolites were associated with incident CVD.

Methods

Cumulative average alcohol consumption (g/day) was derived from the total consumption of beer, wine, and liquor on average of 19 years in 2428 Framingham Heart Study Offspring participants (mean age 56 years, 52% women). We used linear mixed models to investigate the associations of alcohol consumption with 211 log-transformed plasma metabolites, adjusting for age, sex, batch, smoking, diet, physical activity, BMI, and familial relationship. Cox models were used to test the association of alcohol-related metabolite scores with fatal and nonfatal incident CVD (myocardial infarction, coronary heart disease, stroke, and heart failure).

Results

We identified 60 metabolites associated with cumulative average alcohol consumption (p < 0.05/211 ≈ 0.00024). For example, 1 g/day increase of alcohol consumption was associated with higher levels of cholesteryl esters (e.g., CE 16:1, beta = 0.023 ± 0.002, p = 6.3e − 45) and phosphatidylcholine (e.g., PC 32:1, beta = 0.021 ± 0.002, p = 3.1e − 38). Survival analysis identified that 10 alcohol-associated metabolites were also associated with a differential CVD risk after adjusting for age, sex, and batch. Further, we built two alcohol consumption weighted metabolite scores using these 10 metabolites and showed that, with adjustment age, sex, batch, and common CVD risk factors, the two scores had comparable but opposite associations with incident CVD, hazard ratio 1.11 (95% CI = [1.02, 1.21], p = 0.02) vs 0.88 (95% CI = [0.78, 0.98], p = 0.02).

Conclusions

We identified 60 long-term alcohol consumption-associated metabolites. The association analysis with incident CVD suggests a complex metabolic basis between alcohol consumption and CVD.

Peer Review reports

Background

Alcohol drinking is a common lifestyle in many cultures and is a modifiable risk factor associated with over 200 health problems, including cardiovascular disease (CVD) [1], dementia [2], neuropsychiatric conditions [3], liver cirrhosis [4], diabetes [5], and several types of cancer, e.g., gastric cancer [6] and liver cancer [7]. For example, a systematic review and meta-analysis of 23 observational studies found that moderate alcohol consumption was related to a higher cardiovascular risk within 24 h after alcohol intake; however, after 24 h, moderate alcohol consumption seemed to have a protective effect on myocardial infraction [8]. In contrast, heavy alcohol intake had a continued risk for cardiovascular events [8].

Metabolites are small molecules that are intermediates or end-products of metabolism in many cellular processes [9, 10]. Metabolites can be quantified via high-throughput liquid chromatography with tandem mass spectrometry (LC/MS) methods [11]. Association analyses of alcohol consumption and metabolites may help us gain insights into the effect of alcohol consumption on disease pathways. Several European studies investigated associations of total alcohol intake with circulating concentrations of metabolites [12,13,14]. One study investigated the association of 123 metabolites with alcohol intake using a discovery (n = 1983) and replication (n = 991) study design in healthy participants. This study found that 72 metabolites were significantly associated with alcohol consumption in the discovery set, and 34 metabolites remained significance in the replication set [12]. A population-based study investigated 131 metabolites in 2090 individuals and revealed that 40 metabolites in men and 18 metabolites in women significantly differed in their concentrations between moderate-to-heavy and light alcohol drinking [14].

While these metabolomic studies help us to understand the potential molecular basis of alcohol consumption, most studies analyzed the alcohol consumption measured at a single time point, which may not represent the habitual or long-term alcohol consumption. A better understanding of the relationships between long-term alcohol consumption and circulating metabolites, as well as the relationships of alcohol-associated metabolites with CVD risk, may help elucidate the complex relationship of alcohol consumption with etiology and progression of alcohol-associated diseases. To that end, our study used longitudinal data from the Framingham Heart Study (FHS) Offspring cohort to investigate the following three aims (Fig. 1). First, we conducted association analyses of cumulative average alcohol consumption over 20 years with circulating metabolites. Second, we analyzed the specific associations of metabolites with cumulative consumption of three types of alcoholic beverages: beer, wine, and liquor. Third, we conducted association analyses of alcohol-associated metabolites with incident CVD.

Fig. 1
figure 1

The study flowchart

Methods

Population

The Offspring cohort (n = 5124) of the FHS was recruited in 1971, including the children of the Original cohort and spouses of these children [15]. The Offspring participants underwent in-person health examinations every 4 to 8 years to collect comprehensive demographic and clinical measures such as risk factors for cardiovascular and neurodegenerative diseases [16]. This study analyzed Offspring participants who had metabolite measurements at the fifth examination (n = 2526). We excluded participants who attended less than three examinations between examination 1 and examination 5 (n = 25) and those who were younger than 18 years old at examination 1 (n = 73). The final sample size was 2428 for statistical analyses (Fig. 1). All participants provided written informed consent. All FHS study protocols were approved by the Boston University Medical Center Institutional Review Board.

Metabolite quantification

Targeted metabolite profiling was conducted by LC/MS platform and hydrophilic interaction chromatography method as previously described [17,18,19]. The targeted metabolite profiling includes negatively charged polar metabolites (i.e., organic acids, bile acids, and sugars), positively charged polar metabolites (i.e., amino acids, urea cycle intermediates, nucleotides), and lipid metabolite species (i.e., cholesterol esters (CE), diacylglycerols (DAGs), lysophosphatidylcholines (LPCs), lysophosphatidylethanolamines (LPEs), phosphatidylcholines (PCs), sphingomyelins (SMs), triacylglycerols (TAGs)) [17,18,19,20]. Two metabolites, TAG 54:5 and TAG 54:6, had identical values for all participants; we removed TAG 54:6 from the present analysis. Five metabolites were also excluded owing to large proportion of missingness in greater than 1000 participants, resulting in 211 metabolites for the subsequent analyses. Each metabolite measurement was natural log transformed and then standardized to have mean at zero and standard deviation (SD) at 1. Metabolite values beyond ± 4 SD away from the mean of that metabolite were set to be at ± 4 SD [19].

Alcohol consumption

At each examination, alcohol consumption information was collected through the FHS technician-administered questionnaires [21]. The frequency of three types of alcoholic beverages (i.e., standard drinks of beer, wine, or liquor) in a typical week or month were collected. In this study, we used “grams of alcohol consumed per day”, which was derived from the summation of the total standard drinks of beer, wine (red or white), or liquor. One standard drink was defined as one 12 oz. beer, one 4 oz. wine, or one 1.5 oz. 80 proof alcohol spirit, and one drink contains approximately 14 g ethanol [21, 22]. We calculated the grams per day for all types of alcoholic beverages (i.e., total alcohol consumption), as well as for each type (i.e., grams per day of beer, wine, or liquor). To reflect long-term alcohol consumption, we calculated the cumulative average alcohol consumption, i.e., the mean over up to five examinations. We categorized participants into three groups according to their total average alcohol consumption: non-drinkers (0 g/day); moderate drinker (0.1–28 g/day in men and 0.1–14 g/day in women); and heavy drinkers (> 28 g/day in men and > 14 g/day in women).

Disease traits

CVD events were identified through adjudication by a panel of three physicians [23]. The CVD event data were obtained from annual health history updates based on inpatient or outpatient medical records, physical examinations, and mortality registry. CVD was comprised of myocardial infarction (MI), coronary heart disease (CHD), stroke, heart failure (HF), and death to any cardiovascular conditions [24]. We analyzed incident CVD outcome after removal of the prevalent cases at the fifth examination.

Covariates

Demographic information (sex and age), smoking status, dietary intake, time and intensity of physical activities, and medication use were obtained from standard questionnaire [25]. We calculated a physical activity index using the following formula: 1 × sleep hours/day + 1.1 × sedentary hours/day + 1.5 × slight activity hours/day + 2.4 × moderate activity hours/day + 5 × heavy activity hours/day. We used Dietary Approaches to Stop Hypertension (DASH) score to define diet quality. We calculated the diet quality score based on consumptions of eight dietary components: vegetables, fruits, nuts and legumes, whole grains, dairy, red and processed meat, sugar-sweetened beverages, and sodium [26]. Systolic blood pressure (SBP) was measured twice by physicians and the average blood pressure was used in the present analysis. Anthropometry was measured by technicians. Standard assays were used to determine levels of fasting glucose and plasma lipids. We defined diabetes if a participant had fasting glucose ≥ 126 mg/dL or using diabetes medications. All covariates were measured at the fifth examination.

Statistical analysis

Association analysis of metabolites and alcohol consumption

In primary association analyses, we analyzed the relationships of metabolites (dependent variables) with cumulative average alcohol consumption of total alcohol, beer, wine, and liquor as independent variables (Fig. 1). We used linear mixed models to quantify these associations, adjusting for age, sex, (metabolite measurement) batch, smoking status, BMI, physical activity, and diet quality score as fixed effect covariates (all measured at the fifth examination), and familial relationship as a random effect covariate. We used Bonferroni correction (p < 0.05/211 \(\approx\) 0.00024) to account for multiple hypothesis testing with 211 metabolites.

We conducted two secondary analyses as described below. The level of alcohol consumption is usually larger in men compared to women [27]; therefore, we tested sex–alcohol interaction by adding a production term of sex and alcohol consumption. We additionally adjusted for a serum creatinine-based estimated glomerular filtration rate (eGFR) to investigate if kidney function confounded the relationships between the cumulative average alcohol consumption and metabolites. (Fig. 1).

Association analysis of alcohol-associated metabolites and CVD

Metabolites significantly associated with cumulative average alcohol consumption were used as independent variables in testing for associations with incident CVD. A Cox proportional hazards regression model was used, adjusting for sex, age, and batch in the base model. In the multivariable model, we additionally adjusted for BMI, SBP, hypertension treatment status, diabetes, smoking status, and total and high-density lipoprotein cholesterol levels. All covariates included were collected at the fifth examination.

To study an aggregate effect of metabolites with the development of CVD, we constructed weighted composite scores with CVD-associated metabolites that were identified in the base model. We used estimates from the association analyses between cumulative average alcohol consumption and the metabolites to weigh the concentrations of the corresponding metabolites. We then calculated metabolite scores as the linear combination of weighted metabolite concentrations: \({T}_{ij}={\sum }_{j}{W}_{j}{V}_{ij}\). Two metabolite scores were calculated based on the direction of the associations between alcohol consumption, metabolites, and incident CVD. The first score aggregated metabolites with consistent direction, i.e., those were either positively or negatively associated with both alcohol consumption and CVD. The second score aggregated metabolites with opposite direction, i.e., those were positively associated with alcohol consumption and negatively associated with CVD or those were negatively associated with alcohol consumption and positively associated with CVD. The composite scores were standardized with mean at 0 and standard deviation at 1. We performed association analysis of the continuous composite scores with the development of CVD, adjusting for age, sex, batch, BMI, SBP, hypertension treatment status, diabetes, smoking status, and total and high-density lipoprotein cholesterol levels. All association analyses were performed using R software (version 4.0.5).

Pathway analysis for metabolites

Pathway analyses were performed to identify biological pathways related to the alcohol-associated metabolites using MetaboAnalyst 5.0 [28]. Two pathway libraries used in these analyses are Kyoto Encyclopaedia of Genes and Genomes (KEGG) and Small Molecule Pathway Database (SMPDB) for Homo sapiens. We conducted hypergeometric test with default parameters in MetaboAnalyst and reported significant pathways at false discovery rate (FDR) < 0.05.

Results

Participant characteristics

This study included mostly middle-aged to older community dwelling men and women (n = 2428, mean age = 55.9 \(\pm\) 9.3 years, 52% women at fifth examination) (Table 1). The median of the cumulative average of total alcohol consumption was 7.7 g/day (interquartile range 16.8 g/day). Our participants were overweight (mean BMI 27.5 kg/m2) and 18.3% of them were current smokers at the fifth examination. Compared to women, men had greater alcohol consumption at all examinations: median consumption was 13.8 g/day in men versus 4.4 g/day in women at fifth examination. The pairwise Spearman correlation between total alcohol consumption at all five examinations is shown in Additional file 1: Fig. S1.

Table 1 Participant characteristics

Associations of total alcohol consumption and metabolites

We found that 60 metabolites were significantly associated with the cumulative average total alcohol consumption (p < 0.05/211 \(\approx\) 0.00024) adjusting for potential confounders (Table 2, Additional file 8: Table. S1). Of the 60 metabolites, 40 metabolites displayed positive associations with the cumulative average alcohol consumption (Fig. 2, Additional file 8: Table. S1). The most significant metabolite was cholesteryl palmitoleate (CE16:1), a plasma cholesteryl ester that is involved in cholesterol metabolism [29]. One gram per day higher alcohol consumption was associated with higher level of cholesteryl palmitoleate (beta = 0.023, p = 6.3e − 45) in the blood. Several phosphatidylcholine metabolites (e.g., PC 32:1 and PC 34:1) were positively associated with alcohol consumption, for example, 1 g/day higher of alcohol consumption was associated with higher level of PC 32:1 (beta = 0.021, p = 3.1e − 38) in blood. Among the 20 metabolites (Fig. 2) that were negatively associated with alcohol consumption, triacylglycerol 54:4 (TAG 54:4) displayed the most significant association (beta =  − 0.017, p = 6.17e − 22). With additional adjustment for eGFR, the association between alcohol consumption and metabolites remained largely the same (Pearson correlation coefficient = 0.99, Additional file 2: Fig. S2). Compared to the associations with the cumulative average total alcohol consumption, analysis using the current alcohol consumption at the fifth examination yielded similar alcohol–metabolite associations (Pearson r = 0.99; Additional file 3: Fig. S3).

Table 2 The top 20 significant metabolites associated with cumulative average alcohol consumption
Fig. 2
figure 2

Volcano plot for association of cumulative average alcohol consumption and metabolites. Linear mixed model was performed adjusting for age, sex, batch, smoking status, BMI, physical activity index, and diet score as fixed effect, and family relationship as random effect. Blue dots represent significant relationship (p value < 0.00024), grey dots represent non-significant relationship. CE cholesteryl ester, PC phosphatidylcholine, TAG triacylglycerol

Association of each type of alcoholic beverages and metabolites

We found that 19 metabolites were significantly associated with the cumulative average consumption of beer, 30 metabolites were significantly associated with the cumulative average consumption of wine, and 32 metabolites were significantly associated with the cumulative average consumption of liquor (Fig. 3, Additional file 8: Table. S1). Among the significant ones, seven metabolites (CE 16:1, LPC 20:5, PC 32:0, PC 32:1, PC 34:1, PC-B 36:4, and fumarate-malate) were significantly associated with the cumulative consumption of all three types of alcoholic beverages (Additional file 4: Fig. S4).

Fig. 3
figure 3

Comparison of association between different types of alcohol beverages. A Comparison of effect sizes between different types of alcohol beverages. B Comparison of -log 10 (p values) between different types of alcohol beverages. All linear mixed models adjusted for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect, and family relationship as random effect. CE cholesterol esters, DAG diacylglycerols, LPC lysophosphatidylcholines, LPE lysophosphatidylethanolamines, PC phosphatidylcholines, SM sphingomyelins, TAG triacylglycerols. Metabolites shown in panel had significant association with at least one type of alcohol (beer, wine, and liquor) or total alcohol consumption

While the three types of alcoholic beverages displayed largely consistent associations with metabolites (Additional file 5: Fig. S5), beer consumption appeared to have slightly weaker associations compared to consumptions of wine and liquor (Additional file 5: Fig. S5). Among 50 metabolites associated with at least one type of alcohol consumption, we calculated the pair-wise ratio of regression coefficients for each metabolite. When we used an absolute ratio of two as cutoff, our analysis revealed that three metabolites (LPC 16:1, PC 34:3 and PC 38:4) displayed stronger positive association with the cumulative average consumption of beer; one metabolite (3-ureidopropionic acid) showed positive association with liquor consumption liquor; and three metabolites (betaine, 2-hydroxyglutaric acid, xanthurenate) were positively and two metabolites (isoleucine and valine) were negatively associated with wine consumption (Additional file 9: Table. S2).

Sex-specific associations of metabolite with cumulative average total alcohol consumption

We conducted interaction test for sex and alcohol consumption and found that 13 metabolites (betaine, CE 14:0, CE 16:1, PC 32:0, PC 32:1, PC 34:1, TAG 48:0, TAG 52:4, TAG 54:4, TAG 54:5, xanthurenate, dimethylglycine, and kynurenic acid) had a significant interaction term (p < 0.05/211 = 0.00024). Eleven (CE 14:0, CE 16:1, PC 32:0, PC 32:1, PC 34:1, TAG 48:0, TAG 52:4, TAG 54:4, TAG 54:5, xanthurenate, and dimethylglycine) of these 13 metabolites were in 60 metabolites that were significantly associated with the cumulative average total alcohol consumption. We applied sex-specific analysis to examine the associations of cumulative average total alcohol consumption with these 13 metabolites. For women, the associations were significant for all these 13 metabolites (p < 0.05/13 = 0.0038), while for men, 10 metabolites (CE 16:1, PC 32:0, PC 32:1, PC 34:1, TAG 48:0, TAG 52:4, TAG 54:4, TAG 54:5, xanthurenate, and dimethylglycine) had significant associations (Additional file 10: Table. S3). The observed sex-alcohol interaction was due to the stronger associations in women than that in men (Fig. 4). For example, 1 g/day higher alcohol consumption was associated with 0.036 standard deviation increases in the level of CE 16:1 in women (p value = 7.26e − 25), while 0.019 in men (p value = 4.25e − 22).

Fig. 4
figure 4

Comparison of associations of cumulative average alcohol consumption and metabolites between men and women. Linear mixed model was performed in women-only and men-only samples. Covariates included age, batch, smoking status, BMI, physical activity index and diet score as fixed effect, and family relationship as random effect. CE cholesteryl ester, PC phosphatidylcholine, TAG triacylglycerol

Pathway analysis for alcohol-associated metabolites

We conducted pathway analyses of the 60 metabolites that were significantly associated with the cumulative average total alcohol consumption. Using the KEGG library, we found three pathways to demonstrate evidence of enrichment among our significant metabolites including (1) arginine biosynthesis (FDR = 0.002); (2) valine, leucine, and isoleucine biosynthesis (FDR = 0.04); and (3) aminoacyl-tRNA biosynthesis (FDR = 0.007) (Additional file 11: Table. S4). Using the SMPDB library, we did not observe significant pathways at FDR 0.05 level (Additional file 12: Table. S5).

Association with incident CVD

We included 2138 participants without prevalent CVD in association analyses of the 60 alcohol-associated metabolites with incident CVD. A total of 636 participants developed incident CVD with a median follow-up of 18.1 years. With adjustment for sex and age, alcohol consumption was not associated with incident CVD, hazard ratio (HR) = 1, p = 0.87. Additional adjustments for BMI, SBP, hypertension treatment status, diabetes, smoking status, and total and high-density lipoprotein cholesterol levels (i.e., multivariable model) did not substantially alter the association. In addition, we observed no significant pairwise comparisons in nondrinkers, moderate drinkers, and heavy drinkers. In the multivariable model, HR was 1.09 for nondrinkers vs. moderate drinkers (p = 0.62) and 1.15 for heavy drinkers vs. moderate drinkers (p = 0.16).

Ten metabolites (TAG 50:2, TAG 50:1, TAG 48:1, TAG 48:0, isoleucine, leucine, TAG 52:3, TAG 46:0, PC 32:1, and glutamine) were significantly associated with the development of CVD (p < 0.05/60 \(\approx\) 0.0008) in the base model adjusting for age, sex, and batch (Additional file 13: Table. S6). Out of the ten metabolites, higher levels of nine metabolites (TAG 50:2, TAG 50:1, TAG 48:1, TAG 48:0, TAG 46:0, PC 32:1, TAG 52:3, isoleucine and leucine) and lower levels of one metabolite (glutamine) were associated with an increased risk of incident CVD (HR ranges from 1.16 to 1.33, p < 0.05/60≈0.0008). After additionally adjusting for BMI, SBP, hypertension treatment status, diabetes, smoking status, and total and high-density lipoprotein cholesterol level, the association with incident CVD remained nominally significant at p < 0.05 for four metabolites (TAG 50:2, TAG 50:1, TAG 48:1, and glutamine; Additional file 13: Table. S6). We further analyzed the associations of the 10 metabolites with the incidence of four CVD subtypes, including coronary heart disease (CHD), myocardial infarction (MI), heart failure (HF), and stroke. As shown in Additional file 14: Table. S7, with adjustment for age, sex, and batch, we observed that TAG 52:3 was associated with MI (HR = 1.36, p = 1.07e − 4), PC 32:1, TAG 48:0, TAG 50:1,TAG 50:2 and isoleucine were associated with HF (HR ranges from 1.27 to 1.40, p < 0.0008), TAG 48:1, TAG 50:1, TAG 50:2, TAG 52:3, isoleucine and leucine were associated with CHD (HR ranges from 1.23 to 1.32, p < 0.0008). However, for metabolites significantly related to CVD subtypes in the base model, additional adjustment for clinical risk factors attenuated the associations with CVD subtypes to nonsignificant (p > 0.05).

We created two weighted metabolite scores; the first score included six metabolites (TAG 50:2, TAG 50:1, TAG 48:1, TAG 48:0, TAG 46:0 and PC 32:1) that had positive association with both the cumulative average total alcohol drinking and incident CVD in the base models and glutamine, which was negatively associated with both the cumulative average total alcohol drinking and incident CVD in the base model (Additional file 15: Table. S8). The second score included three metabolites (TAG 52:3, isoleucine and leucine), which had inverse association with the cumulative average total alcohol drinking and positive association with incident CVD in the base model (Additional file 15: Table. S8). In multivariable model analysis, per standard deviation increase of the first metabolite score was associated with 11% higher hazard of incident CVD (p = 0.02, 95%CI = [1.02, 1.21]; Fig. 5) and per standard deviation increase of the second metabolite score was associated with 12% lower hazard of incident CVD (p = 0.02, 95%CI = [0.78, 0.98]; Fig. 5).

Fig. 5
figure 5

Forest plot for CVD with metabolite composite score. Seven-met, the weighted score obtained from seven metabolites, TAG 50:2, TAG 50:1, TAG 48:1, TAG 48:0, TAG 46:0, PC 32:1, and glutamine; three-met, the weighted score obtained from three metabolites, TAG 52:3, isoleucine, and leucine. Weights were obtained from the beta coefficients of association test between cumulative average alcohol consumption and metabolites in linear mixed models. In the Cox regression, base model included age, sex, and batch as covariates, and multivariable model additionally adjusted for BMI, SBP, hypertension treatment status, diabetes, smoking status, and total and high-density lipoprotein cholesterol level

Nonlinear association of alcohol consumption with metabolites

We compared metabolite level among nondrinkers, moderate drinkers, heavy drinkers, adjusting for age, sex, batch, BMI, SBP, hypertension treatment status, diabetes, smoking status, and total and high-density lipoprotein cholesterol levels. At p < 0.05, we observed that for two metabolites (LPC 16:1 and histidine), their level in moderate drinkers was either lower or higher than both of that in nondrinkers and heavy drinkers (Additional file 6: Fig. S6). None of the two metabolites were associated with incident CVD.

Discussion

In this present study, we investigated the associations between the cumulative average alcohol consumption and 211 circulating metabolites in the 2428 FHS participants. Of the 211 metabolites, the cumulative average alcohol consumption was associated with 60 metabolites. We found that nine metabolites were more strongly associated with a specific type of alcoholic beverage. We also found that the cumulative average total alcohol consumption displayed stronger associations with 13 metabolites in women than men. Furthermore, we created two alcohol consumption-associated metabolite scores and showed that they had comparable but opposite association with incident CVD. Taken together, with targeted metabolomic profiling, our study identified a series of alcohol consumption-associated circulating metabolites, and via these metabolites, alcohol consumption may have counteractive effects on CVD risk.

Our results showed that higher level of alcohol consumption was associated with higher plasma levels for about two-thirds of the 60 significant metabolites. Among the top positively associated metabolites were cholesteryl esters (e.g., CE16:1 and CE20:5), phosphatidylcholine (e.g., PC 32:1), and lysophosphatidylcholine (e.g., LPC 20:5). In addition, we observed that 14 plasma triacylglycerols (TAGs) were significantly associated with total alcohol consumption. Among alcohol-related TAGs, six (TAG 52:3, TAG 52:4, TAG 52:5, TAG 54:3, TAG 54:4, and TAG 54:5) displayed negative associations with alcohol consumption (i.e., lower alcohol consumption was associated with higher levels of TAGs) while eight (TAG 46:0, TAG 48:0, TAG 48:1, TAG 50:1, TAG 50:2, TAG 58:10, TAG 58:11, and TAG 60:12) displayed positive associations with alcohol consumption. TAGs emerge as biomarkers of a liver-to-β-cell axis that links hepatic β-oxidation to β-cell functional mass and insulin secretion in pancreas [30]. TAGs are broken into glycerol and free fatty acids in the process of lipolysis, and free fatty acids are either processed by beta-oxidation or converged to ketone [31]. In our association analyses between alcohol-associated metabolites and incident CVD, we showed that six TAGs (TAG 50:2, TAG 50:1, TAG 48:1, TAG 48:0, TAG 46:0, and TAG 52:3) were associated with incident CVD using the base models. Specifically, PC 32:1 and TAG 48:0 were positively associated with HF. In addition, TAG 52:3 was positively associated with both MI and CHD. Furthermore, TAG 50:1 and TAG 50:2 were positively associated with both HF and CHD. Among these six TAGs, the association remained significant (p < 0.05) with CVD for TAG 50:2, TAG 50:1, and TAG 48:1 after adjusting for common cardiometabolic CVD risk factors. However, likely due to the reduced number of cases for each CVD subtype (relative to the analysis using the composite incident CVD as the outcome variable), our multivariable analysis may lack sufficient statistical power. Therefore, we observed that additional adjustment for common cardiometabolic CVD risk factors attenuated associations with CVD subtypes (p > 0.05). Overall, these observations are in line with the well-known effects of alcohol intake on lipid metabolism [32].

Alcohol consumption was also associated with several types of circulating metabolites, other than TAGs. For example, we showed that alcohol consumption was associated with reduced levels of dimethylglycine. Dimethylglycine plays an important role in one-carbon metabolism as a methyl donor [33, 34]. This function may be related to our previous observations regarding the strong correlation of alcohol consumption with DNA methylation [22]. Nonetheless, the extent to which alcohol consumption affects the source of one-carbon metabolism and the subsequent impact on the risk of developing CVD warrant further investigation.

Among the metabolites that were negatively associated with alcohol consumption, valine, isoleucine, and leucine are branched-chain amino acids (BCAA). A recent review summarized the complex relationship between impaired BCAA homeostasis to CVD [35]. Several lines of evidence suggest that higher BCAA levels are associated with increased risk of obesity and diabetes, which is consistent with the present observations on the positive association of leucine and isoleucine with incident CVD. Nonetheless, future studies in larger sample size and diverse populations are needed to examine the relationships between alcohol consumption, BCAA, and CVD. Several lines of evidence suggest that higher BCAA levels are associated with increased risk of obesity and diabetes, e.g., Ho’s study demonstrated positive associations between BCAA levels and BMI, fasting glucose, Homeostatic Model Assessment for Insulin Resistance (HOMA-IR) in the same FHS participants [17]. These data are consistent with our previous findings on the inverse associations of alcohol consumption with obesity and type 2 diabetes [21], as well as the observations on the positive association of leucine and isoleucine with incident CVD in the present study. BCAAs are nitrogen donors for hepatic gluconeogenesis [36]. Our observations may support the notion that alcohol consumption, mainly moderate consumption, suppresses gluconeogenesis via lowering BCAA levels and subsequently controlling blood glucose to maintain normoglycemia. Nonetheless, future studies with larger sample sizes and diverse populations are needed to validate our observations, and experimental studies and clinical trials are needed to examine the relationships between alcohol consumption, BCAA, and CVD risk.

Our observations also highlight the complex relationship between alcohol consumption and circulating metabolites, which was demonstrated by the analysis using the two metabolite scores. Our observations suggest that, via circulating metabolites, alcohol drinking may have both positive and negative effects on CVD, and the two effects seemed to cancel each other out in our study samples. However, if certain factors disrupt the balance, it is possible that one effect may prevail over the other effect and leads to either increased or decreased risk of developing CVD. As such, future studies are warranted to understand what factors may modify the association of alcohol consumption and circulating metabolites, as well as their impact on the relationship of alcohol consumption with CVD development.

We observed that wine consumption and liquor had stronger associations with TAG, CE, and SM lipid metabolites, while beer had stronger associations with PC lipid metabolites. We also found that wine and liquor had different associations with amino acids, quinoline carboxylic acids, and hydroxy acids. Liquor consumption was significantly related to higher levels of 3-ureidopropionic acid, whereas wine consumption had stronger association with betaine, 2-hydroxyglutaric acid, xanthurenate, isoleucine, and valine compared to liquor. These observations suggest that consumption of different types of alcoholic beverages are associated with different metabolomic responses. However, the observed associations may also be driven by confounders such as unmeasured components in different alcoholic beverages or dietary and other environmental factors that were not adjusted for in the present analysis. Some of these metabolites such as betaine and isoleucine [35, 37] may play important roles in CVD development. Perhaps due to the short list of metabolites examined, the present study did not support the notion that a certain type of alcohol may bring specific benefits to reduce CVD risk. Future studies including a comprehensive list of metabolites are warranted to investigate this issue.

It was found by van Roekel et al. that total alcohol intake was associated with 34 circulating metabolites, including 3 acylcarnitines, the amino acid citrulline, 4 lysophosphatidylcholines, 13 diacylphosphatidylcholines, 7 acyl-alkylphosphatidylcholines, and 6 sphingomyelins [12] among middle-aged (mean age = 58.3) participants in European Prospective Investigation into Cancer and Nutrition Cohort. Among the 34 metabolites significant in the van Roekel’s study, 13 metabolites were also measured by our metabolite platforms, including LPC 16:0, LPC 16:1, LPC 20:4, PC 32:0, PC 32:1, PC 32:2, PC 34:1, PC 34:3, PC 34:4, PC-B 36:4, PC 38:6, PC 32:1, and SM 24:1. All of the 13 metabolites were significant in our analysis for total alcohol consumption (Additional file 7: Fig. S7, Additional file 16: Table. S9).

In a cross-sectional analysis, Würtz et al. examined the association of alcohol intake with 76 metabolites among young adults (aged 25–45) [38]. In that study, 36 metabolites were considered significantly associated with alcohol consumption. Of these 36 metabolites, 11 metabolites (glutamine, glycine, alanine, isoleucine, leucine, valine, lactate, pyruvate, glycerol, citrate, and creatinine) were also included in the present study, and 3 of 11 metabolites (glutamine, glycerol and leucine) reached significance using our data (Additional file 17: Table. S10). Glutamine and glycerol displayed similar association across the two cohorts. Interestingly, compared to our findings, the Würtz et al. study showed an opposite association between total alcohol consumption and leucine. Leucine is one member of the second metabolite score in the present study, which was inversely associated with incident CVD. As such, this observation may echo our abovementioned hypothesis that certain unknown factor(s) may modify the alcohol–metabolite association in different study samples. In addition, six metabolites (glutamine, glycine, alanine, isoleucine, leucine, valine) were included in all three studies (the Würtz et al. study, the van Roekel’s study, and the present study). Four of the six common metabolites were significant in our study and three of them were significant in Würtz et al.’s study, but none of them were significant in van Roekel’s study (Additional file 17: Table. S10). Together, these observations highlight the need for future studies to comprehensively investigate the heterogeneity with respect to alcohol–metabolite associations in different populations using harmonized metabolite panels.

We observed that, for 13 metabolites, their association strength with alcohol consumption was stronger in women compared to that in men. Women generally have smaller body sizes; consuming the same amount of alcohol would end up with a higher blood alcohol concentration for an average woman compared to an average man. In addition, women may have greater ethanol clearance than men, given the same lean body mass [39]. This may, at least partly, explain our observation with respect to the stronger metabolite response in women. Whether the observed sex–alcohol interaction can be affected by other factors, as well as its impact on CVD and clinical outcomes, may need future investigations.

Strengths and limitations of the study

Our study had several strengths. The most important advantage was that alcohol drinking data from five exams across around 20 years were utilized in the current study. Our study had several limitations. Because of the observational nature of our findings and no experimental validation, causality cannot be inferred. Despite that the sample size in our study was large, the population was primarily white, middle-aged participants. Therefore, the findings from our study may not be generalizable to populations of different races and age groups. As shown in Additional file 1: Fig. S1, cumulative average consumption is a good proxy for long-term alcohol consumption. However, the analysis using total alcohol consumption at exam 5 yielded a similar alcohol–metabolite association. This observation may be driven by the possibility that many of our study participants maintained a low-to-moderate alcohol drinking habit. Future studies are needed to validate our findings in other populations with longitudinal alcohol measurements. Alcohol consumption was measured by questionnaires and calculated based on standard serving size. This approach is cost-effective; however, measurement errors may bias the observed association between alcohol intake and metabolite levels. We acknowledge that the metabolites data analyzed in our study are comprised of a targeted panel. Association between alcohol drinking and untargeted metabolites remained to be studied. Research utilizing metabolite platforms including those from exogenous sources, e.g., ethyl glucuronide and derivates of resveratrol, are needed to better understand the alcohol–metabolite relationship. We observed that sex may modify the association between total alcohol consumption and metabolite levels. Several other factors may also modify the observed associations. For example, the presence or absence of food in the stomach can change the rate of alcohol absorption and metabolism [40, 41] and subsequently affect the association of alcohol consumption with circulating metabolite profiles. In addition, the use of average alcohol consumption in the present analysis may not reflect participants’ diverse drinking patterns. For example, the information on alcohol consumption patterns such as drinking alcohol with or without meals or weekend binge drinking was not collected. Participants’ genetic background is a key factor that may modify alcohol metabolism [42], which may need further studies with larger sample size.

Conclusions

The present study identified a series of alcohol-associated circulating metabolites. Our observations suggest that, via some metabolites, alcohol consumption may have counteractive effects on CVD. Future studies are required to validate our findings and to investigate factors that may modify the associations between alcohol consumption and metabolite levels, particularly for metabolites that potentially contribute to CVD risk, in larger and diverse study samples.

Availability of data and materials

The Framingham Heart Study datasets analyzed in the present study are available at the dbGAP Study Accession: phs000007.v33.p14 (https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000007.v33.p14) via the Controlled Access Data.

Abbreviations

CE:

Cholesterol ester

CHD:

Coronary heart disease

CVD:

Cardiovascular disease

DAG:

Diacylglycerol

DASH:

Dietary Approaches to Stop Hypertension

eGFR:

Estimated glomerular filtration rate

FDR:

False discovery rate

FHS:

Framingham Heart Study

HF:

Heart failure

HR:

Hazard ratio

KEGG:

Kyoto Encyclopaedia of Genes and Genomes

LC/MS:

Liquid chromatography with tandem mass spectrometry

LPC:

Lysophosphatidylcholine

LPE:

Lysophosphatidylethanolamine

MI:

Myocardial infarction

PC:

Phosphatidylcholine

SBP:

Systolic blood pressure

SD:

Standard deviation

SE:

Standard error

SM:

Syphingomyelin

SMPDB:

Small Molecule Pathway Database

TAG:

Triacylglycerol

References

  1. Rosoff DB, Davey Smith G, Mehta N, Clarke TK, Lohoff FW. Evaluating the relationship between alcohol consumption, tobacco use, and cardiovascular disease: a multivariable Mendelian randomization study. PLoS Med. 2020;17(12):e1003410.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Sabia S, Fayosse A, Dumurgier J, Dugravot A, Akbaraly T, Britton A, et al. Alcohol consumption and risk of dementia: 23 year follow-up of Whitehall II cohort study. BMJ. 2018;362:k2927.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Kaddurah-Daouk R, Bogdanov MB, Wikoff WR, Zhu H, Boyle SH, Churchill E, et al. Pharmacometabolomic mapping of early biochemical changes induced by sertraline and placebo. Transl Psychiatry. 2013;3:e223.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Roerecke M, Vafaei A, Hasan OSM, Chrystoja BR, Cruz M, Lee R, et al. Alcohol consumption and risk of liver cirrhosis: a systematic review and meta-analysis. Am J Gastroenterol. 2019;114(10):1574–86.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Engler PA, Ramsey SE, Smith RJ. Alcohol use of diabetes patients: the need for assessment and intervention. Acta Diabetol. 2013;50(2):93–9.

    Article  PubMed  Google Scholar 

  6. D’Avanzo B, La Vecchia C, Franceschi S. Alcohol consumption and the risk of gastric cancer. Nutr Cancer. 1994;22(1):57–64.

    Article  CAS  PubMed  Google Scholar 

  7. Grewal P, Viswanathen VA. Liver cancer and alcohol. Clin Liver Dis. 2012;16(4):839–50.

    Article  PubMed  Google Scholar 

  8. Mostofsky E, Chahal HS, Mukamal KJ, Rimm EB, Mittleman MA. Alcohol and immediate risk of cardiovascular events: a systematic review and dose-response meta-analysis. Circulation. 2016;133(10):979–87.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Suhre K, Shin SY, Petersen AK, Mohney RP, Meredith D, Wagele B, et al. Human metabolic individuality in biomedical and pharmaceutical research. Nature. 2011;477(7362):54–60.

    Article  CAS  PubMed  Google Scholar 

  10. Wang YP, Li JT, Qu J, Yin M, Lei QY. Metabolite sensing and signaling in cancer. J Biol Chem. 2020;295(33):11938–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Nicholson JK, Lindon JC, Holmes E. “Metabonomics”: understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data. Xenobiotica. 1999;29(11):1181–9.

    Article  CAS  PubMed  Google Scholar 

  12. van Roekel EH, Trijsburg L, Assi N, Carayol M, Achaintre D, Murphy N, et al. Circulating metabolites associated with alcohol intake in the european prospective investigation into cancer and nutrition cohort. Nutrients. 2018;10(5).

  13. Langenau J, Boeing H, Bergmann MM, Nothlings U, Oluwagbemigun K. The association between alcohol consumption and serum metabolites and the modifying effect of smoking. Nutrients. 2019;11(10).

  14. Jaremek M, Yu Z, Mangino M, Mittelstrass K, Prehn C, Singmann P, et al. Alcohol-induced metabolomic differences in humans. Transl Psychiatry. 2013;3:e276.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Feinleib M, Kannel WB, Garrison RJ, McNamara PM, Castelli WP. The Framingham Offspring Study. Design and preliminary data Prev Med. 1975;4(4):518–25.

    CAS  PubMed  Google Scholar 

  16. Kannel WB, Feinleib M, McNamara PM, Garrison RJ, Castelli WP. An investigation of coronary heart disease in families. The Framingham offspring study. Am J Epidemiol. 1979;110(3):281–90.

  17. Ho JE, Larson MG, Ghorbani A, Cheng S, Chen MH, Keyes M, et al. Metabolomic profiles of body mass index in the Framingham Heart Study reveal distinct cardiometabolic phenotypes. PLoS ONE. 2016;11(2):e0148361.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Rhee EP, Cheng S, Larson MG, Walford GA, Lewis GD, McCabe E, et al. Lipid profiling identifies a triacylglycerol signature of insulin resistance and improves diabetes prediction in humans. J Clin Invest. 2011;121(4):1402–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, et al. Metabolite profiles and the risk of developing diabetes. Nat Med. 2011;17(4):448–53.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Wang TJ, Ngo D, Psychogios N, Dejam A, Larson MG, Vasan RS, et al. 2-Aminoadipic acid is a biomarker for diabetes risk. J Clin Invest. 2013;123(10):4309–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Sun X, Ho JE, Gao H, Evangelou E, Yao C, Huan T, et al. Associations of alcohol consumption with cardiovascular disease-related proteomic biomarkers: the Framingham Heart Study. J Nutr. 2021;151(9):2574–82.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Liu C, Marioni RE, Hedman AK, Pfeiffer L, Tsai PC, Reynolds LM, et al. A DNA methylation biomarker of alcohol consumption. Mol Psychiatry. 2018;23(2):422–33.

    Article  CAS  PubMed  Google Scholar 

  23. Ho JE, Lyass A, Courchesne P, Chen G, Liu C, Yin X, et al. Protein biomarkers of cardiovascular disease and mortality in the community. J Am Heart Assoc. 2018;7(14).

  24. Ho JE, Larson MG, Ghorbani A, Cheng S, Coglianese EE, Vasan RS, et al. Long-term cardiovascular risks associated with an elevated heart rate: the Framingham Heart Study. J Am Heart Assoc. 2014;3(3):e000668.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Tsao CW, Vasan RS. Cohort Profile: The Framingham Heart Study (FHS): overview of milestones in cardiovascular epidemiology. Int J Epidemiol. 2015;44(6):1800–13.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Kim Y, Huan T, Joehanes R, McKeown NM, Horvath S, Levy D, et al. Higher diet quality relates to decelerated epigenetic aging. Am J Clin Nutr. 2022;115(1):163–70.

    Article  PubMed  Google Scholar 

  27. Wilsnack RW, Vogeltanz ND, Wilsnack SC, Harris TR, Ahlstrom S, Bondy S, et al. Gender differences in alcohol consumption and adverse drinking consequences: cross-cultural patterns. Addiction. 2000;95(2):251–65.

    Article  CAS  PubMed  Google Scholar 

  28. Pang Z, Chong J, Zhou G, de Lima Morais DA, Chang L, Barrette M, et al. MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights. Nucleic Acids Research. 2021;49(W1):W388-W96.

  29. The Metabolomics Innovation Center. Showing metabocard for CE(16:1(9Z)) (HMDB0000658): HMDB; 2005 [Available from: https://hmdb.ca/metabolites/HMDB0000658.

  30. Sanchez-Archidona AR, Cruciani-Guglielmacci C, Roujeau C, Wigger L, Lallement J, Denom J, et al. Plasma triacylglycerols are biomarkers of beta-cell function in mice and humans. Mol Metab. 2021;54:101355.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Edwards M, Mohiuddin SS. Biochemistry, lipolysis. StatPearls. Treasure Island (FL)2023.

  32. Baraona E, Lieber CS. Effects of ethanol on lipid metabolism. J Lipid Res. 1979;20(3):289–315.

    Article  CAS  PubMed  Google Scholar 

  33. Ducker GS, Rabinowitz JD. One-carbon metabolism in health and disease. Cell Metab. 2017;25(1):27–42.

    Article  CAS  PubMed  Google Scholar 

  34. Ueland PM, Holm PI, Hustad S. Betaine: a key modulator of one-carbon metabolism and homocysteine status. Clin Chem Lab Med. 2005;43(10):1069–75.

    Article  CAS  PubMed  Google Scholar 

  35. McGarrah RW, White PJ. Branched-chain amino acids in cardiovascular disease. Nat Rev Cardiol. 2023;20(2):77–89.

    Article  CAS  PubMed  Google Scholar 

  36. Neinast M, Murashige D, Arany Z. Branched chain amino acids. Annu Rev Physiol. 2019;81:139–64.

    Article  CAS  PubMed  Google Scholar 

  37. Lever M, Slow S. The clinical significance of betaine, an osmolyte with a key role in methyl group metabolism. Clin Biochem. 2010;43(9):732–44.

    Article  CAS  PubMed  Google Scholar 

  38. Wurtz P, Cook S, Wang Q, Tiainen M, Tynkkynen T, Kangas AJ, et al. Metabolic profiling of alcohol consumption in 9778 young adults. Int J Epidemiol. 2016;45(5):1493–506.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Kwo PY, Ramchandani VA, O’Connor S, Amann D, Carr LG, Sandrasegaran K, et al. Gender differences in alcohol metabolism: relationship to liver volume and effect of adjusting for body mass. Gastroenterology. 1998;115(6):1552–7.

    Article  CAS  PubMed  Google Scholar 

  40. Cederbaum AI. Alcohol metabolism. Clin Liver Dis. 2012;16(4):667–85.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Wilkinson PK, Sedman AJ, Sakmar E, Kay DR, Wagner JG. Pharmacokinetics of ethanol after oral administration in the fasting state. J Pharmacokinet Biopharm. 1977;5(3):207–24.

    Article  CAS  PubMed  Google Scholar 

  42. Saunders GRB, Wang X, Chen F, Jang S-K, Liu M, Wang C, et al. Genetic diversity fuels gene discovery for tobacco and alcohol use. Nature. 2022;612(7941):720–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank the staff and participants of the Framingham Heart Study for their important contributions.

Funding

YL, JM, and CL are supported by R01AA028263. The FHS was supported by NIH contracts N01-HC-25195, HHSN268201500001I, and 75N92019D00031. JIR was supported in part by the National Center for Advancing Translational Sciences, CTSI grant UL1TR001881, and the National Institute of Diabetes and Digestive and Kidney Disease Diabetes Research Center (DRC) grant DK063491 to the Southern California Diabetes Endocrinology Research Center. Infrastructure for the CHARGE Consortium is supported in part by the National Heart, Lung, and Blood Institute (NHLBI) grant R01HL105756.

Author information

Authors and Affiliations

Authors

Contributions

All authors read and approved the final manuscript. Conceptualization: CL and JM. Methodology: CL, JM, and YL. Investigation: CL, JM, and YL. Visualization: CL, JM, and YL. Funding acquisition: CL and JM. Project administration: CL and JM. Supervision: CL and JM. Writing—original draft: YL, CL, and JM. Writing—review and editing: YL, MW, XL, JR, PEM, RJ, TH, XG, JIR, JS, BY, MN, DL, CL, and JM.

Corresponding authors

Correspondence to Chunyu Liu or Jiantao Ma.

Ethics declarations

Ethics approval and consent to participate

All participants provided written informed consent. All FHS study protocols were approved by the Boston University Medical Center Institutional Review Board.

Consent for publication

All authors approve publication of the final manuscript.

Competing interests

We have no conflicts of interest to disclose, and there has been no financial support for this work that could have influence its outcome.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig.S1.

The Pairwise Spearman Correlation Coefficient between Total Alcohol Consumption at All Five FHS Exams.

Additional file 2: Fig. S2.

Comparison of Association of Alcohol Consumption and Metabolites with or without Adjusting for eGFR. Panel A: comparison of beta; Panel B: comparison of -log10 (p value). Two models both were performed adjusting for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect, and family relationship as random effect. eGFR, estimated glomerular filtration rate.

Additional file 3: Fig. S3.

Comparison with Association of Alcohol Consumption at Exam5 and Metabolites. Panel A: comparison of beta; Panel B: comparison of -log10 (p value). Two models both were performed adjusting for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect, and family relationship as random effect.

Additional file 4: Fig. S4.

Number of Metabolites Significantly Associated with Each Type of Alcohol Consumption. Panel A: three type of alcohol consumption; Panel B: three type of alcohol consumption and total alcohol consumption. All models were performed adjusting for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect, and family relationship as random effect.

Additional file 5: Fig. S5.

Comparison of Association Analyses of Alcohol Consumption with Metabolites. Panel A-C, comparison of effect size from association analyses of cumulative average total alcohol consumption and each type of alcohol consumption with metabolites. Panel D-F, comparison of effect size from association analyses each type of alcohol consumption with metabolites. All models were performed adjusting for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect, and family relationship as random effect. The pairwise Pearson correlation coefficients of regression coefficients was 0.61, 0.64, and 0.88 for beer vs. wine, beer vs. liquor, and wine vs. liquor.

Additional file 6: Fig. S6.

Comparisons between nondrinkers, Moderate drinkers, and heavy drinkers. Values are regression coefficients and 95% confidence interval calculated using moderate drinkers as reference. Models were adjusted for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect, and family relationship as random effect.

Additional file 7: Fig. S7.

Comparison to Roekel’s Study. Panel A: comparison of beta; Panel B: comparison of t value. Only in this analysis, alcohol consumption (g/day) in this study was plus 1 and natural log transformed. The sex and batch-adjusted residual of metabolites from linear mixed model were used as outcome. Then we applied linear mixed model for alcohol drinking and residual of metabolites in linear mixed model, adjusting for age, sex, smoking status, BMI, physical activity, diet score as fixed effect, and familial relationship as random effect. But Roekel’s study used linear model and covariates included age at blood collection, sex, country, fasting status at blood collection, smoking status, BMI, Cambridge physical activity index, and daily intake of energy, meat and meat products, fish, and shellfish.

Additional file 8: Table. S1.

Association between Cumulative Average Alcohol Consumption and Metabolites. Total, total amount of alcohol consumption; beer, beer consumption; wine, wine consumption; liquor, liquor consumption. These alcohol consumptions were average consumption across exams 1 to 5. Linear mixed models were performed adjusting for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect and family relationship as random effect.

Additional file 9: Table. S2.

Comparison between Each Type of Cumulative Average Alcohol Consumption and Metabolites. Total, total amount of alcohol consumption; beer, beer consumption; wine, wine consumption; liquor, liquor consumption. These alcohol consumptions were average consumption across exams 1 to 5. Linear mixed models were performed adjusting for age, sex, batch, smoking status, BMI, physical activity index and diet score as fixed effect and family relationship as random effect.

Additional file 10: Table. S3.

Association between Cumulative Average Alcohol Consumption and Metabolites by Sex. Linear mixed models were performed adjusting for age, batch, smoking status, BMI, physical activity index and diet score as fixed effect and family relationship as random effect.

Additional file 11: Table. S4.

Pathway Analysis Results for 60 Metabolites Significantly Associated with Total Alcohol Consumption under KEGG Library.

Additional file 12: Table. S5.

Pathway Analysis Results for 60 Metabolites Significantly Associated with Total Alcohol Consumption under SMPDB Library.

Additional file 13: Table. S6.

Association between Incident CVD and alcohol-related Metabolites. Base models were performed adjusting for age, sex and batch. Multivariable models were performed additionally adjusting for BMI, SBP, hypertension treatment status, diabetes, smoking status, total and high-density lipoprotein cholesterol level.

Additional file 14: Table. S7.

Association between Incident CVD subtypes and alcohol-related Metabolites. Base models were performed adjusting for age, sex and batch. Multivariable models were performed additionally adjusting for BMI, SBP, hypertension treatment status, diabetes, smoking status, total and high-density lipoprotein cholesterol level.

Additional file 15: Table. S8.

Beta Direction from Association between Cumulative Average Alcohol Consumption and Metabolites, Association between Incident CVD and alcohol-related Metabolites.

Additional file 16: Table. S9.

Comparison to Van Roekel's Study. Only in this analysis, our alcohol consumption (g/day) was plus 1 and natural log-transformed. The sex and batch-adjusted residual of metabolites from linear mixed model were used as outcome. Linear mixed models were performed for alcohol drinking and residual of metabolites in linear mixed model, adjusting for age, sex, smoking status, BMI, physical activity, diet score as fixed effect, and familial relationship as random effect.  Refer to reference 12 for details in Van Roekel's study.

Additional file 17: Table. S10.

Comparison to Würtz’s study. Refer to reference 38 for detail statistical analysis in the Wurtz’s study, and significant threshold in their study is 0.0016. And reference 12 for details in Van Roekel's study, and FDR p value<0.05 was significance level.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Wang, M., Liu, X. et al. Circulating metabolites may illustrate relationship of alcohol consumption with cardiovascular disease. BMC Med 21, 443 (2023). https://doi.org/10.1186/s12916-023-03149-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-023-03149-2

Keywords