Association of Plasma Phospholipid n-3 and n-6 Polyunsaturated Fatty Acids with Type 2 Diabetes: The EPIC-InterAct Case-Cohort Study

Abstract

Background

Whether and how n-3 and n-6 polyunsaturated fatty acids (PUFAs) are related to type 2 diabetes (T2D) is debated. Objectively measured plasma PUFAs can help to clarify these associations.

Methods and Findings

Plasma phospholipid PUFAs were measured by gas chromatography among 12,132 incident T2D cases and 15,919 subcohort participants in the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct study across eight European countries. Country-specific hazard ratios (HRs) were estimated using Prentice-weighted Cox regression and pooled by random-effects meta-analysis. We also systematically reviewed published prospective studies on circulating PUFAs and T2D risk and pooled the quantitative evidence for comparison with results from EPIC-InterAct. In EPIC-InterAct, among long-chain n-3 PUFAs, [alpha] -linolenic acid (ALA) was inversely associated with T2D (HR per standard deviation [SD] 0.93; 95% CI 0.88-0.98), but eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) were not significantly associated. Among n-6 PUFAs, linoleic acid (LA) (0.80; 95% CI 0.77-0.83) and eicosadienoic acid (EDA) (0.89; 95% CI 0.85-0.94) were inversely related, and arachidonic acid (AA) was not significantly associated, while significant positive associations were observed with [gamma]-linolenic acid (GLA), dihomo-GLA, docosatetraenoic acid (DTA), and docosapentaenoic acid (n6-DPA), with HRs between 1.13 to 1.46 per SD. These findings from EPIC-InterAct were broadly similar to comparative findings from summary estimates from up to nine studies including between 71 to 2,499 T2D cases. Limitations included potential residual confounding and the inability to distinguish between dietary and metabolic influences on plasma phospholipid PUFAs.

Conclusions

These large-scale findings suggest an important inverse association of circulating plant-origin n-3 PUFA (ALA) but no convincing association of marine-derived n3 PUFAs (EPA and DHA) with T2D. Moreover, they highlight that the most abundant n6-PUFA (LA) is inversely associated with T2D. The detection of associations with previously less well-investigated PUFAs points to the importance of considering individual fatty acids rather than focusing on fatty acid class.

Background

Whether and how n-3 and n-6 polyunsaturated fatty acids (PUFAs) are related to type 2 diabetes (T2D) is debated. Objectively measured plasma PUFAs can help to clarify these associations.

Methods and Findings

Plasma phospholipid PUFAs were measured by gas chromatography among 12,132 incident T2D cases and 15,919 subcohort participants in the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct study across eight European countries. Country-specific hazard ratios (HRs) were estimated using Prentice-weighted Cox regression and pooled by random-effects meta-analysis. We also systematically reviewed published prospective studies on circulating PUFAs and T2D risk and pooled the quantitative evidence for comparison with results from EPIC-InterAct. In EPIC-InterAct, among long-chain n-3 PUFAs, [alpha] -linolenic acid (ALA) was inversely associated with T2D (HR per standard deviation [SD] 0.93; 95% CI 0.88-0.98), but eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) were not significantly associated. Among n-6 PUFAs, linoleic acid (LA) (0.80; 95% CI 0.77-0.83) and eicosadienoic acid (EDA) (0.89; 95% CI 0.85-0.94) were inversely related, and arachidonic acid (AA) was not significantly associated, while significant positive associations were observed with [gamma]-linolenic acid (GLA), dihomo-GLA, docosatetraenoic acid (DTA), and docosapentaenoic acid (n6-DPA), with HRs between 1.13 to 1.46 per SD. These findings from EPIC-InterAct were broadly similar to comparative findings from summary estimates from up to nine studies including between 71 to 2,499 T2D cases. Limitations included potential residual confounding and the inability to distinguish between dietary and metabolic influences on plasma phospholipid PUFAs.

Conclusions

These large-scale findings suggest an important inverse association of circulating plant-origin n-3 PUFA (ALA) but no convincing association of marine-derived n3 PUFAs (EPA and DHA) with T2D. Moreover, they highlight that the most abundant n6-PUFA (LA) is inversely associated with T2D. The detection of associations with previously less well-investigated PUFAs points to the importance of considering individual fatty acids rather than focusing on fatty acid class.

Author Summary

Why Was This Study Done?

* Most dietary guidelines recommend the consumption of polyunsaturated fatty acids for cardiovascular health, but it is unclear whether or how n-3 and n-6 types of polyunsaturated fatty acids are related to type 2 diabetes.

* Health concerns have been raised previously about a diet high in linoleic acid (n-6 fatty acid), but its association with type 2 diabetes is unclear.

* Major limitations in previous studies have included the error-prone subjective assessment of the habitual consumption of polyunsaturated fatty acids when dietary intakes were evaluated and a small number of type 2 diabetes cases (n = 71 to 673) when objective biomarkers of polyunsaturated fatty acids were measured.

What Did the Researchers Do and Find?

* We measured circulating individual polyunsaturated fatty acids in the blood samples of individuals within a large study from across eight countries of Europe among a reference sample of 15,919 individuals as well as 12,132 individuals who subsequently developed diabetes. Individuals were followed up for an average of approximately 10 y.

* We investigated the association between individual polyunsaturated fatty acids and the risk of future type 2 diabetes using statistical analyses that accounted for factors that could be potential alternative explanations for any observed associations.

* We found that higher levels of blood alpha-linolenic acid, a plant-origin n-3 fatty acid, and n-6 linoleic acid, the most abundant type of polyunsaturated fatty acid, were associated with a lower risk of future type 2 diabetes. In contrast, higher levels of four other minor individual n-6 fatty acids were associated with higher type 2 diabetes risk, while the blood marine-origin n-3 fatty acids were not associated with future diabetes.

What Do These Findings Mean?

* Our findings show that it is important to consider individual circulating polyunsaturated fatty acids for association with type 2 diabetes risk, rather than placing emphasis on the class of circulating polyunsaturated fatty acids.

* We found that blood n-6 linoleic acid, the most abundant polyunsaturated fatty acid, is inversely associated with type 2 diabetes.

* We found no evidence that blood total n-6 polyunsaturated fatty acids may elevate the risk of type 2 diabetes, but several individual minor blood n-6 polyunsaturated fatty acids were associated with increased type 2 diabetes risk, highlighting the importance of polyunsaturated fatty acid metabolism in the development of type 2 diabetes.

Citation: Forouhi NG, Imamura F, Sharp SJ, Koulman A, Schulze MB, Zheng J, et al. (2016) Association of Plasma Phospholipid n-3 and n-6 Polyunsaturated Fatty Acids with Type 2 Diabetes: The EPIC-InterAct Case-Cohort Study. PLoS Med 13(7): e1002094. doi:10.1371/journal.pmed.1002094

Academic Editor: Ronald C. W. Ma, Chinese University of Hong Kong, CHINA

Received: May 6, 2016; Accepted: June 16, 2016; Published: July 19, 2016

Copyright: © 2016 Forouhi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: EPIC-InterAct Study data cannot be deposited publicly as these collaborative data originate from multiple research institutions across eight European countries with different legal frameworks. The authors confirm that researchers seeking the analysis dataset for this work can submit a data request to the EPIC-InterAct study central contact point by emailing [email protected].

Funding: Funding for the InterAct project was provided by the EU FP6 programme (grant number LSHM_CT_2006_037197). In addition, InterAct investigators acknowledge funding from the following sources: Medical Research Council Epidemiology Unit MC_UU_12015/1 and MC_UU_12015/5, and Medical Research Council Human Nutrition Research MC_UP_A090_1006 and Cambridge Lipidomics Biomarker Research Initiative G0800783; TJK: Cancer Research UK; JMH and MJT: Health Research Fund of the Spanish Ministry of Health; Murcia Regional Government (N° 6236); MG: Regional Government of Navarre; -IS, DLvdA, AMWS, YTvdS: Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands; Verification of diabetes cases in EPIC-NL was additionally funded by NL Agency grant IGE05012 and an Incentive Grant from the Board of the UMC Utrecht; PWF: Swedish Research Council, Novo Nordisk, Swedish Diabetes Association, Swedish Heart-Lung Foundation; RK: German Cancer Aid, German Ministry of Research (BMBF); KTK: Medical Research Council UK, Cancer Research UK; PMN: Swedish Research Council; KO and AT: Danish Cancer Society; JRQ: Asturias Regional Government; OR: The Västerboten County Council; RT: AIRE-ONLUS Ragusa, AVIS-Ragusa, Sicilian Regional Government; ER: Imperial College Biomedical Research Centre. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: PWF has received consulting honoraria from Ely Lily & Co and Sanofi Aventis. CL receives a stipend as a specialty consulting editor for PLOS Medicine and serves on the journal's editorial board. NJW served as a guest editor on PLOS Medicine's Diabetes Prevention Special Issue.

Abbreviations: AA, arachidonic acid; ALA, [alpha] -linolenic acid; BMI, body mass index; DGLA, dihomo-[gamma]-linolenic acid; DHA, docosahexaenoic acid; DTA, docosatetraenoic acid; EDA, eicosadienoic acid; EPA, eicosapentaenoic acid; EPIC, European Prospective Investigation into Cancer and Nutrition; FAME, fatty acid methyl esters; FFQ, food frequency questionnaire; GLA, [gamma]-linolenic acid; IGF-1, insulin-like growth factor-1; LA, linoleic acid; HR, hazard ratio; MRC, Medical Research Council; n6-DPA, docosapentaenoic acid; PRISMA, Preferred Reporting Items for Systematic Reviews and Meta-analyses; PUFA, polyunsaturated fatty acid; SD, standard deviation; STROBE, Strengthening the Reporting of Observational Studies in Epidemiology; T2D, type 2 diabetes

Introduction

Recognition of the importance of dietary fat quality rather than a focus on total fat intake has led to the promotion of dietary guidelines that encourage the consumption of polyunsaturated fatty acids (PUFAs) for cardiometabolic health [1]. Dietary n-3 (omega-3) PUFAs, including alpha-linolenic acid (ALA, 18:3n-3) from plant sources and eicosapentaenoic acid (EPA, 20:5n-3), docosahexaenoic acid (DHA, 22:6n-3), and docosapentaenoic acid (DPA, 22:5n-3) from fish or seafood sources, are postulated to be beneficial for the prevention of type 2 diabetes (T2D), but there is inconclusive evidence from both intervention and observational studies. For instance, n-3 PUFA supplementation with EPA, DHA, or fish oil did not significantly improve glycaemic or insulin markers in randomised trials among general populations or among those with T2D [2,3]. Four separate meta-analyses of prospective studies of the association between dietary long-chain n-3 PUFAs or fish intake, a major source of n-3 PUFAs, and incidence of T2D reported heterogeneous findings with no association, higher risk, or lower risk depending on geographical location [4-7]. Dietary n-6 (omega-6) PUFAs, although they are about 5-fold to 7-fold more abundant in the Western diet than n-3 PUFAs [8,9], have not been systematically examined, and the findings are inconsistent for association with T2D [10-13]. There has also been an ongoing debate about whether dietary linoleic acid (LA; 18:2n6), the most abundant n-6 PUFA, has adverse health effects and therefore whether its consumption should be limited in the population [14-16].

To our knowledge, there is no trial assessing effects of long-term exposure to n-3 or n-6 PUFAs on incidence of T2D, while past evidence from observational studies is limited by accuracy in measurement of habitual PUFA consumption from self-report of dietary intake of foods or supplements. Few studies have examined the associations of objectively measured circulating blood PUFAs with incidence of T2D, with limited sample sizes and the inclusion of a variable number of individual PUFAs [9,17-20].

We aimed to investigate the associations between objectively measured individual n-3 and n-6 PUFAs in the plasma phospholipid fraction and incident T2D in the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct case-cohort study [21]. As a secondary objective, to consider the totality of existing evidence, we also conducted a systematic literature review and performed a meta-analysis to compare with the findings from EPIC-InterAct.

Methods

Study Design and Population

The methods of the EPIC-InterAct project have previously been described in detail [21]. From among 340,234 persons with 3.99 million person-years of follow-up (1991-2007) in the eight countries of the EPIC study, we ascertained 12,403 T2D cases and selected a random subcohort of 16,835 individuals with baseline plasma samples. After exclusions (n = 548 prevalent diabetes; n = 133 uncertain diabetes status), the subcohort retained 16,154 individuals and included 778 individuals with incident T2D during follow-up (a feature of the case-cohort design) [21]. From the total 27,779 participants, we excluded 483 without blood fatty acid data, therefore including 12,132 T2D cases and 15,919 subcohort participants in the analysis (with 755 incident cases included by design within the subcohort), with a mean follow-up of 9.8 y. All participants gave written informed consent, and the study was approved by the local ethics committees in the participating countries and the Internal Review Board of the International Agency for Research on Cancer.

T2D Case Adjudication

Incident T2D was ascertained until 31 December 2007 by reviewing multiple sources of evidence, reported previously [21]: self-report, linkage to primary-care registers and secondary-care registers, medication use, hospital admissions, and mortality data. No diabetes cases were ascertained solely by self-report, and we sought further evidence for cases with information on incident T2D from fewer than two independent sources.

Measurement of Plasma Phospholipid Fatty Acids

Fatty acids were profiled at the Medical Research Council Human Nutrition Research, Cambridge, United Kingdom, by analysing plasma samples stored at baseline at -196°C (-150°C in Denmark), a temperature at which fatty acids remain stable [22]. The assay methods were described previously [23]. In short, the plasma phospholipid fraction was obtained using solid phase extraction and hydrolysed and methylated, yielding fatty acid methyl esters (FAME), which were separated by gas chromatography (J&W HP-88, 30 m length) equipped with flame ionisation detection (7890N GC Agilent Technologies, United States). Samples from cases and subcohort participants were processed in a random order by centre, and laboratory staff were blinded to any participant characteristics. Fatty acids were identified by their retention times compared with those of commercial standards and expressed as percent of total phospholipid fatty acids (mol%). Among 27 fatty acids with relative concentrations above 0.05%, we identified 11 individual PUFAs (4 long-chain n-3 PUFAs and 7 n-6 PUFAs) (Fig 1, showing the biosynthesis pathways and major food sources). Human and equine plasma (Sera Laboratories International, UK) were used as quality control samples and included in each batch [23]. The coefficients of variation were less than 8% except for ALA (13%) and [gamma]-linolenic acid (GLA, 21%) [23]. Using the fatty acid measurements, we calculated the sums of total n-3 PUFAs and total n-6 PUFAs. We also estimated the ratios of PUFA variables in the context of biological processes: 18:3n6 / 18:2n6 and 20:4n6 / 20:3n6, indicating estimated activity of [delta]6 and [delta]5 desaturase enzymes, respectively, catalysing the conversion of LA to GLA and of dihomo-[gamma]-linolenic acid (DGLA) to arachidonic acid (AA); 20:3n6 / 18:2n6, indicating the conversion of LA to DGLA; and ratio of total n-6 PUFAs / total n-3 PUFAs, for its clinical interest in balance between the two major PUFA subclasses.

Fig 1. Schematic showing the names of individual n-3 and n-6 PUFAs measured in EPIC-InterAct, their biosynthesis pathways, and major dietary sources.

http://dx.doi.org/10.1371/journal.pmed.1002094.g001

Measurement of Covariables at Baseline

Weight and height were measured by trained professionals to standardised protocols, and body mass index (BMI) was calculated as weight divided by squared height (kg/m2). Lifestyle questionnaires were used to assess demographics, smoking status, medical history, and education level. Physical activity was assessed by a validated questionnaire from which an ordered four-category variable was derived. Glycated haemoglobin (HbA1c) was measured (SHL, Etten-Leur, Netherlands) in the erythrocyte fraction from samples stored at -196°C (-80°C in Umeå), using the Tosoh-G8 (Tosoh Bioscience, Japan).

Statistical Analysis

All analyses were performed using Stata, version 13 (Stata, College Station, Texas). We estimated the association with incident T2D of each of the individual n-3 and n-6 PUFAs, total n-3 PUFAs, total n-6 PUFAs, and predefined ratios of PUFA variables. Country-specific hazard ratios (HRs) per standard deviation (SD; calculated in the subcohort) and 95% confidence intervals (95% CI) were estimated using Prentice-weighted Cox regression [21], which allows for over-representation of cases in a case-cohort design, and pooled using random-effects meta-analysis. Heterogeneity between countries was expressed as I2 values [24]. We adjusted for potential confounders as follows: Model 1 included age (as underlying timescale), study centre, sex, physical activity index (inactive, moderately inactive, moderately active, or active), smoking status (never, former, or current), education level (none, primary school, technical or professional school, secondary school, or higher education) and BMI (continuous). Model 2 further adjusted for total energy intake (continuous), alcohol (none, >0-<6, 6-<12, 12-<24 and ≥24 g/d), and the other potential dietary confounders unlikely to be major sources of PUFAs: meat, fruits, vegetables, dairy products, and soft drinks (continuous). Model 3 additionally adjusted for possible dietary sources of PUFAs (fish and shellfish, nuts and seeds, vegetable oil, olive oil, and margarine; continuous; Fig 1). In a secondary analysis, we also examined the association of quintiles of each of the PUFA exposures with T2D.

In sensitivity analyses, using Model 3, we examined additional confounding or mediating effects of plasma phospholipid saturated fatty acids (even-chain 14:0, 16:0, and 18:0; and odd-chain 15:0 and 17:0; continuous), prevalent comorbidity (including myocardial infarction, stroke, and cancer; yes/no for each), and baseline HbA1c. To minimise the possibility of reverse causality, we repeated analyses after excluding 2,348 individuals with baseline HbA1c ≥6.5% (≥48 mmol/mol) and after excluding T2D cases (n = 1,048) occurring within 2 y after baseline.

Systematic Review and Meta-analysis of Published Studies

Our objective was to do a comparative meta-analysis of the existing evidence in order to provide a systematic comparison with EPIC-InterAct findings rather than providing a narrative comparison. To compare the findings from EPIC-InterAct with the published literature, we identified prospective studies published by 3 November 2015 that reported on the association of circulating n-3 or n-6 PUFAs with T2D using PubMed (S1 Text for further details). We performed a random effects meta-analysis of the association of circulating n-3 or n-6 PUFAs with T2D from published studies identified as eligible from the literature review, using estimates from a model with the greatest degree of adjustment for potential confounders. The overall effect was estimated per study-specific SD [25] and meta-analysed by lipid fraction and across lipid fractions for each fatty acid. Heterogeneity was assessed using Q and I2 statistics [24]. Publication bias was investigated using Egger's and Begg's tests.

Results

Table 1 summarises the baseline characteristics by case status, and S1 and S2 Tables show characteristics by age, sex, BMI, and country. Among subcohort participants, PUFAs comprised a mean 42.7±2.1% of the total phospholipid fatty acids, of which n-6 PUFAs were more abundant (35.9±2.9%) than n-3 PUFAs (6.7±2.0%). PUFAs with relatively high concentrations, e.g., >1%, were LA (22%), DGLA (3.1%), AA (9.2%), EPA (1.2%), and DHA (4.3%) (Figs 2 and 3).

Table 1. The distribution of sociodemographic and dietary factors at baseline by future case and noncase status: EPIC-InterAct study.

http://dx.doi.org/10.1371/journal.pmed.1002094.t001

Fig 2. HRs of T2D and 95% CIs per 1 SD increase in total and individual n-3 PUFAs (ALA, EPA, DPA, and DHA): EPIC-InterAct study.

Model 1: displayed as diamonds. Age as underlying time variable, and adjusted for centre, sex, physical activity (inactive, moderately inactive, moderately active, or active), smoking (never, former, or current), education level (none, primary school, technical or professional school, secondary school, or higher education), and BMI (continuous, kg/m2). Model 2: displayed as circles. Adjusted as in Model 1 + total energy intake (continuous, kcal/d), alcohol (none, >0-<6, 6-<12, 12-<24 and ≥24 g/d), and (continuous, g/d intake of) meat, fruit and vegetables, dairy products, and soft drinks. Model 3: displayed as squares. Adjusted as in Model 2 + (continuous, g/d intake of) fish and shellfish, nuts and seeds, vegetable oil, olive oil, and margarine.

http://dx.doi.org/10.1371/journal.pmed.1002094.g002

Fig 3. HRs of T2D and 95% CIs per 1 SD increase in total and individual n-6 PUFAs (LA, GLA, EDA, DGLA, AA, DTA, and n-6 DPA): EPIC-InterAct study.

Model 1: displayed as diamonds. Age as underlying time variable, and adjusted for centre, sex, physical activity (inactive, moderately inactive, moderately active, or active), smoking (never, former, or current), education level (none, primary school, technical or professional school, secondary school, or higher education), and BMI (continuous, kg/m2). Model 2: displayed as circles. Adjusted as in Model 1 + total energy intake (continuous, kcal/d), alcohol (none, >0-<6, 6-<12, 12-<24 and ≥24 g/d), and (continuous, g/d intake of) meat, fruit and vegetables, dairy products, and soft drinks. Model 3: Displayed as squares. Adjusted as in Model 2 + (continuous, g/d intake of) fish and shellfish, nuts and seeds, vegetable oil, olive oil, and margarine.

http://dx.doi.org/10.1371/journal.pmed.1002094.g003

The associations we observed did not change substantially by covariables examined (Model 1, 2, and 3) (Fig 2 and S7 Table). There was no evidence of an association between total n-3 PUFAs and T2D, but an inverse association was evident with ALA (18:3n3), with a per 1 SD HR of 0.93 (95% CI 0.88-0.98; Model 3). There was no association with EPA (20:5n3) and DHA (22:6n3), and a borderline inverse association with DPA (22:5n3). Heterogeneity of associations by country appeared to be weak to moderate (I2 = 0% to 50%) (S1 Fig).

Sensitivity analyses were not markedly different for the majority of associations (S3 Table). As an exception, n-3 PUFAs were found to be sensitive to additional adjustment, showing a positive association of EPA with T2D [HR 1.09; 95% CI 1.04-1.14] when adjusted for odd-chain saturated fatty acids and an inverse association of DHA (HR 0.92; 95% CI 0.87-0.98] when adjusted for even-chain saturated fatty acids. Secondary analysis of quintiles of n-3 PUFAs showed similar overall patterns (S4 Table).

Associations with n-6 PUFAs are shown in Fig 3. LA (18:2n6) was strongly inversely associated with T2D (per 1 SD HR 0.80; [95% CI 0.77-0.83], Model 3]. AA (20:4n6) was not associated, while four of the individual n-6 PUFAs were positively associated with T2D: Model 3, GLA (18:3n6) (HR 1.19; 95% CI 1.14-1.24], DGLA (20:3n6) (HR 1.46; 95% CI 1.34-1.59), docosatetraenoic acid (DTA, 22:4n6) (HR 1.13; 95% CI 1.06-1.21), and n-6 DPA (22:5n6) (HR 1.14; 95% CI 1.05-1.24). In sensitivity analyses, results were similar, except adjustment for saturated fatty acids altered findings for 22:4n6 and 22:5n6 to the null (S3 Table). Associations by country are shown in S1 Fig. The findings using quintile distributions of n-6 PUFAs were consistent (S4 Table). For instance, comparing the top and bottom quintiles, there was an inverse association with T2D for LA (HR 0.54; 95% CI 0.47-0.61, p-trend < 0.001).

The ratio of GLA (18:3n6) to LA (18:2n6), representing estimated [delta]6 desaturase activity, was positively associated with T2D (HR 1.21, 95% CI 1.16-1.26; Model 3), while the ratio of AA (20:4n6) to DGLA (20:3n6), representing [delta]5 desaturase activity, was inversely associated (HR 0.73; 95% CI 0.67-0.80). The ratio of DGLA to LA was positively associated with T2D (HR 1.44; 95% CI 1.33-1.57), but there was no association of total n6 to n3 ratio with T2D (HR 0.98; 95% CI 0.93-1.04).

The flow diagram for the systematic literature search and study profiles is shown in S2 Fig and S5 Table. The combined findings from the previously published evidence were based on between 71 to a maximum of 2,499 T2D cases from between one to nine studies, depending on individual fatty acids that were measured (Table 2). Compared with the findings in EPIC-InterAct (12,131 T2D cases), the key differences were that each of circulating total n-6 PUFAs, LA, eicosadienoic acid (EDA), DTA and n-6 DPA, as well as the estimated [delta]6 desaturase activity, were not significantly associated with T2D. There was no evidence of publication bias for any fatty acid variables (pBegg's and pEgger's > 0.1), except for the 20:3n6/18:2n6 ratio (pBegg = 0.04 and pEgger = 0.003) for which exclusion of a potential outlying estimate retained the significant positive association (S3 Fig). Despite differences in lipid fraction (e.g., plasma or erythrocyte phospholipids or total sera), laboratory setting, and population characteristics, there was little heterogeneity for ALA and GLA (I2 = 0) and moderate heterogeneity for major individual PUFAs (EPA, DPA, DHA, and LA) (I2 = 30% to 53%), but substantial heterogeneity for total n-3, total n-6, AA, and the other minor PUFAs. The observed heterogeneity was not explained to a significant degree by lipid fraction in metaregression analysis, except for total n-3 and AA where I2 as the measure of heterogeneity became 0.0% after adjustment for lipid fractions (S6 Table).

Table 2. Results for associations of circulating n-3 and n-6 PUFAs and fatty acid ratios with incidence of T2D: EPIC-InterAct study and comparative meta-analysis of the published prospective studies*.

http://dx.doi.org/10.1371/journal.pmed.1002094.t002

Discussion

EPIC-InterAct provides the largest and most extensive evaluation to date of the association between 11 individual plasma phospholipid n-3 and n-6 PUFAs and incident T2D. These findings indicate that there is no significant association between total n-3 PUFAs and T2D, but that the plant-derived n-3 PUFA, ALA, is robustly inversely associated with T2D incidence. In contrast, the largely marine-derived n-3 PUFAs were not associated, except for a modest inverse association of DPA with T2D. Total plasma phospholipid n-6 PUFAs were inversely associated, largely driven by the most abundant circulating PUFA, LA. Our findings also highlight that four n-6 PUFAs (GLA, DGLA, DTA, and n-6DPA) were associated with higher T2D incidence, while AA (20:4n6) was not associated with T2D. There was no association of the total n6-to-n3 PUFA ratio, but estimated [delta]6 desaturase activity was positively--and [delta]5 desaturase activity inversely--associated with T2D.

Findings in Context: n-3 Fatty Acids

The overall null associations of circulating EPA and DHA with T2D in EPIC-InterAct are consistent with our concurrent comparative meta-analysis of the published literature. Moreover, our findings for circulating EPA and DHA are consistent with the lack of evidence of association between dietary EPA+DHA intakes and T2D risk in prospective cohorts [6] or of net benefit or harm of marine-origin n-3 supplementation on glycaemia or insulin sensitivity in randomised trials [2,3,26]. However, our sensitivity analyses showed that the null findings for EPA and DHA should be interpreted with caution because there was a positive association of T2D with EPA after further adjustment for the odd-chain saturated fatty acids (15:0 and 17:0) and, conversely, an inverse association with DHA after adjustment for even-chain saturated fatty acids (sum of 14:0, 16:0, and 18:0). This indicates that the observed findings may be influenced by complex inter-relationships and metabolic effects of different fatty acids [27], but the totality of current evidence does not support benefits of EPA or DHA for T2D prevention.

On the other hand, we found a modest inverse association between DPA and T2D in EPIC-InterAct as well as in our comparative meta-analysis, though it has been less well investigated previously [6,28]. n-3 DPA is an intermediate between EPA and DHA and appears to have unique biological functions--for example, having a greater potency than EPA and DHA as a precursor of eicosanoids [28]--but little is known of DPA's glycaemic effects [29,30]. In genome-wide association analyses, gene variants associated with circulating n-3 PUFA were found to explain 11.4% of the variance of DPA, much higher than that for EPA or DHA (<1%) [31], highlighting that DPA levels are related to fatty acid metabolism. Our current findings suggest that there may be a potential beneficial (inverse) association between DPA and T2D, but further biochemical and clinical research is needed to examine a causal role of DPA on glycaemic outcomes and regulation of DPA synthesis.

Our findings provide convincing evidence for an inverse association between the plant food-based circulating ALA and T2D. This is in line with supportive evidence from randomised trials that ALA or flaxseed oil (rich in ALA) intervention could improve glucose homeostasis [32-35] and also in line with prior weak evidence of dietary ALA showing nonsignificant inverse association with T2D with relative risk 0.93 (95% CI 0.83-1.04) [6]. Conversion of ALA (an essential fatty acid from the diet) to long-chain n-3 PUFAs is not likely to explain the observation, because the conversion is small (0.2%-8% of ALA) [36] and because neither EPA nor DHA was associated with T2D. The precise mechanisms of how ALA may influence diabetes risk are unknown, but there is evidence suggesting that ALA may induce insulin secretion through direct actions on G-protein receptors expressed in pancreatic [Beta]-cells and through stimulating enteroendocrine L-cells, as well as by enhancing insulin sensitivity through hepatic insulin-like growth factor-1 (IGF-1) related pathways [37].

Findings in Context: n-6 Fatty Acids

EPIC-InterAct findings provide convincing evidence for a strong inverse association between circulating LA (18:2n6), the most abundant of all PUFAs, and T2D, with a 20% lower incidence per SD higher LA level. A prevailing concern relates to LA being the metabolic precursor to proinflammatory metabolites, particularly AA-derived eicosanoids, though the deleterious effects are not definitively proven and may depend on a complex interplay of metabolic factors [15,38]. Our observation in EPIC-InterAct that AA (20:4n6) is not associated--but intermediate metabolites GLA (18:3n6), DGLA (20:3n6), and estimated [delta]6 desaturase activity are adversely (positively) associated--is in line with a previous Mendelian randomisation study suggesting a potential causal adverse association of estimated [delta]6 desaturase activity (desaturation of LA[right arrow] GLA) with T2D, but not for [delta]5 desaturase activity (DGLA[right arrow]AA) [9]. Despite the causal inference for estimated [delta]6 desaturase activity on T2D risk, the biological mechanisms are speculated but not yet proven [9].

There is little other evidence to compare the EPIC-InterAct findings of the three relatively low-abundance n-6 PUFAs that were associated with lower risk (EDA, 20:2n6) or higher risk (DTA, 22:4n6, and n6-DPA, 22:5n6) of T2D.

Strengths and Limitations

Key strengths of EPIC-InterAct are its large sample size, prospective study design, and long follow-up. Moreover, the inclusion of EPIC-InterAct populations from eight European countries with diverse dietary intakes, comprehensive T2D case ascertainment and verification, the adjustment for many relevant potential confounders, and a series of sensitivity analyses enabled us to report robust findings, though we cannot rule out residual confounding. Our meta-analysis of other available studies enabled us to provide a systematic comparison with EPIC-InterAct findings and to place them in context of the wider evidence. We did not seek to provide pooled summary estimates that included findings from EPIC-InterAct as its large size with n = 12,132 T2D cases far outweighs the incident cases contributed by the other studies (ranging between 71 to 2,499 cases).

Several limitations merit attention. Our measurements of PUFAs might have involved misclassification because of measurement only at baseline and long-term storage of plasma samples. However, samples stored at -196°C are likely to be stable over time [22], and any errors were likely to be at random with respect to case status. We could only examine relative and not absolute concentrations of fatty acids, but this is a valid approach adopted commonly in epidemiological research and tends to provide better interpretation of metabolic inter-relationships [22]. Plasma phospholipid PUFAs may reflect both dietary and metabolic influences, but we could not distinguish these. We were unable to examine correlations between dietary and circulating individual PUFAs since in EPIC-InterAct we did not have dietary data on individual PUFAs. However, in general the correlations with circulating levels are low to modest because of measurement error issues of self-report and the linkage with food composition databases, as well as the different time frames of assessment between recall of habitual intake in a food frequency questionnaire (FFQ) versus circulating blood fatty acids. In contrast, there is evidence from dietary intervention trials that essential PUFAs in the blood (n-3 PUFA ALA, and n-6 PUFA LA) reflect dietary intakes of PUFAs [22]. However, we also need to note that circulating PUFAs may be exchanged across tissues and metabolised to diverse molecules [28]. Our clinical case ascertainment could be limited by potential misclassification due to false positive or negative diagnoses. However, we minimised false positives by applying rigorous verification criteria, while potential false negatives due to undiagnosed incident diabetes can be assumed to be nondifferential with regard to the exposure and hence any potential bias would be unlikely to alter our conclusions based on a relative scale (HR) in our analysis. Being based in European populations, our findings may have limited generalisability to other populations. Our comparative meta-analysis may be limited by using risk estimates per study-specific SD because the biological meaning of one SD may vary by studies. Thus, the estimates should be interpreted cautiously as risks per relative ranking on average, though the per-SD estimates are also advantageous given differences in lipid fractions, the number of fatty acids assessed, and laboratory methods. These differences in PUFA measurements across different cohorts point to the advantage of the standardised laboratory methods for measuring PUFA in EPIC-InterAct.

Implications

The lack of association of circulating EPA and DHA with T2D in EPIC-InterAct differs from the evidence for their inverse association with cardiovascular disease [39]. Whether or not dietary long-chain n-3 PUFAs reduce cardiovascular events is currently unresolved [39,40], but their overall beneficial effects on lipids (e.g., lowering of triglycerides), together with proposed anti-inflammatory, antiarrhythmic, and antithrombotic effects, and effects on improved endothelial function and blood pressure are plausible mechanisms [26]. In contrast, our current report of a lack of association between circulating EPA or DHA and incident T2D is supported by findings from randomised trials that found no evidence of net benefit or harm of long-chain n-3 supplementation on glycaemia or insulin sensitivity [2,3,26]. Therefore, our current findings advance scientific understanding of the possible differential links between circulating n-3 fatty acids and different health endpoints. Moreover, our findings for DPA and for ALA, both being inversely associated with incident T2D, deserve further follow-up to understand mechanisms and causality.

There has been an ongoing debate about the pros and cons of a diet high in the n6-PUFA LA (18:2n6), mainly for cardiovascular health [14,39,41]. This is important to resolve as LA is ubiquitous, being present in a wide variety of foods [38]. The strong evidence that dietary LA influences circulating LA has been reviewed previously [22], and in that context our finding of the inverse association between circulating LA and T2D incidence may have potential implications for the development of dietary advice. However, future work is warranted to build a body of evidence using diverse study designs that compare dietary LA, circulating LA, and their metabolites for the prevention of T2D.

Conclusion

We provide robust large-scale evidence that circulating plasma phospholipid n-3 PUFA ALA and n-6 PUFA LA concentrations are associated with lower incidence of T2D. In contrast four n-6 PUFAs (GLA, DGLA, DTA, and n-6DPA) are associated with higher T2D incidence in EPIC-InterAct. These findings (i) highlight that it is important to consider individual PUFAs rather than focusing on overall circulating n-3 or n-6 PUFA groups; (ii) provide robust evidence that circulating LA, the most abundant PUFA, is inversely associated with T2D; and (iii) meaningfully advance our understanding that some fatty acids (EPA and DHA) may be related differently with T2D or cardiovascular disease, while stimulating a scientific debate about the potential role of other, previously less well studied individual fatty acids.

Supporting Information

S1 Fig. Forest plots showing the country-specific HRs and 95% CIs for the association between per 1 SD of (A) total and individual n-3 PUFAs (ALA, EPA, DPA, and DHA) and T2D and (B) total and individual n-6 PUFAs (LA, GLA, EDA, DGLA, AA, DTA, and n6-DPA), as well as (C) fatty acid ratios and T2D from random-effects meta-analysis across the eight countries of EPIC-InterAct.

(PDF)

S2 Fig. Flow diagram of systematic literature search and identification of studies for the prospective associations of circulating concentrations of n-3 and n-6 PUFAs with incident T2D.

(PDF)

S3 Fig. Funnel plot for the prospective association of 20:3n6/18:2n6 with incident T2D.

(PDF)

S1 PRISMA. Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) checklist

(DOC)

S1 STROBE. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) checklist.

(DOC)

S1 Table. The distribution of individual and total PUFAs in the study subcohort by age, sex, BMI, and by case or noncase status--EPIC-InterAct study.

(DOC)

S2 Table. The distribution of individual and total PUFAs in the study subcohort by country--EPIC-InterAct study.

(DOC)

S3 Table. Sensitivity analyses for the association between individual and total PUFAs and fatty acid ratios and incident T2D: EPIC-InterAct Study.

(DOC)

S4 Table. Pooled HR and 95% CI for associations between quintiles of individual and total PUFAs and fatty acid ratios and incident T2D: EPIC-InterAct Study.

(DOC)

S5 Table. Study characteristics of prospective studies examining the association of n-3 and n-6 PUFAs with incident T2D.

(DOC)

S6 Table. Meta-analysis results by lipid fraction for the prospective associations of circulating fatty acids with incidence of T2D.

(DOC)

S7 Table. Associations of plasma phospholipid n-3 and n-6 PUFAs with incidence of T2D for Models 1, 2 and 3.

(DOC)

S1 Text. Methods for systematic review and meta-analysis of published literature.

(DOC)

S2 Text. Analysis plan (August 2014).

(PDF)

Acknowledgments

We thank all EPIC participants and staff for their contribution to the study. We thank the lab teams at the Medical Research Council (MRC) Epidemiology Unit for sample management and at MRC Human Nutrition Research for biochemical analysis and Nicola Kerrison for data management.

Author Contributions

Conceived and designed the experiments: NGF FI SJS AK MBS JZ ZY IS MG JMH JK LYW KS JLG EJMF AA PA HB CD GF PWF CG RK TJK KTK TK LMM PMN KO VP DP SP JRQ MRB OR CS AS NS AMWS AT MJT RT DLvA YTvdS CL ER NJW. Performed the experiments: NGF FI SJS AK MBS JZ ZY IS MG JMH JK LYW KS JLG EJMF AA PA HB CD GF PWF CG RK TJK KTK TK LMM PMN KO VP DP SP JRQ MRB OR CS AS NS AMWS AT MJT RT DLvA YTvdS CL ER NJW. Analyzed the data: SJS FI. Contributed reagents/materials/analysis tools: NGF FI SJS AK MBS JZ ZY IS MG JMH JK LYW KS JLG EJMF AA PA HB CD GF PWF CG RK TJK KTK TK LMM PMN KO VP DP SP JRQ MRB OR CS AS NS AMWS AT MJT RT DLvA YTvdS CL ER NJW. Wrote the first draft of the manuscript: NGF. Contributed to the writing of the manuscript: NGF FI SJS AK MBS JZ ZY IS MG JMH JK CL NJW. Enrolled patients: NGF FI SJS MBS JZ ZY IS MG JMH JK EJMF AA PA HB CD GF PWF CG RK TJK KTK TK LMM PMN KO VP DP SP JRQ MRB OR CS AMWS AT MJT RT DLvA YTvdS CL ER NJW. Agree with the manuscript's results and conclusions: NGF FI SJS AK MBS JZ ZY IS MG JMH JK LYW KS JLG EJMF AA PA HB CD GF PWF CG RK TJK KTK TK LMM PMN KO VP DP SP JRQ MRB OR CS AS NS AMWS AT MJT RT DLvA YTvdS CL ER NJW. All authors have read, and confirm that they meet, ICMJE criteria for authorship.

References

Eyre H, Kahn R, Robertson RM, Clark NG, Doyle C, Hong Y, et al. (2004) Preventing cancer, cardiovascular disease, and diabetes: a common agenda for the American Cancer Society, the American Diabetes Association, and the American Heart Association. Circulation 109: 3244-3255. pmid:15198946 doi: 10.1161/01.cir.0000133321.00456.00

Akinkuolie AO, Ngwa JS, Meigs JB, Djousse L (2011) Omega-3 polyunsaturated fatty acid and insulin sensitivity: a meta-analysis of randomized controlled trials. Clin Nutr 30: 702-707. doi: 10.1016/j.clnu.2011.08.013. pmid:21959352

Hartweg J, Perera R, Montori V, Dinneen S, Neil HA, Farmer A (2008) Omega-3 polyunsaturated fatty acids (PUFA) for type 2 diabetes mellitus. Cochrane Database Syst Rev: CD003205. doi: 10.1002/14651858.CD003205.pub2. pmid:18254017

Xun P, He K (2012) Fish Consumption and Incidence of Diabetes: meta-analysis of data from 438,000 individuals in 12 independent prospective cohorts with an average 11-year follow-up. Diabetes Care 35: 930-938. doi: 10.2337/dc11-1869. pmid:22442398

Wallin A, Di Giuseppe D, Orsini N, Patel PS, Forouhi NG, Wolk A (2012) Fish Consumption, Dietary Long-chain n-3 Fatty Acids, and Risk of Type 2 Diabetes: Systematic review and meta-analysis of prospective studies. Diabetes Care 35: 918-929. doi: 10.2337/dc11-1631. pmid:22442397

Wu JH, Micha R, Imamura F, Pan A, Biggs ML, Ajaz O, et al. (2012) Omega-3 fatty acids and incident type 2 diabetes: a systematic review and meta-analysis. Br J Nutr 107 Suppl 2: S214-S227. doi: 10.1017/S0007114512001602. pmid:22591895

Zheng JS, Huang T, Yang J, Fu YQ, Li D (2012) Marine N-3 polyunsaturated fatty acids are inversely associated with risk of type 2 diabetes in Asians: a systematic review and meta-analysis. PLoS ONE 7: e44525. doi: 10.1371/journal.pone.0044525. pmid:22984522

National Diet and Nutrition Survey. Results from Years 1, 2, 3 and 4 (combined) of the Rolling Programme (2008/2009-2011/2012). London. 2014. https://www.gov.uk/government/uploads/system/uploads/attachment_data/file/310995/NDNS_Y1_to_4_UK_report.pdf

Kroger J, Zietemann V, Enzenbach C, Weikert C, Jansen EH, Doring F, et al. (2011) Erythrocyte membrane phospholipid fatty acids, desaturase activity, and dietary fatty acids in relation to risk of type 2 diabetes in the European Prospective Investigation into Cancer and Nutrition (EPIC)-Potsdam Study. AmJ ClinNutr 93: 127-142. doi: 10.3945/ajcn.110.005447

Alhazmi A, Stojanovski E, McEvoy M, Garg ML (2014) Macronutrient intake and type 2 diabetes risk in middle-aged Australian women. Results from the Australian Longitudinal Study on Women's Health. Public Health Nutr 17: 1587-1594. doi: 10.1017/S1368980013001870. pmid:23866795

Song Y, Manson JE, Buring JE, Liu S (2004) A prospective study of red meat consumption and type 2 diabetes in middle-aged and elderly women: the women's health study. Diabetes Care 27: 2108-2115. pmid:15333470 doi: 10.2337/diacare.27.9.2108

van Dam RM, Willett WC, Rimm EB, Stampfer MJ, Hu FB (2002) Dietary fat and meat intake in relation to risk of type 2 diabetes in men. Diabetes Care 25: 417-424. pmid:11874924 doi: 10.2337/diacare.25.3.417

Salmeron J, Hu FB, Manson JE, Stampfer MJ, Colditz GA, Rimm EB, et al. (2001) Dietary fat intake and risk of type 2 diabetes in women. AmJClinNutr 73: 1019-1026

Ramsden CE, Zamora D, Leelarthaepin B, Majchrzak-Hong SF, Faurot KR, Suchindran CM, et al. (2013) Use of dietary linoleic acid for secondary prevention of coronary heart disease and death: evaluation of recovered data from the Sydney Diet Heart Study and updated meta-analysis. BMJ 346: e8707. doi: 10.1136/bmj.e8707. pmid:23386268

Choque B, Catheline D, Rioux V, Legrand P (2014) Linoleic acid: between doubts and certainties. Biochimie 96: 14-21. doi: 10.1016/j.biochi.2013.07.012. pmid:23900039

Ramsden CE, Zamora D, Majchrzak-Hong S, Faurot KR, Broste SK, Frantz RP, et al. (2016) Re-evaluation of the traditional diet-heart hypothesis: analysis of recovered data from Minnesota Coronary Experiment (1968-73). BMJ 353: i1246. doi: 10.1136/bmj.i1246. pmid:27071971

Hodge AM, English DR, O'Dea K, Sinclair AJ, Makrides M, Gibson RA, et al. (2007) Plasma phospholipid and dietary fatty acids as predictors of type 2 diabetes: interpreting the role of linoleic acid. AmJ ClinNutr 86: 189-197

Patel PS, Sharp SJ, Jansen E, Luben RN, Khaw KT, Wareham NJ, et al. (2010) Fatty acids measured in plasma and erythrocyte-membrane phospholipids and derived by food-frequency questionnaire and the risk of new-onset type 2 diabetes: a pilot study in the European Prospective Investigation into Cancer and Nutrition (EPIC)-Norfolk cohort. AmJ ClinNutr 92: 1214-1222. doi: 10.3945/ajcn.2010.29182

Wang L, Folsom AR, Zheng ZJ, Pankow JS, Eckfeldt JH (2003) Plasma fatty acid composition and incidence of diabetes in middle-aged adults: the Atherosclerosis Risk in Communities (ARIC) Study. AmJ ClinNutr 78: 91-98. doi: 10.1016/s0939-4753(03)80029-7

Krachler B, Norberg M, Eriksson JW, Hallmans G, Johansson I, Vessby B, et al. (2008) Fatty acid profile of the erythrocyte membrane preceding development of Type 2 diabetes mellitus. NutrMetab CardiovascDis 18: 503-510. doi: 10.1016/j.numecd.2007.04.005

Langenberg C, Sharp S, Forouhi NG, Franks PW, Schulze MB, Kerrison N, et al. (2011) Design and cohort description of the InterAct Project: an examination of the interaction of genetic and lifestyle factors on the incidence of type 2 diabetes in the EPIC Study. Diabetologia 54: 2272-2282. doi: 10.1007/s00125-011-2182-9. pmid:21717116

Hodson L, Skeaff CM, Fielding BA (2008) Fatty acid composition of adipose tissue and blood in humans and its use as a biomarker of dietary intake. Prog Lipid Res 47: 348-380. doi: 10.1016/j.plipres.2008.03.003. pmid:18435934

Wang LY, Summerhill K, Rodriguez-Canas C, Mather I, Patel P, Eiden MY S., et al. (2013) Development and validation of a robust automated analysis of plasma phospholipid fatty acids for metabolic phenotyping of large epidemiological studies. Genome Med 5: 39. doi: 10.1186/gm443. pmid:23618465

Higgins JP, Thompson SG, Deeks JJ, Altman DG (2003) Measuring inconsistency in meta-analyses. BMJ 327: 557-560. pmid:12958120 doi: 10.1136/bmj.327.7414.557

Greenland S, Longnecker MP (1992) Methods for trend estimation from summarized dose-response data, with applications to meta-analysis. Am J Epidemiol 135: 1301-1309. pmid:1626547

Balk EM, Lichtenstein AH, Chung M, Kupelnick B, Chew P, Lau J (2006) Effects of omega-3 fatty acids on serum markers of cardiovascular disease risk: a systematic review. Atherosclerosis 189: 19-30. pmid:16530201 doi: 10.1016/j.atherosclerosis.2006.02.012

Forouhi NG, Koulman A, Sharp SJ, Imamura F, Kroger J, Schulze MB, et al. (2014) Differences in the prospective association between individual plasma phospholipid saturated fatty acids and incident type 2 diabetes: the EPIC-InterAct case-cohort study. Lancet Diabetes Endocrinol 2: 810-818. doi: 10.1016/S2213-8587(14)70146-9. pmid:25107467

Kaur G, Cameron-Smith D, Garg M, Sinclair AJ (2011) Docosapentaenoic acid (22:5n-3): a review of its biological effects. Prog Lipid Res 50: 28-34. doi: 10.1016/j.plipres.2010.07.004. pmid:20655949

Mahendran Y, Agren J, Uusitupa M, Cederberg H, Vangipurapu J, Stancakova A, et al. (2014) Association of erythrocyte membrane fatty acids with changes in glycemia and risk of type 2 diabetes. AmJClinNutr 99: 79-85. doi: 10.3945/ajcn.113.069740. pmid:24153340

Lankinen MA, Stancakova A, Uusitupa M, Agren J, Pihlajamaki J, Kuusisto J, et al. (2015) Plasma fatty acids as predictors of glycaemia and type 2 diabetes. Diabetologia 58: 2533-2544. doi: 10.1007/s00125-015-3730-5. pmid:26277381

Lemaitre RN, Tanaka T, Tang W, Manichaikul A, Foy M, Kabagambe EK, et al. (2011) Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium. PLoS Genet 7: e1002193. doi: 10.1371/journal.pgen.1002193. pmid:21829377

Bloedon LT, Balikai S, Chittams J, Cunnane SC, Berlin JA, Rader DJ, et al. (2008) Flaxseed and cardiovascular risk factors: results from a double blind, randomized, controlled clinical trial. J Am Coll Nutr 27: 65-74. pmid:18460483 doi: 10.1080/07315724.2008.10719676

Rhee Y, Brunt A (2011) Flaxseed supplementation improved insulin resistance in obese glucose intolerant people: a randomized crossover design. Nutr J 10: 44. doi: 10.1186/1475-2891-10-44. pmid:21554710

Hutchins AM, Brown BD, Cunnane SC, Domitrovich SG, Adams ER, Bobowiec CE (2013) Daily flaxseed consumption improves glycemic control in obese men and women with pre-diabetes: a randomized study. Nutr Res 33: 367-375. doi: 10.1016/j.nutres.2013.02.012. pmid:23684438

Gomes PM, Hollanda-Miranda WR, Beraldo RA, Castro AV, Geloneze B, Foss MC, et al. (2015) Supplementation of alpha-linolenic acid improves serum adiponectin levels and insulin sensitivity in patients with type 2 diabetes. Nutrition 31: 853-857. doi: 10.1016/j.nut.2014.12.028. pmid:25933493

Burdge G (2004) Alpha-linolenic acid metabolism in men and women: nutritional and biological implications. Curr Opin Clin Nutr Metab Care 7: 137-144. pmid:15075703 doi: 10.1097/00075197-200403000-00006

Bhaswant M, Poudyal H, Brown L (2015) Mechanisms of enhanced insulin secretion and sensitivity with n-3 unsaturated fatty acids. J Nutr Biochem 26: 571-584. doi: 10.1016/j.jnutbio.2015.02.001. pmid:25841249

Whelan J (2008) The health implications of changing linoleic acid intakes. Prostaglandins Leukot Essent Fatty Acids 79: 165-167. doi: 10.1016/j.plefa.2008.09.013. pmid:18990554

Chowdhury R, Warnakula S, Kunutsor S, Crowe F, Ward HA, Johnson L, et al. (2014) Association of dietary, circulating, and supplement fatty acids with coronary risk: a systematic review and meta-analysis. Ann Intern Med 160: 398-406. doi: 10.7326/M13-1788. pmid:24723079

Rizos EC, Ntzani EE, Bika E, Kostapanos MS, Elisaf MS (2012) Association between omega-3 fatty acid supplementation and risk of major cardiovascular disease events: a systematic review and meta-analysis. JAMA 308: 1024-1033. doi: 10.1001/2012.jama.11374. pmid:22968891

Al-Khudairy L, Hartley L, Clar C, Flowers N, Hooper L, Rees K (2015) Omega 6 fatty acids for the primary prevention of cardiovascular disease. Cochrane Database Syst Rev: CD011094. doi: 10.1002/14651858.CD011094.pub2. pmid:26571451