Estimates were pooled using random effects meta-analysis. See eFigure 1 in the Supplement for results using inverse-variance weights. CHS indicates Cardiovascular Health Study20; EPIC-Norfolk, European Prospective Investigation of Cancer (Norfolk)21; HPFS, Health Professionals Follow-up Study23; KIHD, Kuopio Ischaemic Heart Disease Risk Factor Study25; MCCS, Melbourne Collaborative Cohort Study26; MESA, Multi-Ethnic Study of Atherosclerosis27; NHS, Nurses’ Health Study30; NSHDS, Northern Sweden Health and Disease Study28,29; PHS, Physician’s Health Study31; SCHS, Singapore Chinese Health Study33; and ULSAM, Uppsala Longitudinal Study of Adult Men.35
Estimates were pooled using random effects meta-analysis. See eFigure 1 in the Supplement for results using inverse-variance weights. CHS indicates Cardiovascular Health Study20; EPIC-Norfolk, European Prospective Investigation of Cancer (Norfolk)21; HPFS, Health Professionals Follow-up Study23; KIHD, Kuopio Ischaemic Heart Disease Risk Factor Study25; MCCS, Melbourne Collaborative Cohort Study26; MESA, Multi-Ethnic Study of Atherosclerosis27; NHS, Nurses’ Health Study30; NSHDS, Northern Sweden Health and Disease Study28,29; PHS, Physician’s Health Study31; SCHS, Singapore Chinese Health Study33; SHHEC, Scottish Heart Health Extended Chart32; and ULSAM, Uppsala Longitudinal Study of Adult Men.35
Customize your JAMA Network experience by selecting one or more topics from the list below.
Del Gobbo LC, Imamura F, Aslibekyan S, et al. ω-3 Polyunsaturated Fatty Acid Biomarkers and Coronary Heart Disease: Pooling Project of 19 Cohort Studies. JAMA Intern Med. 2016;176(8):1155–1166. doi:10.1001/jamainternmed.2016.2925
Are seafood and plant-derived ω-3 fats related to first coronary heart disease (CHD) events?
In a consortium of 19 studies, biomarkers of seafood and plant-derived ω-3 fats were associated with a significantly lower risk of fatal CHD. In contrast, associations with nonfatal myocardial infarction were generally less robust.
Blood biomarkers of seafood and plant-derived omega-3 fats are linked to a modestly lower risk of dying from heart attacks. Consumption of ω-3–rich foods should be encouraged.
The role of ω-3 polyunsaturated fatty acids for primary prevention of coronary heart disease (CHD) remains controversial. Most prior longitudinal studies evaluated self-reported consumption rather than biomarkers.
To evaluate biomarkers of seafood-derived eicosapentaenoic acid (EPA; 20:5ω-3), docosapentaenoic acid (DPA; 22:5ω-3), and docosahexaenoic acid (DHA; 22:6ω-3) and plant-derived α-linolenic acid (ALA; 18:3ω-3) for incident CHD.
A global consortium of 19 studies identified by November 2014.
Available prospective (cohort, nested case-control) or retrospective studies with circulating or tissue ω-3 biomarkers and ascertained CHD.
Data Extraction and Synthesis
Each study conducted standardized, individual-level analysis using harmonized models, exposures, outcomes, and covariates. Findings were centrally pooled using random-effects meta-analysis. Heterogeneity was examined by age, sex, race, diabetes, statins, aspirin, ω-6 levels, and FADS desaturase genes.
Main Outcomes and Measures
Incident total CHD, fatal CHD, and nonfatal myocardial infarction (MI).
The 19 studies comprised 16 countries, 45 637 unique individuals, and 7973 total CHD, 2781 fatal CHD, and 7157 nonfatal MI events, with ω-3 measures in total plasma, phospholipids, cholesterol esters, and adipose tissue. Median age at baseline was 59 years (range, 18-97 years), and 28 660 (62.8%) were male. In continuous (per 1-SD increase) multivariable-adjusted analyses, the ω-3 biomarkers ALA, DPA, and DHA were associated with a lower risk of fatal CHD, with relative risks (RRs) of 0.91 (95% CI, 0.84-0.98) for ALA, 0.90 (95% CI, 0.85-0.96) for DPA, and 0.90 (95% CI, 0.84-0.96) for DHA. Although DPA was associated with a lower risk of total CHD (RR, 0.94; 95% CI, 0.90-0.99), ALA (RR, 1.00; 95% CI, 0.95-1.05), EPA (RR, 0.94; 95% CI, 0.87-1.02), and DHA (RR, 0.95; 95% CI, 0.91-1.00) were not. Significant associations with nonfatal MI were not evident per 1 SD. Across quintiles, lower risk of nonfatal MI was evident with EPA (RR, 0.71; 95% CI, 0.56-0.90) and ALA (RR, 0.87; 95% CI, 0.78-0.97), and lower risk of fatal CHD was evident with DPA (RR, 0.76; 95% CI, 0.65-0.90) and DHA (RR, 0.77; 95% CI, 0.64-0.89). Associations appeared generally stronger in phospholipids and total plasma. Restricted cubic splines did not identify evidence of nonlinearity in dose responses.
Conclusions and Relevance
On the basis of available studies of free-living populations globally, biomarker concentrations of seafood and plant-derived ω-3 fatty acids are associated with a lower incidence of fatal CHD.
A review1 of experimental studies and randomized clinical trials revealed a protective effect of ω-3 polyunsaturated fatty acids (PUFAs) on coronary heart disease (CHD) risk pathways and clinical risk factors. However, key controversies remain. First, randomized clinical trials using fish oil supplements have found mixed effects on CHD events.2-8 However, most trials provided supplements for a few years or less and were conducted in patients with preexisting CHD or at high risk for CHD and taking multiple cardiovascular drugs. Furthermore, the background dietary intake in these trials was not usually measured, and many participants may have had adequate intake from diet alone. Thus, their generalizability and relevance for long-term effects of dietary ω-3 PUFAs for primary CHD prevention are uncertain. Second, although a review9 of several prior observational studies found inverse associations between seafood-derived ω-3 PUFAs and CHD death, potential effects on nonfatal myocardial infarction (MI) or total CHD are less clear. This review9 additionally found that most studies also relied on self-reported dietary questionnaires, which may produce errors or bias in recall. A previous review1 also found that a handful of studies have measured objective circulating or tissue ω-3 PUFA levels, but such findings could be limited by publication bias. In addition, prior individual biomarker studies were generally underpowered to explore potentially relevant differences in effects depending on underlying participant characteristics, medication use, or genetic variation. Thus, uncertainties remain about the effects of ω-3 PUFAs on CHD.
Most studies10,11 have evaluated combined intakes or biomarker levels of long-chain ω-3 PUFAs. However, individual ω-3 PUFAs, including eicosapentaenoic acid (EPA; 20:5ω-3), docosapentaenoic acid (DPA; 22:5ω-3), and docosahexaenoic acid (DHA; 22:6ω-3), may have shared and complementary effects.10 Similarly, most prior research11 has focused on seafood-derived long-chain ω-3 PUFAs, and potential CHD effects of plant-derived α-linolenic acid (ALA; 18:3ω-3) are far less understood.
To address each of these important questions, we developed a global consortium of available studies with circulating or tissue ω-3 biomarkers and ascertained incident CHD. Our aims were to evaluate the associations of individual seafood- and plant-derived ω-3 PUFAs with incident total CHD, fatal CHD, and nonfatal MI. We also explored dose response and potential heterogeneity in effects according to key underlying participant characteristics, including medication use. Our primary hypothesis was that biomarkers of long-chain ω-3 PUFAs would be associated with a lower risk of incident fatal CHD but not nonfatal MI.
The Fatty Acids and Outcomes Research Consortium (FORCe) is an extension of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Fatty Acids Working Group,12,13 originally developed to investigate effects of genetic variation on biomarker fatty acid levels. Using direct expert contacts and by reviewing the literature, we identified large studies that had measured circulating or tissue biomarkers of ω-3 PUFAs in general populations and had ascertained incident CHD. We focused on prospective (cohort, nested case-control) studies; retrospective studies were included if fatty acids were measured in adipose tissue at the time of the first event, given the stability of adipose measures. Studies were asked to join the consortium and participate in standardized pooled analyses of biomarker fatty acids and clinical events. Of 20 identified studies by November 2014, only 1 study14 declined to join the consortium. All cohort participants gave written informed consent according to their local guidelines. All studies received approval from ethical oversight committees. Data were deidentified.
A standardized analysis protocol was developed and provided to each participating study, including harmonized definitions of populations, exposures, outcomes, covariates, effect modifiers, and methods for pooling across studies. These specific methods are detailed below. An experienced analyst from each study performed the analysis for that study using individual-level data and provided all results to the lead investigator (L.C.D.G.) in standardized electronic forms.
To understand the potential effects on primary prevention and minimize reverse causation, participants from each study were excluded if they had a prior history of MI, angina, coronary revascularization, or stroke. All remaining participants with measured ω-3 PUFA biomarker levels were included in the analyses.
Fatty acid biomarker concentrations were assessed in study-specific lipid compartments, including total plasma, erythrocyte and plasma phospholipids, cholesterol esters, triglycerides, and adipose tissue, with all reported as weight percentage of total fatty acids. Details of fatty acid measurement methods for each study are provided in the eMethods in the Supplement.
We evaluated total incident CHD, defined as fatal or nonfatal MI, CHD death, or sudden cardiac death; nonfatal MI, generally defined as chest pain with abnormal cardiac enzyme concentrations or serial electrocardiogram changes; and fatal CHD, defined as fatal MI, CHD death, or sudden cardiac death. Details of the methods for assessing and defining these outcomes in each study are provided in eTable 1 in the Supplement. We excluded soft CHD end points (eg, angina, revascularization) to minimize bias in ascertainment of these events.
Covariates were standardized across studies based on prespecified, harmonized definitions and categorizations, including age (continuous), sex (male, female), race (white, black, and other study-specific subgroups, if available), clinical center (study-specific categories, if applicable), body mass index (BMI [calculated as the weight in kilograms divided by height in meters squared], continuous), educational level (less than high school graduate, high school graduate, some college or vocational school, or college graduate), smoking (current, former, or never or current vs not current, if former not assessed), physical activity (quartiles of metabolic equivalents or, if metabolic equivalents unavailable, 4 categories of physical or leisure activity as defined in the study), alcohol intake (none, 1-6 drinks per week, 1-2 drinks per day, or >2 drinks per day; with 14 g of alcohol equaling 1 standard drink), diabetes mellitus (yes or no, defined as treatment with oral hypoglycemic agents, insulin, or by using fasting glucose levels, or study-specific definitions otherwise), treated hypertension (yes or no, defined as treatment with antihypertensive drugs or as diagnosed hypertension otherwise), treated hypercholesterolemia (yes or no, defined as treatment with lipid-lowering drugs or as diagnosed hypercholesterolemia otherwise), regular aspirin use (yes or no, categorized as ≥3 times per week), and biomarker concentrations of ω-6 PUFA linoleic acid (LA; 18:2ω-6), arachidonic acid (AA; 20:4ω-6), and total trans fatty acids (each continuous). For categorical variables, missing covariate data were included as an indicator category.
Standardized study-specific analyses were performed for each study using individual-level data. Because of variability in the total number of fatty acids assessed across the different assays (12-58 fatty acids) and to permit comparison and pooling of findings across the different biomarker compartments, each fatty acid was analyzed continuously per 1-SD increase and per study-specific quintiles, performing meta-analysis within each quintile and comparing each quintile-specific pooled result with the first quintile.
For prospective cohorts, Cox proportional hazards regression models estimated the hazard ratio for incident total CHD, nonfatal MI, and fatal CHD, with follow-up from the date of biomarker measurement to date of event, end of follow-up, loss to follow-up, or death, whichever occurred first. For prospective and retrospective matched case-control studies, conditional logistic regression was used to estimate the hazard ratio (if matched on time to event) or odds ratio. All analyses used robust SEs.
For pooled meta-analyses, hazard ratios and odds ratios were considered to approximate relative risks (RRs). Our prespecified primary analysis used inverse-variance weights. Owing to observed moderate heterogeneity in some analyses, we report findings using random effects models using the method of DerSimonian and Laird.15 Results were first pooled separately for each biomarker compartment and then pooled together. If more than 1 biomarker had been measured in a given study, main results included the biomarker that might best represent diet: for long-chain ω-3 fatty acids, total plasma or serum measures were used, and for ALA, adipose tissue measures were used in the main results.16
To statistically test potential nonlinear dose responses, we performed multivariate meta-regression, modeling restricted cubic splines.17 These analyses evaluated total plasma and phospholipids; cholesterol esters and adipose tissue fatty acids were not evaluated using splines because of fewer numbers of studies in these compartments.
Overall heterogeneity was assessed using the I2 statistic, with heterogeneity considered low if the I2 was less than 35%. Heterogeneity by prespecified subgroups was evaluated by pooling individual study results based on defined strata, including age (<60 or ≥60 years), sex (male or female), race/ethnicity (white, black, Hispanic, or Chinese), LA and AA biomarker concentrations (less than or greater than or equal to the median value in each study), type 2 diabetes (yes or no), statin use (yes or no), regular aspirin use (yes or no), year of biomarker sampling (before 2000 or 2000 or later), and (for ALA) EPA, DPA, and DHA biomarker concentration (less than or greater than or equal to the median value in each study). Effect estimates in each study-specific strata were pooled, and the statistical significance of differences between subgroups of potential sources of heterogeneity was assessed using meta-regression. In a meta-analysis13 of 7 studies with available genetic data, we also examined potential interaction by variants in FADS (OMIM 208150) desaturase genes (rs174546, rs968567), pooling study-specific estimates using an additive genetic model (eMethods in the Supplement). Within each study, interaction terms for each single-nucleotide polymorphism (SNP) were constructed as a cross-product term of the ω-3 PUFA biomarker (continuous) by the SNP (ordinal: 0, 1, or 2 T alleles), including the main effects in the model; these interaction terms were pooled using meta-analysis.
Meta-analyses were performed using STATA statistical software, version 12 (StataCorp). We considered 2-tailed P < .05 to be statistically significant.
Several sensitivity analyses were performed. We also performed continuous meta-analysis using absolute (percentage of fatty acids), rather than study-specific, units. To minimize potential reverse causation attributable to preexisting subclinical disease, we excluded cases identified in the first 2 years after biomarker sampling. To minimize exposure misclassification attributable to within-person changes in biomarker levels over time, we censored participants after the first 6 years of follow-up. We performed additional analyses limited to prospective studies only and those without self-reported CHD events.
The 19 studies18-35 included 45 637 unique participants from 16 countries, including the United States, Australia, Costa Rica, Finland, France, Germany, Ireland, Israel, Italy, Norway, Singapore, the Soviet Union, Spain, Sweden, Switzerland, and the United Kingdom (Table 1 and eMethods in the Supplement). Most were prospective cohort18,20,25-27,32,34,35 (n = 10 studies) or prospective nested case-control21,23,24,28-31,33 (n = 7) studies; 2 were retrospective case-control studies19,22 that used adipose tissue biomarkers. Biomarker types included phospholipids (plasma or erythrocyte) (n = 10), total plasma (n = 6), adipose tissue (n = 4), and cholesterol esters (n = 1); 2 studies (Nurses’ Health Study [NHS]30 I and Health Professionals Follow-up Study [HPFS]23) measured total plasma and erythrocyte phospholipids. As expected, ω-3 PUFA concentrations and distributions varied across biomarker types and study assays (eTable 2 and eTable 3 in the Supplement), with coefficients of variation less than 10% for most fatty acids and biomarkers (eMethods in the Supplement). For example, median ALA concentrations were generally higher in adipose (0.72% of total fatty acids) than in circulating biomarkers (0.20% of total fatty acids), whereas EPA, DPA, and DHA concentrations were generally higher in phospholipids than in other compartments, consistent with prior reports.16,36
Median age at baseline was 59 years (range, 18-97 years), and 28 660 participants (62.8%) were male (eTable 4 in the Supplement). During a median 10 years of follow-up (range, 1.3-42 years), 7973 incident CHD events, 7157 nonfatal MIs, and 2781 fatal CHD events occurred. Because not all studies ascertained every exposure and outcome (Table 1), cases of nonfatal MI and fatal CHD do not sum to the number of total CHD cases. Medication use varied across studies; by design, all participants in the analysis were free of prevalent CHD or stroke (eTable 5 in the Supplement). Most participants were white, but the Cardiovascular Health Study (CHS)20 and Multi-Ethnic Study of Atherosclerosis (MESA)27 included relatively higher proportions of African American individuals (12.3% and 24.7%, respectively); MESA and the Costa-Rican Case-Control Study of Myocardial Infarction (Costa Rican adults study)19 included Hispanic individuals (24.9% and 100%, respectively), and MESA27 and the Singapore Chinese Health Study (SCHS)33 included Chinese participants (25.0% and 100%, respectively). Across all studies, the median BMI was 26, and generally up to 30% of participants were current smokers, except for higher rates in the Scottish Heart Health Extended Cohort32 (44.6%) and Uppsala Longitudinal Study of Adult Men35 50 (51.0%). Alcohol intake was modest, with most participants consuming up to 1 drink per day. Fish oil supplementation was infrequently assessed across studies; use was low (0%-4% of participants) in 5 of 6 studies with these data (eTable 6 in the Supplement); only the European Prospective Investigation of Cancer (Norfolk) (EPIC-Norfolk)21 had a higher prevalence (33%).
In continuous (per 1-SD increase) multivariable-adjusted analyses, each ω-3 PUFA was associated with an approximately 9% lower risk of fatal CHD, with RRs of 0.91 (95% CI, 0.82-1.00) for EPA, 0.90 (95% CI, 0.85-0.96) for DPA, and 0.90 (95% CI, 0.84-0.96) for DHA (Figure 1 and Figure 2). There was moderate heterogeneity in this association for EPA and DHA (I2 = 37% for both) and low heterogeneity for DPA and DHA (I2 = 0% and 31%, respectively) (Figure 1 and Figure 2). The sum of EPA, DPA, and DHA was associated with an 11% lower risk of fatal CHD (eFigure 2 in the Supplement). By contrast, none of the long-chain ω-3 PUFAs were significantly associated with nonfatal MI (Table 2). We found that DPA, but not EPA or DHA, was associated with a significantly lower risk of total CHD, with RRs of 0.94 (95% CI, 0.90-0.99), 0.94 (95% CI, 0.87-1.02), and 0.95 (95% CI, 0.91-1.00), respectively (Table 2).
Some differences were observed when EPA, DPA, and DHA were evaluated across quintiles (eTable 7 in the Supplement). Upper quintiles of EPA and DHA were associated with a lower risk of nonfatal MI (quintile 5 vs quintile 1 comparison: RR, 0.71 [95% CI, 0.56-0.90] vs 0.87 [95% CI, 0.78-0.97]). Top quintiles of DPA and DHA were associated with a lower risk of fatal CHD (quintile 5 vs quintile 1 comparison: RR, 0.76 [95% CI, 0.65-0.90] vs 0.77 [95% CI, 0.64-0.89]) (eTable 7 in the Supplement). Across biomarker types, inverse associations were generally stronger in the phospholipid and total plasma compartments (Table 2); heterogeneity decreased when adipose tissue and cholesterol ester estimates were excluded from meta-analyses. All findings were generally similar or stronger when using inverse-variance weights (eTable 7 and eFigures 1 and 3-5 in the Supplement).
In continuous analysis, ALA was associated with a 9% lower risk of fatal CHD (RR, 0.91 [95% CI, 0.84-0.98]) (Figure 1) but not total CHD or nonfatal MI. For ALA, no single biomarker compartment consistently had stronger associations across outcomes.
Restricted cubic splines did not identify evidence of nonlinear associations of any of the ω-3 PUFA biomarkers and outcomes in total plasma or phospholipids (P for nonlinearity > .05 for each) (eFigure 6 in the Supplement).
No significant differences in associations of ω-3 PUFA biomarkers with incident CHD events were observed by age, sex, ω-6 PUFA (LA or AA) concentrations, type 2 diabetes status, statin use, regular aspirin use, year of biomarker sampling, or (for ALA) EPA, DPA, and DHA concentrations (P for heterogeneity > .05 for each) (eTable 8 in the Supplement). In a post hoc meta-regression, we observed no significant differences in associations by median length of follow-up (<10 or ≥10 years) (P for heterogeneity > .05 for each). Effect modification was suggested by ethnicity. Compared with white individuals, ALA was associated with a significantly lower risk of nonfatal MI among African American individuals (P for heterogeneity = .001). For EPA, ALA was associated with a significantly lower risk of total CHD and nonfatal MI among Chinese participants (P for heterogeneity = .02 and .01, respectively). Among 7 studies10,17,19,33-35 with SNP data (eTable 9 in the Supplement), no significant interaction was identified by FADS desaturase gene variants (eTable 10 in the Supplement).
Compared with the main findings, no appreciable differences were observed after excluding cases identified in the first 2 years after biomarker sampling, censoring participants at the first 6 years of follow-up, excluding retrospective studies, or excluding studies with self-reported events (eTable 11 in the Supplement). For 2 studies23,30 that analyzed phospholipids in addition to total plasma (NHS I and HPFS), results were similar when phospholipid results from these studies were used in place of total plasma in the meta-analysis.
In this consortium which pooled individual-level harmonized analyses from 19 studies, including 45 637 unique participants and nearly 8000 first CHD events in 16 countries, the ω-3 biomarkers ALA, DPA, and DHA were associated with a modestly lower risk of fatal CHD. The magnitude of observed effect sizes for fatal CHD (approximately 9% per 1-SD increase) are consistent with findings for cardiac death from a meta-analysis8 of trials. In contrast, associations with nonfatal MI were generally less robust. Across these diverse studies, findings were consistent by age, sex, ω-6 PUFA LA and AA levels, year of biomarker sampling, presence or absence of diabetes, statin use, and aspirin use. By contrast, effect modification by ethnicity was observed for select exposure-outcome pairings. Separate analyses by biomarker compartment provided further insights, with generally stronger inverse associations in phospholipids and total plasma for EPA, DPA, and DHA. This investigation provides the most comprehensive estimates to date of the associations between seafood and plant-based ω-3 PUFAs, assessed using biomarkers, and primary incidence of CHD in generally healthy, free-living populations around the world.
In randomized clinical trials, long-chain ω-3 PUFAs benefit multiple major cardiovascular risk factors, including triglyceride levels, blood pressure, heart rate, heart rate variability, endothelial function, and myocardial oxygen demand.1,10 Compared with many other single nutrients, these demonstrated physiologic benefits provide strong biologic plausibility to support an effect on clinical events. Both EPA and DHA are also precursors to bioactive lipid metabolites, including specialized proresolving mediators36 and cytochrome P450–generated monoepoxides,37 that could contribute to lower CHD risk. Our findings are consistent with prior experimental evidence that long-chain ω-3 PUFAs may have membrane stabilizing actions in the setting of ischemia-induced ventricular fibrillation1,38 and observational evidence indicating that benefits of fish consumption are most related to arrhythmic events39 and fatal CHD.40
Use of biomarkers allowed separate investigation of each ω-3 PUFA. Of interest, EPA, DPA, and DHA were similarly associated with lower risk of fatal CHD. Biomarker levels of these fatty acids are only moderately interrelated (eg, r = 0.43, 0.51, and 0.13 for EPA and DHA, EPA and DPA, and DPA and DHA, respectively).41 Thus, observed associations for any 1 of these fatty acids are unlikely to be fully explained by the others. Compared with EPA and DHA, comparatively less is known about the molecular and physiologic effects of DPA; a previous study10 suggests that DPA may inhibit ex vivo collagen-stimulated platelet aggregation, thromboxane production, and cyclooxygenase 1 activity. In addition, whereas circulating and tissue levels of EPA and DHA are strongly influenced by dietary seafood consumption, DPA concentrations appear to be mainly derived from endogenous elongation of EPA; interconversion between DPA and DHA is very limited.10
Given diverse global dietary sources, the effects of ALA on CHD are of particular interest. However, a meta-analysis11 of prior reports of dietary and biomarker ALA and CHD risk revealed inconsistency, perhaps owing to methodologic and analytic differences in these investigations. A key strength of our analysis was the prespecified, harmonized analytic plan using individual-level data within each study, which allowed consistent assessment of how ALA relates to CHD. In addition, the inclusion of 19 studies minimized the effect of publication bias, wherein studies with positive findings are more likely reported. Compared with prior work,42 we also took advantage of a far larger number of events (eg, >4-fold larger for fatal CHD) and observed a similarly lower risk of fatal CHD for ALA as for long-chain ω-3 PUFAs. Mechanistically, these findings are supported by effects of ALA on thrombosis, inflammation, arrhythmia, and endothelial function.43-47 Our findings, combined with relative affordability, global accessibility, and sustainability of ALA,48 support the potential importance of ALA for improving global cardiovascular health.
The large number of cases in our analysis and inclusion of multiple ethnic subgroups allowed us to explore potential effect modification in ω-3 PUFAs and CHD outcome associations by race/ethnicity. The stronger effects among African American and Chinese individuals for select ω-3–outcome pairings could be attributable to differences in consumption patterns, true biologic diversity in pathophysiologic pathways, or chance. For instance, preparation methods among Chinese populations could differ from that among white populations, with raw or steamed fish in the former population vs more deep fried fish in the latter. Differential pathophysiologic pathways for cardiometabolic risk have also been documented among Asian vs white populations; for example, Asian individuals develop type 2 diabetes at a much lower BMI than white individuals,49,50 and at least 3 meta-analyses51-53 have reported a protective effect of seafood-derived ω-3 fatty acids on incident diabetes among Asian individuals, but not white individuals. Finally, our findings of effect modification by race could be attributable to chance; few studies were available with African American (n = 2 studies: CHS and MESA) and Chinese subgroups (n = 2 studies: SCHS and MESA).27,33 Our results highlight the need for further work to better understand differences in associations by race for ω-3 fatty acids and incident CHD phenotypes.
Randomized clinical trials of fish oil supplements have found mixed effects, although overall pooled findings indicate a benefit in lowering risk of cardiac death,8 consistent with our findings. Our investigation tests unique and separate questions from these trials attributable to differences in the ω-3 PUFA source (dietary in our study vs supplements in trials), population (primary prevention vs secondary prevention or high risk), and duration (habitual intake vs short-term supplementation). In cohort-harmonized, stratified analyses, we found little evidence of effect modification by statin or aspirin use, which have been hypothesized to reduce benefits of ω-3 fatty acids.54 We also found little evidence of effect modification by circulating ω-6 PUFAs (LA or AA), which have been hypothesized but never found in humans to reduce cardiovascular benefits of ω-3 PUFAs.55 Altogether, our findings suggest that seafood and plant-derived ω-3 PUFAs are beneficial for fatal CHD prevention across diverse population subgroups.
Our stratification by biomarker type provides new insights on which lipid compartments may most influence CHD. Adipose tissue has been suggested for measurement of fatty acids because it reflects long-term dietary intake (approximately 1 year).56 However, EPA and DHA are more highly concentrated in phospholipids than other compartments, and phospholipids also respond rapidly to dietary changes57,58 and may best reflect effects on membrane receptors.1 The best compartment for assessing ALA has also been unclear, with possible advantages to adipose tissue16 but relatively little prior evaluation of other compartments. Our findings suggest generally stronger inverse associations with CHD for EPA, DPA, and DHA in phospholipids and total plasma, although fewer estimates and cases were available for adipose and cholesterol esters compartments, decreasing precision. Our results support use of phospholipids or total plasma for long-chain ω-3 PUFA exposure and a need for further investigation of optimal biomarkers for assessing ALA.
Our analysis has several important strengths. Use of biomarkers provided measures of exposure free of recall error and allowed separate evaluation of different ω-3 PUFAs and different circulating and tissue lipid compartments. We used standardized definitions and modeling for the populations, exposures, outcomes, covariates, effect modifiers, and analysis, reducing heterogeneity and potential investigator bias. Most studies18-25,27-35 used centralized adjudication processes or registry linkage rather than self-report alone in ascertaining events, reducing information bias. Findings were similar in several sensitivity analyses, suggesting the robustness of results to varying assumptions. Our inclusion of all available studies substantially reduces the potential for publication bias. The studies included both sexes, multiple races, and a range of socioeconomic statuses across 16 countries, increasing generalizability.
Potential limitations should be considered. Relatively few studies19,22,32,35 were available for some lipid compartments (eg, adipose, cholesterol esters), limiting inference for these specific compartments. All cohorts assessed ω-3 PUFA exposure once at baseline, and changes over time would attenuate findings toward the null, causing underestimation of associations. Reduced statistical power in quintile analyses and spline analyses made it difficult to ascertain the specific shape of dose responses. Because we analyzed observational studies, biomarkers and risk factors were not evenly distributed by randomization of study participants; hence, unmeasured or residual confounding cannot be excluded, and ω-3 biomarker measures may in part be markers of a long-term healthy lifestyle. However, our findings were independent of a range of major cardiovascular risk factors, and the relative specificity for fatal CHD and known biologic effects of ω-3 PUFAs also argue against residual confounding as the sole explanation for our findings. Because our findings reflect habitual ω-3 exposure for primary CHD prevention, our results do not imply that fish oil supplementation will necessarily reduce CHD events, particularly in the context of short-term treatment or in those with preexisting CHD. Awareness of ω-3 PUFAs, laboratory ascertainment methods, and medical treatment of hypertension, hypercholesterolemia, diabetes, angina, and MI have changed over time and could contribute to heterogeneity of results; however, we observed no significant differences in the results for studies assaying ω-3 PUFA levels before vs after the year 2000.
Our pooled collaboration of global studies using biomarkers of ω-3 PUFAs in participants without prevalent CHD reveals that habitual consumption of seafood and plant-based ω-3 PUFAs is associated with a lower risk of fatal CHD.
Correction: This article was corrected on February 11, 2019, to correct the affiliation for Mohammad Y. Yakoob, PhD, and on October 3, 2016, for numerous errors in the Abstract and in reporting of analyses in the text and Supplement.
Accepted for Publication: May 3, 2016.
Corresponding Author: Liana C. Del Gobbo, PhD, Department of Medicine, Division of Cardiovascular Medicine, Stanford University School of Medicine, Stanford, CA 94305 (email@example.com).
Published Online: June 27, 2016. doi:10.1001/jamainternmed.2016.2925
Author Contributions: Dr Del Gobbo 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.
Study concept and design: All authors.
Acquisition, analysis, or interpretation of data: Del Gobbo conducted meta-analyses and data analysis for CHS. Aslibekyan performed analysis for the Costa-Rican adults study. Marklund performed analysis for ULSAM 50 and 70. Virtanen performed analysis for KIHD. Wennberg performed analysis for NSHDS I and II. Yakoob performed analysis for NHS I. Chiuve performed analysis for HPFS. dela Cruz performed analysis for MCCS. Frazier-Wood performed analysis for MESA. Guallar performed analysis for EURAMIC. Matsumoto performed analysis for PHS. Prem performed analysis for SCHS. Samieri performed analysis for the 3C study. Tanaka performed analysis for InCHIANTI. Wu performed analysis for SHHECS. Zhou performed analysis for ARIC. Imamura performed spline analyses. All authors interpreted data.
Drafting of the manuscript: Del Gobbo, Virtanen, Aslibekyan, Chaves, Fretts, Djousse, Mozaffarian.
Critical revision of the manuscript for important intellectual content: All authors.
Obtained funding: Helmer, Ingelsson, Yuan, Barberger-Gateau, Campos, Chaves, Djoussé, Giles, Gómez-Aracena, Hodge, Hu, Jansson, Johansson, Khaw, Koh, Lemaitre, Lind, Luben, Mozaffarian, Rimm, Risérus, Samieri, Franks, Siscovick, Stampfer, L. M. Steffen, B. T. Steffen, Tsai, van Dam, Voutilainen, Willett, Woodward.
Administrative, technical, or material support: All authors.
Study supervision: Mozaffarian.
Conflict of Interest Disclosures: Dr Barberger-Gateau reported receiving grants and nonfinancial support from Danone Research, Vifor Pharma, and Groupe Lipides et Nutrition; personal fees and nonfinancial support from Nutricia, ILSI Europe, and PiLeJe; and a grant from the Centre National Interprofessionnel de l'Industrie Laitière. Dr Del Gobbo reported receiving ad hoc consulting fees from the Life Sciences Research Organization. Dr Djoussé reported receiving investigator-initiated grants from Merck, Amarin Pharma Inc, California Walnut Commission, and the National Institutes of Health and serving as ad hoc consultant for Amarin Pharma Inc. Dr Helmer reported receiving fees for a conference from Novartis. Dr Mozaffarian reported receiving ad hoc honoraria from Bunge, Pollock Institute, and Quaker Oats; ad hoc consulting for Foodminds, Life Sciences Research Organization, Nutrition Impact, Amarin, AstraZeneca, Winston, and Strawn LLP; membership in Unilever North America Scientific Advisory Board; and chapter royalties from UpToDate. Dr Woodward reported receiving consulting fees from Amgen and Novartis and research funding from Sanofi. No other disclosures were reported.
Funding/Support: The Atherosclerosis Risk in Communities study was performed as a collaborative study supported by contracts HHSN268201100005C, HHSN268201100006C, HHSN268201100007C, HHSN268201100008C, HHSN268201100009C, HHSN268201100010C, HHSN268201100011C, HHSN268201100012C, R01HL087641, R01HL59367, and R01HL086694 from the National Heart, Lung, and Blood Institute; contract U01HG004402 from the National Human Genome Research Institute; and contract HHSN268200625226C from the National Institutes of Health. We thank the staff and participants of the Atherosclerosis Risk in Communities study for their important contributions. Infrastructure was partly supported by grant UL1RR025005, a component of the National Institutes of Health and National Institutes of Health Roadmap for Medical Research. The Cardiovascular Health Study (CHS) was supported by contracts HHSN268201200036C, HHSN268200800007C, N01HC55222, N01HC85079, N01HC85080, N01HC85081, N01HC85082, N01HC85083, and N01HC85086 and grant U01HL080295 from the National Heart, Lung, and Blood Institute, with additional contribution from the National Institute of Neurological Disorders and Stroke. Additional support was provided by grant R01AG023629 from the National Institute on Aging. A full list of principal CHS investigators and institutions can be found at CHS-NHLBI.org. The Costa Rican adult study was supported by grant R01HL081549 from the National Institutes of Health. The European Study on Antioxidants, Myocardial Infarction and Cancer was supported by the Commission of the European Communities, as a Concerted Action within Directorate General-XII, with additional support from Directorate General-V Europe against Cancer. The national studies were financed by the Dutch Ministry of Health. Ulster Cancer Foundation and Milk Intervention Board, grant AKT76 from Cancer Research Switzerland, grant 32-9257-87 from the Swiss National Science Foundation, The Spanish Fondo de Investigaciones Sanitarias and Ministry of Science and Education, and German Federal Health Office. EPIC-Norfolk was funded by grants from Medical Research Council and Cancer Research UK. Dr Imamura received support from grant MC_UU_12015/5 from the Medical Research Council Epidemiology Unit. The Health Professionals Follow-up Study (HPFS) was supported by grants UM1 CA167552, R01 HL35464, AA11181, HL35464, CA55075, HL60712, and P30 DK46200 from the National Institutes of Health. The InChianti study was supported by grant ICS 110.1\RS97.71 as a targeted project by the Italian Ministry of Health and in part by contracts N01-AG-916413, N01-AG-821336, 263 MD 9164 13, and 263 MD 821336 from the Intramural Research Program of the National Institutes of Health. The Kuopio Ischaemic Heart Disease Risk Factor Study was supported by grants 41471 and 1041086 from the Academy of Finland, Helsinki, Finland. The Melbourne Collaborative Cohort Study recruitment was funded by VicHealth and Cancer Council Victoria and was further supported by grants 209057, 251553, and 504711 from Australia’s National Health and Medical Research Council and by infrastructure provided by Cancer Council Victoria. Cases and their vital status were ascertained through the Victorian Cancer Registry and the Australian Institute of Health and Welfare, including the National Death Index and the Australian Cancer Database. The Multi-Ethnic Study of Atherosclerosis (MESA) and the MESA SHARe project are conducted and supported by the National Heart, Lung, and Blood Institute in collaboration with MESA investigators. Support for MESA is provided by contracts N01-HC-95159, N01-HC-95160, N01-HC-95161, N01-HC-95162, N01-HC-95163, N01-HC-95164, N01-HC-95165, N01-HC-95166, N01-HC-95167, N01-HC-95168, N01-HC-95169, UL1-TR-001079, and UL1-TR-000040 from the National Heart, Lung, and Blood Institute. Funding for SHARe genotyping was provided by contract N02-HL-64278 from the National Heart, Lung, and Blood Institute. Genotyping was performed at Affymetrix and the Broad Institute of Harvard and MIT using the Affymetrix Genome-Wide Human SNP Array 6.0. The Northern Sweden Health and Disease Studies I and II were supported by the Swedish Cancer Society and the Swedish Research Council. The Nurses Health Study was supported by research grants UM1 CA186107, R01 CA49449, R01 HL034594, P01CA87969, R01HL034594, and R01HL088521 from the National Institutes of Health. The Physician’s Health Study was supported by grants R21 HL088081, CA-34944, CA-40360, and CA-097193 from the National Cancer Institute and grants HL-26490 and HL-34595 from the National Heart, Lung, and Blood Institute. The Three-City study was conducted under a partnership agreement between the Institut National de la Santé et de la Recherche Médicale, the University Bordeaux 2 Victor Segalen, and Sanofi. The Fondation pour la Recherche Médicale funded the preparation and initiation of the study. The Three-City study was also supported by the Caisse Nationale Maladie des Travailleurs Salariés, Direction Générale de la Santé, MGEN, Institut de la Longévité, Conseils Régionaux d’Aquitaine et Bourgogne, Fondation de France, Ministry of Research–Institut National de la Santé et de la Recherche Médicale Programme Cohortes et collections de données biologiques, grant COGINUT ANR-06-PNRA-005 from the Agence Nationale de la Recherche, grant FCS 2009-2012 from the Fondation Plan Alzheimer, and the Caisse Nationale pour la Solidarité et l’Autonomie. Dr Samieri was sponsored by a grant from the Fondation Plan Alzheimer. The Scottish Heart Health Extended Cohort study was funded by the Scottish Health Department Chief Scientist Organization, British Heart Foundation, and FP Fleming Trust. The Singapore Chinese Health Study was supported by grant NMRC 1270/2010 from the Singapore National Medical Research Council and grants R01CA 144034 and UM1 CA182876 from the National Institutes of Health. The Uppsala Longitudinal Studies of Adult Men 50 and 70 were funded by the Swedish Research Council for Health, Working Life and Welfare, Uppsala City Council, and Swedish Research Council.
Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and the decision to submit the manuscript for publication.
Disclaimer: The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Additional Contributions: We acknowledge Roger Tavendale, PhD, for his work with the Scottish Heart Health Study.
Group Information: The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Fatty acids and Outcomes Research Consortium (FORCe) members are as follows: Liana C. Del Gobbo, PhD, Fumiaki Imamura, PhD, Stella Aslibekyan, PhD, Matti Marklund, PhD, Jyrki K. Virtanen, PhD, Maria Wennberg, PhD, Mohammad Y. Yakoob, PhD, Stephanie E. Chiuve, ScD, Luicito dela Cruz, PhD, Alexis C. Frazier-Wood, PhD, Amanda M. Fretts, MPH, PhD, Eliseo Guallar, PhD, Chisa Matsumoto, PhD, MD, Kiesha Prem, MSc, Tosh Tanaka, PhD, Jason H. Y. Wu, PhD, Xia Zhou, PhD, Catherine Helmer, MD, PhD, Erik Ingelsson, MD, PhD, Jian-Min Yuan, MD, PhD, Pascale Barberger-Gateau, PhD, Hannia Campos, PhD, Paulo H. M. Chaves, MD, PhD, Luc Djoussé, MD, ScD, Graham G. Giles, PhD, Jose Gómez-Aracena, PhD, Allison M. Hodge, PhD, Frank B. Hu, PhD, MD, MPH, Jan-Håkan Jansson, MD, PhD, Ingegerd Johansson, PhD, Kay-Tee Khaw, PhD, MD, Woon-Puay Koh, PhD, Rozenn N. Lemaitre, PhD, MPH, Lars Lind, PhD, Robert N. Luben, PhD, Eric B. Rimm, ScD, Ulf Risérus, PhD, MD, Cecilia Samieri, PhD, Paul W. Franks, PhD, David S. Siscovick, MPH, MD, Meir Stampfer, DrPH, MD, Lyn M. Steffen, PhD, MPH, Brian T. Steffen, PhD, Michael Y. Tsai, PhD, Rob M. van Dam, PhD, Sari Voutilainen, PhD, Walter C. Willett, DrPH, MD, Mark Woodward, PhD, and Dariush Mozaffarian, MD, DrPH.