Intended for healthcare professionals

CCBY Open access
Research

Association of omega 3 polyunsaturated fatty acids with incident chronic kidney disease: pooled analysis of 19 cohorts

BMJ 2023; 380 doi: https://doi.org/10.1136/bmj-2022-072909 (Published 18 January 2023) Cite this as: BMJ 2023;380:e072909
  1. Kwok Leung Ong, senior research fellow1,
  2. Matti Marklund, senior research fellow234,
  3. Liping Huang, senior data analyst2,
  4. Kerry-Anne Rye, professor1,
  5. Nicholas Hui, resident physician1,
  6. Xiong-Fei Pan, research professor256,
  7. Casey M Rebholz, associate professor7,
  8. Hyunju Kim, assistant scientist7,
  9. Lyn M Steffen, professor8,
  10. Anniek C van Westing, doctoral student9,
  11. Johanna M Geleijnse, professor9,
  12. Ellen K Hoogeveen, attending physician10111213,
  13. Yun-Yu Chen, research scientist,
  14. Kuo-Liong Chien, professor11,
  15. Amanda M Fretts, associate professor14,
  16. Rozenn N Lemaitre, research professor15,
  17. Fumiaki Imamura, senior investigator scientist16,
  18. Nita G Forouhi, professor16,
  19. Nicholas J Wareham, professor16,
  20. Anna Birukov, postdoctoral fellow1718,
  21. Susanne Jäger, research scientist1718,
  22. Olga Kuxhaus, research scientist1718,
  23. Matthias B Schulze, professor171819,
  24. Vanessa Derenji de Mello, university senior researcher20,
  25. Jaakko Tuomilehto, professor emeritus212223,
  26. Matti Uusitupa, professor emeritus20,
  27. Jaana Lindström, research manager21,
  28. Nathan Tintle, associate scientist2425,
  29. William S Harris, professor242627,
  30. Keisuke Yamasaki, attending physician27,
  31. Yoichiro Hirakawa, assistant professor28,
  32. Toshiharu Ninomiya, professor27,
  33. Toshiko Tanaka, staff scientist29,
  34. Luigi Ferrucci, scientific director29,
  35. Stefania Bandinelli, principal investigator30,
  36. Jyrki K Virtanen, associate professor20,
  37. Ari Voutilainen, data manager20,
  38. Tharusha Jayasena, postdoctoral fellow31,
  39. Anbupalam Thalamuthu, senior research fellow31,
  40. Anne Poljak, senior research scientist32,
  41. Sonia Bustamante, research scientist3231,
  42. Perminder S Sachdev, scientia professor,
  43. Mackenzie K Senn II, research coordinator33,
  44. Stephen S Rich, professor3435,
  45. Michael Y Tsai, professor36,
  46. Alexis C Wood, assistant professor33,
  47. Markku Laakso, professor3738,
  48. Maria Lankinen, postdoctoral fellow20,
  49. Xiaowei Yang, doctoral student39,
  50. Liang Sun, associate professor39,
  51. Huaixing Li, professor39,
  52. Xu Lin, professor3940,
  53. Christoph Nowak, assistant professor41,
  54. Johan Ärnlöv, professor41,
  55. Ulf Risérus, professor4,
  56. Lars Lind, professor42,
  57. Mélanie Le Goff, statistician43,
  58. Cécilia Samieri, research director43,
  59. Catherine Helmer, research director43,
  60. Frank Qian, resident physician4445,
  61. Renata Micha, associate professor4647,
  62. Adrienne Tin, associate professor748,
  63. Anna Köttgen, professor749,
  64. Ian H de Boer, professor505152,
  65. David S Siscovick, senior research scientist53,
  66. Dariush Mozaffarian, professor47,
  67. Jason HY Wu, professor254
  1. 1Lipid Research Group, School of Biomedical Sciences, University of New South Wales, Sydney, NSW, Australia
  2. 2The George Institute for Global Health, Faculty of Medicine and Health, University of New South Wales, Sydney, NSW, Australia
  3. 3The Welch Center for Prevention, Epidemiology, and Clinical Research, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA
  4. 4Department of Public Health and Caring Sciences, Clinical Nutrition and Metabolism, Uppsala University, Uppsala, Sweden
  5. 5Key Laboratory of Birth Defects and Related Diseases of Women and Children (Sichuan University), Ministry of Education, West China Second University Hospital, Sichuan University, Chengdu, Sichuan, China
  6. 6Shuangliu Institute of Women's and Children's Health, Shuangliu Maternal and Child Health Hospital, Chengdu, Sichuan, China
  7. 7Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA
  8. 8University of Minnesota School of Public Health, Minneapolis, MN, USA
  9. 9Division of Human Nutrition and Health, Wageningen University, Wageningen, The Netherlands
  10. 10Department of Nephrology, Jeroen Bosch Hospital, Den Bosch, The Netherlands
  11. 11Institute of Epidemiology and Preventive Medicine College of Public Health, National Taiwan University, Taipei, Taiwan
  12. 12Department of Medical Research, Taichung Veterans General Hospital, Taichung, Taiwan
  13. 13Heart Rhythm Center, Division of Cardiology, Department of Medicine, Taipei Veterans General Hospital, Taipei, Taiwan
  14. 14Department of Epidemiology, University of Washington, Seattle, WA, USA
  15. 15Department of Medicine, University of Washington, Seattle, WA, USA
  16. 16MRC Epidemiology Unit, University of Cambridge School of Clinical Medicine, Cambridge, UK
  17. 17Department of Molecular Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, Nuthetal, Germany
  18. 18German Center for Diabetes Research (DZD), München-Neuherberg, Germany
  19. 19Institute of Nutritional Science, University of Potsdam, Potsdam, Germany
  20. 20Institute of Public Health and Clinical Nutrition, School of Medicine, University of Eastern Finland, Kuopio, Finland
  21. 21Population Health Unit, Finnish Institute for Health and Welfare, Helsinki, Finland
  22. 22Department of Public Health, University of Helsinki, Helsinki, Finland
  23. 23Saudi Diabetes Research Group, King Abdulaziz University, Jeddah, Saudi Arabia
  24. 24The Fatty Acid Research Institute, Sioux Falls, SD, USA
  25. 25Department of Population Health Nursing Science, College of Nursing, University of Illinois—Chicago, Chicago, IL, USA
  26. 26Department of Internal Medicine, Sanford School of Medicine, University of South Dakota, Sioux Falls, SD, USA
  27. 27Department of Epidemiology and Public Health, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan
  28. 28Department of Medicine and Clinical Science, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan
  29. 29Translational Gerontology Branch, National Institute on Aging, National Institute of Health, Baltimore, MD, USA
  30. 30Geriatric Unit, Azienda USL Toscana Centro, Florence, Italy
  31. 31Centre for Healthy Brain Ageing, School of Psychiatry, University of New South Wales, Sydney, NSW, Australia
  32. 32Mark Wainwright Analytical Centre, University of New South Wales, Sydney, NSW, Australia
  33. 33USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine, Houston, TX, USA
  34. 34Department of Public Health Sciences, University of Virginia, Charlottesville, VA, USA
  35. 35Center for Public Health Genomics, University of Virginia, Charlottesville, VA, USA
  36. 36Department of Laboratory Medicine and Pathology, University of Minnesota, Minneapolis, MN, USA
  37. 37Institute of Clinical Medicine, Internal Medicine, University of Eastern Finland, Kuopio, Finland
  38. 38Department of Medicine, Kuopio University Hospital, Kuopio, Finland
  39. 39Shanghai Institute of Nutrition and Health, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai, China
  40. 40Key Laboratory of Systems Health Science of Zhejiang Province, Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Hangzhou, China.
  41. 41Department of Neurobiology, Care Sciences and Society, Karolinska Institutet, Sweden
  42. 42Department of Medical Sciences, Uppsala University, Uppsala, Sweden
  43. 43Bordeaux Population Health Research Centre, INSERM, UMR 1219, University of Bordeaux, Bordeaux, France
  44. 44Department of Nutrition, Harvard T.H. Chan School of Public Health, Harvard Medical School, Boston, MA, USA
  45. 45Department of Medicine, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, USA
  46. 46Department of Food Science and Nutrition, University of Thessaly, Karditsa, Greece
  47. 47The Friedman School of Nutrition Science and Policy, Tufts University, Boston, MA, USA
  48. 48Department of Medicine, University of Mississippi Medical Center, Jackson, MS, USA
  49. 49Institute of Genetic Epidemiology, Department of Data Driven Medicine, Faculty of Medicine and Medical Centre, University of Freiburg, Freiburg, Germany
  50. 50Division of Nephrology, Department of Medicine, University of Washington, Seattle, WA, USA
  51. 51Kidney Research Institute, University of Washington, Seattle, WA, USA
  52. 52Puget Sound VA Healthcare System, Seattle, WA, USA
  53. 53New York Academy of Medicine, New York, NY, USA
  54. 54School of Population Health, University of New South Wales, Sydney, NSW, Australia
  1. Correspondence to: K L Ong oklws{at}yahoo.com.hk
  • Accepted 12 December 2022

Abstract

Objective To assess the prospective associations of circulating levels of omega 3 polyunsaturated fatty acid (n-3 PUFA) biomarkers (including plant derived α linolenic acid and seafood derived eicosapentaenoic acid, docosapentaenoic acid, and docosahexaenoic acid) with incident chronic kidney disease (CKD).

Design Pooled analysis.

Data sources A consortium of 19 studies from 12 countries identified up to May 2020.

Study selection Prospective studies with measured n-3 PUFA biomarker data and incident CKD based on estimated glomerular filtration rate.

Data extraction and synthesis Each participating cohort conducted de novo analysis with prespecified and consistent exposures, outcomes, covariates, and models. The results were pooled across cohorts using inverse variance weighted meta-analysis.

Main outcome measures Primary outcome of incident CKD was defined as new onset estimated glomerular filtration rate <60 mL/min/1.73 m2. In a sensitivity analysis, incident CKD was defined as new onset estimated glomerular filtration rate <60 mL/min/1.73 m2 and <75% of baseline rate.

Results 25 570 participants were included in the primary outcome analysis and 4944 (19.3%) developed incident CKD during follow-up (weighted median 11.3 years). In multivariable adjusted models, higher levels of total seafood n-3 PUFAs were associated with a lower incident CKD risk (relative risk per interquintile range 0.92, 95% confidence interval 0.86 to 0.98; P=0.009, I2=9.9%). In categorical analyses, participants with total seafood n-3 PUFA level in the highest fifth had 13% lower risk of incident CKD compared with those in the lowest fifth (0.87, 0.80 to 0.96; P=0.005, I2=0.0%). Plant derived α linolenic acid levels were not associated with incident CKD (1.00, 0.94 to 1.06; P=0.94, I2=5.8%). Similar results were obtained in the sensitivity analysis. The association appeared consistent across subgroups by age (≥60 v <60 years), estimated glomerular filtration rate (60-89 v ≥90 mL/min/1.73 m2), hypertension, diabetes, and coronary heart disease at baseline.

Conclusions Higher seafood derived n-3 PUFA levels were associated with lower risk of incident CKD, although this association was not found for plant derived n-3 PUFAs. These results support a favourable role for seafood derived n-3 PUFAs in preventing CKD.

Introduction

Chronic kidney disease (CKD) affects about 700 million people worldwide, with an estimated global prevalence of one in 11 in the general population.12 Patients with CKD are at higher risk of cardiovascular disease and death34 because the condition could eventually progress to kidney failure that severely impacts health and quality of life.56 Therefore, there is a need to identify factors that might prevent the onset and progression of CKD.

The consumption of long chain omega 3 polyunsaturated fatty acids (n-3 PUFAs) that are mostly obtained from seafood (including eicosapentaenoic acid (EPA), docosapentaenoic acid (DPA), and docosahexaenoic acid (DHA)) confers cardiometabolic benefits.7 Meta-analyses of randomised controlled trials have shown that increased n-3 PUFA intake improved arterial stiffness,8 lowered blood pressure,9 and reduced plasma triglycerides.1011 Because endothelial dysfunction, hypertension, and dyslipidaemia are CKD risk factors,121314 n-3 PUFAs could protect against the development of CKD.

However, there are important limitations to the observational research addressing this important theory. Firstly, most previous studies assessed self-reported intake of n-3 PUFAs using dietary questionnaires and evaluated associations with CKD in cross sectional analyses, with inconsistent associations with prevalent CKD.15161718 Only two longitudinal studies assessed self-reported intake of n-3 PUFAs and risk of incident CKD, and also reported inconsistent findings.1819 Secondly, the use of self-reported n-3 PUFA intake is subject to errors related to misreporting and inaccuracy of food composition databases.20 Circulating n-3 PUFA levels are valid and reliable biomarkers of intake as shown by several randomised controlled trials.2021222324 Moreover, biomarkers reflect the dietary intake and also the underlying metabolism, and therefore indicate the bioavailable n-3 PUFA intake; this might explain why n-3 PUFA biomarkers have been reported to be superior to self-reported n-3 PUFA intake in predicting outcomes such as cardiovascular disease and all cause mortality.25 Thirdly, self-reported intake could not distinguish accurately between individual n-3 PUFAs, which might differ in their biological properties and effects on CKD.26 Finally, studies evaluating the association between α linolenic acid (ALA) biomarker level and incident CKD are lacking, even though ALA is the major dietary n-3 PUFA obtained from specific plant foods and oils that has been associated with lower cardiovascular outcome risk.2728

We pooled cohort studies conducted by different research teams to assess the association of n-3 PUFA biomarkers with incident CKD. This strategy allowed us to more precisely quantify the association with a larger sample size, distinguish between individual n-3 PUFAs, and avoid measurement errors related to self-reported dietary intake.

Methods

Study population

The study was conducted within the Fatty Acids and Outcomes Research Consortium (FORCE) (http://force.nutrition.tufts.edu), a consortium of studies with circulating or adipose tissue fatty acid biomarker measurements and ascertained chronic disease events.2829303132 Studies were identified and study investigators were asked to participate based on previous FORCE projects, expert contacts, and literature searches. All prospective studies of adult participants with n-3 PUFA biomarkers measured in blood fractions or adipose tissue and with data on the primary outcome of incident CKD were eligible. A total of 50 potentially relevant studies were identified in May 2020, but 22 (44%) were not eligible because of lack of exposure or outcome data, and 5 (10%) studies were not included because of a limited number of patients with incident CKD (supplementary table S1). Of the remaining 23 cohorts, 19 (83%) participated in this project, while 4 (17%) declined. The study proposal and analytical plan were approved by FORCE, and all participating cohorts have ethics approval with appropriate data sharing agreements, which enabled the current project to be performed. Detailed cohort information is provided in the supplementary methods.

Measurement of n-3 PUFA biomarkers

Biomarker levels of n-3 PUFAs (ALA, EPA, DPA, and DHA) were measured by each cohort in plasma phospholipids, erythrocyte phospholipids, total plasma, total serum, cholesterol esters, or several lipid compartments using gas chromatography (supplementary methods). Individual fatty acid levels were expressed as the proportion of total fatty acids.

Outcome measurement

All outcomes were determined using estimated glomerular filtration rate (eGFR) at baseline and follow-up visits, with the same definitions in all the 19 participating cohorts. The primary outcome was incident CKD using the same definition as used by the Chronic Kidney Disease Genetics (CKDGen) Consortium.33 Briefly, incident CKD was defined as an eGFR that decreased to <60 mL/min/1.73 m2 during follow-up among participants with baseline eGFR≥60 mL/min/1.73 m2. eGFR was estimated using the 2009 Chronic Kidney Disease Epidemiology Collaboration equation with calibrated serum or plasma creatinine.34 If creatinine levels were not calibrated and were obtained with a Jaffé assay before 2009, then data were corrected by multiplying creatinine levels by 0.95 before eGFR calculation, as described previously.35 A more stringent definition of incident CKD was used in a sensitivity analysis: both eGFR<60 mL/min/1.73 m2 and <75% of baseline eGFR during follow-up among participants with baseline eGFR≥60 mL/min/1.73 m2.36

We also assessed two secondary outcomes: annual absolute change in eGFR from baseline to the last follow-up as a continuous variable; and ≥40% decrease in eGFR from baseline to the last follow-up.37 Unlike the primary outcomes, analysis of the secondary outcomes comprised all participants, including those with baseline eGFR<60 mL/min/1.73 m2.

Cohort analyses

A prespecified and harmonised protocol with the same participant inclusion and exclusion criteria, exposures, outcomes, covariates, and analytical methods were sent to all cohorts to perform de novo analyses of individual data. For the outcomes of incident CKD and ≥40% decrease in eGFR, Cox proportional hazards models were used to estimate hazard ratios, considered to approximate relative risk, for cohorts with eGFR available in two or more follow-up visits. In cohorts with data available for only one follow-up eGFR, Poisson regression models were used to estimate relative risks. For the outcome of annual absolute change in eGFR from baseline to the last follow-up, multivariable linear regression was used. All analyses estimated robust standard errors.

Because n-3 PUFA biomarker levels were measured in different lipid compartments in the participating cohorts, they were evaluated as continuous variables using study specific interquintile range (difference between midpoint of lowest fifth (10th centile) and highest fifth (90th centile)).29303132 To assess potential nonlinear associations, n-3 PUFA levels were also evaluated as variables (in fifths) using study specific cut-off points. In the primary multivariable model (model 1), analyses were adjusted for age, sex, race, clinical centre or field site, education, occupation, body mass index, smoking, alcohol intake, physical activity, prevalent coronary heart disease, and use of lipid lowering drugs at baseline, when applicable. In model 2, additional covariates that could confound or mediate the association of n-3 PUFAs with incident CKD were considered, including prevalent diabetes mellitus, urine albumin-creatinine ratio, systolic blood pressure, and use of antihypertensive drugs at baseline (see supplementary methods for definitions of all covariates). Participants with missing data on n-3 PUFAs were excluded from the analysis. Missing data on covariates were handled according to the usual practice of each cohort and study investigators (see supplementary methods for cohort specific details).

In each study, a statistical interaction was evaluated by fitting the n-3 PUFA biomarker (continuous), each potential effect modifier, and their cross product. The prespecified effect modifiers included age (≥60 v <60 years), race, hypertension, diabetes, coronary heart disease, and eGFR (60-89 v ≥90 mL/min/1.73 m2) at baseline. These modifiers were selected by considering their potential biological relevance.

Data pooling and meta-analysis

Effect estimates from each participating cohort were provided to the lead author in standardised electronic forms and pooled using inverse variance weighted meta-analysis. Results from participating cohorts were pooled overall and within each specific type of lipid compartment. For studies with measures in multiple compartments, we prioritised the compartments in the overall pooled analysis on the basis of which best reflect long term intake, as prespecified in the following order: erythrocyte phospholipids, plasma phospholipids, cholesterol esters, and total plasma or serum.22 In sensitivity analyses, the overall pooled analyses were performed after excluding data from one cohort at a time or using alternative biomarker compartments. To assess potential nonlinear relations, exploratory analysis was conducted to model the association within each compartment using meta-regression with restricted cubic splines constructed from variables (in fifths) using study specific cut-off points.282930313238 Interaction terms were pooled using inverse variance weighted meta-regression with Bonferroni correction for multiple testing, and such analysis was considered exploratory. All meta-analyses were performed using Stata (version 16.0) or R (version 4.1.0), with two tailed P<0.05 considered statistically significant.

Patient and public involvement

Patients and members of the public were not directly involved in this research owing to lack of funding, staff, and infrastructure to facilitate their involvement in this pooled analysis of data from 19 participating cohorts.

Results

Study cohorts

The primary outcome analysis included 25 570 participants without prevalent CKD from 19 cohorts in a total of 12 countries (table 1). Among the 19 cohorts, the mean age of participants ranged from 49 to 77 years and the mean body mass index from 23.2 to 28.3. The weighted median follow-up duration was 11.3 years. Sixteen cohorts recruited men and women, and most (n=15) recruited predominantly white participants (supplementary table S2). The mean baseline eGFR ranged from 76.1 to 99.8 mL/min/1.73 m2 (supplementary table S3), and the mean annual change in eGFR ranged from −6.2 to 0.3 mL/min/1.73 m2/year. In total, 4944 participants (19.3%) developed the primary outcome of incident CKD during follow-up (supplementary table S4).

Table 1

Baseline characteristics of 19 prospective cohorts that participated in analysis of association of n-3 PUFA biomarkers with incident chronic kidney disease (as part of FORCE Consortium)

View this table:

Most of the participating cohorts measured fatty acid levels in erythrocyte or plasma phospholipids (n=11), followed by total plasma or serum (n=7), and cholesterol esters (n=3; table 1). The proportion of total fatty acids for each n-3 PUFA varied between different lipid compartments (supplementary table S5; supplementary fig S1).

Seafood n-3 PUFAs and incident CKD

Higher levels of total seafood derived n-3 PUFAs were associated with 8% lower risk of incident CKD per interquintile range with low heterogeneity in the observed associations between cohorts (relative risk 0.92, 95% confidence interval 0.86 to 0.98; P=0.009, I2=9.9%; table 2 and fig 1). Findings were not appreciably altered by further adjustment for covariates that could be potential confounders or mediators (0.91, 0.85 to 0.97; P=0.006, I2=17.2%). Similar results were obtained when the n-3 PUFAs were assessed categorically—participants with total seafood n-3 PUFA level in the highest fifth had 13% lower risk of incident CKD than those in the lowest fifth (0.87, 0.80 to 0.96; P=0.005, I2=0.0%; table 2 and supplementary fig S2). Protective associations were found for the individual seafood n-3 PUFAs (table 2), and associations appeared consistent across lipid compartments (fig 1; supplementary table S6). Little evidence could be found of nonlinear associations between total seafood n-3 PUFAs and incident CKD in any of the biomarker compartments (supplementary fig S3).

Table 2

Association of seafood n-3 PUFA biomarkers with primary outcome of incident CKD

View this table:
Fig 1
Fig 1

Association of total seafood n-3 PUFAs (eicosapentaenoic acid, docosapentaenoic acid, and docosahexaenoic acid) with incident CKD. Incident CKD was defined as an eGFR<60 mL/min/1.73 m2 during follow-up among participants with baseline eGFR≥60 mL/min/1.73 m2. Analyses were adjusted for age, sex, race, field centre if applicable, education, occupation, body mass index, smoking, physical activity, alcohol intake, prevalent coronary heart disease, and use of lipid lowering drugs. Study specific estimates per interquintile range (difference between midpoint of lowest fifth and highest fifth) of total n-3 PUFAs were pooled separately for different lipid compartments and overall. AOC=Alpha Omega Cohort; ARIC=Atherosclerosis Risk in Communities Study; CCCC=Chin-Shan Community Cardiovascular Cohort Study; CHS=Cardiovascular Health Study; CKD=chronic kidney disease; eGFR=estimated glomerular filtration rate; EPIC-Norfolk=European Prospective Investigation into Cancer and Nutrition (Norfolk); EPIC-Potsdam=European Prospective Investigation into Cancer and Nutrition (Potsdam); FDPS=Finnish Diabetes Prevention Study; FHS=Framingham Heart Study; InCHIANTI=Invecchiare in Chianti Study; KIHD=Kuopio Ischaemic Heart Disease Risk Factor Study; MAS=Memory and Ageing Study; MESA=Multi-Ethnic Study of Atherosclerosis; METSIM=Metabolic Syndrome in Men; n-3 PUFA=omega 3 polyunsaturated fatty acid; NHAPC=Nutrition and Health of Ageing Populations in China; PIVUS=Prospective Investigation of the Vasculature in Uppsala Seniors; WHIMS=Women’s Health Initiative Memory Study; ULSAM=Uppsala Longitudinal Study of Adult Men; 3C=Three City Study

When incident CKD was defined using a more stringent definition (new onset of eGFR<60 mL/min/1.73 m2 and <75% of baseline eGFR) in a sensitivity analysis, the results were largely consistent with the primary outcome analysis. Higher levels of total seafood n-3 PUFAs were associated with a lower risk of incident CKD (0.91, 0.84 to 0.99 per interquintile range; P=0.009, I2=27.3%; supplementary table S7). In additional sensitivity analyses, the association of total seafood n-3 PUFAs with incident CKD remained significant after excluding data from one cohort at a time, or using alternative biomarker compartments in the overall pooled results (supplementary table S8).

Seafood n-3 PUFAs and secondary outcomes

Total seafood n-3 PUFA level was not associated with the outcome of ≥40% decrease in eGFR in continuous analyses, but was associated with a lower risk (about 15%) for those in the highest fifth versus those in the lowest fifth (relative risk 0.85, 95% confidence interval 0.74 to 0.98; P=0.03, I2=44.1%; table 3).

Table 3

Association of seafood n-3 PUFA biomarkers with secondary outcome of ≥40% decrease in eGFR from baseline

View this table:

Higher levels of total seafood n-3 PUFAs, especially DHA, were associated with a slower annual decline in eGFR (table 4). For instance, the annual decline in eGFR was 0.07 mL/min/1.73 m2 lower (95% confidence interval 0.02 to 0.13; P=0.007, I2=42.2%) for people with total seafood n-3 PUFA level in the highest fifth than those in the lowest fifth.

Table 4

Association of seafood n-3 PUFA biomarkers with secondary outcome of annual change in eGFR

View this table:

ALA and incident CKD

ALA was not associated with incident CKD, with low heterogeneity in the observed associations between cohorts (relative risk 1.00 per interquintile range, 95% confidence interval 0.94 to 1.06; P=0.94, I2=5.8%; table 2 and supplementary fig S4). No significant association was observed in analyses using a more stringent definition of incident CKD (supplementary table S7). No significant association was observed in the secondary outcomes of ≥40% decrease in eGFR (0.98, 0.89 to 1.07; P=0.08, I2=35.8%; table 3) and annual change in eGFR (mean difference −0.03 mL/min/1.73 m2 per interquintile range, 95% confidence interval −0.06 to 0.00; P=0.16, I2=24.1%; table 4).

Effect modification analyses

The association of n-3 PUFAs with the primary outcome of incident CKD was not modified by key demographic and clinical characteristics at baseline, including age, race, hypertension, diabetes, coronary heart disease, and eGFR (P for pooled interaction terms >0.05). No difference in association was observed based on study location, follow-up duration, and sample size (supplementary tables S9-11).

Discussion

Principal findings

This de novo pooled analysis of 19 cohort studies gathered data from more than 25 000 patients as part of the FORCE consortium. We found that higher levels of seafood n-3 PUFA biomarkers were associated with a modestly lower risk of incident CKD and slower decline in renal function, whereas these associations were not found with higher levels of plant n-3 PUFAs (ALA). The results were consistent across a range of sensitivity analyses, highlighting the robustness of the findings. Similar trends were observed using secondary outcomes of >40% decrease in eGFR and annual changes in eGFR. These findings support a favourable role for increased intake of seafood n-3 PUFAs for the primary prevention of CKD.

Comparison with previous studies

Randomised controlled trials found increased intake of seafood n-3 PUFAs reduced blood pressure, a key risk factor for CKD development.9 Conversely, ALA generally exhibits weaker to no effect on metabolic risk factors such as lipid, glucose, and inflammatory profile in randomised controlled trials compared with EPA and DHA.39 Furthermore, findings from previous experimental animal studies provide strong biological plausibility similar to our results of the beneficial effects of seafood n-3 PUFAs on CKD.404142 For example, in a mouse model of nephrectomy induced chronic renal failure, supplementation with EPA and DHA for 12 weeks mitigated tubulointerstitial injury through reductions in fibrosis, inflammation, and oxidative stress.40 Increased endogenous production of DHA by a transgenic technique also protected against the development of kidney fibrosis and inflammation in mice with unilateral ureter obstruction induced nephropathy.41 Renoprotective mechanisms might be mediated through the endogenous conversion of n-3 PUFAs to specialised inflammation resolving mediators such as resolvins.4344 Most previous experiments focused on EPA and DHA, but less is known about DPA, although limited evidence also suggests DPA has favourable properties on cardiometabolic profile.4546

Although our findings do not prove a causal relation between seafood n-3 PUFAs and CKD risk, they are supportive and consistent with current clinical guidelines that recommend adequate intake of seafood as part of healthy dietary patterns, especially when seafood replaces the intake of less healthy foods.474849 Although ALA could be enzymatically transformed to the longer chain n-3 PUFAs like EPA and DHA, such endogenous conversion only occurs at a low rate27 and our findings suggest intake of this plant derived n-3 PUFA alone might not maintain renal health. However, although ALA is an essential fatty acid, circulating ALA generally has a weaker correlation with ALA intake compared with the seafood derived n-3 PUFAs,225051 potentially explained by a higher oxidation rate of ALA compared with other fatty acids.52 Therefore, if ALA intake were associated with incident CKD, we might not be able to detect such an association using an ALA biomarker. In our study, the levels of EPA and DHA probably mostly reflected variations in dietary intake from seafood, rather than n-3 PUFA supplements, because relatively small proportions of the general population used n-3 PUFA supplements, especially for studies that were conducted before the early 2000s when fish oil supplements were infrequently used (most of the included cohorts).53 Current dietary guidelines generally recommend at least two servings per week of oily fish to provide around 250 mg/day of long chain n-3 PUFAs for the general population.2749

Given the renoprotective effects of seafood n-3 PUFAs observed in experimental animal studies, human trials have also assessed the use of n-3 PUFA supplementation to prevent decline in kidney function among patients with existing CKD or kidney failure. A recent systematic review of randomised controlled trials found n-3 PUFA supplementation could improve lipid profile and reduce oxidative stress, but not blood pressure in patients with CKD.54 Another recent systematic review of randomised controlled trials found n-3 PUFA supplementation could prevent the progression to kidney failure in patients with CKD who were not receiving renal replacement therapy, although the certainty of evidence is low.55 Therefore, current clinical guidelines do not recommend the use of n-3 PUFA supplementation to prevent further decline in renal function for those with existing CKD, although using n-3 PUFA supplementation to treat hypertriglyceridaemia is recommended.56 Our findings highlight the need for large randomised controlled trials to assess increased intake of seafood n-3 PUFAs for the primary prevention of CKD and decline in renal function. One such randomised controlled trial of 1312 patients with type 2 diabetes found that fish oil supplementation (daily dose of 1 g, containing 465 mg EPA and 375 mg DHA) did not affect the change in eGFR over five years.57 Several other randomised controlled trials have also reported null or small beneficial effects of n-3 PUFA supplementation on kidney function, but these were secondary outcomes or post hoc analyses and were not powered to detect incident CKD as a primary outcome.5859 Further, randomised controlled trials of n-3 PUFA supplementation often included participants with sufficient seafood intake, in which additional n-3 PUFA supplementation might not confer additional cardiometabolic benefits.60 Therefore, future randomised controlled trials could target increased seafood intake, instead of n-3 PUFA supplementation, because consumption of seafood can replace the intake of less healthy foods. Finally, all these trials tested n-3 PUFA supplementation mainly in high risk populations, limiting the applicability of the study findings to the general population.

Strengths and limitations of this study

Our study has several strengths. Our literature search only identified two high quality prospective observational analyses of the association of n-3 PUFA biomarkers with incident CKD or decline in kidney. These two studies were included in our pooled analysis with updated patient numbers and harmonised analytical protocols.6162 Our findings expand upon these previous individual cohort based analyses with a sevenfold larger sample size and enhanced statistical power to more precisely quantify the association of n-3 PUFAs with CKD. The large sample size also enabled a detailed evaluation of potential effect modifiers of the association between n-3 PUFAs and CKD. By collaborating with participating cohorts using a standardised approach to define exposures, covariates, and outcome variables, we reduced potential heterogeneity that is inherent in a publication based meta-analysis. Our collaborative approach enabled the inclusion of most of the cohorts with available exposure and outcome data, which reduced the likelihood of publication bias.

The use of n-3 PUFA biomarkers allowed us to distinguish between individual n-3 PUFAs that could have differing biological effects,26 and avoid measurement errors related to self-reported dietary intake. Conducting analyses across cohorts from diverse demographic, dietary, and medical backgrounds with different incident CKD rates enhanced the generalisability of our findings, and the use of several related secondary outcomes and different sensitivity analyses enhanced the robustness of the findings. However, more studies are needed to further understand the generalisability of our findings, especially in countries experiencing high incidence of CKD such as China and India.

Our study also has some limitations. The n-3 PUFA biomarkers and covariates were only measured once at baseline and this might lead to increasing random misclassification over time, which would tend to bias the associations towards the null. However, n-3 PUFA biomarkers have good long term reproducibility as determined in a previous study with repeated within person measurements across 13 years.63 While we standardised our methods across cohorts as much as possible, fatty acid assays, determination of outcomes, and covariate measurements varied by cohort and might contribute to the heterogeneity in our findings. Among different prespecified covariates, many participating cohorts did not have data for urine albumin-creatinine ratio. However, in cohorts that measured urine albumin-creatinine ratio, association of n-3 PUFA levels with incident CKD was similar whether it was adjusted for or not (results not shown). Although, the possibility of residual confounding by other unmeasured or imprecisely measured covariates could not be excluded, the differential association of plant derived ALA and seafood n-3 PUFAs with incident CKD suggests residual confounding was unlikely to entirely account for our findings. In this analysis, data were not adjusted for n-6 PUFA biomarker levels because experimental evidence on their relation with kidney function was limited. Because n-6 PUFA levels are associated with a lower CVD risk,30 future studies assessing their association with incident CKD are warranted, pending more experimental work assessing their possible impact on CKD.

Conclusion

In this large de novo pooled analysis across 19 cohorts with more than 25 000 patients, higher seafood n-3 PUFA levels were associated with a lower incident CKD risk and a slower decline in renal function. Although the magnitude of these associations was modest, our findings suggest adequate consumption of seafood and oily fish should be part of healthy dietary patterns. Additionally, further randomised controlled trials are warranted to assess the potential beneficial role of seafood n-3 PUFAs in preventing and managing CKD.

What is already known on this topic

  • Animal studies suggest omega 3 polyunsaturated fatty acids (n-3 PUFAs) have beneficial effects on kidney function, but evidence from human studies is limited

  • Most previous studies on chronic kidney disease (CKD) assessed self-reported intake of n-3 PUFAs using dietary questionnaires, which are subject to errors related to misreporting and inaccuracy of food composition databases

What this study adds

  • Higher seafood n-3 PUFA levels (eicosapentaenoic acid, docosapentaenoic acid, and docosahexaenoic acid) were associated with a lower incident CKD risk and a slower decline in renal function

  • Plant derived n-3 PUFA level (α linolenic acid) was not associated with incident CKD

  • Findings are consistent with dietary guidelines recommending seafood and oily fish consumption as part of healthy dietary patterns, and provide strong evidence for further trials analysing seafood derived n-3 PUFAs in CKD prevention

Ethics statements

Ethical approval

All the 19 participating cohorts have ethics approval from their respective institutions with appropriate data sharing agreements that enabled the current project of pooled data analysis (Alpha Omega Cohort—AOC: Medical Ethics Committee South-West Holland at the Haga Hospital and local ethics committees of all participating hospitals (approval No L01.049); Atherosclerosis Risk in Communities Study—ARIC: Johns Hopkins Bloomberg School of Public Health Institutional Review Board (approval No 11012/CR1014); CCCC—Chin-Shan Community Cardiovascular Cohort Study: Institutional Review Board of the National Taiwan University Hospital (approval No 201003001R); CHS—Cardiovascular Health Study: University of Washington Human Subject Committee (approval No STUDY00000109); EPIC-Norfolk—European Prospective Investigation into Cancer and Nutrition (Norfolk): Norwich Local Ethics Committee (approval No 98CN01); EPIC-Potsdam—European Prospective Investigation into Cancer and Nutrition (Potsdam): Ethics Committee of the State of Brandenburg, Germany (approval No AS 29/93 and S 9/2002); FDPS—Finnish Diabetes Prevention Study: the Ethics Committee of the National Public Health Institute of Helsinki (dated 6 March 1992, approval No n/a); FHS—Framingham Heart Study, and WHIMS—Women’s Health Initiative Memory Study: University of South Dakota Institutional Review Board (approval No IRB-21-136); Hisayama: The Kyushu University Institutional Review Board for Clinical Research (approval No 2021-457); InCHIANTI—Invecchiare in Chianti Study: Internal Review Board of the Intramural Research Program of the National Institutes of Health (approval No AG001050; exemption 11796); KIHD—Kuopio Ischaemic Heart Disease Risk Factor Study: Ethics Working Committee of Kuopio University (approval No 01/12/1983), and Ethics Committee of Kuopio University and Kuopio University Hospital (approval No 143/97); MAS—Memory and Ageing Study: Human Research Ethics Committee of the University of New South Wales (approval no. HC200506); MESA—Multi-Ethnic Study of Atherosclerosis: Intuitional Review Board Human Subjects Committee at The University of Minnesota (approval No 9805M00034), Johns Hopkins Medicine Intuitional Review Board (approval No 99-11-10-06), University of California Los Angeles Intuitional Review Board (approval No 99-057-22A), North western University Intuitional Review Board (approval No 200806-3515), Wake Forest University Intuitional Review Board (approval No BG00-035), Columbia University Medical Center Intuitional Review Board (approval No AAAA7791), and Baylor College of Medicine Intuitional Review Board (approved protocol No H-49021); METSIM—Metabolic Syndrome in Men: Ethics Committee of the Kuopio University Hospital (approval No 174/2004); NHAPC—Nutrition and Health of Ageing Populations in China: Institutional Review Board of the Institute for Nutritional Sciences, Chinese Academy of Sciences (approval No E-2005-01 for baseline survey and E-2009-01 for follow-up survey); PIVUS—Prospective Investigation of the Vasculature in Uppsala Seniors, and ULSAM—Uppsala Longitudinal Study of Adult Men: Ethics Committee of Uppsala University (approval No 2011/045, 00-419, and 251/90); 3C—Three City Study: Consultative Committee for the Protection of Persons participating in Biomedical Research at Kremlin-Bicêtre University Hospital (approval No 99-28).

Data availability statement

Individual participant data are owned by individual participating cohorts and are available to researchers on consent from participating cohorts. For further queries or requests, please contact force@tufts.edu. Further details are available at the FORCE website: http://force.nutrition.tufts.edu.

Footnotes

  • Contributors: KLO and JHYW designed the study and drafted the manuscript. KLO, MM, KAR, NH, XFP, RM, AT, AK, IHdB, DSS, DM, and JHY contributed to the development of study proposal and protocol. MM, CMR, HK, LMS, ACvW, JMG, EKH, YYC, KLC, AMF, RNL, FI, NGF, NJW, AB, SJ, OK, MBS, VDdM, JT, MU, JL, NT, WSH, KY, YH, TN, TT, LF, SBa, JKV, AV, TJ, AT, AP, SBu, PSS, MKS, SSR, MYT, ACW, MLaa, MLan, XY, LS, HL, XL, CN, JA, UR, LL, MLG, CS, and CH contributed to cohort specific data (principal investigators, analysts, associates, and trainees). KLO and LH performed the pooled analysis of cohort data. KLO, MM, KAR, XFP, WSH, FQ, AT, AK, IhdB, DSS, DM, and JHYW contributed to the primary writing team. All authors contributed to the data interpretation, critically revised the manuscript and approved the final version for publication. KLO is the guarantor of this manuscript and accepts full responsibility for the work and the conduct of the study, had access to the data, and controlled the decision to publish. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

  • Funding: Cohort specific funding is outlined in the supplementary methods. KLO was supported by the Australian National Health and Medical Research Council Career Development Fellowship (1122854) and the University of New South Wales Safety Net Fellowship. The funders had no role in the study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the article for publication in the analysis and interpretation of data, in the writing of the report, and in the decision to submit the article for publication. We operated independently from the funders.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at https://www.icmje.org/disclosure-of-interest and declare: support from Australian National Health and Medical Research Council Career Development Fellowship and the University of New South Wales Safety Net Fellowship for the submitted work; MM reports research funding from Resolve to Save Lives, World Health Organisation and North western University, and support as invited speaker in the Nordic Dairy Congress 2022. CMR and LMS report research funding from the US National Institute of Health (NIH). CMR reports participation on the Data Safety Monitoring Boards of the SUPER and ADEPT trials and leadership role as associate editor of Diabetes Care. ACvW reports research funding from Jaap Schouten Foundation. JMG reports research funding from Jaap Schouten Foundation, EU Horizon 2020, and Ministry of Health, Welfare and Sports, The Netherlands, and leadership role as the Vice President of the Dutch Health Council. EKH reports research funding from Dutch Kidney Foundation. FI and NGF receive the MRC Epidemiology Unit core support. NGF reports research fundings from the NIHR Cambridge Biomedical Research Centre Theme on Nutrition, Diet and Lifestyle. JT declares conference support from the University of Antioquia, Colombia and possession of stocks from Orion Pharma. MU declares possession of stocks from Orion Pharma. WSH declares possession of stock in OmegaQuant Analytics, LLC (a laboratory that offers blood fatty acid testing to healthcare providers, researchers and consumers). PSS reports research funding from the National Health and Medical Research Council of Australia and honoraria from Biogen Australia and Roche Australia. ACW reports research funding from the US Department of Agriculture/Agricultural Research Service, NIH, National Cattlemen’s Beef Association and Hass Avocado Board Avocado Nutrition Research Center. JA has received research honoraria for lectures from AstraZeneca and Novartis, and has participated in the advisory board for AstraZeneca and Boerhinger Ingelheim, unrelated to the present study. UR reports research funding from Swedish Research Council Forma and Swedish Diabetes Foundation. RM reports research funding from the US NIH, Gates Foundation, Nestle and Danone, and consulting fees from Development Initiatives with leadership role as chair of the Independent Expert Group, Global Nutrition Report. AT declares as co-chair of ClinGen Gout Genetic Curation Panel and reports research funding from the US National Institute of Health; AK reports research funding from the German Research Foundation. DM reports research funding from the US NIH, the Gates Foundation, The Rockefeller Foundation, Vail Innovative Global Research, and the Kaiser Permanente Fund at East Bay Community Foundation; personal fees from Acasti Pharma and Barilla; scientific advisory board, Beren Therapeutics, Brightseed, Calibrate, Elysium Health, Filtricine, HumanCo, Instacart, January Inc., Perfect Day, Tiny Organics, and (ended) Day Two, Discern Dx, and Season Health; stock ownership in Calibrate and HumanCo; and chapter royalties from UpToDate. No other relationships or activities that could appear to have influenced the submitted work were reported.

  • The lead author (KLO) affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as originally planned in the protocol have been explained.

  • Dissemination to participants and related patient and public communities: We plan to share the results of this study through press releases, presentations, institutional newsletters, and multiple social media platforms, including Twitter and LinkedIn.

  • Provenance and peer review: Not commissioned; externally peer reviewed.

http://creativecommons.org/licenses/by/4.0/

This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/.

References