Next Article in Journal
DNA Damage Stress Response and Follicle Activation: Signaling Routes of Mammalian Ovarian Reserve
Next Article in Special Issue
Modifying Effect of the Interleukin-18 Level on the Association between BDNF Methylation and Long-Term Cardiovascular Outcomes in Patients with Acute Coronary Syndrome
Previous Article in Journal
E-Cadherin Expression in Relation to Clinicopathological Parameters and Survival of Patients with Epithelial Ovarian Cancer
Previous Article in Special Issue
Genetic and Functional Effects of Adiponectin in Type 2 Diabetes Mellitus Development
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 22
10.3390/ijms232214380
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mendelian Randomization Indicates a Causal Role for Omega-3 Fatty Acids in Inflammatory Bowel Disease

by
Courtney Astore
,
Sini Nagpal
and
Greg Gibson
*
Center for Integrative Genomics, School of Biological Sciences, Georgia Institute of Technology, Krone EBB1 Building, 950 Atlantic Drive, Atlanta, GA 30332, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(22), 14380; https://doi.org/10.3390/ijms232214380
Submission received: 21 October 2022 / Revised: 16 November 2022 / Accepted: 17 November 2022 / Published: 19 November 2022
(This article belongs to the Special Issue Lipid Metabolism and Genes)
Browse Figures
Review Reports Versions Notes

Abstract

:
Inflammatory bowel disease (IBD) is characterized by chronic inflammation of the gastrointestinal system. Omega-3 (ω3) fatty acids are polyunsaturated fatty acids (PUFAs) that are largely obtained from diet and have been speculated to decrease the inflammatory response that is involved in IBD; however, the causality of this association has not been established. A two-sample Mendelian randomization (MR) was used to assess genetic associations between 249 circulating metabolites measured in the UK Biobank as exposures and IBD as the outcome. The genome-wide association study summary level data for metabolite measurements and IBD were derived from large European ancestry cohorts. We observed ω3 fatty acids as a significant protective association with IBD, with multiple modes of MR evidence replicated in three IBD summary genetic datasets. The instrumental variables that were involved in the causal association of ω3 fatty acids with IBD highlighted an intronic SNP, rs174564, in FADS2, a protein engaged in the first step of alpha-linolenic acid desaturation leading to anti-inflammatory EPA and thence DHA production. A low ratio of ω3 to ω6 fatty acids was observed to be a causal risk factor, particularly for Crohn’s disease. ω3 fatty acid supplementation may provide anti-inflammatory responses that are required to attenuate inflammation that is involved in IBD.

1. Introduction

The human metabolome gauges environmental factors and genetic influences that can be indicative of abnormal biochemical and/or physiological health [1]. While IBD has been previously thought to only impact the Western world, the occurrence is emerging in developing countries, possibly due to modern environmental exposures [2]. Approximately 1.4 million Americans are impacted by IBD, which may be the result of environmental exposures such as diet and smoking status [3]. Identification of environment-metabolite-disease associations may aid in the discovery of diagnostic biomarkers or therapeutic targets, including simple dietary interventions for conditions such as inflammatory bowel disease (IBD), which is a complex gastrointestinal disease that is characterized by immune dysregulation, microbiota alterations, and intestinal permeability. Metabolomic profiles have been previously used to discriminate IBD patients from controls. It has also been shown that the gut microbiome can impact metabolite levels, especially those of individuals with IBD [3,4]. IBD has a complex pathogenesis with multiple genetically and phenotypically distinct subtypes, including Crohn’s disease (CD) and ulcerative colitis (UC). The complexity of IBD pathogenesis also contributes to the multiplicity of possible treatment options [5]. There are currently four classes of approved drugs with hundreds of molecules undergoing various phases of the drug discovery process for IBD [6]. It has also been realized that there are modifiable environmental factors such as diet that can benefit the treatment and management of IBD [7,8].
Omega-3 (ω3) fatty acids are a subtype of polyunsaturated fatty acids (PUFA) which are fat molecules having more than one unsaturated carbon bond in the backbone of the structure. There are three notable forms of ω3 fatty acids, alpha-linolenic acid (ALA), eicosapentaenoic acid (EPA), and docasahexaenoic acid (DHA). Each type of ω3 fatty acid differs slightly in its chemical structure and source; ALA is obtained mostly from plant oils while EPA and DHA are obtained primarily from fish/seafood. ALA is not synthesized in the human body and can only be obtained through diet. However, the body can convert small amounts of ALA into EPA and then to DHA, in a series of steps catalyzed by the FADS2 and FADS1 desaturase enzymes [9]. An extensive literature has investigated the association of ω3 fatty acids with health benefits [10,11] including a likely reduction of the inflammatory responses that are seen in individuals with IBD, although the genetics and causal relationship remain unclear [12,13].
To assess the possible causal associations between ω3 fatty acids and IBD, we herein employ Mendelian randomization (MR) [14]. This approach takes advantage of the invariant nature of genotypes in individuals which can be used as instruments to assess whether an exposure (such as a metabolite) mediates the condition. Essentially, we regress genotypic disease effects on all highly significant metabolite effects that were obtained from genome-wide association studies. MR is advantageous as there are multiple methods for adjusting for confounding influences, and hence assessing causal associations, and one can also assess the direction of effect of the relationship (protective or risk) [15]. A recent MR study assessed the causal association between ω6 fatty acids and IBD [16], concluding that anti-inflammatory PUFAs downstream of FADS1 may be protective. Both pro- and anti-inflammatory mechanisms have been linked with increased levels of various PUFAs; however, such mechanisms remain challenging to unveil [17,18]. Our study focuses on the causal relationship between ω3 fatty acids and IBD, with implications regarding their effect on other autoimmune diseases such as rheumatoid arthritis [19,20].

2. Results

2.1. Causal Association between Circulating ω3 Fatty Acids and IBD

To evaluate whether ω3 fatty acids are likely to be causally involved in the pathogenesis of IBD, we identified 31 IBD-associated SNPs in the Sample one IBD-Genetics Consortium (GC) GWAS (Table S1), 38 SNPs in the Sample two IBD-GC GWAS (Table S2), and subsequently 43 SNPs in the Sample three IBD-FinnGen GWAS (Table S3) as genetic instruments. Using the 31 sample one SNPs, which had an average F-statistic of 10, implying strong confidence as instrumental variables overall, it was observed that increased circulating ω3 fatty acids are likely to mediate decreased risk for IBD (IVW OR 0.87; 95% CI: 0.82–0.92; p = 6.89 × 10–6). The p-values for the other MR methods (i.e., MR-Egger, weighted median, simple mode, and weighted mode), were all significant with p < 10−3 and with ORs less than 1.0, also consistent with a protective role of ω3 fatty acids. Replicating this result for the earlier release of the sample two IBD-GC GWAS weights as genetic instruments, which had an average F-statistic of 7.4, it was further observed that ω3 fatty acids resulted in a decreased risk for IBD (IVW OR 0.79; 95% CI: 0.75–0.89; p = 7.42 × 10−6). The p-values for the other MR methods (i.e., MR-Egger, weighted median, simple mode, and weighted mode), were consistently less than 0.01, and the ORs were less than 1.0. In addition, the independent replication Sample 3 IBD GWAS from FinnGen, for which the SNPs had an average F-statistic of 6.0, it was also observed that ω3 fatty acids resulted in a decreased risk for IBD (IVW OR 0.86; 95% CI: 0.77–0.96; p = 6.9 × 10−3). The ORs were less than one across all methods with p-values < 0.01, except for the simple mode MR method. The slopes representing the causal association of each MR method are shown in Figure 1 for ω3 fatty acids in the Sample one and Sample two IBD-GC GWAS datasets.
When removing the SNPs with an F-statistic < 10 in Sample one, the p-values remained in the range of 10−5 < p < 10−2 with ORs less than 1.0. By Fisher’s method for combining p-values from all three samples, all MR modes were significant at p < 0.0001, except for the simple mode (p < 0.01). Full summary statistics of each MR method are shown in Table 1.
Given the heterogeneity of ω3 fatty acid effect sizes, with one third of the SNPs showing nominal effects in one of the samples in the opposite direction to the overall trend (that is, with positive effects on disease) we performed multiple sensitivity analyses. However, neither Cochran’s Q test p-value for the MR-Egger and IVW methods nor the MR-Egger intercept p-value resulted in p > 0.05, indicating no strong evidence for horizontal pleiotropy involving the genotypes that were used for either IBD GWAS dataset (Table 1). Similarly, leave-one-out analysis indicated the same direction of effect of the causal association despite removal of each genetic instrument for all three IBD GWAS datasets (Figure S1 showing the result for Sample one IBD GWAS). These results imply that the inferred protective effect of ω3 fatty acids for IBD is attributable to multiple variants, not just the major-effect SNP in FADS2.

2.2. Regulatory Effects of the Genetic Instruments in Blood and Colon

To identify which tissue the effect may be mediated in, we examined the eQTL profiles for all SNPs used as ω3 fatty acid genetic IVs for the IBD GWAS. The regulatory effects on gene expression, namely eQTLs, were extracted from the GTEx database for sigmoid colon, transverse colon, and whole blood. There were some instances where a SNP had multiple corresponding genes, but just nine of the SNPs that did not have an eQTL in any of the specified tissues. Comparing the effect sizes of the eQTLs between the sigmoid and transverse colon to whole blood, we observed consistent positive correlation (Figure 2), implying that the genotypes are likely to have similar influences on ω3 fatty acid abundance in all three tissues.
Notably, the largest eQTL effect sizes were consistently observed for FADS2, which encodes a fatty acid desaturase, and showed the largest GWAS association with IBD. Since FADS1 has also been implicated in IBD risk [16], we further contrasted the eQTL effects of the lead GWAS variant at the FADS1/2 locus, rs174564, on both genes, finding a much stronger association with FADS2 in blood and terminal ileum, and similar effect sizes but in opposite directions in the colon (Table S4). Intriguingly, the sign of the regulatory effect is also reversed in the liver, where FADS1 is more strongly influenced. Other consistently large associations were observed for DOCK6 and DOCK7 (which are located on different chromosomes), CATSPER2, INCENP, and ZMIZ2, the latter in the opposite direction with respect to the effect of the alternate allele on expression. LPL (encoding lipoprotein lipase) appears to be a blood-specific factor, whereas MAU2, DAGLA (encoding diacylglycerol lipase alpha), and SIK3 may have opposite signs in the sigmoid and transverse colon. These results suggest that further investigation of the role of some of these genes in fatty acid metabolism may illuminate the mechanisms of association with IBD.

2.3. Causal Association of Additional Metabolites for IBD

MR screening of 249 metabolite measures as exposures revealed a total of five additional metabolite measures that have evidence for a causal association with IBD according to the primary MR method, IVW. These were concordant at the multiple testing threshold of 10−4 in both the Sample one IBD GWAS and Sample two IBD GWAS datasets as shown in the volcano plots in Figure 3A,B, respectively.
The same measures were also the most significant in the FinnGen replication Sample three, although only at the nominal p < 0.05 level. Of the five metabolite measures, four were protective (OR < 1.0) and one was a risk factor (OR > 1.0). Notably, the concordant protective associations all involve ω3 fatty acids in some capacity: Docosahexaenoic acid (DHA, the most unsaturated of the lipids), the ratio of DHA to total fatty acids, and the ratio of ω3 fatty acids to total fatty acids. The risk factor that is concordant with both the IBD discovery and replication was the ratio of ω6 fatty acids to ω3 fatty acids, corroborating the protective role of ω3 fatty acids. One source of DHA is ghee, which has been previously speculated to aid in the prevention of some diseases [21].
We also considered whether ω3 fatty acids interact with genetic risk to influence IBD prevalence in the UK Biobank. While ω3 fatty acid, and particularly DHA, levels are significantly reduced in IBD cases relative to the controls, the difference is similar across the spectrum of the IBD polygenic risk score (Figure 4). Similarly, we did not observe any statistical interaction between rs174564 and DHA or ω3 fatty acids and IBD prevalence, implying that genetic and dietary influences have additive effects on risk. It is noteworthy that smoking status, type of bread consumption, and alcohol intake all do strongly interact with polygenic risk, as does fresh fruit intake, but fish consumption does not [22].

3. Discussion

The positive causal association between ω3 fatty acids and IBD was investigated. From this, ω3 fatty acids were identified as one of the most significant protective associations with IBD, with MR evidence in the range of 10−6 < p < 10−2 across all five MR methods in three IBD GWAS samples. This study enhances the findings in [16] regarding PUFA involvement by establishing that it is ω3 rather than ω6 fatty acids that are likely to be protective, and implicating FADS2 over FADS1 as the key enzyme in the small bowel specifically influencing Crohn’s disease risk. We did not find any evidence of a significant bidirectional effect of IBD on ω3 fatty acids.
Multiple lines of evidence support the robustness of our findings. First, the MR results are replicated in two versions of the IBD-GC’s large-scale GWAS that produced different SNP effect weights, and in the independent FinnGen IBD study. Second, five different implementations of two-sample MR that each filter different SNPs based on certain criteria, all support the inference that the ω3 fatty acid-associated genetic instruments jointly exert protective effects. Third, the leave-one-out approach shows that the result is not driven solely by the very strong effect of rs174564, since the exclusion of that variant reduces the significance of the result but does not eliminate it. Also, the core result remains significant when just the five SNPs with F-statistics greater than 10 are used. Fourth, individuals with IBD have lower levels of ω3 fatty acid levels and DHA levels than that of controls in the UK BioBank.
The mechanism by which ω3 fatty acids protect against IBD remains to be clarified. They are generally thought to be anti-inflammatory, but some studies also indicate a pro-inflammatory function, for example in rheumatoid arthritis joints [19]. It is thus important to experimentally evaluate the effect of PUFAs on diverse immune cell subsets in the context of tissue residency. Previous studies have shown that ω3 fatty acids inhibit cyclooxygenase (COX), which produce prostaglandin hormones that have pro-inflammatory effects [23]. This also indicates that ω3 fatty acid-mediated COX inhibition resembles the mechanism of action of aspirin [23]. Modification of hormone and cytokine production is also thought to polarize lymphoid and myeloid subsets, which could have a variety of local and systemic effects. The major genetic risk factor in African Americans [24] 2 PTGER4, the receptor for prostaglandin e2, raising the possibility that dietary effects may be strongly mediated through PUFAs to influence disease in this population.
This study implicates several IBD-associated loci as being involved in fatty acid metabolism, but the molecular mechanisms remain to be elucidated. Of particular interest are the pair of guanine nucleotide exchange factor-encoding genes, DOCK6 and DOCK7, which are better known for their roles in the regulation of cytokinesis in neurite outgrowth, skin, and bone development. The genes are unlinked and hence are clearly independent associations operating in similar directions with respect to ω3 fatty acid abundance and IBD risk. SIK3 encodes a kinase that is implicated in mTOR signaling, and ZMIZ2 is a transcription factor that enhances androgen receptor signaling, but MAU2 is engaged in sister chromatid cohesion; these may not be the causal genes at the loci. Associations with multiple lipoproteins also deserve further investigation.
One gene whose role in ω3 fatty acid production is more direct is the one that is adjacent to the strongest genetic instrument, the intronic SNP, rs174564, in FADS2, a desaturase protein engaged in the first desaturation step of alpha-linolenic acid metabolism from C18:3 to C18:4. Decreased FADS2 expression has been previously linked to the pro-inflammatory response in IBD [11]. A recent MR study [18] showed that this SNP associates with risk of rheumatoid arthritis and is mostly responsible for the inferred causal role of ω3 fatty acids in promoting that autoimmune condition since it drove those MR results. We also observed that FADS2 expression is elevated in both the blood and colon in individuals with the IBD-protective variant, consistent with a key role for that enzyme in both tissue compartments. Notably, rs174564 was identified as an eQTL for FADS1 as well, although with different directions of effect in comparison to FADS2 for some tissues. In addition, we found that the effect allele for rs174564 was attributed with decreased levels of ω3 fatty acid and DHA levels in the UKBB as shown in Figure S3. Further studies are needed to elucidate the reasons for the inverted effects on IBD and RA, and to establish whether these also explain some of the reported negative health effects of taking an excess of ω3 fatty acid supplements.
Despite the compelling evidence for a causal role of ω3 fatty acid in protection against IBD, our study has several limitations. First, although each of the genetic instruments is genome-wide significant, and hence robustly associated with ω3 fatty acid levels in the UK Biobank study, it would be desirable to confirm these effects in more diverse cohorts, particularly with extensive African, Asian, or Hispanic ancestry. Second, although we see a strong correlation in eQTL effects between blood and colon, the study does not establish which tissue type is the more important mediator. Since the protection against ulcerative colitis appears to be much weaker than for Crohn’s disease, a specific colonic role is implicated, but note that the metabolites were measured circulating in serum. The MR results for ω3 fatty acids and Crohn’s disease were significant prior to genetic instrument outlier detection via IVW-radial, but not formally significant with outlier detection, and also not significant for ulcerative colitis as shown in Supplementary Figure S2. Third, although two-sample Mendelian randomization is a powerful discovery tool, one sample MR would confirm the findings and allow for direct evaluation of confounders and covariates that may contribute to pleiotropic influences. The inclusion of ω3 fatty acid measurement in ongoing IBD GWAS would facilitate such studies.
However, these limitations should not detract from the conclusion that dietary supplementation of ω3 fatty acids ought to be reconsidered as a therapeutic option. Our results only directly inform the prevalence of the disease and should not be taken to imply that ω3 fatty acid supplementation will be curative. It may, for example, need to be administered for an extended period, or at a particular age, and likely only benefits a subset of patients. Clear differences in FADS2 haplotypes across ancestries have been proposed to explain variable results in clinical trials involving ω3 fatty acid supplements [25], and dietary shifts may feasibly contribute to population differences in the rates of increase in IBD prevalence [26].

4. Materials and Methods

4.1. Study Design

A two-sample MR was employed to assess the causal associations between genetically predicted ω3 fatty acid levels as exposures and IBD as the outcome [27]. MR was performed to evaluate the hypothesis that the selected genetic instruments, which predict the modifiable exposure ω3 fatty acids, are causally associated with the exposure-related outcome [28]. Stringent criteria for genetic instrument selection, as described in the Supplemental Materials, were implemented in the open source 2SampleMR package in R [29], also supplemented with additional controls. The study design is shown in Figure 5 and includes two releases of the IBD-GC GWAS summary statistics [30,31] and an independent FinnGen GWAS [32] for IBD associations, as well as metabolome GWAS summary statistics from the UK Biobank.
We followed the STROBE-MR guidelines for performance and reporting of MR studies [33,34]. Details on the summary statistics are provided in the Supplementary Methods online, along with the criteria that was used for genetic instrument selection.

4.2. Statistical Analysis

There were three methods that were performed to assess the potential impact of pleiotropy involving the selected ω3 fatty acid genetic instruments [35]: Cochran’s Q test [36], MR-Egger intercept test [37], and leave-one-out analysis [38], also explained in more detail in the Supplement [39]. A total of five MR methods were used to assess the causal relationship between metabolites and IBD. The IVW method was utilized as the primary method as it assumes the validity of all genetic instruments, while the other four (MR-Egger, weighted median, simple mode, and weighted mode) serve as sensitivity checks [30]. The IVW-radial method was used for outlier detection and removal [40]. These ensure that the causal association is replicated across multiple MR methods despite their differences in SNP selection for generating the regression [41].

4.3. Gut vs. Blood eQTL Analysis

To further explore the possible functional impact of the genetic instruments determined by MR to mediate the causal relationship between ω3 fatty acids and IBD, whole blood, transverse colon, and sigmoid colon eQTL normalized effect sizes and p-values were extracted from the Genotype-Tissue Expression (GTEx) database for each of the gene-SNP pairs using the GTEx eQTL calculator [42]. The nearby genes for each SNP were annotated via ANNOVAR [43] Notably, some of the SNPs had an association with the expression of multiple genes, in which case every possible pair was included to minimize bias. Whole blood, gut, and liver eQTL normalized effect sizes were compared to assess the tissue specificity of the lead FADS1/2 SNP.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232214380/s1.

Author Contributions

Conceptualization, C.A. and G.G.; formal analysis, C.A.; investigation, S.N.; writing—original draft preparation, C.A.; writing—review and editing, G.G.; visualization, C.A.; supervision, G.G.; funding acquisition, G.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by R01 DK119991 awarded to GG as an ancillary study for the NIH NIDDK IBD-Genetics Consortium.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The GWAS summary level data can be obtained via MR-Base. The eQTL summary level data can be obtained via GTEx. UK Biobank data must be requested via www.ukbiobank.ac.uk accessed on 30 May 2022.

Acknowledgments

We thank Judy Cho and Ted Denson for comments on the manuscript. We are indebted to the UK Biobank study participants study data for whom was obtained for this study under protocol 17984 registered to GG.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aldars-García, L.; Gisbert, J.P.; Chaparro, M. Metabolomics insights into inflammatory bowel disease: A comprehensive review. Pharmaceuticals 2021, 14, 1190. [Google Scholar] [CrossRef] [PubMed]
  2. Kaplan, G.G.; Windsor, J.W. The four epidemiological stages in the global evolution of inflammatory bowel disease. Nat. Rev. Gastroenterol. Hepatol. 2021, 18, 56–66. [Google Scholar] [CrossRef]
  3. Abraham, C.; Cho, J.H. Inflammatory bowel disease. N. Engl. J. Med. 2009, 361, 2066–2078. [Google Scholar] [CrossRef] [PubMed]
  4. Cai, Z.; Wang, S.; Li, J. Treatment of inflammatory bowel disease: A comprehensive review. Front. Med. 2021, 8, 765474. [Google Scholar] [CrossRef] [PubMed]
  5. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef]
  6. Carreras-Torres, R.; Ibáñez-Sanz, G.; Obón-Santacana, M.; Duell, E.J.; Moreno, V. Identifying environmental risk factors for inflammatory bowel diseases: A Mendelian randomization study. Sci. Rep. 2020, 10, 19273. [Google Scholar] [CrossRef]
  7. Rizzello, F.; Spisni, E.; Giovanardi, E.; Imbesi, V.; Salice, M.; Alvisi, P.; Valerii, M.C.; Gionchetti, P. Implications of the westernized diet in the onset and progression of IBD. Nutrients 2019, 11, 1033. [Google Scholar] [CrossRef] [Green Version]
  8. Lattka, E.; Illig, T.; Koletzko, B.; Heinrich, J. Genetic variants of the FADS1 FADS2 gene cluster as related to essential fatty acid metabolism. Curr. Opin. Lipidol. 2010, 21, 64–69. [Google Scholar] [CrossRef]
  9. Shahidi, F.; Ambigaipalan, P. Omega-3 polyunsaturated fatty acids and their health benefits. Annu. Rev. Food Sci. Technol. 2018, 9, 345–381. [Google Scholar] [CrossRef]
  10. Merino, D.M.; Ma, D.W.; Mutch, D.M. Genetic variation in lipid desaturases and its impact on the development of human disease. Lipids Health Dis. 2010, 9, 63. [Google Scholar] [CrossRef]
  11. Marton, L.T.; Goulart, R.D.A.; Carvalho, A.C.A.D.; Barbalho, S.M. Omega fatty acids and inflammatory bowel diseases: An overview. Int. J. Mol. Sci. 2019, 20, 4851. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Barbalho, S.M.; de Alvares Goulart, R.; Quesada, K.; Bechara, M.D.; de Carvalho, A.D.C.A. Inflammatory bowel disease: Can omega-3 fatty acids really help? Ann. Gastroenterol. 2016, 29, 37–43. [Google Scholar] [PubMed]
  13. Sanderson, E.; Glymour, M.M.; Holmes, M.V.; Kang, H.; Morrison, J.; Munafò, M.R.; Palmer, T.; Schooling, C.M.; Wallace, C.; Zhao, Q.; et al. Mendelian randomization. Nat. Rev. Methods Prim. 2022, 2, 6. [Google Scholar] [CrossRef]
  14. Evans, D.M.; Davey Smith, G. Mendelian randomization: New applications in the coming age of hypothesis-free causality. Annu. Rev. Genom. Hum. Genet. 2015, 16, 327–350. [Google Scholar] [CrossRef]
  15. Di’Narzo, A.F.; Houten, S.M.; Kosoy, R.; Huang, R.; Vaz, F.M.; Hou, R.; Wei, G.; Wang, W.; Comella, P.H.; Dodatko, T.; et al. Integrative analysis of the inflammatory bowel disease serum metabolome improve our understanding of genetic etiology and points to novel putative therapeutic targets. Gastroenterology 2022, 162, 828–843. [Google Scholar] [CrossRef]
  16. Giacobbe, J.; Benoiton, B.; Zunszain, P.; Pariante, C.M.; Borsini, A. The anti-inflammatory role of omega-3 polyunsaturated fatty acids metabolites in pre-clinical models of psychiatric, neurodegenerative, and neurological disorders. Front. Psych. 2020, 11, 122. [Google Scholar] [CrossRef] [Green Version]
  17. Natto, Z.S.; Yaghmoor, W.; Alshaeri, H.K.; Van Dyke, T.E. Omega-3 fatty acids effects on inflammatory biomarkers and lipid profiles among diabetic and cardiovascular disease patients: A systematic review and meta-analysis. Sci. Rep. 2019, 9, 18867. [Google Scholar] [CrossRef] [Green Version]
  18. Zhu, G.; Zhou, S.; Xu, Y.; Gao, R.; Li, H.; Zhai, B.; Liu, X.; He, Y.; Wang, X.; Han, G.; et al. Mendelian randomization study on the causal effects of omega-3 fatty acids on rheumatoid arthritis. Clin. Rheumatol. 2022, 41, 1305–1312. [Google Scholar] [CrossRef]
  19. Kostoglou-Athanassiou, I.; Athanassiou, L.; Athanassiou, P. The effect of omega-3 fatty acids on rheumatoid arthritis. Mediterr. J. Rheumatol. 2020, 31, 190–194. [Google Scholar] [CrossRef]
  20. Joshi, K.S. Docosahexaenoic acid content is significantly higher in ghrita prepared by traditional Ayurvedic method. J. Ayurveda Integr. Med. 2014, 5, 85–88. [Google Scholar] [CrossRef]
  21. Nagpal, S.; Tandon, R.; Gibson, G. Canalization of the polygenic risk for common diseases and traits in the UK Biobank cohort. Mol. Biol. Evol. 2022, 39, msac053. [Google Scholar] [CrossRef] [PubMed]
  22. Massaro, M.; Habib, A.; Lubrano, L.; Turco, S.D.; Lazzerini, G.; Bourcier, T.; Weksler, B.B.; Caterina, R.D. The omega-3 fatty acid docosahexaenoate attenuates endothelial cyclooxygenase-2 induction through both NADP(H) oxidase and PKC epsilon inhibition. Proc. Natl. Acad. Sci. USA 2006, 103, 15184–15189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Vane, J.R.; Botting, R.M. The mechanism of action of aspirin. Thromb. Res. 2003, 110, 255–258. [Google Scholar] [CrossRef]
  24. Somineni, H.K.; Nagpal, S.; Venkateswaran, S.; Cutler, D.J.; Okou, D.T.; Haritunians, T.; Simpson, C.L.; Begum, F.; Datta, L.W.; Quiros, A.J.; et al. Whole-genome sequencing of African Americans implicates differential genetic architecture in inflammatory bowel disease. Am. J. Hum. Genet. 2021, 108, 431–445. [Google Scholar] [CrossRef]
  25. Chilton, F.H.; Manichaikul, A.; Yang, C.; O’Connor, T.D.; Johnstone, L.M.; Blomquist, S.; Schembre, S.M.; Sergeant, S.; Zec, M.; Tsai, M.Y.; et al. Interpreting clinical trials with omega-3 supplements in the context of ancestry and FADS genetic variation. Front Nutr. 2021, 8, 808054. [Google Scholar] [CrossRef]
  26. Windsor, J.W.; Kaplan, G.G. Evolving epidemiology of IBD. Curr. Gastroenterol. Rep. 2019, 21, 40. [Google Scholar] [CrossRef]
  27. Hemani, G.; Tilling, K.; Davey Smith, G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017, 13, e1007081. [Google Scholar]
  28. Eilbeck, K.; Quinlan, A.; Yandell, M. Settling the score: Variant prioritization and Mendelian disease. Nat. Rev. Genet. 2017, 18, 599–612. [Google Scholar] [CrossRef]
  29. Hemani, G.; Zheng, J.; Elsworth, B.; Wade, K.H.; Haberland, V.; Baird, D.; Laurin, C.; Burgess, S.; Bowden, J.; Langdon, R.; et al. The MR-Base platform supports systematic causal inference across the human phenome. eLife 2018, 7, e34408. [Google Scholar] [CrossRef]
  30. De Lange, K.M.; Moutsianas, L.; Lee, J.C.; Lamb, C.A.; Luo, Y.; Kennedy, N.A.; Jostins, L.; Rice, D.L.; Gutierrez-Achury, J.; Ji, S.-G.; et al. Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease. Nat. Genet. 2017, 49, 256–261. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, J.Z.; Van Sommeren, S.; Huang, H.; Ng, S.C.; Alberts, R.; Takahashi, A.; Ripke, S.; Lee, J.C.; Jostins, L.; Shah, T.; et al. Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations. Nat. Genet. 2015, 47, 979–986. [Google Scholar] [CrossRef] [PubMed]
  32. Kurki, M.I.; Karjalainen, J.; Palta, P.; Sipilä, T.P.; Kristiansson, K.; Donner, K.; Reeve, M.P.; Laivuori, H.; Aavikko, M.; Kaunisto, M.A.; et al. FinnGen: Unique genetic insights from combining isolated population and national health register data. medRxiv 2022. [Google Scholar] [CrossRef]
  33. Skrivankova, V.W.; Richmond, R.C.; Woolf, B.A.; Yarmolinsky, J.; Davies, N.M.; Swanson, S.A.; VanderWeele, T.J.; Higgins, J.P.T.; Timpson, N.J.; Dimou, N.; et al. Strengthening the reporting of observational studies in epidemiology using Mendelian randomization: The STROBE-MR statement. JAMA 2021, 326, 1614–1621. [Google Scholar] [CrossRef] [PubMed]
  34. Bowden, J.; Del Greco, M.F.; Minelli, C.; Davey Smith, G.; Sheehan, N.; Thompson, J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat. Med. 2017, 36, 1783–1802. [Google Scholar] [CrossRef] [Green Version]
  35. Bowden, J.; Del Greco, M.F.; Minelli, C.; Zhao, Q.; Lawlor, D.A.; Sheehan, N.A. Improving the accuracy of two-sample summary-data Mendelian randomization: Moving beyond the NOME assumption. Int. J. Epidemiol. 2019, 48, 728–742. [Google Scholar] [CrossRef] [Green Version]
  36. Hemani, G.; Bowden, J.; Davey Smith, G. Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum. Mol. Genet. 2018, 27, R195–R208. [Google Scholar] [CrossRef] [Green Version]
  37. Bowden, J.; Davey Smith, G.; Haycock, P.C.; Burgess, S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 2016, 40, 304–314. [Google Scholar] [CrossRef] [Green Version]
  38. Zhu, X. Mendelian randomization and pleiotropy analysis. Quant. Biol. 2021, 9, 122–132. [Google Scholar] [CrossRef]
  39. Bowden, J.; Spiller, W.; Del Greco, M.F.; Sheehan, N.; Thompson, J.; Minelli, C.; Davey Smith, G. Improving the visualization, interpretation and analysis of two-sample summary data Mendelian randomization via the Radial plot and Radial regression. Int. J. Epidemiol. 2018, 47, 1264–1278. [Google Scholar] [CrossRef] [Green Version]
  40. Yang, J.; He, X.; Qian, L.; Zhao, B.; Fan, Y.; Gao, F.; Yan, B.; Zhu, F.; Ma, X. Association between plasma proteome and childhood neurodevelopmental disorders: A two-sample Mendelian randomization analysis. EBioMedicine 2022, 78, 103948. [Google Scholar] [CrossRef]
  41. Burgess, S.; Butterworth, A.; Thompson, S.G. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 2013, 37, 658–665. [Google Scholar] [CrossRef] [PubMed]
  42. GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 2020, 369, 1318–1330. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucl. Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 14380 g001 550
Figure 1. (A) Slopes representing the causal association for each MR method for ω3 fatty acids on the Sample one IBD GWAS. (B) Slopes representing the causal association for each MR method for ω3 fatty acids on the Sample two IBD GWAS. Each point is labeled with the SNP and nearby gene(s).
Figure 1. (A) Slopes representing the causal association for each MR method for ω3 fatty acids on the Sample one IBD GWAS. (B) Slopes representing the causal association for each MR method for ω3 fatty acids on the Sample two IBD GWAS. Each point is labeled with the SNP and nearby gene(s).
Ijms 23 14380 g001
Ijms 23 14380 g002 550
Figure 2. (A) eQTL effect size comparisons between the sigmoid colon vs. whole blood. The color of the points represent the −log10(whole blood p-value (eQTL)) and size of the points represent −log10(sigmoid colon p-value (eQTL)). (B) eQTL effect size comparisons between the transverse colon vs. whole blood. The color of the points represent the −log10(whole blood p-value (eQTL)) and size of the points represent −log10(transverse colon p-value (eQTL)). The grey boxes indicate SNPs that generally have eQTL effects with p > 0.05 in both tissues.
Figure 2. (A) eQTL effect size comparisons between the sigmoid colon vs. whole blood. The color of the points represent the −log10(whole blood p-value (eQTL)) and size of the points represent −log10(sigmoid colon p-value (eQTL)). (B) eQTL effect size comparisons between the transverse colon vs. whole blood. The color of the points represent the −log10(whole blood p-value (eQTL)) and size of the points represent −log10(transverse colon p-value (eQTL)). The grey boxes indicate SNPs that generally have eQTL effects with p > 0.05 in both tissues.
Ijms 23 14380 g002
Ijms 23 14380 g003 550
Figure 3. (A) IVW method MR results of screening of 249 metabolites as exposures on IBD (Sample one). (B) IVW method MR results of screening of 249 metabolites as exposures on IBD (Sample two). Each point represents the causal relationship between a metabolite and the outcome. Red points indicate that the metabolite is a risk factor for IBD, and blue points indicate that the metabolite is protective for IBD. The dashed grey line represents the multiple testing p-value threshold of 1.0 × 10−4 (0.05/2 × 249).
Figure 3. (A) IVW method MR results of screening of 249 metabolites as exposures on IBD (Sample one). (B) IVW method MR results of screening of 249 metabolites as exposures on IBD (Sample two). Each point represents the causal relationship between a metabolite and the outcome. Red points indicate that the metabolite is a risk factor for IBD, and blue points indicate that the metabolite is protective for IBD. The dashed grey line represents the multiple testing p-value threshold of 1.0 × 10−4 (0.05/2 × 249).
Ijms 23 14380 g003
Ijms 23 14380 g004 550
Figure 4. (A) Density plot comparing IBD case vs. control ω3 fatty acid levels. (B) Box plot comparing IBD case vs. control ω3 fatty acid levels. The number above the boxes represents the Welch’s t-test p-value. (C) The prevalence vs. risk of IBD for high vs. low ω3 fatty acid levels. (D) Density plot comparing IBD case vs. control DHA levels. (E) Box plot comparing IBD case vs. control DHA levels. The number above the boxes represents the Welch’s t-test p-value. (F) The prevalence vs. risk of IBD for high vs. low DHA levels. (G) Density plot comparing IBD case vs. control ω6 fatty acid to ω3 fatty acid ratio. (H) Box plot comparing IBD case vs. control ω6 fatty acid to ω3 fatty acid ratio. The number above the boxes represents the Welch’s t-test p-value. (I) The prevalence vs. risk of IBD for high vs. low ω6 fatty acid to ω3 fatty acid ratio.
Figure 4. (A) Density plot comparing IBD case vs. control ω3 fatty acid levels. (B) Box plot comparing IBD case vs. control ω3 fatty acid levels. The number above the boxes represents the Welch’s t-test p-value. (C) The prevalence vs. risk of IBD for high vs. low ω3 fatty acid levels. (D) Density plot comparing IBD case vs. control DHA levels. (E) Box plot comparing IBD case vs. control DHA levels. The number above the boxes represents the Welch’s t-test p-value. (F) The prevalence vs. risk of IBD for high vs. low DHA levels. (G) Density plot comparing IBD case vs. control ω6 fatty acid to ω3 fatty acid ratio. (H) Box plot comparing IBD case vs. control ω6 fatty acid to ω3 fatty acid ratio. The number above the boxes represents the Welch’s t-test p-value. (I) The prevalence vs. risk of IBD for high vs. low ω6 fatty acid to ω3 fatty acid ratio.
Ijms 23 14380 g004
Ijms 23 14380 g005 550
Figure 5. Methodology depicting the MR study design for selecting ω3 fatty acid genetic instruments and assessing the causal relationship between circulating ω3 fatty acids on IBD.
Figure 5. Methodology depicting the MR study design for selecting ω3 fatty acid genetic instruments and assessing the causal relationship between circulating ω3 fatty acids on IBD.
Ijms 23 14380 g005
Table
Table 1. MR results for ω3 fatty acid IVs on IBD (Samples one, two, and three).
Table 1. MR results for ω3 fatty acid IVs on IBD (Samples one, two, and three).
MR MethodIBD GWASNumber of SNPsbSEp-ValueOR (95% CI)Cochran’s QHeterogeneity Test p-ValueMR-Egger Intercept p-Value
MR EggerSample one31−0.170.040.00050.85
(0.78, 0.92)		30.870.370.37
Sample two38−0.240.060.00030.79
(0.70, 0.89)		42.020.230.34
Sample three43−0.240.080.0030.79
(0.68, 0.91)		33.000.800.09
Weighted medianSample one31−0.160.041.4 × 10−50.85
(0.79, 0.92)		......
Sample two38−0.220.051.6 × 10−50.80
(0.72, 0.89)		......
Sample three43−0.190.070.0080.83
(0.72, 0.95)		......
Inverse variance weightedSample one31−0.140.036.9 × 10−60.87
(0.82, 0.92)		31.760.38..
Sample two38−0.200.047.4 × 10−60.82
(0.75, 0.89)		43.110.23..
Sample three43−0.150.060.0070.86
(0.77, 0.96)		36.000.72..
Simple modeSample one31−0.310.110.0080.73
(0.59, 0.91)		......
Sample two38−0.280.140.050.76
(0.58, 0.99)		......
Sample three43−0.110.160.50.90
(0.66, 1.23)		......
Weighted modeSample one31−0.150.040.00020.86
(0.80, 0.92)		......
Sample two38−0.210.042.7 × 10−50.81
(0.74, 0.88)		......
Sample three43−0.190.070.0070.83
(0.73, 0.94)		......
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Astore, C.; Nagpal, S.; Gibson, G. Mendelian Randomization Indicates a Causal Role for Omega-3 Fatty Acids in Inflammatory Bowel Disease. Int. J. Mol. Sci. 2022, 23, 14380. https://doi.org/10.3390/ijms232214380

AMA Style

Astore C, Nagpal S, Gibson G. Mendelian Randomization Indicates a Causal Role for Omega-3 Fatty Acids in Inflammatory Bowel Disease. International Journal of Molecular Sciences. 2022; 23(22):14380. https://doi.org/10.3390/ijms232214380

Chicago/Turabian Style

Astore, Courtney, Sini Nagpal, and Greg Gibson. 2022. "Mendelian Randomization Indicates a Causal Role for Omega-3 Fatty Acids in Inflammatory Bowel Disease" International Journal of Molecular Sciences 23, no. 22: 14380. https://doi.org/10.3390/ijms232214380

APA Style

Astore, C., Nagpal, S., & Gibson, G. (2022). Mendelian Randomization Indicates a Causal Role for Omega-3 Fatty Acids in Inflammatory Bowel Disease. International Journal of Molecular Sciences, 23(22), 14380. https://doi.org/10.3390/ijms232214380

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop