- 1Department of Endocrinology, The First Affiliated Hospital of Yangtze University, Jingzhou, Hubei, China
- 2Health Science Center, Yangtze University, Jingzhou, Hubei, China
- 3Southern Theater General Hospital, Guangzhou, Guangdong, China
- 4Department of Nursing, Renmin Hospital of Wuhan University, Wuhan, Hubei, China
- 5Department of Neurology, The Second Affiliated Hospital of Zhengzhou University, Zhengzhou, China
Background and aims: Previous research has underscored the association between oily fish intake and type 2 diabetes (T2DM), yet the causality remains elusive.
Methods: A bidirectional univariable Mendelian Randomization (MR) analysis was employed to evaluate the causal effects of oily fish and non-oily fish intake on T2DM. Replication analysis and meta-analysis were conducted to ensure robust results. Multivariable MR analysis was utilized to assess confounders, and further mediation MR analysis discerned mediating effects. Linkage Disequilibrium Score (LDSC) analysis was undertaken to compute genetic correlations. Inverse variance weighted (IVW) was the primary method, complemented by a series of sensitivity analyses.
Results: The LDSC analysis unveiled a significant genetic correlation between oily fish intake and T2DM (Genetic correlation: -0.102, p = 4.43 × 10−4). For each standard deviation (SD) increase in genetically predicted oily fish intake, the risk of T2DM was reduced by 38.6% (OR = 0.614, 95% CI 0.504 ~ 0.748, p = 1.24 × 10−6, False Discovery Rate (FDR) = 3.72 × 10−6). The meta-analysis across three data sources highlighted a persistent causal association (OR = 0.728, 95% CI 0.593 ~ 0.895, p = 0.003). No other causal effects were identified (all p > 0.5, FDR > 0.5). The main outcomes remained consistent in most sensitivity analyses. Both MVMR and mediation MR analyses emphasized the mediating roles of triglycerides (TG), body mass index (BMI), and 25-hydroxyvitamin D (25OHD) levels.
Conclusion: To encapsulate, there’s an inverse association between oily fish intake and T2DM risk, suggesting potential benefits of oily fish intake in T2DM prevention.
Introduction
Diabetes Mellitus (DM), a leading metabolic disorder marked by chronic hyperglycemia, is on the rise globally (1). Type 2 diabetes mellitus (T2DM) accounts for approximately 90% of all diabetes cases worldwide (2). Current estimates from the World Health Organization (WHO) indicate that over 422 million individuals suffer from diabetes globally, with projections reaching up to 629 million by 2045 (3, 4). The prevalence of diabetes, especially in developing nations, has been on an upward trajectory. This surge poses a considerable direct and indirect financial burden on society (5). Consequently, identifying novel modifiable risk factors for T2DM is crucial to inform clinical management strategies and mitigate the disease’s onset and progression.
Dietary habits emerge as a pivotal modifiable risk determinant, with burgeoning evidence suggesting that functional foods harbor bioactive components correlated with physiological health benefits pertinent to the prevention and management of T2DM (6). Notably, the consumption of oily fish rich in omega-3 fatty acids has consistently piqued interest. Omega-3 fatty acids, particularly eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) are renowned for their anti-inflammatory and cardioprotective attributes (7). Their potential roles in glucose metabolism, anti-inflammatory actions, lipid profile amelioration, and insulin sensitivity could shed light on mechanisms underpinning T2DM prevention (7). A population-based prospective study has indicated an association between oily fish consumption and reduced risk of T2DM (8). Yet, a cross-sectional study by Rashmi D Sahay and colleagues posits a direct correlation between oily fish intake and diabetes prevalence in the general populace (9). Furthermore, a meta-analysis amalgamating 13 cohort studies has furnished evidence of a tepid positive relationship between fish and n-3 fatty acid intake and T2DM risk (10). The discordance in research outcomes, coupled with intrinsic limitations of observational studies such as residual confounding and potential reverse causality, obscures the definitive causal relationship between oily fish intake and T2DM pathophysiology.
Mendelian Randomization (MR) operates as an analytical tool, harnessing genetic variants as instrumental variables (IVs) extracted from genome-wide association studies (GWAS) to examine the causal connection between a risk factor (or exposure) and an associated outcome (11). When set against classic observational methods, MR provides an enhanced perspective of the causal impact of the exposure on the outcome by lessening the impact of confounding elements (12). LD Score Regression (LDSC) functions as a strategy to gauge a trait’s heritability, essentially determining the extent of variance within a trait due to genetic influences (13). Additionally, it evaluates the genetic overlap among varied traits based on summary statistics from GWAS (14). Against this backdrop, our study is designed to explore the genetic tendencies linked to oily fish intake concerning T2DM. Moreover, we delve deeper into the potential roles of nine confounders through Multivariable MR (MVMR) and mediation MR frameworks.
Materials and methods
Study design
The foundational data for our study was retrieved from open-access, summary-tier datasets stemming from GWAS. We deployed a suite of techniques including Univariate MR (UVMR), MVMR, genetic correlation, mediation MR, and meta-analysis to delve deeply into the causal nexus between the intake of oily and non-oily fish and T2DM. We adopted three pivotal principles in selecting IVs for the exposure: (i) the pinpointed genetic variant, acting as the instrumental variable, demonstrates a robust affiliation with the exposure; (ii) this genetic variant maintains no ties with probable confounders; and (iii) the influence of such variants on the outcome operates exclusively via the exposure, negating the possibility of tangential channels (15). The intricacies of the MR framework can be gleaned from Figure 1. Pertinent details regarding summary statistics data sources are cataloged in Table 1 and the Supplementary File.
Figure 1. Overview of research design and analysis strategy. Overview of the research design. Exposures come from UKB, with outcomes including T2DM. The MR framework is based on three fundamental MR assumptions, with MVMR analyses adjusting for eight mediating factors for positive results. MVMR: Multivariate Mendelian Randomization; UVMR: Univariate Mendelian Randomization; BMI: Body Mass Index; SNP, Single Nucleotide Polymorphism; MR-PRESSO, MR Pleiotropy Residual Sum and Outlier; HbA1c, Glycated Hemoglobin A1c; LDL-C, Low Density Lipoprotein Cholesterol; HDL-C, High Density Lipoprotein Cholesterol; TG, Triglyceride; T2DM, Type 2 Diabetes Mellitus; CRP, C-Reactive protein; IR, insulin resistance; 25OHD, 25-hydroxyvitamin D; DIAGRAM, DIAbetes Genetics Replication And Meta-analysis; UKB, UK Biobank.
All the encompassed GWAS researchers have secured endorsements from their respective institutional review boards. Given that our endeavor pivots on the secondary scrutiny of accessible data, there wasn’t a necessity for further ethical clearance.
Selection of genetic instrumental variables
The data for oily fish (460,443 individuals) and non-oily fish (460,880 individuals) intake is sourced from the UK Biobank (UKB) Consortium (16). To ensure the accuracy of the MR estimates, these SNPs were required to satisfy specific criteria:
1. All selected Single Nucleotide Polymorphism (SNP)s, when considered as IVs, demonstrated an association with the specific exposure at a genome-wide significance threshold (p < 5 × 10−8).
2. For the prevention of biases from linkage disequilibrium, the SNPs were ensured not to be associated with any potential confounders and to be independent of each other (r2 < 0.001, clumping distance = 10,000 kb).
3. The strength of instrumental variables was evaluated using F-statistics, defined by the formula: F = beta2/se2 (where beta represents the SNP-exposure association and variance is denoted by se). A higher F-statistic signified enhanced instrumental strength, with a mandatory requirement for all SNPs to have an F-statistic exceeding 10 (17).
4. MR-Steiger filtering was integrated to reinforce the validity of our results, thereby omitting variants that had more substantial associations with outcomes than with exposures (18).
5. If the outcome dataset lacks SNP, this study will not use proxy SNPs (r2 > 0.8) to ensure the accuracy of the outcomes.
6. It was essential that the effect of an SNP on exposure and its effect on the outcome both aligned with the identical allele. SNPs not meeting this alignment were excluded.
Source of outcome phenotypes
The summary-level GWAS meta-analysis for T2DM (80,154 cases and 853,816 controls) integrated 32 cohorts, sourced from the AMP-T2D Knowledge Portal and the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) consortium (19). Replication analysis and meta-analysis data for T2DM were sourced from the FinnGen (R9) consortium (20, 57,698 cases and 308,252 controls) (20) and the GWAS by Xue et al. (21). (62,892 cases and 596,424 controls), respectively.
Data sources for possible mediators
We further obtained genetic associations for Body Mass Index (BMI) from the Genetic Investigation of Anthropometric Traits (GIANT) consortium (22), Glycated Hemoglobin A1c (HbA1c) from Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) (23), triglycerides (TG) from Richardson, Tom et al. (24), Low-Density Lipoprotein Cholesterol (LDL-C) and High-Density Lipoprotein Cholesterol (HDL-C) from Global Lipids Genetics Consortium (GLGC) (25), C-Reactive protein (CRP) from Ligthart et al. (26), insulin resistance (IR) from Dupuis et al. (27), 25-hydroxyvitamin D (25OHD) levels from Manousaki et al. (28).
Statistical analyses
Primary MR analysis
In the UVMR analysis, we applied the Wald ratio test to individual IVs. For assessing multiple IVs (≥2), we employed the multiplicative random-effects inverse-variance-weighted (IVW) method, supplemented by MR-Egger and weighted median methods. The weighting in IVW is directly proportional to the Wald ratio estimate and inversely proportional to its variance for each SNP (29). While IVW provides reliable estimates when all genetic variations are valid, the weighted median is more suitable when over half are invalid, and MR-Egger when all are invalid (30). We adjusted for multiple comparisons using the False Discovery Rate (FDR), with a p-value <0.05 post-adjustment denoting a significant causal link. If an unadjusted p-value was <0.05 but the FDR-adjusted p-value was >0.05, it was only suggestive.
For the Mediation MR & MVMR analysis, given potential confounding by factors like BMI, HbA1c, TG, LDL-C, HDL-C, CRP, IR, and 25OHD in the exposure-outcome pathway, we proceeded with MVMR to gauge the direct causal effect. The first MVMR assumption pertains to genetic variation’s link with one or more exposures, while subsequent assumptions mirror those of UVMR (31). We also assessed the proportion of mediated effects. Initially, MR effect estimates for exposure-associated (Oily fish intake) outcomes were derived using IVW. This was followed by MVMR to measure the impact of the eight mediating factors on the outcome. Multiplying the two estimates for each outcome gave the indirect effect of the exposure. By dividing the mediated effect by the total effect, we could determine the contribution of each mediator to the overall effect.
Genetic correlation analysis
Linkage Disequilibrium Score (LDSC) regression, optimized for summary-level GWAS datasets, provides an advanced methodology to analyze genetic correlations in multifaceted diseases or phenotypic traits. This method adeptly differentiates true polygenic associations from confounders like cryptic relatedness or population stratification (32). If a genetic correlation is both statistically and quantitatively significant, it suggests that the overall phenotype correlation is not solely driven by environmental confounding factors (32). The LDSC platform, accessible at https://github.com/bulik/ldsc, serves as a tool for investigating genetic correlations between exposures and their associated outcome phenotypes.
Sensitivity analysis
In the UVMR framework, a battery of tests was conducted to ensure analytical rigor. Cochran’s Q test was employed to gauge the heterogeneity across the selected genetic markers, with a p-value below 0.05 indicating significant discrepancies among the investigated SNPs (33). The potential for directional pleiotropy within the MR schema was assessed via MR-Egger regression (34). An MR-Egger’s intercept p-value under 0.05 suggests notable directional pleiotropy, though this approach may have relative limitations in precision (35). To identify outliers and evaluate horizontal pleiotropy, the MR Pleiotropy Residual Sum and Outlier (MR-PRESSO) technique was adopted, with significance inferred when the global p-value is less than 0.05 (36). By excluding outliers, this method aids in refining our corrections. The leave-one-out technique was also applied to determine the impact of singular SNPs on the cumulative findings (37).
The R2 value was determined using the equation 2 × MAF×(1-MAF) × beta2, where MAF stands for the minor allele frequency for each instrumental SNP. The resulting values were then aggregated to obtain the coefficient essential for the power calculation (38). The statistical power was calculated based on the mRnd website (39).1
Results
Genetic instrument selection and genetic correlation between phenotypes
The study reports F statistics exceeding 40 for all instrumental variants, signifying a robust reduction of bias from weak instruments. The SNPs selected as IVs ranged from 9 to 177, accounting for an explained variance of 0.05 to 27.29%, and the power statistics obtained ranged from 14 to 100% (Supplementary Table S1).
LDSC analysis revealed a significant genetic correlation between oily fish intake and T2DM (rg: -1.102, p = 4.43 × 10−4), but no genetic association was found with non-oily fish intake (rg: 0.024, p = 0.477). The SNP-based liability-scale heritability (h2) spanned from 2.49 to 5.76% (Supplementary Table S2).
Association of genetically predicted oily fish and non-oily fish intake with T2DM
The scatter plot provides a clear visualization of all positive results (Figure 2). After applying the FDR correction for multiple comparisons, our primary IVW analysis demonstrated strong causal evidence (Figure 3). Specifically, for every standard deviation (SD) increase in genetically predicted oily fish intake, the risk of T2DM reduced by 38.6% (OR = 0.614, 95% CI 0.504 ~ 0.748, p = 1.24 × 10−6, FDR = 3.72 × 10−6). Additionally, when the OR was 0.614, we possessed a 100% Power-statistic to detect a correlation between oily fish intake and T2DM. No causal relationship was found between non-oily fish intake and T2DM (OR = 0.902, 95% CI 0.540 ~ 1.505, p = 0.692, FDR = 0.695). In further reverse MR analysis, no causal relationship was observed as well (all p > 0.05, FDR >0.05) (Supplementary Table S3).
Figure 2. Genetic associations with oily fish intake with T2DM at a genome-wide level of significance. (A) Scatter plot. (B) Funnel plot. (C) Forest plots. (D) Leave-one-out plot. Horizontal and vertical lines represent 95% confidence intervals for the genetic associations. T2DM, T2DM, Type 2 Diabetes Mellitus.
Figure 3. Summary of estimates of genetically predicted oily fish and non-oily fish intake for T2DM. FDR, False Discovery Rate; OR odds ratio; CI, confidence interval; T2DM, Type 2 diabetes mellitus; MR, Mendelian randomization.
A series of sensitivity analyses confirmed the robustness of the UVMR results (Supplementary Table S4). MR-Egger detected no horizontal pleiotropy (p > 0.05). The leave-one-out analysis further validated that the causal relationship wasn’t influenced by any single SNP, and the funnel plot showed symmetry (Figure 2; Supplementary Figures S1–S3). The Steiger test indicated that all SNPs passed the test, and the direction of causality remained unchanged, further solidifying the results.
Meta-analysis of oily fish intake’s effect on T2DM across three data sources
To assess the consistency of the causal association between oily fish intake and T2DM across different data sources, primary analyses were replicated based on the FinnGen consortium and the GWAS by Xue et al. Using the GWAS data from Xue et al., the causal effect of oily fish intake on T2DM remained evident (OR = 0.735, 95% CI 0.583 ~ 0.927, p = 0.009); however, the causal effect was no longer significant in the FinnGen consortium data. In the replication analyses of FinnGen Consortium and Xue et al., we had 45 and 95% Power-statistic detection of an association between oily fish intake and T2DM, with ORs of 0.877 and 0.735, respectively (Supplementary Table S1). A series of sensitivity analyses confirmed the robustness of the results. Ultimately, a meta-analysis aggregating three distinct data sources provided consistent evidence of a causal link between oily fish intake and T2DM (OR = 0.728, 95% CI 0.593 ~ 0.895, p = 0.003) (Figure 4).
Figure 4. Mendelian randomization association of genetically predicted oily fish intake with T2DM in meta-analysis. Odds ratios are scaled per SD increase in the genetically predicted oily fish intake. OR, odds ratio; CI, confidence interval; T2DM, type 2 diabetes.
Mediation & MVMR analysis
In the MVMR analysis, only outcome phenotypes with strong associations were selected for multivariable analysis. The MR estimates for oily fish intake and T2DM lost statistical significance after adjusting for BMI, HbA1c, CRP, IR, and 250HD levels, indicating the influence of confounding factors (Table 2). Further mediation MR analysis revealed that BMI mediated the primary mediating effect in the causal association between oily fish and T2DM (mediation effect: 33.57%). Additionally, TG (mediation effect: 8.55%) and 25OHD levels (mediation effect: 11.65%) played roles in partial mediation (Table 3). The causal relationships among the three conformed to the principles of mediation MR.
Table 2. Multivariable Mendelian randomization associations of oily fish intake with T2DM adjusting for confounding factors.
Table 3. Mediation analysis of the mediation effect of oily fish intake on T2DM via eight confounding factors
Discussion
This study utilized a comprehensive MR analysis to investigate the relationship between genetic susceptibility related to oily and non-oily fish intake and T2DM. MR results support previous epidemiological studies (8), suggesting a causal relationship between oily fish consumption and a decreased T2DM risk. Such associations were consistent across replication and meta-analyses. MVMR analysis highlighted the impact of various confounders. Mediation MR revealed that BMI mediated a third of the effect, while TG and 250HD contributed further. The study will delve into the significance of these findings from various perspectives, including the roles of inflammatory and oxidative stress biomarkers, homocysteine metabolism, Omega-3 fatty acids, and mediating factors.
Oily fish, a rich source of omega-3 fatty acids like EPA and DHA, plays a significant role in anti-inflammatory processes. EPA and DHA mitigate inflammation by inhibiting the production of pro-inflammatory eicosanoids derived from arachidonic acid. A double-blind study by Jisun So and colleagues highlights that while DHA exhibits stronger anti-inflammatory effects, EPA notably balances pro-inflammatory and anti-inflammatory proteins (40). Additionally, both EPA and DHA downregulate the expression of inflammatory cytokines. In vitro and animal studies demonstrate their effectiveness in reducing the release of tumor necrosis factor-alpha (TNF-α), interleukin-6 (IL-6), and interferon gamma-induced protein 10 (IP-10), cytokines often overexpressed in T2DM patients (41, 42). Extensive research indicates that chronic inflammation accelerates oxidative stress, leading to cellular and tissue damage, including in pancreatic β-cells (43, 44). Omega-3 fatty acids can lower reactive oxygen species (ROS) levels and enhance the activity of antioxidant enzymes such as superoxide dismutase (SOD) and glutathione peroxidase, pivotal in neutralizing harmful free radicals and maintaining cellular balance. This reduction in oxidative stress damage, in line with Marilena Lepretti’s findings (45), underscores the potential mechanism behind oily fish consumption’s protective role against T2DM.
Elevated homocysteine levels are considered a risk factor for T2DM, with Tao Huang’s MR and meta-analysis establishing a causal relationship (46). Homocysteine is metabolized through two primary pathways: remethylation and transsulfuration. Remethylation involves converting homocysteine back to methionine, a process critically dependent on vitamins B6, B12, and folate. Jie Zhu’s research suggests that a combination of omega-3 fatty acids, folate, and vitamins B6 and B12 is more effective in reducing plasma homocysteine levels (47), with omega-3 fatty acids indirectly supporting this pathway by maintaining healthy levels of these vitamins. Conversely, transsulfuration converts homocysteine to cysteine, which is further used in glutathione synthesis, an effective antioxidant (48). Omega-3 fatty acids enhance the transsulfuration pathway, reducing homocysteine levels and bolstering antioxidant defenses. Additionally, Sohrab A. Shaikh’s review links elevated homocysteine levels with IR (49), where omega-3’s anti-inflammatory and metabolic regulatory functions can mitigate resistance. Carlos O Mendivil’s study also highlights the role of oily fish consumption in regulating homocysteine levels (50). Thus, given the association between homocysteine levels and T2DM, emphasizing oily fish intake is crucial in T2DM management.
While observational studies have inherent limitations, such as potential confounding factors and the challenge of distinguishing causation from correlation, the use of MR analysis can help overcome these biases, aiming to identify genuine causal associations. The omega-3 fatty acids in oily fish have gained significant attention for their possible protective effects against diabetes. Marilena Lepretti’s research suggests that the presence of EPA and DHA in cell membranes might enhance insulin receptor function (45). This is believed to be due to these fatty acids increasing membrane fluidity, which facilitates the insulin receptor’s movement and its interaction with insulin. Silvia V. Conde’s study further supports this, demonstrating that a higher omega-3 intake can enhance insulin signaling and sensitivity (51). Apart from omega-3 s, other nutrients in oily fish, like proteins and fats, have essential roles in blood glucose stabilization. A study from the Joslin Diabetes Center found that a combination of proteins and fats could slow carbohydrate digestion and its absorption, maintaining steady glucose levels (52). Regular fluctuations in blood glucose can lead to increased insulin secretion, which, over time, might result in IR.
Omega-3 fatty acids are also known for their anti-inflammatory properties. They can reduce pro-inflammatory cytokine production, such as tumor necrosis factor-α and interleukin-6, potentially reducing IR (53). Afiat Berbudi’s comprehensive review highlights the importance of these anti-inflammatory effects in decreasing T2DM risk (54). In the context of lipid metabolism, A Jaca’s systematic review indicates that omega-3 s can significantly lower triglyceride levels, possibly by inhibiting Very low-density lipoproteins (VLDL) synthesis and secretion (55). Minglei Ma’s research further clarifies this, suggesting that reduced triglyceride levels are linked to increased insulin sensitivity (56). This provides more context to the mediator MR’s findings about the effects of TG and BMI. Oily fish stands as a substantial source of vitamin D; increased consumption can potentially raise plasma levels of 25-hydroxyvitamin D. On the role of vitamin D, a pivotal randomized controlled trial (RCT) by Anastassios G Pittas and his team delineated the consequential role of vitamin D in the prevention of T2DM. It was found that vitamin D can augment the function of pancreatic β-cells, thereby increasing insulin secretion. Additionally, it can potentiate the action of insulin, resulting in a reduction of blood glucose levels (57).
The study has several strengths. First, our MR analysis novelly establishes the causal relationship between oily fish consumption and T2DM, which is further supported by consistent meta-analysis evidence. Second, the exclusive identification of SNPs as IVs within European populations reduces potential population stratification biases. Third, our rigorous methodologies, paired with F-statistics exceeding 10, minimize potential biases from weak instruments. Fourth, the influence of confounders was examined through MVMR, and mediation MR further elaborated the roles of BMI, TG, and 250HD. Fifth, we included varied statistical model-based sensitivity analyses and “leave-one-out” techniques to bolster our findings’ reliability. However, our study is not without limitations. Notably, while MR-PRESSO analysis attempted to eliminate horizontal pleiotropy, we could not fully exclude it, but supplementary sensitivity evaluations reaffirm our results’ robustness. Also, our dependence on summary-level GWAS data means subgroup analyses were not possible.
Conclusion
To encapsulate, there’s an inverse association between oily fish intake and T2DM risk, suggesting the potential benefits of oily fish intake in T2DM prevention.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
Ethical approval was not required for the studies involving humans because all the encompassed GWAS researchers have secured endorsements from their respective institutional review boards. Given that our endeavor pivots on the secondary scrutiny of accessible data, there wasn’t a necessity for further ethical clearance. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
YZ: Writing – original draft, Writing – review & editing. ER: Investigation, Supervision, Writing – review & editing. CZ: Validation, Writing – review & editing. YW: Writing – review & editing. XC: Writing – review & editing. LL: Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We thank all GWAS participants and investigators for publicly making the summary statistics data available.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2024.1288886/full#supplementary-material
Footnotes
References
1. Oh, MY, Kim, SS, Kim, IJ, Lee, IK, Baek, HS, Lee, HW, et al. Clinical characteristics of diabetic patients transferred to Korean referral hospitals. Diabetes Metab J. (2014) 38:388–94. doi: 10.4093/dmj.2014.38.5.388
2. Dugani, SB, Wood-Wentz, CM, Mielke, MM, Bailey, KR, and Vella, A. Assessment of disparities in diabetes mortality in adults in US rural vs nonrural counties, 1999-2018. JAMA Netw Open. (2022) 5:e2232318. doi: 10.1001/jamanetworkopen.2022.32318
3. Saeed, A, Lambojon, K, Saeed, H, Saleem, Z, Anwer, N, Aziz, MM, et al. Access to insulin products in Pakistan: a National Scale Cross-Sectional Survey on prices, availability, and affordability. Front Pharmacol. (2022) 13:820621. doi: 10.3389/fphar.2022.820621
4. Muralidhara, S, Lucieri, A, Dengel, A, and Ahmed, S. Holistic multi-class classification & grading of diabetic foot ulcerations from plantar thermal images using deep learning. Health Inf Sci Syst. (2022) 10:21. doi: 10.1007/s13755-022-00194-8
5. Low, B-H, Lin, Y-D, Huang, B-W, Chia, T, Bau, J-G, and Huang, H-Y. Impaired microvascular response to muscle stretching in chronic smokers with type 2 diabetes. Front Bioeng Biotechnol. (2020) 8:602. doi: 10.3389/fbioe.2020.00602
6. Alkhatib, A, Tsang, C, Tiss, A, Bahorun, T, Arefanian, H, Barake, R, et al. Functional foods and lifestyle approaches for diabetes prevention and management. Nutrients. (2017) 9:1310. doi: 10.3390/nu9121310
7. Alfaddagh, A, Elajami, TK, Saleh, M, Mohebali, D, Bistrian, BR, and Welty, FK. An omega-3 fatty acid plasma index ≥4% prevents progression of coronary artery plaque in patients with coronary artery disease on statin treatment. Atherosclerosis. (2019) 285:153–62. doi: 10.1016/j.atherosclerosis.2019.04.213
8. Chen, G-C, Arthur, R, Qin, L-Q, Chen, L-H, Mei, Z, Zheng, Y, et al. Association of Oily and Nonoily Fish Consumption and fish oil supplements with incident type 2 diabetes: a large population-based prospective study. Diabetes Care. (2021) 44:672–80. doi: 10.2337/dc20-2328
9. Sahay, RD, Ollberding, NJ, Missoni, S, Novokmet, N, Sarac, J, Sarić, T, et al. Fish and shellfish intake and diabetes in a costal population of the Adriatic. Coll Antropol. (2015) 39:401–9.
10. Zhou, Y, Tian, C, and Jia, C. Association of fish and n-3 fatty acid intake with the risk of type 2 diabetes: a meta-analysis of prospective studies. Br J Nutr. (2012) 108:408–17. doi: 10.1017/S0007114512002036
11. Tin, A, and Kottgen, A. Mendelian randomization analysis as a tool to gain insights into causes of diseases: a primer. J Am Soc Nephrol. (2021) 32:2400–7. doi: 10.1681/ASN.2020121760
12. Davies, NM, and Holmes GDS, MV. Reading Mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ. (2018) 362:k601. doi: 10.1136/bmj.k601
13. Bulik-Sullivan, BK, and Neale, BM. LD score regression distinguishes confounding from Polygenicity in GWAS. Nat Genet. (2015) 47:291–5. doi: 10.1038/ng.3211
14. Tashman, KC, Cui, R, O’Connor, LJ, Neale, BM, and Finucane, HK. Significance testing for small annotations in stratified LD-score regression. medRxiv. (2021). doi: 10.1101/2021.03.13.21249938
15. Lawlor, DA, Harbord, RM, Sterne, JAC, Timpson, N, and Davey, SG. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. (2008) 27:1133–63. doi: 10.1002/sim.3034
16. Cole, JB, Florez, JC, and Hirschhorn, JN. Comprehensive genomic analysis of dietary habits in UK biobank identifies hundreds of genetic associations. Nat Commun. (2020) 11:1467. doi: 10.1038/s41467-020-15193-0
17. Bowden, J, Del Greco, MF, Minelli, C, Davey Smith, G, Sheehan, NA, and Thompson, JR. Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-egger regression: the role of the I2 statistic. Int J Epidemiol. (2016) 45:dyw220–dyw1974. doi: 10.1093/ije/dyw220
18. Hemani, G, Tilling, K, and Davey, SG. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. (2017) 13:e1007081. doi: 10.1371/journal.pgen.1007081
19. Mahajan, A, Spracklen, CN, Zhang, W, Ng, MCY, Petty, LE, Kitajima, H, et al. Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation. Nat Genet. (2022) 54:560–72. doi: 10.1038/s41588-022-01058-3
20. Kurki, MI, Karjalainen, J, Palta, P, Sipilä, TP, Kristiansson, K, Donner, KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. (2023) 613:508–18. doi: 10.1038/s41586-022-05473-8
21. Xue, A, Wu, Y, Zhu, Z, Zhang, F, Kemper, KE, Zheng, Z, et al. Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes. Nat Commun. (2018) 9:2941. doi: 10.1038/s41467-018-04951-w
22. Yengo, L, Sidorenko, J, Kemper, KE, Zheng, Z, Wood, AR, Weedon, MN, et al. Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry. Hum Mol Genet. (2018) 27:3641–9. doi: 10.1093/hmg/ddy271
23. Soranzo, N, Sanna, S, Wheeler, E, Gieger, C, Radke, D, Dupuis, J, et al. Common variants at 10 genomic loci influence hemoglobin a₁(C) levels via glycemic and nonglycemic pathways. Diabetes. (2010) 59:3229–39. doi: 10.2337/db10-0502
24. Richardson, TG, Sanderson, E, Palmer, TM, Ala-Korpela, M, Ference, BA, Davey Smith, G, et al. Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: a multivariable Mendelian randomisation analysis. PLoS Med. (2020) 17:e1003062. doi: 10.1371/journal.pmed.1003062
25. Willer, CJ, Schmidt, EM, Sengupta, S, Peloso, GM, Gustafsson, S, Kanoni, S, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet. (2013) 45:1274–83. doi: 10.1038/ng.2797
26. Ligthart, S, Vaez, A, Võsa, U, Stathopoulou, MG, de Vries, PS, Prins, BP, et al. Genome analyses of >200,000 individuals identify 58 loci for chronic inflammation and highlight pathways that link inflammation and complex disorders. Am J Hum Genet. (2018) 103:691–706. doi: 10.1016/j.ajhg.2018.09.009
27. Dupuis, J, Langenberg, C, Prokopenko, I, Saxena, R, Soranzo, N, Jackson, AU, et al. New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet. (2010) 42:105–16. doi: 10.1038/ng.520
28. Manousaki, D, Mitchell, R, Dudding, T, Haworth, S, Harroud, A, Forgetta, V, et al. Genome-wide association study for vitamin D levels reveals 69 independent loci. Am J Hum Genet. (2020) 106:327–37. doi: 10.1016/j.ajhg.2020.01.017
29. Hemani, G, Zheng, J, Elsworth, B, Wade, KH, Haberland, V, Baird, D, et al. The MR-base platform supports systematic causal inference across the human phenome. eLife. (2018) 7:e34408. doi: 10.7554/eLife.34408
30. Bowden, J, Davey Smith, G, Haycock, PC, and Burgess, S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
31. Burgess, S, and Thompson, SG. Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. (2015) 181:251–60. doi: 10.1093/aje/kwu283
32. Tobin, MD, Minelli, C, Burton, PR, and Thompson, JR. Commentary: development of Mendelian randomization: from hypothesis test to Mendelian deconfounding. Int J Epidemiol. (2004) 33:26–9. doi: 10.1093/ije/dyh016
33. Kulinskaya, E, Dollinger, MB, and Bjørkestøl, K. On the moments of Cochran’s Q statistic under the null hypothesis, with application to the meta-analysis of risk difference. Res Synth Methods. (2020) 11:920. doi: 10.1002/jrsm.1446
34. Burgess, S, and Thompson, SG. Interpreting findings from Mendelian randomization using the MR-egger method. Eur J Epidemiol. (2017) 32:377–89. doi: 10.1007/s10654-017-0255-x
35. Wu, F, Huang, Y, Hu, J, and Shao, Z. Mendelian randomization study of inflammatory bowel disease and bone mineral density. BMC Med. (2020) 18:312. doi: 10.1186/s12916-020-01778-5
36. Verbanck, M, Chen, C-Y, Neale, B, and Do, R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
37. Cheng, H, Garrick, DJ, and Fernando, RL. Efficient strategies for leave-one-out cross validation for genomic best linear unbiased prediction. J Anim Sci Biotechnol. (2017) 8:38. doi: 10.1186/s40104-017-0164-6
38. Guan, W, Steffen, BT, Lemaitre, RN, Wu, JHY, Tanaka, T, Manichaikul, A, et al. Genome-wide association study of plasma N6 polyunsaturated fatty acids within the CHARGE consortium. Circ Cardiovasc Genet. (2014) 7:321–31. doi: 10.1161/CIRCGENETICS.113.000208
39. Brion, M-JA, Shakhbazov, K, and Visscher, PM. Calculating statistical power in Mendelian randomization studies. Int J Epidemiol. (2013) 42:1497–501. doi: 10.1093/ije/dyt179
40. So, J, Wu, D, Lichtenstein, AH, Tai, AK, Matthan, NR, Maddipati, KR, et al. EPA and DHA differentially modulate monocyte inflammatory response in subjects with chronic inflammation in part via plasma specialized pro-resolving lipid mediators: a randomized, double-blind, crossover study. Atherosclerosis. (2021) 316:90–8. doi: 10.1016/j.atherosclerosis.2020.11.018
41. Szabó, Z, Marosvölgyi, T, Szabó, É, Bai, P, Figler, M, and Verzár, Z. The potential beneficial effect of EPA and DHA supplementation managing cytokine storm in coronavirus disease. Front Physiol. (2020) 11:752. doi: 10.3389/fphys.2020.00752
42. Cheng, M, Zhang, S, Ning, C, and Huo, Q. Omega-3 fatty acids supplementation improve nutritional status and inflammatory response in patients with lung Cancer: a randomized clinical trial. Front Nutr. (2021) 8:686752. doi: 10.3389/fnut.2021.686752
43. Alarifi, S, Ali, D, Alkahtani, S, and Almeer, RS. ROS-mediated apoptosis and genotoxicity induced by palladium nanoparticles in human skin malignant melanoma cells. Oxidative Med Cell Longev. (2017) 2017:8439098–10. doi: 10.1155/2017/8439098
44. Chang, Y-Q, Zhou, G-J, Wen, H-M, He, D-Q, Xu, C-L, Chen, Y-R, et al. Treatment of radiation-induced brain injury with bisdemethoxycurcumin. Neural Regen Res. (2023) 18:416–21. doi: 10.4103/1673-5374.346549
45. Lepretti, M, Martucciello, S, Burgos Aceves, MA, Putti, R, and Lionetti, L. Omega-3 fatty acids and insulin resistance: focus on the regulation of mitochondria and endoplasmic reticulum stress. Nutrients. (2018) 10:350. doi: 10.3390/nu10030350
46. Huang, T, Ren, J, Huang, J, and Li, D. Association of homocysteine with type 2 diabetes: a meta-analysis implementing Mendelian randomization approach. BMC Genomics. (2013) 14:867. doi: 10.1186/1471-2164-14-867
47. Zhu, J, Xun, P-C, Kolencik, M, Yang, K-F, Fly, AD, and Kahe, K. Do B vitamins enhance the effect of Omega-3 polyunsaturated fatty acids on cardiovascular diseases? A systematic review of clinical trials. Nutrients. (2022) 14:1608. doi: 10.3390/nu14081608
48. Miran, M, Salami, M, and Emam-Djomeh, Z. Potential role of functional foods and antioxidants in relation to oxidative stress and Hyperhomocysteinemia In: MI Waly, editor. Nutritional management and metabolic aspects of Hyperhomocysteinemia. Cham: Springer International Publishing (2021). 177–97.
49. Shaikh, SA, Varatharajan, R, and Muthuraman, A. Palm oil derived Tocotrienol-Rich fraction attenuates vascular dementia in type 2 diabetic rats. Int J Mol Sci. (2022) 23:13531. doi: 10.3390/ijms232113531
50. Mendivil, CO. Fish consumption: a review of its effects on metabolic and hormonal health. Nutr Metab Insights. (2021) 14:11786388211022378. doi: 10.1177/11786388211022378
51. Tucker, LA. Macronutrient intake and insulin resistance in 5665 randomly selected, non-diabetic U.S. Adults Nutrients. (2022) 14:918. doi: 10.3390/nu14050918
52. Joslin Education Team, Carbs, Protein and Fats – Their effect on glucose levels. Joslin Diabetes Center (2021) United States
53. Molfino, A, Amabile, MI, Monti, M, and Muscaritoli, M. Omega-3 polyunsaturated fatty acids in critical illness: anti-inflammatory, Proresolving, or both? Oxidative Med Cell Longev. (2017) 2017:5987082. doi: 10.1155/2017/5987082
54. Berbudi, A, Rahmadika, N, Tjahjadi, AI, and Ruslami, R. Type 2 diabetes and its impact on the immune system. Curr Diabetes Rev. (2020) 16:442–9. doi: 10.2174/1573399815666191024085838
55. Jaca, A, Durão, S, and Harbron, J. Omega-3 fatty acids for the primary and secondary prevention of cardiovascular disease. S Afr Med J. (2020) 110:1158–9. doi: 10.7196/SAMJ.2020.v110i12.14730
56. Ma, M, Liu, H, Yu, J, He, S, Li, P, Ma, C, et al. Triglyceride is independently correlated with insulin resistance and islet beta cell function: a study in population with different glucose and lipid metabolism states. Lipids Health Dis. (2020) 19:121. doi: 10.1186/s12944-020-01303-w
Keywords: Mendelian randomization, oily fish, type 2 diabetes, causal relationship, Linkage Disequilibrium Score
Citation: Zhang Y, Ren E, Zhang C, Wang Y, Chen X and Li L (2024) The protective role of oily fish intake against type 2 diabetes: insights from a genetic correlation and Mendelian randomization study. Front. Nutr. 11:1288886. doi: 10.3389/fnut.2024.1288886
Edited by:
Marcello Pinti, University of Modena and Reggio Emilia, ItalyReviewed by:
Marcos Edgar Herkenhoff, University of São Paulo, BrazilVanessa Fernandes Rodrigues, University of São Paulo, Brazil
Cosmin Mihai Vesa, University of Oradea, Romania
Copyright © 2024 Zhang, Ren, Zhang, Wang, Chen and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lin Li, 2022721046@yangtzeu.edu.cn