Plasma Metabolites From Choline Pathway and Risk of Cardiovascular Disease in the PREDIMED (Prevention With Mediterranean Diet) Study
Background The relationship between plasma concentrations of betaine and choline metabolism and major cardiovascular disease (CVD) end points remains unclear. We have evaluated the association between metabolites from the choline pathway and risk of incident CVD and the potential modifying effect of Mediterranean diet interventions.
Methods and Results We designed a case‐cohort study nested within the PREDIMED (Prevention With Mediterranean Diet) trial, including 229 incident CVD cases and 751 randomly selected participants at baseline, followed up for 4.8 years. We used liquid chromatography–tandem mass spectrometry to measure, at baseline and at 1 year of follow‐up, plasma concentrations of 5 metabolites in the choline pathway: trimethylamine N‐oxide, betaine, choline, phosphocholine, and α‐glycerophosphocholine. We have calculated a choline metabolite score using a weighted sum of these 5 metabolites. We used weighted Cox regression models to estimate CVD risk. The multivariable hazard ratios (95% confidence intervals) per 1‐SD increase in choline and α‐glycerophosphocholine metabolites were 1.24 (1.05–1.46) and 1.24 (1.03–1.50), respectively. The baseline betaine/choline ratio was inversely associated with CVD. The baseline choline metabolite score was associated with a 2.21‐fold higher risk of CVD across extreme quartiles (95% confidence interval, 1.36–3.59; P<0.001 for trend) and a 2.27‐fold higher risk of stroke (95% confidence interval, 1.24–4.16; P<0.001 for trend). Participants in the higher quartiles of the score who were randomly assigned to the control group had a higher risk of CVD compared with participants in the lower quartile and assigned to the Mediterranean diet groups (P=0.05 for interaction). No significant associations were observed for 1‐year changes in individual plasma metabolites and CVD.
Conclusions A metabolite score combining plasma metabolites from the choline pathway was associated with an increased risk of CVD in a Mediterranean population at high cardiovascular risk.
Clinical Trial Registration URL: http://www.controlled-trials.com. Unique identifier: ISRCTN35739639.
What Is New?
Baseline plasma concentrations of choline pathway metabolites (choline, phosphocholine, and α‐glycerophosphocholine) and a metabolite score combining choline pathway metabolites were associated with the risk of major cardiovascular events after 4.8 years of follow‐up in a Mediterranean population at high cardiovascular risk.
What Are the Clinical Implications?
Our findings suggest that metabolites from choline pathways may play a role in the development of cardiovascular disease.
The present work may have important public health implications for focusing on dietary interventions, such as adhering to a Mediterranean diet and reducing the intake of animal products, to improve metabolite profiles and, consequently, the risk of cardiovascular disease.
However, further research is needed to replicate these results in other populations and investigate the potential mechanisms underlying the associations between metabolite profiles and cardiovascular disease.
In the past few years, interest has emerged on the role of plasma metabolites in the development of cardiometabolic diseases and on how diet can influence disease risk through inducing changes in the metabolome profile. Lately, metabolites from the choline pathway have received increasing attention for its link with the gut microbiota and the potential mechanistic relevance of them.1 In particular, these studies have focused on metabolites associated with betaine and choline metabolism and the gut microbiota–dependent metabolite, trimethylamine N‐oxide (TMAO), which have been associated with atherosclerosis development and cardiovascular disease (CVD).1, 2 Choline can be oxidized to betaine in humans, and the betaine/choline ratio has been inversely associated with nonalcoholic fatty liver disease,3 but its association with CVD has not been previously described.
Dietary lipid phosphatidylcholine is the primary source of choline, and catabolism of betaine and choline by intestinal microbes leads to TMAO production.4 Foods, such as meat, egg yolks, and high‐fat dairy products, which are rich in phosphatidylcholine, choline, carnitine, and TMA, serve as dietary precursors for TMAO generation in mice and humans.2 Dietary supplementation with these metabolites in mice promoted atherosclerosis, and suppression of the microbiota with antibiotics prevented the effect of dietary choline in enhancing atherosclerosis.4 This evidence suggests that diet could modulate intestinal microbiota composition and, subsequently, the ability to synthesize trimethylamine and TMAO, thus providing a mechanistic link between diet, gut microbiome, and atherosclerosis.1, 2, 5, 6 However, not all human studies have confirmed the association between TMAO and cardiovascular events.7, 8, 9, 10 Indeed, the relation between plasma concentration of choline pathway metabolites and the incidence of hard CVD end points in a population‐based level and the role that diet may play in these associations remain to be elucidated.
We hypothesized that higher concentrations at baseline and 1‐year changes in choline pathway metabolites (TMAO, betaine, choline, phosphocholine, and α‐glycerophosphocholine) would be associated with the risk of cardiovascular events. The Mediterranean diet (MedDiet), which is high in plant‐based foods, has been linked with a reduced risk of cardiovascular risk factors and CVD.11 Because previous evidence suggests that excessive consumption of phosphatidylcholine and choline could be reduced by consuming a plant‐based or high‐fiber diet,4 such as the MedDiet, we hypothesized that the associations between choline metabolites and CVD will be modified by MedDiet interventions. Using a case‐cohort study conducted in participants at high cardiovascular risk from the PREDIMED (Prevention With Mediterranean Diet) Study, we aimed to address the following: (1) whether baseline choline pathway metabolites (and their 1‐year changes) predict future risk of CVD and (2) whether these associations are mitigated by a MedDiet intervention.
We designed a case‐cohort study in the framework of the PREDIMED Study. The design and protocol of the PREDIMED Study (http://www.predimed.es) have been described in detail elsewhere.11, 12 Briefly, the PREDIMED Study was a large, multicenter, parallel‐group, randomized, controlled trial aimed at evaluating the effect of the MedDiet on the primary prevention of CVD. A total of 7447 participants, aged 55 to 80 years and at high cardiovascular risk but free of CVD at baseline, were randomly assigned to receive 1 of the following 3 interventions: MedDiet supplemented with extravirgin olive oil, MedDiet supplemented with mixed nuts, or an equivalent education program with advice to reduce the intake of all types of fat (control group). The primary end point of the PREDIMED trial was a composite of cardiovascular events (myocardial infarction, stroke, or death from cardiovascular causes); 288 incident CVD cases occurred during 4.8 years of follow‐up.
In a case‐cohort study design, all incident cases and a randomly selected sample of the original cohort were selected as participants. In the present study, we included all the incident cases of CVD diagnosed during the follow‐up and a random sample of 10% of the original PREDIMED trial (referred to as the “subcohort,” which by its random selection may include some cases) with available EDTA plasma samples at baseline. With this design, we maintained the original randomization scheme of the trial. The present analysis included a total of 980 participants, 751 were noncases and 229 were cases (there were 37 overlapping cases between the subcohort and the total cases). Of these participants, 923 of the 980 had available samples after 1 year of follow‐up and were included in the 1‐year change analyses (Figure S1). We defined cases as the participants who developed a major cardiovascular event (stroke, myocardial infarction, or cardiovascular death) during follow‐up. We defined subcohort as the random sample selected from the full roster of the PREDIMED Study (including 37 incident CVD cases, which were defined as internal cases). The 37 overlapping cases were treated as cases in the analysis (total CVD cases, 229). We defined external cases as the cases that were not randomly included in the subcohort. All participants provided written informed consent, according to a protocol approved by the institutional review boards before inclusion in the study.
At baseline and at yearly follow‐up visits, trained nurses collected fasting blood samples from the PREDIMED participants. After an overnight fast, tubes for EDTA plasma were collected and aliquots were coded and kept refrigerated until they were stored at −80°C in freezers. Pairs of samples (baseline and first‐year visit) were randomly ordered and shipped on dry ice to the Broad Institute (Cambridge, MA) for the metabolomics analysis.
The present work was a hypothesis‐driven analysis of 5 metabolites from the choline pathway. Polar plasma metabolites, including TMAO, choline, betaine, phosphocholine, and α‐glycerophosphocholine, were profiled using liquid chromatography–tandem mass spectrometry on a system composed of a Shimadzu Nexera X2 U‐HPLC coupled to a Q Exactive hybrid quadrupole orbitrap mass spectrometer.13 Metabolite extracts were prepared from plasma samples (10 μL) via protein precipitation with the addition of 9 volumes of 74.9:24.9:0.2 v/v/v acetonitrile/methanol/formic acid containing stable isotope‐labelled internal standards (valine‐d8 and phenylalanine‐d8). The samples were centrifuged (10 minutes, 9000g, 4°C), and the supernatants were injected directly onto a 150×2‐mm, 3‐μm Atlantis HILIC column. The column was eluted isocratically at a flow rate of 250 μL/min with 5% mobile phase A (10 mmol/L ammonium formate and 0.1% formic acid in water) for 0.5 minutes, followed by a linear gradient to 40% mobile phase B (acetonitrile with 0.1% formic acid) over 10 minutes. Mass spectrometry analyses were performed using electrospray ionization in the positive ion mode using full‐scan analysis over 70 to 800 m/z at 70 000 resolutions and a 3‐Hz data acquisition rate. Other mass spectrometry settings were as follows: sheath gas, 40; sweep gas, 2; spray voltage, 3.5 kV; capillary temperature, 350°C; S‐lens, RF 40; heater temperature, 300°C; microscans, 1; automatic gain control target, 1e6; and maximum ion time, 250 milliseconds. Metabolite identities were confirmed using authentic reference standards. Raw data were processed using TraceFinder software and Progenesis QI. Internal standard peak areas were monitored for quality control and to ensure system performance throughout analyses. Pooled plasma reference samples were also inserted every 20 samples as an additional quality control.
Ascertainment of CVD Cases
For the present analysis, the primary end point was a composite of cardiovascular events (myocardial infarction, stroke, or death from cardiovascular causes). As a secondary end point, we separately analyzed incident stroke, because this was the most common element included in the definition of the composite primary end point in the PREDIMED Study. Information on major cardiovascular events was updated on a yearly basis and obtained from the continuous contact with participants and their families during the trial, contact with general practitioners, the yearly comprehensive review of medical records, and consultation of the National Death Index. The end point was determined by review of the Endpoint Adjudication Committee that was blinded to the intervention group. Only confirmed cases were included in the analyses.
At baseline and at yearly follow‐up visits, a questionnaire about lifestyle variables, educational achievement, history of diseases, medication use, and family history of disease was administered. Physical activity was assessed using the validated Spanish version of the Minnesota Leisure‐Time Physical Activity questionnaire.14 Participants were considered to have diabetes mellitus, hypercholesterolemia, or hypertension if they had previously been diagnosed and/or they were being treated with antidiabetic, cholesterol‐lowering, or antihypertensive agents, respectively. Trained dietitians completed a 137‐item validated semiquantitative food frequency questionnaire in face‐to‐face interviews with participants.15 We used Spanish food composition tables to estimate energy and nutrient intake.16 Trained personnel took anthropometric and blood pressure measurements.
We applied a rank‐based inverse normal transformation to approximate a normal distribution of metabolite levels.17 Baseline characteristics are presented per case status and for quartiles of the choline score as the mean (SD) for quantitative traits and number (percentage) for categorical variables. Baseline characteristics were compared between cases and noncases using t tests for continuous variables and χ2 tests for categorical variables.
We used Cox proportional hazard models, with Barlow weights (to account for oversampling of cases in the study design18), to estimate hazard ratios (HRs) and their 95% confidence intervals (CIs) for the primary combined end point of CVD and, separately, for nonfatal stroke. In all Cox models, Barlow weights were used to account for oversampling of cases in the study design. Noncases in the subcohort were weighted inversely proportional to the sampling fraction (ie, in our study, the weight was 10). For noncases, observation time began on the date of their randomization and ended on their date of death (from noncardiovascular causes), the study end date (December 1, 2010), or the date of the last medical visit (if they had stopped attending scheduled intervention appointments), whichever came first. Cases had a weight of 1 at the instant the individual experienced an end point, and observation time was the time when the participant developed CVD. Incident internal cases were treated as noncases (observation time began on the date that they were randomized and ended the instant the individual died) until they became cases, and at such instant their weight changed.
Cox models were adjusted for age (years), sex (men/women), family history of premature heart disease (yes/no), smoking (never, former, or current), and body mass index (kg/m2), and were stratified by intervention group (both MedDiet interventions and low‐fat control group) (model 1). Model 2 was additionally adjusted for physical activity (metabolic equivalent tasks in min/d), baseline hypertension (yes/no), dyslipidemia (yes/no), and diabetes mellitus (yes/no). Baseline individual plasma metabolites were analyzed as both continuous variables (1‐SD increment in transformed concentration of metabolites) and using quartiles (using cut points defined among noncases). To test for the linear trend across quartiles, the median of each quartile was assigned and analyzed as a continuous variable. To account for multiple testing, we adjusted P values of the multivariable‐adjusted associations between 1‐SD increment in metabolite concentration and CVD risk using the Benjamin‐Hochberg procedure.
In addition, we calculated a choline metabolite score as the weighted sum of concentrations of 5 metabolites from the choline pathway and modeled the score as a main exposure variable in the Cox model. The weight for each metabolite was the regression coefficient for a 1‐SD increment in the plasma concentration estimated from the multivariable Cox regression model. We also calculated the ratio between betaine and choline (by dividing the raw values and then applying inverse normal transformations). We conducted stratified analyses by sex, age (<65 versus ≥65 years), obesity (<30 versus ≥30 kg/m2), smoking status (current/former versus never), family history of premature coronary heart disease, baseline type 2 diabetes mellitus, baseline hypertension, and baseline dyslipidemia. The interactions between these variables of stratification and the choline score were tested by adding multiplicative terms into the multivariable Cox models; a likelihood ratio test was used for testing statistical significance of the interaction term.
We also examined the associations of 1‐year changes in the individual plasma metabolites with CVD. We used the same Cox regression models as in the baseline analyses, but we further adjusted for baseline levels and an interaction term between these metabolites at baseline and 1‐year change (both as continuous variables). We repeated the analysis using incident stroke as the outcome, both including and excluding the nonstroke CVD cases (ie, treating the 113 nonstroke CVD cases as noncases or removing them from analyses), and found similar results. Therefore, we present stroke results only after removing the nonstroke CVD cases.
To examine the combined associations of the MedDiet interventions and choline pathway metabolites with the incidence of CVD, we categorized the participants into joint subgroups, according to the assignment to the intervention group (MedDiet merged versus control group), quartiles of the choline metabolite score, and betaine/choline ratio. We used participants assigned to the control group and the first quartile of the choline score as the reference group; we also used participants in the MedDiet group and the first quartile of the betaine/choline ratio as the reference group. We estimated HRs and their 95% CIs for CVD, according to this joint categorization, using the same models as in the previous analysis. To test the interaction between the MedDiet and the metabolite score and betaine/choline ratio, we added a multiplicative term (1 df) between intervention group (MedDiet merged versus control group) and the continuous score into the multivariable Cox models. In secondary analysis, we further adjusted multivariable model 2 for previous metabolites that predicted CVD in our population, including a weighted score of branched‐chain amino acids (isoleucine, leucine, and valine)19; short‐chain (C2carnitine to C7carnitine), medium‐chain (C8carnitine to C14:2carnitine), and long‐chain (C16carnitine to C26carnitine) acylcarnitines score20; and glutamine/glutamate ratio.21 In addition, we further adjusted the models for medication use (dyslipidemia, antihypertensive, and antidiabetic medications).
All statistical analyses were performed using SAS version 9.4 and R version 2.13.0. P<0.05 was considered statistically significant.
The baseline characteristics of the 980 participants are shown in Table 1, according to whether they had a major CVD event during the 4.8 years of follow‐up and quartiles of the choline metabolite score. Participants in the higher quartile of the score had higher risk profiles at baseline than those in the lower quartile of the score, including older age and higher prevalence of hypertension, dyslipidemia, and diabetes mellitus; they were also more likely to smoke. The characteristics of the metabolites included in the analysis and the pathways in which they are involved are described in Table S1. A diagram of the choline pathway is presented in Figure S2. A heat map of Spearman correlation coefficients of the plasma metabolites analyzed in the present study is shown in Figure S3. Means and SD of plasma metabolite concentrations are presented in Table S2 for the total population and in Table S3 stratified by the intervention group.
Plasma Concentrations of Choline Pathway Metabolites and CVD
The associations between individual plasma metabolites and the choline metabolite score with the risk of CVD and stroke are presented in Table 2. Compared with participants in the lowest quartile of choline, those in the highest quartile had a significantly higher risk of CVD after adjusting for cardiovascular risk factors (HR [95% CI], 1.72 [1.05–2.81]; P=0.01 for trend). A 1‐SD increment in baseline concentrations of choline and α‐glycerophosphocholine was associated with higher risk of CVD incidence. Similar, but weakened, associations were observed when stroke was used as the end point (Table 2). In the multivariable adjusted models, the choline score was associated with a 2.21‐fold higher risk of CVD and a 2.27‐fold higher risk of stroke across quartiles (HR [95% CI], 2.21 [1.36–3.59] [P<0.001 for trend] and 2.27 [1.24–4.16] [P<0.001 for trend], respectively). The betaine/choline ratio was inversely associated with the incidence of CVD (HR [95% CI], 0.57 [0.35–0.92]; P=0.04 for trend) but not with the incidence of stroke alone. The HRs associated with a 1‐SD increment in the choline metabolite score and the betaine/choline ratio were 2.27 (95% CI, 1.36–3.80) and 0.80 (95% CI, 0.68–0.94), respectively. Subgroup analysis for the association between the choline score and the risk of CVD is shown in Table 3. We observed similar associations between the score and CVD after adjusting for baseline plasma concentrations of branched chain amino acids, acylcarnitines, and glutamine/glutamate ratio: the HR per SD was 2.27 (95% CI, 1.27–4.08). Similarly, the HR per SD for betaine/choline ratio after adjustment for these metabolites was 0.77 (95% CI, 0.64–0.92). The associations between 1‐SD increment in choline and α‐glycerophosphocholine remained significant after multiple‐comparison adjustment. The adjusted P values were 0.03 and 0.04, respectively. Further adjusting the models for medication use (dyslipidemia, antihypertensive, and antidiabetic medications) did not change the results. For the overall score, the HR per SD was 2.36 (95% CI, 1.39–4.02); and for the betaine/choline ratio, the HR per SD was 0.82 (95% CI, 0.69–0.96).
Choline Pathway Metabolites and MedDiet Interventions
As illustrated in the Figure, we observed a higher risk of CVD events for those participants who had a higher choline score (quartiles 2–4) and were assigned to the control group compared with those with lower levels in the score and assigned to the MedDiet intervention groups (HR, 2.37; 95% CI, 1.34–4.18). Participants with a higher betaine/choline ratio and assigned to the MedDiet intervention group had a lower risk of CVD compared with those with a lower ratio and in the control group (HR, 0.36; 95% CI, 0.21–0.63) (Figure). The P values for interaction between the continuous choline score, the betaine/choline ratio, and the intervention group (MedDiet versus control group) and CVD were nonsignificant at P=0.05 and P=0.93, respectively (Figure). Likewise, participants with a higher score and in the control group also had a higher risk of stroke (HR, 2.66; 95% CI, 1.22–5.80) (Table S4).
One‐Year Changes in Plasma Concentrations of Metabolites and CVD
No significant associations were observed for quartiles of 1‐year changes in plasma concentrations of individual metabolites and the risk of CVD and stroke. In multivariable model 1, a 1‐SD increment of betaine after 1 year was associated with lower risk of CVD (HR [95% CI], 0.80 [0.66–0.98]), but the association was attenuated in multivariable model 2 and when analyzed as quartiles of the change (Table S5).
In this case‐cohort study nested within the PREDIMED Study, we observed that baseline plasma concentrations of choline pathway metabolites (choline, phosphocholine, and α‐glycerophosphocholine) and a metabolite score combining choline pathway metabolites were associated with the risk of major cardiovascular events after 4.8 years of follow‐up. Higher concentrations of a choline metabolite score were associated with higher risk of CVD and stroke, whereas the ratio of betaine/choline was inversely associated with CVD but not with stroke alone. Therefore, our findings suggest that metabolites from choline pathways may play a role in the development of CVD.
Plasma metabolites in the choline pathway, which are derivatives from dietary sources and may also reflect alterations in the gut microbiota, have been previously related to multifactorial diseases, such as obesity, type 2 diabetes mellitus, and CVD1, 22 Several mechanisms are probably involved, including modulation of host energy metabolism, gut epithelial permeability, gut peptide hormone secretion, and an increase in metabolic endotoxemia and inflammatory status.23, 24, 25 One of the mechanisms explaining the potential associations between some of these metabolites and CVD is the potential proatherogenic effects.6 It has been shown that TMAO can inhibit reverse cholesterol transport, alter sterol metabolism,26 and affect platelet activation and thrombosis risk.27 However, the associations between plasma metabolites within the choline pathway, identified using novel high‐throughput metabolomics techniques, and the incidence of CVD at a population‐based level have only recently been reported.
To our knowledge, this is the first case‐cohort study in a clinical trial setting showing that higher concentrations of plasma choline and α‐glycerophosphocholine are associated with an increased risk of CVD; higher concentrations of choline and phosphocholine are associated with a higher risk of stroke. The associations between dietary and plasma choline and other metabolites from this pathway and CVD events in humans are still unclear. Epidemiologic studies have reported inconsistent associations between systemic levels of choline and betaine and cardiovascular events.6, 7, 28, 29 Recently, it has been reported that higher plasma choline and betaine levels were associated with 1.9‐ and 1.4‐fold increased risks, respectively, of major adverse cardiac events after 3 years of follow‐up in patients undergoing elective diagnostic coronary angiography.29 However, these associations were attenuated and no longer significant when TMAO was included in the models.29 In contrast, betaine supplementation attenuated atherosclerotic lesions in apolipoprotein E−/− mice.30 In our study, the betaine/choline ratio was inversely associated with CVD in individuals at high cardiovascular risk, suggesting that the balance between systemic concentrations of betaine and choline could be of importance for the prevention of cardiovascular events. However, we found no significant associations between 1‐year changes in choline pathway metabolites and the incidence of CVD. It is possible that 1 year was a short period to observe these associations, or other mechanisms and pathways are implicated in the pathogenesis of CVD in our population.
Interestingly, in a cross‐sectional subset of the Nutrition, Aging and Memory in Elders cohort,31 higher plasma choline concentrations were associated with an unfavorable cardiometabolic risk factor profile (lower high‐density lipoprotein cholesterol and higher body mass index). They were also associated with greater odds of large‐vessel cerebral vascular disease and history of CVD but lower odds of small‐vessel cerebral vascular disease. On the contrary, plasma betaine concentrations were associated with a favorable cardiometabolic risk factor profile (lower low‐density lipoprotein cholesterol and triglycerides) and lower odds of diabetes mellitus. In the same study, higher plasma phosphatidylcholine was associated with characteristics of both favorable (higher high‐density lipoprotein cholesterol and lower body mass index) and unfavorable (higher low‐density lipoprotein cholesterol and triglycerides) cardiometabolic risk. These findings suggest divergent roles of choline metabolites in the pathogenesis of cardiometabolic risk factors and cerebrovascular disease.31 Similarly, in our study, we have observed that baseline choline concentrations are marginally associated with the prevalence of being overweight and obesity (odds ratios [95% CIs] per 1‐SD increase, 1.13 [1.00–1.29] and 1.15 [1.00–1.34], respectively). Significant positive correlations, although they were weak, were observed for choline and triglycerides (Spearman correlation coefficient=0.08; P=0.04). We have also observed cross‐sectional inverse associations between plasma betaine and prevalence of diabetes mellitus and being overweight (odds ratios [95% CIs] per 1‐SD increase, 0.64 [0.55–0.74] and 0.88 [0.76–1.06], respectively). Plasma betaine was inversely correlated with plasma glucose and body mass index (Spearman correlation coefficient=−0.18 and −0.07 [P<0.01 and P=0.01], respectively). Altogether, these findings suggested a relationship between these metabolites and cardiovascular risk factors.
The associations between dietary choline and betaine on the risk of CVD have been evaluated in several prospective studies.32, 33, 34 Compared with the lowest quartile of intake, incident coronary heart disease after 14 years of follow‐up in participants from the ARIC (Atherosclerosis Risk in Communities) study was nonsignificantly higher in the highest quartile of choline and choline plus betaine (HRs [95% CIs], 1.22 [0.91–1.64] and 1.14 [0.85–1.53], respectively), compared with the lowest quartile.32 In women from the Jackson Heart Study, greater choline intake was associated with lower left ventricular mass (P=0.0006 for trend across choline quartiles) and with abdominal aortic calcium score.34 Similarly, higher dietary intakes of betaine and choline were not associated with CVD risk in postmenopausal Dutch women followed up for 97 months.33 On the other hand, findings from the Nurses' Health Study and Health Professionals Follow‐Up Study, including 80 978 women and 39 434 men, respectively, followed up to 32 years, suggested that higher phosphatidylcholine intake was associated with an increased risk of CVD mortality after adjustment for potential confounders.35
Previous published literature on this topic has mostly focused on the formation of the gut microbiota–dependent metabolite, TMAO.2, 4, 6 One of the most compelling studies linking these metabolites with CVD complications was reported by Wang and collaborators a few years ago.6 The authors identified a novel pathway linking dietary phosphatidylcholine/choline, which can be converted by the intestinal microbiota into trimethylamine and subsequently converted into TMAO by hepatic flavin‐containing mono‐oxygenases, with CVD pathogenesis.13 These findings are consistent with later in vivo and in vitro studies, suggesting that choline metabolites may increase the risk of atherosclerosis and coronary heart disease.4, 29, 36, 37 In addition, the ability of oral broad‐spectrum antibiotics to temporarily suppress the production of TMAO suggests that intestinal microorganisms may play an important role in the production of TMAO from phosphatidylcholine in humans.4 Nevertheless, our data did not support positive significant associations between plasma concentrations of TMAO and the risk of CVD in individuals at high cardiovascular risk after 4.8 years of follow‐up. It is plausible that TMAO may be more relevant for predicting recurrent CVD or CVD survival, being a marker of established disease. Although TMAO was associated with CVD mortality in participants undergoing incident dialysis in a nested case‐control study, when they replicated the findings in an independent study, no differences were found in TMAO concentrations between cases of cardiovascular death and controls.7 One of the potential reasons for these associations is that elevated TMAO levels are strongly related with renal function in patients with chronic kidney disease, and they normalize after renal transplantation.38 Therefore, because kidney disease is a risk factor for CVD, TMAO was probably associated with CVD mortality because individuals undergoing dialysis have more risk of developing CVD.
Findings from a recent study, including 817 participants of the CARDIA (Coronary Artery Risk Development in Young Adults), followed up for 10 years, showed that TMAO was not associated with coronary artery calcium incidence or other measures of CVD risk, including insulin resistance, inflammatory markers, and lipid profile. The authors suggest that it is possible that the adverse effect of TMAO may be more relevant when concentrations of TMAO are higher than in their sample or in later stages of the disease process.10 In another study comparing the concentrations of TMAO, carnitine, and choline in 34 obese individuals undergoing bariatric surgery, TMAO was not elevated in obese patients before the surgery but increased ≈2‐fold after bariatric surgery.39 The authors reported that these results were unexpected and suggested that one explanation could be adaptive shifts in the gut microbiota with increased ability to metabolize dietary choline and carnitine to TMAO precursors.39 In contrast, carnitine and choline, which are abundant in nutrients, such as in red meat and eggs, and not microbiota dependent, were reduced after lifestyle interventions and rebounded after bariatric surgery.39
Along these lines, choline and its metabolites can be obtained from food sources or synthesized de novo. In a randomized, controlled, crossover trial, it was observed that dimethylamine and TMAO were significantly higher after consuming a diet high in fish.40 The authors reported that plasma and urinary concentrations of TMAO were associated with recent fish intake; therefore, high urinary TMAO concentrations can be associated with healthy diets that are rich in fish.40 Our data also showed that plasma TMAO is correlated with the intake of fish and total meat (Spearman correlation coefficients, 0.12 and 0.08, respectively; P<0.01). However, gut bacteria can synthesize TMAO from choline and, hence, high urinary and plasma TMAO concentrations can also originate from red meat consumption, which is generally associated with adverse health outcomes. Thus, the global pattern of metabolites, which reflects the totality of the diet, might be more important in indicating dietary patterns than are individual biomarkers.40 Findings from Tang et al suggest that excessive consumption of phosphatidylcholine and choline should be avoided by consuming a vegetarian or high‐fiber diet,4 which has similar characteristics as the MedDiet intervention in the PREDIMED Study. It is established that foods such as meat, egg yolks, and high‐fat dairy products, which are high in phosphatidylcholine, choline, carnitine, and TMA, can serve as precursors of TMAO production.6
Another point that has not been raised before and deserves consideration is the possible attenuation of detrimental effects of metabolite profiles on CVD risk by adhering to a MedDiet. Multiple lines of evidence could account for the observed benefits of the MedDiet on CVD, including the reduction of low‐grade inflammation,41 enhanced endothelial function,42 lower oxidative stress,43 and lower levels of oxidized low‐density lipoprotein.44 Integrating diet information with metabolomics data holds great potential for further understanding biomarkers and mechanisms that relate diet, health, and disease. The present work may have important public health implications for focusing on dietary interventions, such as adhering to the MedDiet and reducing the intake of animal products, to improve metabolite profiles and, consequently, the risk of CVD.
Several limitations and strengths of the present study warrant mentioning. First, because the participants were elderly Mediterranean individuals at high cardiovascular risk, the results may not be extrapolated to other populations. Notwithstanding, because of their increased risk of CVD, they are an important therapeutic target. In addition, the case‐cohort design maximized the efficiency of the high‐throughput metabolomics profiling, which allowed us to extend the results to all PREDIMED participants. Second, although we adjusted for several potential confounders, residual confounding by other unknown or unmeasured variables cannot be ruled out. Third, liquid chromatography–tandem mass spectrometry–based metabolite measurements may not have a direct clinical translation for each metabolite trait. Finally, in the present analysis, we have only included choline pathway metabolites, as subrogates of diet and gut microbiota–related metabolites, that have been identified in our samples with the use of the metabolomics platform from the present project. However, in the coming years, new metabolites will probably need to be included in this group. The strengths of the present work include the prospective design, the ability to control for potential confounders because of the well‐characterized comprehensive data of the randomized PREDIMED trial, and the accurate and blind assessment of incident CVD cases.
In conclusion, our findings demonstrated, for the first time, that metabolites from the choline pathway were associated with an increased risk of major adverse cardiovascular events, independent of traditional risk factors, in a Mediterranean population at high cardiovascular risk.
Guasch‐Ferré, Hu, Martínez‐González, and Salas‐Salvadó conceived and designed research; Guasch‐Ferré and Liang performed statistical analysis; Hu, Martínez‐González, and Salas‐Salvadó handled funding and supervision; Guasch‐Ferré, Hu, Ruiz‐Canela, Bulló, Toledo, Wang, Corella, Gómez‐Gracia, Fiol, Estruch, Lapetra, Fitó, Arós, Majem, Ros, Dennis, Liang, Clish, Martínez‐González, and Salas‐Salvadó acquired the data; Guasch‐Ferré and Salas‐Salvadó drafted the article Hu, Ruiz‐Canela, Bulló, Toledo, Wang, Corella, Gómez‐Gracia, Fiol, Estruch, Lapetra, Fitó, Arós, Majem, Ros, Liang, Clish, Martínez‐González, and Salas‐Salvadó critically revised the article for key intellectual content.
Sources of Funding
This work was supported by NIH research grant HL118264. The PREDIMED trial was supported by the official funding agency for biomedical research of the Spanish Government, Instituto de Salud Carlos III, through grants provided to research networks specifically developed for the trial (RTIC G03/140 to Estruch; RTIC RD 06/0045 to Martínez‐González) and through Centro de Investigación Biomédica en Red de Fisiopatología de la Obesidad y Nutrición; and by grants from Centro Nacional de Investigaciones Cardiovasculares (06/2007), Fondo de Investigación Sanitaria–Fondo Europeo de Desarrollo Regional (PI04–2239, PI05/2584, CP06/00100, PI07/0240, PI07/1138, PI07/0954, PI 07/0473, PI10/01407, PI10/02658, PI11/01647, P11/02505, and PI13/00462), Ministerio de Ciencia e Innovación (AGL‐2009–13906‐C02 and AGL2010–22319‐C03), Fundación Mapfre 2010, Consejería de Salud de la Junta de Andalucía (PI0105/2007), Public Health Division of the Department of Health of the Autonomous Government of Catalonia, Generalitat Valenciana (ACOMP06109, GVA‐COMP2010–181, GVACOMP2011–151, CS2010‐AP‐111, and CS2011‐AP‐042), and Regional Government of Navarra (P27/2011). Guasch‐Ferré was supported by a postdoctoral fellowship granted by the Agency for Administration of University and Research Grants of the Autonomous Government of Catalonia (2014‐BP‐A 00017). None of the funding sources played a role in the design, collection, analysis, or interpretation of the data or in the decision to submit the article for publication.
Hu and Ros received grants from the California Walnut Commission. Martínez‐González and Salas‐Salvadó received grants from the International Nut Council. None of the other authors has a relevant conflict of interest to report related to this research.
Table S1. Characteristics of the Choline Pathway Metabolites
Table S2. Means and SD of Baseline Plasma Concentration of Metabolites (μmol/L) in the Total Population
Table S3. Means and SD of Plasma Concentration of Metabolites (μmol/L) Per Intervention Group at Baseline and 1 Year of Intervention
Table S4. Joint Associations of the Scores and Intervention Group With Risk of Cardiovascular Disease and Stroke
Table S5. Risk of Cardiovascular Disease and Stroke by 1‐Year Changes in Choline Pathway Metabolites in the PREDIMED Study
Figure S1. Flow chart of study participants.
Figure S2. Diagram of choline pathway.
Figure S3. Correlation matrix for plasma metabolite levels.
We thank all the participants for their collaboration, all the PREDIMED personnel for their assistance, and all the personnel of affiliated primary care centers for making the study possible.
- ↵Griffin JL, Wang X, Stanley E. Does our gut microbiome predict cardiovascular risk? A review of the evidence from metabolomics. Circ Cardiovasc Genet. 2015;8:187–191.
- ↵Wang Z, Roberts AB, Buffa JA, Levison BS, Zhu W, Org E, Gu X, Huang Y, Zamanian‐Daryoush M, Culley MK, DiDonato AJ, Fu X, Hazen JE, Krajcik D, DiDonato JA, Lusis AJ, Hazen SL. Non‐lethal inhibition of gut microbial trimethylamine production for the treatment of atherosclerosis. Cell. 2015;163:1585–1595.
- ↵Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, Feldstein AE, Britt EB, Fu X, Chung Y‐M, Wu Y, Schauer P, Smith JD, Allayee H, Tang WHW, DiDonato JA, Lusis AJ, Hazen SL. Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. 2011;472:57–63.
- ↵Kalim S, Clish CB, Wenger J, Elmariah S, Yeh RW, Deferio JJ, Pierce K, Deik A, Gerszten RE, Thadhani R, Rhee EP. A plasma long‐chain acylcarnitine predicts cardiovascular mortality in incident dialysis patients. J Am Heart Assoc. 2013;2:e000542. DOI: 10.1161/JAHA.113.000542.
- ↵Vaarhorst AAM, Verhoeven A, Weller CM, Böhringer S, Göraler S, Meissner A, Deelder AM, Henneman P, Gorgels APM, van den Brandt PA, Schouten LJ, van Greevenbroek MM, Merry AHH, Verschuren WMM, van den Maagdenberg AMJM, van Dijk KW, Isaacs A, Boomsma D, Oostra BA, van Duijn CM, Jukema JW, Boer JMA, Feskens E, Heijmans BT, Slagboom PE. A metabolomic profile is associated with the risk of incident coronary heart disease. Am Heart J. 2014;168:45–52.e7.
- ↵Yin J, Liao S‐X, He Y, Wang S, Xia G‐H, Liu F‐T, Zhu J‐J, You C, Chen Q, Zhou L, Pan S‐Y, Zhou H‐W. Dysbiosis of gut microbiota with reduced trimethylamine‐N‐oxide level in patients with large‐artery atherosclerotic stroke or transient ischemic attack. J Am Heart Assoc. 2015;4:e002699. DOI: 10.1161/JAHA.115.002699.
- ↵Meyer KA, Benton TZ, Bennett BJ, Jacobs DR, Lloyd‐Jones DM, Gross MD, Carr JJ, Gordon‐Larsen P, Zeisel SH. Microbiota‐dependent metabolite trimethylamine N‐oxide and coronary artery calcium in the coronary artery risk development in young adults study (CARDIA). J Am Heart Assoc. 2016;5:e003970. DOI: 10.1161/JAHA.116.003970.
- ↵Estruch R, Ros E, Salas‐Salvado J, Covas MI, Corella D, Aros F, Gomez‐Gracia E, Ruiz‐Gutierrez V, Fiol M, Lapetra J, Lamuela‐Raventos RM, Serra‐Majem L, Pinto X, Basora J, Munoz MA, Sorli JV, Martinez JA, Martinez‐Gonzalez MA; PREDIMED Study Investigators . Primary prevention of cardiovascular disease with a Mediterranean diet. N Engl J Med. 2013;368:1279–1290.
- ↵Martinez‐Gonzalez MA, Corella D, Salas‐Salvado J, Ros E, Covas MI, Fiol M, Warnberg J, Aros F, Ruiz‐Gutierrez V, Lamuela‐Raventos RM, Lapetra J, Munoz MA, Martinez JA, Saez G, Serra‐Majem L, Pinto X, Mitjavila MT, Tur JA, Portillo Mdel P, Estruch R; PREDIMED Study Investigators . Cohort profile: design and methods of the PREDIMED study. Int J Epidemiol. 2012;41:377–385.
- ↵Fernandez‐Ballart JD, Pinol JL, Zazpe I, Corella D, Carrasco P, Toledo E, Perez‐Bauer M, Martinez‐Gonzalez MA, Salas‐Salvado J, Martin‐Moreno JM. Relative validity of a semi‐quantitative food‐frequency questionnaire in an elderly Mediterranean population of Spain. Br J Nutr. 2010;103:1808–1816.
- ↵Mataix J. Tablas de Composición de Alimentos. 4th ed. Granada: Universidad de Granada; 2003.
- ↵Ruiz‐Canela M, Toledo E, Clish CB, Hruby A, Liang L, Salas‐Salvado J, Razquin C, Corella D, Estruch R, Ros E, Fito M, Gomez‐Gracia E, Aros F, Fiol M, Lapetra J, Serra‐Majem L, Martinez‐Gonzalez MA, Hu FB. Plasma branched‐chain amino acids and incident cardiovascular disease in the PREDIMED trial. Clin Chem. 2016;62:582–592.
- ↵Guasch‐Ferré M, Zheng Y, Ruiz‐Canela M, Hruby A, Martínez‐González MA, Clish CB, Corella D, Estruch R, Ros E, Fitó M, Dennis C, Morales‐Gil IM, Arós F, Fiol M, Lapetra J, Serra‐Majem L, Hu FB, Salas‐Salvadó J. Plasma acylcarnitines and risk of cardiovascular disease: effect of Mediterranean diet interventions. Am J Clin Nutr. 2016;103:1408–1416.
- ↵Zheng Y, Hu FB, Ruiz‐Canela M, Clish CB, Dennis C, Salas‐Salvado J, Hruby A, Liang L, Toledo E, Corella D, Ros E, Fitó M, Gómez‐Gracia E, Arós F, Fiol M, Lapetra J, Serra‐Majem L, Estruch R, Martínez‐González MA. Metabolites of glutamate metabolism are associated with incident cardiovascular events in the PREDIMED PREvención con DIeta MEDiterránea (PREDIMED) Trial. J Am Heart Assoc. 2016;5:e003755. DOI: 10.1161/JAHA.116.003755.
- ↵Bäckhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI. Host‐bacterial mutualism in the human intestine. Science. 2005;307:1915–1920.
- ↵Samuel BS, Shaito A, Motoike T, Rey FE, Backhed F, Manchester JK, Hammer RE, Williams SC, Crowley J, Yanagisawa M, Gordon JI. Effects of the gut microbiota on host adiposity are modulated by the short‐chain fatty‐acid binding G protein‐coupled receptor, Gpr41. Proc Natl Acad Sci U S A. 2008;105:16767–16772.
- ↵Koeth RA, Wang Z, Levison BS, Buffa JA, Org E, Sheehy BT, Britt EB, Fu X, Wu Y, Li L, Smith JD, DiDonato JA, Chen J, Li H, Wu GD, Lewis JD, Warrier M, Brown JM, Krauss RM, Tang WHW, Bushman FD, Lusis AJ, Hazen SL. Intestinal microbiota metabolism of L‐carnitine, a nutrient in red meat, promotes atherosclerosis. Nat Med. 2013;19:576–585.
- ↵Zhu W, Gregory JC, Org E, Buffa JA, Gupta N, Wang Z, Li L, Fu X, Wu Y, Mehrabian M, Sartor RB, McIntyre TM, Silverstein RL, Tang WHW, DiDonato JA, Brown JM, Lusis AJ, Hazen SL. Gut microbial metabolite TMAO enhances platelet hyperreactivity and thrombosis risk. Cell. 2016;165:111–124.
- ↵Roe AJ, Zhang S, Bhadelia RA, Johnson EJ, Lichtenstein AH, Rogers GT, Rosenberg IH, Smith CE, Zeisel SH, Scott TM. Choline and its metabolites are differently associated with cardiometabolic risk factors, history of cardiovascular disease, and MRI‐documented cerebrovascular disease in older adults. Am J Clin Nutr. 2017;105:1283–1290. DOI: 10.3945/ajcn.116.137158.
- ↵Millard HR, Musani SK, Dibaba DT, Talegawkar SA, Taylor HA, Tucker KL, Bidulescu A. Dietary choline and betaine; associations with subclinical markers of cardiovascular disease risk and incidence of CVD, coronary heart disease and stroke: the Jackson Heart Study. Eur J Nutr. 2016. Available at: https://jhu.pure.elsevier.com/en/publications/dietary-choline-and-betaine-associations-with-subclinical-markers. Accessed October 6, 2017.
- ↵Zheng Y, Li Y, Rimm EB, Hu FB, Albert CM, Rexrode KM, Manson JE, Qi L. Dietary phosphatidylcholine and risk of all‐cause and cardiovascular‐specific mortality among US women and men. Am J Clin Nutr. 2016;104:173–180.
- ↵Senthong V, Wang Z, Li XS, Fan Y, Wu Y, Tang WHW, Hazen SL. Intestinal microbiota‐generated metabolite trimethylamine‐N‐oxide and 5‐year mortality risk in stable coronary artery disease: the contributory role of intestinal microbiota in a COURAGE‐like patient cohort. J Am Heart Assoc. 2016;5:e002816. DOI: 10.1161/JAHA.115.002816.
- ↵Senthong V, Li XS, Hudec T, Coughlin J, Wu Y, Levison B, Wang Z, Hazen SL, Tang WHW. Plasma trimethylamine N‐oxide, a gut microbe‐generated phosphatidylcholine metabolite, is associated with atherosclerotic burden. J Am Coll Cardiol. 2016;67:2620–2628.
- ↵Trøseid M, Hov JR, Nestvold TK, Thoresen H, Berge RK, Svardal A, Lappegård KT. Major increase in microbiota‐dependent proatherogenic metabolite TMAO one year after bariatric surgery. Metab Syndr Relat Disord. 2016;14:197–201.
- ↵Garcia‐Perez I, Posma JM, Gibson R, Chambers ES, Hansen TH, Vestergaard H, Hansen T, Beckmann M, Pedersen O, Elliott P, Stamler J, Nicholson JK, Draper J, Mathers JC, Holmes E, Frost G. Objective assessment of dietary patterns by use of metabolic phenotyping: a randomised, controlled, crossover trial. Lancet Diabetes Endocrinol. 2017;5:184–195.
- ↵Esposito K, Marfella R, Ciotola M, Di Palo C, Giugliano F, Giugliano G, D'Armiento M, D'Andrea F, Giugliano D. Effect of a Mediterranean‐style diet on endothelial dysfunction and markers of vascular inflammation in the metabolic syndrome: a randomized trial. JAMA. 2004;292:1440–1446.
- ↵Fitó M, Guxens M, Corella D, Sáez G, Estruch R, de la Torre R, Francés F, Cabezas C, López‐Sabater MDC, Marrugat J, García‐Arellano A, Arós F, Ruiz‐Gutierrez V, Ros E, Salas‐Salvadó J, Fiol M, Solá R, Covas M‐I; PREDIMED Study Investigators . Effect of a traditional Mediterranean diet on lipoprotein oxidation: a randomized controlled trial. Arch Intern Med. 2007;167:1195–1203.