Skip to main content
  • More from ADA
    • Diabetes
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care
  • Subscribe
  • Log in
  • My Cart
  • Follow ada on Twitter
  • RSS
  • Visit ada on Facebook
Diabetes Care

Advanced Search

Main menu

  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • Special Article Collections
    • ADA Standards of Medical Care
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • Special Article Collections
    • ADA Standards of Medical Care
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Journal Policies
    • Instructions for Authors
    • ADA Peer Review
  • More from ADA
    • Diabetes
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care

User menu

  • Subscribe
  • Log in
  • My Cart

Search

  • Advanced search
Diabetes Care
  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • Special Article Collections
    • ADA Standards of Medical Care
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • Special Article Collections
    • ADA Standards of Medical Care
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Journal Policies
    • Instructions for Authors
    • ADA Peer Review
Cardiovascular and Metabolic Risk

Plasma Lipidomic Profiling and Risk of Type 2 Diabetes in the PREDIMED Trial

  1. Cristina Razquin1,2,3,
  2. Estefanía Toledo1,2,3,
  3. Clary B. Clish4,
  4. Miguel Ruiz-Canela1,2,3,
  5. Courtney Dennis4,
  6. Dolores Corella3,5,
  7. Christopher Papandreou3,6,
  8. Emilio Ros3,7,
  9. Ramon Estruch3,8,
  10. Marta Guasch-Ferré3,6,9,
  11. Enrique Gómez-Gracia3,10,
  12. Montserrat Fitó3,11,
  13. Edward Yu9,
  14. José Lapetra3,12,
  15. Dong Wang9,
  16. Dora Romaguera3,13,
  17. Liming Liang14,15,
  18. Angel Alonso-Gómez3,16,
  19. Amy Deik4,
  20. Mónica Bullo3,6,
  21. Lluis Serra-Majem3,17,
  22. Jordi Salas-Salvadó3,6,
  23. Frank B. Hu9,18 and
  24. Miguel A. Martínez-González1,2,3,9⇑
  1. 1Department of Preventive Medicine and Public Health, University of Navarra, Pamplona, Spain
  2. 2IdiSNA, Navarra Institute for Health Research, Pamplona, Spain
  3. 3CIBER Fisiopatología de la Obesidad y Nutrición (CIBERObn), Instituto de Salud Carlos III, Madrid, Spain
  4. 4Broad Institute of MIT and Harvard University, Cambridge, MA
  5. 5Department of Preventive Medicine, University of Valencia, Valencia, Spain
  6. 6Human Nutrition Unit, Faculty of Medicine and Health Sciences, Institut d’Investigació Sanitària Pere Virgili, Rovira i Virgili University, Reus, Spain
  7. 7Lipid Clinic, Department of Endocrinology and Nutrition, Institut d’Investigacions Biomediques August Pi Sunyer (IDI-BAPS), Hospital Clinic, University of Barcelona, Barcelona, Spain
  8. 8Department of Internal Medicine, Institut d’Investigacions Biomediques August Pi Sunyer (IDI-BAPS), Barcelona, Spain
  9. 9Department of Nutrition, Harvard T.H. Chan School of Public Health, Boston, MA
  10. 10Department of Preventive Medicine, University of Malaga, Malaga, Spain
  11. 11Cardiovascular and Nutrition Research Group, Institut de Recerca Hospital del Mar, Barcelona, Spain
  12. 12Research Unit, Department of Family Medicine, Distrito Sanitario Atención Primaria Sevilla, Seville, Spain
  13. 13Instituto de Investigación Sanitaria Illes Balears (IdISBa), University Hospital of Son Espases, Palma de Mallorca, Spain
  14. 14Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA
  15. 15Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
  16. 16Department of Cardiology, University Hospital of Alava, Vitoria, Spain
  17. 17Research Institute of Biomedical and Health Sciences (IUIBS), University of Las Palmas de Gran Canaria, and Service of Preventive Medicine, Complejo Hospitalario Universitario Insular Materno Infantil (CHUIMI), Canary Health Service, Las Palmas de Gran Canaria, Spain
  18. 18Channing Division for Network Medicine, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA
  1. Corresponding author: Miguel A. Martínez-González, mamartinez{at}unav.es.
Diabetes Care 2018 Dec; 41(12): 2617-2624. https://doi.org/10.2337/dc18-0840
PreviousNext
  • Article
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF
Loading

Abstract

OBJECTIVE Specific lipid molecular changes leading to type 2 diabetes (T2D) are largely unknown. We assessed lipidome factors associated with future occurrence of T2D in a population at high cardiovascular risk.

RESEARCH DESIGN AND METHODS We conducted a case-cohort study nested within the PREDIMED trial, with 250 incident T2D cases diagnosed during 3.8 years of median follow-up, and a random sample of 692 participants (639 noncases and 53 overlapping cases) without T2D at baseline. We repeatedly measured 207 plasma known lipid metabolites at baseline and after 1 year of follow-up. We built combined factors of lipid species using principal component analysis and assessed the association between these lipid factors (or their 1-year changes) and T2D incidence.

RESULTS Baseline lysophosphatidylcholines and lysophosphatidylethanolamines (lysophospholipids [LPs]), phosphatidylcholine-plasmalogens (PC-PLs), sphingomyelins (SMs), and cholesterol esters (CEs) were inversely associated with risk of T2D (multivariable-adjusted P for linear trend ≤0.001 for all). Baseline triacylglycerols (TAGs), diacylglycerols (DAGs), and phosphatidylethanolamines (PEs) were positively associated with T2D risk (multivariable-adjusted P for linear trend <0.001 for all). One-year changes in these lipids showed associations in similar directions but were not significant after adjustment for baseline levels. TAGs with odd-chain fatty acids showed inverse associations with T2D after adjusting for total TAGs.

CONCLUSIONS Two plasma lipid profiles made up of different lipid classes were found to be associated with T2D in participants at high cardiovascular risk. A profile including LPs, PC-PLs, SMs, and CEs was associated with lower T2D risk. Another profile composed of TAGs, DAGs, and PEs was associated with higher T2D risk.

Introduction

Type 2 diabetes (T2D), which results from insulin resistance and inadequate compensatory insulin secretion (1,2), currently has a global prevalence of 9%, which is projected to rise to 10.4% (642 million cases) by 2040 (3).

Dyslipidemia, characterized by a high plasma triglyceride concentration, low HDL cholesterol concentration, and increased concentration of small dense LDL cholesterol particles, is usually present in T2D (4,5). Triglycerides encompass a large number of individual molecular species, whereas lipoproteins include many different lipid classes containing multiple molecular species (6). The role that these individual molecular lipid species play in T2D development remains unclear. Hypercaloric and low-quality diets also contribute to T2D by leading to an excess of fat depositions, which is enhanced by insulin resistance, resulting in lipotoxicity (7). Lipidomics may help to clarify the biological mechanisms underlying the link between dyslipidemia, nutrition, and T2D.

The PREDIMED trial (Prevención con Dieta Mediterránea), which assessed a Mediterranean diet (MedDiet) intervention, provides an opportunity to identify lipidome profiles associated with T2D and to discern if the intervention changed the lipidome determining the risk of T2D.

The aims of the current study were to 1) assess lipidome patterns associated with subsequent risk of incident T2D, 2) analyze if 1-year changes in these lipid patterns induced by the dietary intervention were associated with subsequent T2D risk, and 3) evaluate if the protective effects of the intervention on T2D were partially explained by changes in the lipidome.

Research Design and Methods

We used an unstratified case-cohort study nested within the PREDIMED trial (www.predimed.es), a primary cardiovascular prevention trial testing Mediterranean diets, as described elsewhere (8,9). In brief, 7,447 participants (men aged 55–80 years and women aged 60–80 years), initially free of cardiovascular disease (CVD) but at high cardiovascular risk, were allocated to three dietary interventions: 1) a Mediterranean diet supplemented with extra virgin olive oil (MedDiet + EVOO), 2) a Mediterranean diet supplemented with mixed nuts (MedDiet + nuts), or 3) a control diet (low-fat diet). In the full PREDIMED cohort, 3,541 participants did not have T2D at baseline. Among these participants, there were 273 incident cases of T2D observed during follow-up. Participants randomized to MedDiet groups, and especially to the MedDiet + EVOO group, had a significantly lower risk of T2D than the control group (10).

The present case-cohort study comprises a random selection of 694 participants (∼20%) from the eligible subjects of the PREDIMED cohort (those with available EDTA plasma samples and who did not have T2D at baseline), together with all incident cases of T2D that occurred during a median follow-up of 3.8 years of intervention. Lipid metabolites were measured for 889 participants in the full PREDIMED cohort. The subcohort used in this study included 639 noncases and 53 overlapping cases. There were an additional 197 T2D cases, yielding a total of 250 incident cases (Supplementary Fig. 1). In addition, 658 participants (501 noncases and 157 cases that occurred after 1 year of follow-up) had 1-year follow-up samples and were included in the 1-year change analyses (Supplementary Fig. 1).

The institutional review boards of the recruitment centers approved the study protocol, and participants provided written informed consent.

Covariate Assessment

At baseline and at yearly follow-up visits, participants completed a questionnaire collecting lifestyle information, educational achievement, history of illnesses, medication use, and family history of disease. Physical activity was assessed using the validated Spanish version of the Minnesota Leisure-Time Physical Activity Questionnaire (11).

Study Samples and Metabolite Profiling

Fasting blood samples were collected at baseline and after 1 year of follow-up. After an overnight fast, plasma EDTA tubes were collected and aliquots were coded and kept refrigerated until they were stored at −80°C. In June 2015, pairs of samples (baseline and 1st-year visits from each participant) were randomly ordered and shipped on dry ice to the Broad Institute for the metabolomics analyses. Specifically, 207 plasma polar and nonpolar lipids were profiled using a Nexera x2 U-HPLC system (Shimadzu Scientific Instruments, Marlborough, MA) coupled to an Exactive Plus orbitrap mass spectrometer (Thermo Fisher Scientific, Waltham, MA). Lipids were extracted from plasma (10 µL) using 190 µL of isopropanol containing 1,2-didodecanoyl-sn-glycero-3-phosphocholine as an internal standard (Avanti Polar Lipids, Alabaster, AL). After centrifugation (10 min, 9,000g, ambient temperature), supernatants (2 µL) were injected directly onto a 100 × 2.1 mm ACQUITY BEH C8 column (1.7 µm) (Waters, Milford, MA). The column was eluted at a flow rate of 450 µL/min isocratically for 1 min at 80% mobile phase A (95:5:0.1 volume for volume for volume [v/v/v] 10 mmol/L ammonium acetate/methanol/acetic acid), followed by a linear gradient to 80% mobile phase B (99.9:0.1 v/v methanol/acetic acid) over 2 min, a linear gradient to 100% mobile phase B over 7 min, and then 3 min at 100% mobile phase B. Mass spectometry analyses were performed using electrospray ionization in the positive ion mode using full scan analysis over m/z 200–1,100 at 70,000 resolution and 3 Hz data acquisition rate. Additional mass spectrometry settings were as follows: ion spray voltage, 3.0 kV; capillary temperature, 300°C; probe heater temperature, 300°C; sheath gas, 50; auxiliary gas, 15; and S-lens RF level, 60. Raw data were processed using Progenesis QI software (NonLinear Dynamics) for feature alignment, nontargeted signal detection, and signal integration. Targeted processing of a subset of lipids was conducted using TraceFinder software (version 3.2; Thermo Fisher Scientific). Lipids were denoted by headgroup and total acyl carbon content and total acyl double bond content.

Clinical Assessment

The PREDIMED protocol included T2D as a prespecified secondary end point. The adjudication of new diagnoses of T2D during follow-up was conducted by the Clinical End point Committee (blinded to the intervention group) as described elsewhere (9,10). The American Diabetes Association criteria (1), namely, two confirmations of fasting plasma glucose ≥7.0 mmol/L or 2-h plasma glucose ≥11.1 mmol/L after a 75-g oral glucose load, were used to adjudicate cases.

Statistical Analysis

Missing values for 26 lipid metabolites (four with >5% of values missing and 22 with <1% missing) were replaced by the half of the minimum detectable value, assuming that they were missing because they were at lower concentrations than the detectable threshold.

Baseline individual lipid values were normalized and scaled in multiples of 1 SD with Blom inverse normal transformation (12). Changes in lipid values (1-year value minus the baseline value) were calculated, and the resulting difference was also normalized and scaled.

The statistical assessment of the association between lipid patterns and T2D was conducted in three steps.

Factor Analysis: Lipid Factors

The first step was an exploratory principal component analysis (PCA) performed considering the 207 lipid metabolites as candidates to be included in the obtained factors. Those factors with an eigenvalue higher than 2 were retained. Fifteen factors (not correlated) were extracted explaining 84% of the total variance. An orthogonal rotation (varimax) was used to better interpret the results. Individual metabolites with absolute loadings >0.40 were considered relevant components of the identified factors (Supplementary Table 1), as previously done based on convention (13). To analyze the association of each extracted factor with T2D, Cox regression models with Barlow weights (14) were fitted. Each factor was introduced in the model either as a continuous variable or categorized in quartiles and was adjusted for age, sex, intervention group, and the rest of the PCA-identified factors. Quartile cutoff points were generated considering only the subcohort, and thereafter cases were categorized according to the same cutoff points.

Similar models were used to evaluate the linear trend among factors considering the median value of each quartile as a quantitative variable.

Grouping by Lipid Families: Lipid Scores

After identifying PCA factors associated with T2D risk, our second step was to evaluate the association between the main lipids represented in those PCA factors and T2D risk. In this second step, lipid molecular species were summed into individual scores based on their lipid class (according to their chemical structure) to clarify potential biological mechanisms. Unlike the lipid factors (obtained only through the data-driven PCA), in the lipid scores, both the known chemical structure and the data-driven result obtained with PCA were accounted for.

Three parallel Cox regression models for lipid scores were designed. Model 1 (M1) was adjusted for age, sex, and intervention group; model 2 (M2) was additionally adjusted for BMI, smoking, leisure-time physical activity, hypertension, and dyslipidemia; and model 3 (M3) was additionally adjusted for baseline glucose (continuous and quadratic term).

To analyze the effects of each lipid score, we calculated hazard ratios (HRs) for incident T2D per 1-SD increase in baseline individual lipid concentrations with weighted Cox regression models using the M3 adjustment (see above). The HRs for individual lipids and their P values were plotted, according to the previously defined lipid scores, in a two-dimensional graph defined by the number of carbon atoms (x-axis) and the number of double bonds (y-axis) in the acyl chain, as we previously reported for CVD (15). Lipids with the same number of carbon atoms and double bonds were slightly pulled apart horizontally to visualize both results. We included an additional graph to plot the residual of each triacylglycerol (TAG) over the total content of the considered TAG, due to the fact that hypertriglyceridemia is a known risk factor for T2D (16).

Areas under the receiver operating curves (AUROCs) were estimated to assess the predictive ability of each score beyond known predictors of T2D: age, sex, BMI, smoking, hypertension, dyslipidemia, leisure-time physical activity, intervention group, and baseline glucose concentrations.

One-Year Changes in Lipid Scores

Our third step was to study the effects of changes in these lipid scores after 1 year of the intervention. Changes for each lipid score were used as the main exposure variable. After excluding T2D cases that occurred during the 1st year of intervention, each change of score was introduced (as a continuous variable or in quartiles) as the main independent variable in Cox models adjusted for their respective baseline scores. The models were the same as those used above to analyze the effects of the scores on T2D risk, plus an additional adjustment for baseline lipid concentration. For 1-year changes in lipids, we also plotted HRs and P values according to number of carbon atoms and double bonds.

An additional analysis was conducted to observe if adjustment for 1-year changes in lipids attenuated the association between the nutritional intervention and T2D using weighted Cox models with robust SE to account for intracluster correlations (9). The models were initially adjusted for age, BMI, smoking, hypertension, dyslipidemia, and baseline glucose (linear and quadratic term) and propensity scores that used 30 baseline variables to estimate the probability of assignment to each of the intervention groups and stratified by center, sex, and educational level (9). In a second step, we additionally adjusted for 1-year changes in lipids to observe if the HRs were attenuated, which would suggest that the lipid changes had a mediating effect. Statistical significance was set a priori at <0.05.

Results

Baseline characteristics by diabetes incident status are shown in Table 1. Subjects who developed T2D during follow-up were more likely to be men and smokers. They showed a mean fasting glucose concentration of 117 ± 18 mg/dL at baseline, suggesting that many may have had prediabetes at baseline.

View this table:
  • View inline
  • View popup
Table 1

Baseline characteristics of study participants according to outcome status

Factor Analysis

Fifteen factors with eigenvalues ≥2 were extracted from the PCA analysis conducted on 207 candidate baseline lipid metabolites (Supplementary Table 1). Four of them were directly associated and three inversely associated with T2D incidence (Table 2 and Supplementary Table 2). Lysophospholipids (LPs), cholesterol esters (CEs), sphingomyelins (SMs), and phosphatidylcholine-plasmalogens (PC-PLs) were widely represented among factors associated with lower T2D risk. TAGs, diacylglycerols (DAGs), and phosphatidylethanolamines (PEs) were preponderantly associated with higher risk.

View this table:
  • View inline
  • View popup
Table 2

Association between baseline lipid factors (PCA extracted) and T2D risk (adjusted for age, sex, and intervention group)

Baseline Lipid Scores

Based on these lipid patterns, seven classes or families of lipids (according to their common chemical structures) were identified. The identified metabolites belonging to each lipid class were summed to build the following scores: 1) LP score, grouping lysophosphatidylcholines (LPCs) and lysophosphatidylethanolamines (LPEs) (n of metabolites = 18); 2) PC-PL score (n of metabolites = 15); 3) SM score (n of metabolites = 11); 4) CE score (n of metabolites = 13); 5) TAG score, including only those TAGs with ≤56 carbon atoms (C) and ≤3 double bonds (n of metabolites = 40); 6) DAG score (n of metabolites = 14); and 7) PE score (n of metabolites = 12). Baseline levels of the scores by intervention arm of the trial can be found in Supplementary Table 3.

Higher LP, PC-PL, SM, and CE baseline scores presented a significant inverse association with T2D (P for linear trend ≤0.001 for all, adjusted for sex, age, and intervention group). These associations were maintained after additional adjustment for BMI, smoking, leisure-time physical activity, hypertension, and dyslipidemia (Table 3). When fasting baseline glucose levels were introduced into the models, the inverse association was maintained for the LP, PC-PL, SM, and CE scores (P for linear trend = 0.040, 0.001, <0.001, and 0.002, respectively). On the contrary, higher TAG, DAG, and PE baselines scores presented a significant direct association with T2D (P for trend <0.001, <0.001, and 0.001, respectively, in fully adjusted models). Further adjustments for baseline glucose showed that the association between the DAG score and T2D was robust (P for linear trend = 0.011). The association between TAGs and incident T2D was attenuated but remained significant (P = 0.044).

View this table:
  • View inline
  • View popup
Table 3

Association between baseline lipid scores and T2D risk

When we assessed associations for each individual lipid by number of carbon atoms and number of double bonds (Fig. 1A), LPs, SMs, and CEs were the most homogenous lipid groups regarding their individual inverse associations with T2D. A clear direct association of TAG, DAG, and PE scores with T2D was evident (Fig. 1B). We did not find any clear pattern of associations with T2D incidence by number of carbon atoms or double bonds. However, for TAGs, we observed that the direct association with T2D was strongly attenuated for odd-chain TAGs. In fact, we observed that odd-chain TAGs adjusted for total TAGs presented an inverse association with T2D as depicted in Fig. 1C, which plots the residual of each individual metabolite beyond the sum of all the considered TAGs (≤56 C and ≤3 double bonds).

  • Download figure
  • Open in new tab
  • Download powerpoint
  • Download figure
  • Open in new tab
  • Download powerpoint
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

A: HRs per 1-SD increase in baseline lipid concentration for lipid groups inversely associated with T2D. Lipid species were inverse normally transformed, and HRs were calculated from weighted Cox models adjusted for age, sex, intervention group, BMI, smoking, hypertension, dyslipidemia, and baseline glucose (linear and quadratic term). B: HRs per 1-SD increase in baseline lipid concentration for lipid groups directly associated with T2D. Lipid species were inverse normally transformed, and HRs were calculated from weighted Cox models adjusted for age, sex, intervention group, BMI, smoking, hypertension, dyslipidemia, and baseline glucose (linear and quadratic term). C: HRs for T2D per 1 SD for the residual of each TAG over the total content of the considered TAG. Lipid species were inverse normally transformed before calculating the residual, and HRs were calculated from weighted Cox models adjusted for age, sex, intervention group, BMI, smoking, hypertension, dyslipidemia, and baseline glucose (linear and quadratic term).

We observed that the sum of all the lipid scores that were inversely associated with T2D incidence (LP, PC-PL, SM, and CE scores) was able to significantly improve the prediction of T2D beyond conventional risk factors, although the size of this improvement was small (AUROC excluding lipid scores = 0.83 [95% CI 0.81–0.86], AUROC including LP, PC-PL, SM, and CE scores = 0.84 [95% CI 0.82–0.87]; P = 0.036 for the improvement).

One-Year Change in Lipid Scores

In the analysis assessing 1-year changes, the number of incident cases (only those occurring after the 1st year and with available plasma sample) was reduced from 251 to 121 and the statistical power was considerably lower. We additionally adjusted for baseline scores to assess the association of 1-year changes beyond baseline predictions. Point estimates for 1-year changes suggested associations similar to those observed at baseline (inverse associations for LP, PC-PL, SM, and CE scores and direct associations for TAG, DAG, and PE scores). However, all of these associations did not remain statistically significant (Supplementary Table 4). We found a significant independent direct association per each SD increase in 1-year changes in PE score (HR 1.25 [95% CI 1.01–1.56]).

The assessment to determine whether 1-year changes in lipid scores mediated the effect of the MedDiet on T2D found that 1-year changes in short TAGs partially mediated both the intervention with MedDiet + EVOO (HR 0.39 [95% CI 0.19–0.80] without TAG change; HR 0.45 [95% CI 0.22–0.93] when adjusting for TAG change) and MedDiet + nuts (HR 0.49 [95% CI 0.25–0.96] without TAG change; HR 0.53 [95% CI 0.26–1.06] after adjusting for TAG change) (Supplementary Table 5). In fact, the MedDiet + nuts intervention was marginally associated with reduced TAG plasma levels after 1 year (B coefficient = −4.81, P = 0.062; data not shown) compared with the control group in a linear regression model adjusted for the same confounders as the Cox models. Moreover, we observed that 1-year changes in DAGs and in PEs in part explained the effects of the MedDiet + EVOO intervention. The HRs for T2D were 0.38 (95% CI 0.19–0.80) before adjustment for changes in DAGs vs. 0.43 (95% CI 0.21–0.89) after adjustment for the changes and 0.36 (95% CI 0.18–0.73) before adjustment for changes in PEs vs. 0.40 (95% CI 0.18–0.87) after additional adjustment for PE change (Supplementary Table 5). Changes in LP, PC-PL, SM, or CE scores showed no apparent mediating effects.

Conclusions

We have identified several individual molecular species and some lipid classes prospectively associated with T2D risk. Baseline LP (LPC and LPE), PC-PL, SM, and CE scores were inversely associated with the risk of T2D, whereas baseline TAG, DAG, and PE scores were directly associated with T2D incidence. For 1-year changes in these scores, associations beyond baseline levels were mainly nonsignificant, although the point estimates remained in the same direction. However, these 1-year change analyses had suboptimal power.

We found that both LPCs and LPEs, grouped as LPs, were associated with reduced risk of T2D. Previous studies found that LPC levels were lower in individuals with obesity, insulin resistance, and T2D (17–19). In fact, increased levels of LPCs have been defined as indicators of metabolic health in obesity, as LPCs appear to have glucose-lowering and anti-inflammatory effects (20). Similarly, LPC and LPE levels were reported to be lower in patients with T2D, and in patients with diabetes, lower levels of these lipids were associated with higher risk of CVD (21).

Similarly to LP, we also found that PC-PLs were inversely related to T2D risk. Plasmalogens have been widely investigated because of their role as endogenous antioxidants, limiting the oxidation of other lipids (7,22). They may also decrease the risk of T2D through other beneficial mechanisms, such as antiapoptotic and anti-inflammatory functions (7).

A few studies (23,24), including the EPIC-Potsdam study (25), have reported an inverse association of SMs with T2D, consistent with our results. A study of a large cohort of patients with prediabetes and diabetes also reported inverse associations between plasma odd-chain SMs and T2D (6). The knockout model of SM synthase results in mitochondrial dysfunction and impaired glucose-stimulated insulin secretion, which provides mechanistic support for our findings (26).

Unexpectedly, we found an inverse association of CEs with T2D. Contrary to our findings, previous cross-sectional studies have reported strong direct associations with T2D (6,27). However, we found similar inverse associations between CEs and CVD in the PREDIMED trial (28). It is possible that we could have detected the defined “atherogenic lipoprotein phenotype” in subjects at high T2D risk (high plasma levels of TAGs, low levels of HDL, and atypically dense LDL particles). In this situation, LDL particles are loaded with TAGs instead of CEs, and after the hydrolysis of TAGs in the liver, lipid-depleted LDL particles (small and dense) are released (29). By losing their lipid core, these particles also lose antioxidant vitamins and become dense and oxidatively damaged, which may trigger foam-cell formation and therefore atherosclerosis. This lipoprotein phenotype has also been suggested for insulin resistance and eventually T2D (30).

We found that DAG, TAG, and PE scores were strongly associated with a higher risk of T2D. Higher circulating levels of DAG and short TAGs have been previously associated with T2D (6,31). Thus, our findings confirm this positive association, highlighting the adverse role of short and saturated/low unsaturated species (31). Interestingly, after adjusting each individual TAG for the total TAG score, we observed that odd-chain TAGs were inversely associated with T2D. Odd-chain fatty acids, especially C15:0 and C17:0, have been described as biomarkers of dairy product intake (32) and have been reported to be inversely associated with T2D (6,33) and CVD (34). Thus, it seems important to consider the specific fatty acid content of TAGs to establish plasma risk profiles for T2D.

Our results confirm previous studies that have found that PEs are associated with high fasting glucose and T2D (6). PEs are a minor species in plasma, but they are important structural lipids in membranes. An increase in PEs and an imbalance between PC and PE have been related to obesity and nonalcoholic fatty liver disease (35–37), which are both related to T2D (6).

Supplementary Table 6 displays the description of the studies used to compare and discuss the lipidome profiles associated with T2D in our population.

The main strengths of our study include the case-cohort design nested within the PREDIMED trial, which enables the extension of the identified lipid patterns to all PREDIMED participants. Additionally, the analyses considering the complete lipidome allowed us to observe the effect of each lipid metabolite and each lipid group in the context of coexisting and interacting with the other plasma lipids.

We acknowledge that despite our extensive adjustments, residual confounding cannot be ruled out. Also, our results may not be generalizable to other populations because all study participants lived in a Mediterranean country and were at high cardiovascular risk. Additionally, many of the participants who developed diabetes during the trial were people at high risk at baseline. At baseline, case subjects with T2D presented high mean levels of fasting glucose (117 ± 18 mg/dL), which could be because many had prediabetes at baseline. Thus, the baseline lipid biomarkers could reflect an established prediabetic profile rather than a nondiabetic risk profile. In this context, a recent study reported that plasma lipid profiles were similar in subjects with both prediabetes and diabetes (6), which indicates that our identified lipid patterns may be discriminating progressors versus nonprogressors rather than healthy subjects versus subjects with T2D.

Our results have important implications in helping to clarify the biological mechanisms underlying the link between dyslipidemia and T2D. They also suggest that in subjects at high T2D risk, a plasma lipid profile characterized by high levels of DAGs, short TAGs, and PEs and low levels of LPs, PC-PLs, SMs, and CEs could be identified before T2D onset, which could enable early intervention.

Article Information

Acknowledgments. The authors thank Elena Hemler, Harvard T.H. Chan School of Public Health, Boston, MA, for her helpful comments on the manuscript.

Funding. This study was supported by the National Institutes of Health (research grant R01-DK-102896) and a grant from Instituto de Salud Carlos III (RD 06/0045) to M.A.M.-G. C.P. was supported by a postdoctoral fellowship granted by the Autonomous Government of Catalonia (PERIS 2016–2020 Incorporació de Científics I Tecnòlegs, SLT002/0016/00428). M.G.-F. was supported by a fellowship granted by the Daniel Bravo Andreu Private Foundation (Spain) and by the American Diabetes Association (grant 1-18-PMF-029). E.Y. was supported by grant F31-DK-114938-01.

Duality of Interest. E.R. and R.E. received grants from the California Walnut Commission and personal fees from pharmaceutical and beverage companies. J.S.-S. has received grants for research through his institution from the International Nut and Dried Fruit Council and is a nonpaid member of its Scientific Advisory Committee. No other potential conflicts of interest relevant to this article were reported.

Author Contributions. C.R. and M.A.M.-G. conceived and designed the work, acquired or analyzed and interpreted data, conducted the statistical analysis, and drafted and critically revised the manuscript for important intellectual content and approved the version to be published. E.T., C.B.C., M.R.-C., J.S.-S., and F.B.H. conceived and designed the work, acquired or analyzed and interpreted data, and critically revised the manuscript for important intellectual content and approved the version to be published. C.D., D.C., C.P., E.R., R.E., M.G.-F., E.G.-G., M.F., E.Y., J.L., D.W., D.R., L.L., A.A.-G., A.D., M.B., and L.S.-M. acquired or analyzed and interpreted data and critically revised the manuscript for important intellectual content and approved the version to be published. M.A.M.-G. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Prior Presentation. Parts of this study were presented in abstract form at the 52nd Annual Scientific Meeting of the European Society for Clinical Investigation, Barcelona, Spain, 30 May–1 June 2018.

Footnotes

  • Clinical trial reg. no. ISRCTN35739639, www.isrctn.org

  • This article contains Supplementary Data online at http://care.diabetesjournals.org/lookup/suppl/doi:10.2337/dc18-0840/-/DC1.

  • Received April 17, 2018.
  • Accepted September 7, 2018.
  • © 2018 by the American Diabetes Association.
http://www.diabetesjournals.org/content/license

Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. More information is available at http://www.diabetesjournals.org/content/license.

References

  1. ↵
    1. American Diabetes Association
    . 2. Classification and diagnosis of diabetes. Diabetes Care 2017;40(Suppl. 1):S11–S24pmid:27979889
    OpenUrlFREE Full Text
  2. ↵
    1. American Diabetes Association
    . Diagnosis and classification of diabetes mellitus. Diabetes Care 2008;31(Suppl. 1):S55–S60pmid:18165338
    OpenUrlFREE Full Text
  3. ↵
    1. Ogurtsova K,
    2. da Rocha Fernandes JD,
    3. Huang Y, et al
    . IDF Diabetes Atlas: global estimates for the prevalence of diabetes for 2015 and 2040. Diabetes Res Clin Pract 2017;128:40–50pmid:28437734
    OpenUrlCrossRefPubMed
  4. ↵
    1. Mooradian AD
    . Dyslipidemia in type 2 diabetes mellitus. Nat Clin Pract Endocrinol Metab 2009;5:150–159pmid:19229235
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    1. Chehade JM,
    2. Gladysz M,
    3. Mooradian AD
    . Dyslipidemia in type 2 diabetes: prevalence, pathophysiology, and management. Drugs 2013;73:327–339pmid:23479408
    OpenUrlCrossRefPubMedWeb of Science
  6. ↵
    1. Meikle PJ,
    2. Wong G,
    3. Barlow CK, et al
    . Plasma lipid profiling shows similar associations with prediabetes and type 2 diabetes. PLoS One 2013;8:e74341pmid:24086336
    OpenUrlCrossRefPubMed
  7. ↵
    1. Huynh K,
    2. Martins RN,
    3. Meikle PJ
    . Lipidomic profiles in diabetes and dementia. J Alzheimers Dis 2017;59:433–444pmid:28582856
    OpenUrlPubMed
  8. ↵
    1. Martínez-González MA,
    2. Corella D,
    3. Salas-Salvadó J, et al.; PREDIMED Study Investigators
    . Cohort profile: design and methods of the PREDIMED study. Int J Epidemiol 2012;41:377–385pmid:21172932
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    1. Estruch R,
    2. Ros E,
    3. Salas-Salvadó J, et al.; PREDIMED Study Investigators
    . Primary prevention of cardiovascular disease with a Mediterranean Diet supplemented with extra-virgin olive oil or nuts. N Engl J Med 2018;378:e34pmid:29897866
    OpenUrlCrossRefPubMed
  10. ↵
    1. Salas-Salvadó J,
    2. Bulló M,
    3. Estruch R, et al
    . Prevention of diabetes with Mediterranean diets: a subgroup analysis of a randomized trial. Ann Intern Med 2014;160:1–10pmid:24573661
    OpenUrlCrossRefPubMedWeb of Science
  11. ↵
    1. Elosua R,
    2. Garcia M,
    3. Aguilar A,
    4. Molina L,
    5. Covas MI,
    6. Marrugat J; Investigators of the MARATDON Group
    . Validation of the Minnesota leisure time physical activity questionnaire in Spanish women. Med Sci Sports Exerc 2000;32:1431–1437pmid:10949009
    OpenUrlCrossRefPubMedWeb of Science
  12. ↵
    1. Blom G
    . Statistical Estimates and Transformed Beta-Variables. New York, John Wiley & Sons A/S, 1958
  13. ↵
    1. Hunter WG,
    2. Kelly JP,
    3. McGarrah RW III,
    4. Kraus WE,
    5. Shah SH
    . Metabolic dysfunction in heart failure: diagnostic, prognostic, and pathophysiologic insights from metabolomic profiling. Curr Heart Fail Rep 2016;13:119–131pmid:27216948
    OpenUrlPubMed
  14. ↵
    1. Barlow WE,
    2. Ichikawa L,
    3. Rosner D,
    4. Izumi S
    . Analysis of case-cohort designs. J Clin Epidemiol 1999;52:1165–1172pmid:10580779
    OpenUrlCrossRefPubMedWeb of Science
  15. ↵
    1. Toledo E,
    2. Wang DD,
    3. Ruiz-Canela M, et al
    . Plasma lipidomic profiles and cardiovascular events in a randomized intervention trial with the Mediterranean diet. Am J Clin Nutr 2017;106:973–983pmid:28814398
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Marott SC,
    2. Nordestgaard BG,
    3. Tybjærg-Hansen A,
    4. Benn M
    . Components of the metabolic syndrome and risk of type 2 diabetes. J Clin Endocrinol Metab 2016;101:3212–3221pmid:27285293
    OpenUrlPubMed
  17. ↵
    1. Barber MN,
    2. Risis S,
    3. Yang C, et al
    . Plasma lysophosphatidylcholine levels are reduced in obesity and type 2 diabetes. PLoS One 2012;7:e41456
    OpenUrlCrossRefPubMed
    1. Zhao X,
    2. Fritsche J,
    3. Wang J, et al
    . Metabonomic fingerprints of fasting plasma and spot urine reveal human pre-diabetic metabolic traits. Metabolomics 2010;6:362–374
    OpenUrlCrossRefPubMedWeb of Science
  18. ↵
    1. Tonks KT,
    2. Coster AC,
    3. Christopher MJ, et al
    . Skeletal muscle and plasma lipidomic signatures of insulin resistance and overweight/obesity in humans. Obesity (Silver Spring) 2016;24:908–916
    OpenUrlCrossRef
  19. ↵
    1. Lehmann R,
    2. Franken H,
    3. Dammeier S, et al
    . Circulating lysophosphatidylcholines are markers of a metabolically benign nonalcoholic fatty liver. Diabetes Care 2013;36:2331–2338pmid:23514731
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. García-Fontana B,
    2. Morales-Santana S,
    3. Díaz Navarro C, et al
    . Metabolomic profile related to cardiovascular disease in patients with type 2 diabetes mellitus: a pilot study. Talanta 2016;148:135–143pmid:26653434
    OpenUrlCrossRefPubMed
  21. ↵
    1. Braverman NE,
    2. Moser AB
    . Functions of plasmalogen lipids in health and disease. Biochim Biophys Acta 2012;1822:1442–1452
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    1. Morris JK,
    2. Piccolo BD,
    3. Shankar K,
    4. Thyfault JP,
    5. Adams SH
    . The serum metabolomics signature of type 2 diabetes is obscured in Alzheimer’s disease. Am J Physiol Metab 2018;314:E584–E596
    OpenUrl
  23. ↵
    1. Xu F,
    2. Tavintharan S,
    3. Sum CF,
    4. Woon K,
    5. Lim SC,
    6. Ong CN
    . Metabolic signature shift in type 2 diabetes mellitus revealed by mass spectrometry-based metabolomics. J Clin Endocrinol Metab 2013;98:E1060–E1065pmid:23633210
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    1. Floegel A,
    2. Stefan N,
    3. Yu Z, et al
    . Identification of serum metabolites associated with risk of type 2 diabetes using a targeted metabolomic approach. Diabetes 2013;62:639–648pmid:23043162
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Yano M,
    2. Watanabe K,
    3. Yamamoto T, et al
    . Mitochondrial dysfunction and increased reactive oxygen species impair insulin secretion in sphingomyelin synthase 1-null mice. J Biol Chem 2011;286:3992–4002pmid:21115496
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Kopprasch S,
    2. Dheban S,
    3. Schuhmann K, et al
    . Detection of independent associations of plasma lipidomic parameters with insulin sensitivity indices using data mining methodology. PLoS One 2016;11:e0164173
    OpenUrl
  27. ↵
    1. Razquin C,
    2. Liang L,
    3. Toledo E, et al
    . Plasma lipidome patterns associated with cardiovascular risk in the PREDIMED trial: a case-cohort study. Int J Cardiol 2018;253:126–132pmid:29306452
    OpenUrlPubMed
  28. ↵
    1. Gardner CD,
    2. Fortmann SP,
    3. Krauss RM
    . Association of small low-density lipoprotein particles with the incidence of coronary artery disease in men and women. JAMA 1996;276:875–881pmid:8782636
    OpenUrlCrossRefPubMedWeb of Science
  29. ↵
    1. Feher MD,
    2. Caslake M,
    3. Foxton J,
    4. Cox A,
    5. Packard CJ
    . Atherogenic lipoprotein phenotype in type 2 diabetes: reversal with micronised fenofibrate. Diabetes Metab Res Rev 1999;15:395–399pmid:10634964
    OpenUrlCrossRefPubMedWeb of Science
  30. ↵
    1. Rhee EP,
    2. Cheng S,
    3. Larson MG, et al
    . Lipid profiling identifies a triacylglycerol signature of insulin resistance and improves diabetes prediction in humans. J Clin Invest 2011;121:1402–1411pmid:21403394
    OpenUrlCrossRefPubMedWeb of Science
  31. ↵
    1. Brevik A,
    2. Veierød MB,
    3. Drevon CA,
    4. Andersen LF
    . Evaluation of the odd fatty acids 15:0 and 17:0 in serum and adipose tissue as markers of intake of milk and dairy fat. Eur J Clin Nutr 2005;59:1417–1422pmid:16118654
    OpenUrlCrossRefPubMed
  32. ↵
    1. Imamura F,
    2. Sharp SJ,
    3. Koulman A, et al
    . A combination of plasma phospholipid fatty acids and its association with incidence of type 2 diabetes: the EPIC-InterAct case-cohort study. PLoS Med 2017;14:e1002409
  33. ↵
    1. Khaw K-T,
    2. Friesen MD,
    3. Riboli E,
    4. Luben R,
    5. Wareham N
    . Plasma phospholipid fatty acid concentration and incident coronary heart disease in men and women: the EPIC-Norfolk Prospective Study. PLoS Med 2012;9:e1001255
  34. ↵
    1. Arendt BM,
    2. Ma DWL,
    3. Simons B, et al
    . Nonalcoholic fatty liver disease is associated with lower hepatic and erythrocyte ratios of phosphatidylcholine to phosphatidylethanolamine. Appl Physiol Nutr Metab 2013;38:334–340pmid:23537027
    OpenUrlCrossRefPubMed
    1. Li Z,
    2. Agellon LB,
    3. Allen TM, et al
    . The ratio of phosphatidylcholine to phosphatidylethanolamine influences membrane integrity and steatohepatitis. Cell Metab 2006;3:321–331pmid:16679290
    OpenUrlCrossRefPubMedWeb of Science
  35. ↵
    1. Weir JM,
    2. Wong G,
    3. Barlow CK, et al
    . Plasma lipid profiling in a large population-based cohort. J Lipid Res 2013;54:2898–2908pmid:23868910
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Diabetes Care: 41 (12)

In this Issue

December 2018, 41(12)
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by Author
  • Masthead (PDF)
Sign up to receive current issue alerts
View Selected Citations (0)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word about Diabetes Care.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Plasma Lipidomic Profiling and Risk of Type 2 Diabetes in the PREDIMED Trial
(Your Name) has forwarded a page to you from Diabetes Care
(Your Name) thought you would like to see this page from the Diabetes Care web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Plasma Lipidomic Profiling and Risk of Type 2 Diabetes in the PREDIMED Trial
Cristina Razquin, Estefanía Toledo, Clary B. Clish, Miguel Ruiz-Canela, Courtney Dennis, Dolores Corella, Christopher Papandreou, Emilio Ros, Ramon Estruch, Marta Guasch-Ferré, Enrique Gómez-Gracia, Montserrat Fitó, Edward Yu, José Lapetra, Dong Wang, Dora Romaguera, Liming Liang, Angel Alonso-Gómez, Amy Deik, Mónica Bullo, Lluis Serra-Majem, Jordi Salas-Salvadó, Frank B. Hu, Miguel A. Martínez-González
Diabetes Care Dec 2018, 41 (12) 2617-2624; DOI: 10.2337/dc18-0840

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Add to Selected Citations
Share

Plasma Lipidomic Profiling and Risk of Type 2 Diabetes in the PREDIMED Trial
Cristina Razquin, Estefanía Toledo, Clary B. Clish, Miguel Ruiz-Canela, Courtney Dennis, Dolores Corella, Christopher Papandreou, Emilio Ros, Ramon Estruch, Marta Guasch-Ferré, Enrique Gómez-Gracia, Montserrat Fitó, Edward Yu, José Lapetra, Dong Wang, Dora Romaguera, Liming Liang, Angel Alonso-Gómez, Amy Deik, Mónica Bullo, Lluis Serra-Majem, Jordi Salas-Salvadó, Frank B. Hu, Miguel A. Martínez-González
Diabetes Care Dec 2018, 41 (12) 2617-2624; DOI: 10.2337/dc18-0840
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Research Design and Methods
    • Results
    • Conclusions
    • Article Information
    • Footnotes
    • References
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Impact of Obesity on Measures of Cardiovascular and Kidney Health in Youth With Type 1 Diabetes as Compared With Youth With Type 2 Diabetes
  • Serum Galectin-3 and Subsequent Risk of Coronary Heart Disease in Subjects With Childhood-Onset Type 1 Diabetes: A Cohort Study
  • OGTT Glucose Response Curves, Insulin Sensitivity, and β-Cell Function in RISE: Comparison Between Youth and Adults at Randomization and in Response to Interventions to Preserve β-Cell Function
Show more Cardiovascular and Metabolic Risk

Similar Articles

Navigate

  • Current Issue
  • Standards of Care Guidelines
  • Online Ahead of Print
  • Archives
  • Submit
  • Subscribe
  • Email Alerts
  • RSS Feeds

More Information

  • About the Journal
  • Instructions for Authors
  • Journal Policies
  • Reprints and Permissions
  • Advertising
  • Privacy Policy: ADA Journals
  • Copyright Notice/Public Access Policy
  • Contact Us

Other ADA Resources

  • Diabetes
  • Clinical Diabetes
  • Diabetes Spectrum
  • Scientific Sessions Abstracts
  • Standards of Medical Care in Diabetes
  • BMJ Open - Diabetes Research & Care
  • Professional Books
  • Diabetes Forecast

 

  • DiabetesJournals.org
  • Diabetes Core Update
  • ADA's DiabetesPro
  • ADA Member Directory
  • Diabetes.org

© 2021 by the American Diabetes Association. Diabetes Care Print ISSN: 0149-5992, Online ISSN: 1935-5548.