- Research
- Open access
- Published:
Genetic overlap between inflammatory bowel disease and iridocyclitis: insights from a genome-wide association study in a European population
BMC Genomic Data volume 25, Article number: 92 (2024)
Abstract
Background
Inflammatory bowel disease (IBD) is occasionally associated with ophthalmic diseases, including iridocyclitis (IC). The co-occurrence of IBD and IC has been increasingly observed, possibly due to shared genetic structures.
Methods
A three-part analysis was executed utilizing genome-wide association study (GWAS) data on IBD and IC. First, the overall genetic correlation between the two traits was observed using linkage disequilibrium score regression (LDSC). Subsequent to this, a local genetic correlation analysis was conducted utilizing the heritability estimation from summary statistics (HESS) methodology. Finally, the conditional/conjunctional false discovery rate (cond/conjFDR) statistical framework was utilized to ascertain the degree of genetic overlap between the two traits.
Results
Positive overall correlations were observed among IBD, ulcerative colitis (UC), and IC, encompassing both acute/subacute and chronic IC presentations. While a significant correlation was identified between Crohn’s disease (CD) and IC, it was not evident for acute/subacute or chronic IC (P > 0.05). Notably, IBD (encompassing CD and UC) demonstrated local genetic correlations with IC and acute/subacute IC, with pronounced enrichment notably on chromosomes 1 and 6, though such correlations were not observed with chronic IC. The conjFDR analysis confirmed the genetic overlap between the two diseases. The shared genes overlapping between IBD (encompassing CD and UC) and IC were IL23R, GPR35, and ERAP1. For acute/subacute IC and chronic IC, there were six overlapping genes (GPR35, RPL23AP12, IL23R, SNAPC4, ERAP1, and INAVA) and one overlapping gene (INAVA), respectively.
Conclusion
This study confirms the existence of a shared genetic structure between IBD and IC, providing a biological basis for their comorbidity. Additionally, this finding has significant implications for preventing and treating these two diseases.
Introduction
Inflammatory bowel disease (IBD) represents an immune-mediated inflammatory disease primarily impacting the gastrointestinal tract. Clinically, it manifests in two primary forms: Crohn’s disease (CD) and ulcerative colitis (UC) [1, 2]. A recent survey underscores a persistent escalation in the global prevalence of IBD, inevitably imposing substantial economic strains on both society and public healthcare systems [3]. Throughout the diagnostic and therapeutic phases, numerous clinicians have noted that IBD often presents with ophthalmic conditions alongside typical intestinal symptoms [4]. Among them, iridocyclitis (IC) (also known as anterior uveitis) warrants particular attention [5]. As the prevailing form of uveitis, IC exhibits a prevalence of 17% among individuals with IBD [6]. Common clinical manifestations of IC include eye pain, photophobia, redness, swelling, and blurred vision [6]. Although IC does not pose a significant threat to human life in most cases, it can precipitate serious complications, including cataracts and secondary glaucoma, and even carries the risk of blindness [4]. As an extraintestinal manifestation of IBD, IC may exhibit a correlation with the activity of the intestinal disease, potentially increasing in frequency as the disease advances, thereby exerting a notable impact on the quality of life for individuals with IBD [7].
Two recent epidemiological studies have suggested a potential association between IBD and IC [8, 9]. However, these observational studies may be influenced by unmeasured confounding factors, reporting bias from participants, and reverse causality, which could lead to biased results. This has prompted researchers to explore the relationship between the two from a genetic perspective, gradually recognizing that genetic factors play an important role in the development of both IBD and IC. A population-based familial study showed that the risk of first-degree relatives of CD patients developing the disease was 10 times higher than that of first-degree relatives of healthy controls, while the risk for UC was 8 times higher [10]. In addition, twin studies also demonstrated a high heritability for IBD, with 17 out of 20 pairs of twins concordant for CD, and similar findings observed for UC [11]. Similarly, in genetic studies of IC, the incidence rate among first-degree relatives of patients was 13%, compared to only 1% among healthy controls [12]. These studies indicate that both IBD and IC are highly heritable diseases. However, no study has systematically explored the genetic relationship between IBD and IC from a genetic perspective. Therefore, identifying whether there is a shared genetic basis between the two diseases could provide important insights into the underlying pathological mechanisms. If we can identify their genetic overlap at the single-locus level and further determine the significant genetic risk loci, it may lead to the discovery of new therapeutic targets and promote the development of precision medicine.
Based on this, our study conducted a three-part analysis of the genome-wide association study (GWAS) datasets for IBD (including CD and UC) and IC (acute/subacute and chronic IC) to progressively elucidate the genetic basis of these two diseases. First, we quantified the overall genetic correlation between IBD and IC using linkage disequilibrium score regression (LDSC) [13]. Second, utilizing the heritability estimation method from summary statistics (HESS), we further analyzed their local genetic correlation [14]. Finally, using the conditional/conjunction false discovery rate (cond/conjFDR) method, we identified shared genetic variants between the two diseases and conducted an in-depth analysis of the polygenic overlap [15]. Through these approaches, we aim to reveal the potential genetic links between IBD and IC, providing new insights into the pathological mechanisms of these two complex diseases.
Methodology
GWAS data
The GWAS datasets for IBD and its subtypes (CD, UC) were retrieved from the IEU GWAS database (https://gwas.mrcieu.ac.uk/), with corresponding ID numbers ebi-a-GCST004131, ebi-a-GCST004132, and ebi-a-GCST004133. The data for IC and its subtypes (acute/subacute and chronic IC) were selected from the FinnGen database (https://r10.finngen.fi/) [16]. Detailed information about the above GWAS data can be found in Table 1. The flow chart of this study is shown in Fig. 1.
Genetic correlation analysis
The LDSC analysis provides an estimated value of the genetic correlation (Rg) between these two traits [17]. Rg reflects their true overall genetic effect, ranging from − 1 to + 1. A genetic correlation score of + 1 signifies a completely positive genetic association between the two traits, denoting that shared genetic variations exert influence on both traits in the same direction. Conversely, a score of -1 indicates a completely negative genetic correlation, suggesting that shared genetic variations manifest opposite effects on the two traits [17]. The analysis process utilized default parameters related to munge_sumstats.py, -rg, -ref-ld-chr, and -w-ld-chr, which are accessible for download from the website (https://alkesgroup.broadinstitute.org/LDSCORE/). Participants of European ancestry were selected for the LDSC analysis to reduce potential biases caused by changes in population genetic structure.
Local genetic correlation analysis
We employed the HESS to estimate the local genetic correlation in order to examine whether IBD has genetic correlations shared with IC in the local independent regions of the genome.
ρ-HESS calculates local genetic correlation through the following steps [14]:
-
1.
Partitioning the genome: HESS begins by dividing the genome into 1703 independent linkage disequilibrium (LD) regions, each around 1.6 Mb in size. This ensures that SNPs within each region are in LD with one another.
-
2.
Using GWAS summary statistics: HESS relies on summary statistics from genome-wide association studies (GWAS), including the effect sizes and standard errors for each SNP. These data are used to estimate the correlation between traits in specific genomic regions.
-
3.
Calculating local genetic covariance: ρ-HESS computes local genetic covariance using the following formula:
$${\rm{Local}}\,{\rm{genetic}}\,{\rm{covariance}} = \widehat {\beta \:}_f^T{V^{ - 1}}{\widehat {\beta \:}_g}$$(1)where \(\:{\widehat{\beta\:}}_{f}\) and \(\:{\widehat{\beta\:}}_{g}\) represent the effect sizes for two traits, and \(\:V\) is the LD matrix. This computation accounts for sample overlap and LD structure within each region.
-
4.
Regularization: To mitigate noise in the LD matrix, which may arise due to small sample sizes, HESS applies truncated singular value decomposition (truncated-SVD) to regularize the LD matrix, enhancing the accuracy of estimates and reducing standard errors.
-
5.
Standardizing genetic covariance: To make results comparable across different traits and regions, HESS standardizes the local genetic covariance to local genetic correlation using the formula:
$$\begin{array}{l}\:{\rm{Local}}\:{\rm{genetic}}\:{\rm{correlation}} = \\\frac{{{\rm{Local}}\:{\rm{genetic}}\:{\rm{covariance}}}}{{\sqrt { - 1} }}\:\overline {\sqrt {h_{f,{\rm{local}}}^2} } \sqrt {h_{g,{\rm{local}}}^2} \end{array}$$where \(\:{h}_{f,\text{l}\text{o}\text{c}\text{a}\text{l}}^{2}\) and \(\:{h}_{g},{\text{l}\text{o}\text{c}\text{a}\text{l}}^{2}\) represent the local SNP heritability for each trait in that region.
By following these steps, HESS identifies genomic regions with significant local genetic correlations, even when no significant correlation is observed at the genome-wide level. The results require Bonferroni correction(p < 0.05/1703 = 2.94E-05).
Conditional quantile–quantile plots
The enrichment of multi-epitope genes can be elucidated through the development of conditional quantile-quantile (Q-Q) plots. This visualization method illustrates enrichment when the proportion of SNPs linked with the primary phenotype (e.g., IBD) rises in correspondence with the strength of association with the secondary phenotype (e.g., IC) [18]. All Q-Q plots were generated using the precimed/mixer software package in Python 3.11 (https://github.com/precimed/mixer).
CondFDR/ConjFDR analysis
The condFDR/conjFDR methodology offers significant advantages in mining comorbid genes by precisely pinpointing specific shared loci that surpass the significance threshold, which is impossible in traditional GWAS analysis [19]. The steps for performing a conditional false discovery rate (condFDR) analysis are as follows:
-
1.
Compute FDR: The false discovery rate (FDR) is computed using empirical cumulative distribution functions (cdfs). SNP enrichment between the two traits is analyzed to adjust the test statistics for SNPs in the primary phenotype.
-
2.
Build a 2D Look-Up Table: The condFDR value for each SNP, conditional on the secondary phenotype, is calculated and organized into a two-dimensional look-up table. This table helps identify SNPs that are significant in both traits.
-
3.
Control for Multiple Testing: The FDR framework inherently controls for multiple testing, avoiding the loss of power associated with traditional methods such as Bonferroni correction.
-
4.
Conjunctional FDR: By repeating the above steps and switching the roles of the primary and secondary phenotypes, conjunctional FDR (conjFDR) is calculated to identify SNPs that are significantly associated with both traits simultaneously.
This procedure enhances the discovery power of genome-wide association studies (GWAS) and identifies shared genomic loci between traits. We assess the effect size of the shared genetic loci through z score. A z score of “–” indicates that the gene exhibits a protective role in the disease. A z score of “+” indicates that the gene presents a risky role in the disease. The magnitude of the value suggests the intensity of the effect.
Results
Overall genetic correlation
In the LDSC analysis of IBD (encompassing CD and UC) and IC (Table 2), IC was found to be positively correlated with IBD (Rg = 0.4344, P = 0.0012), CD (Rg = 0.3088, P = 0.0132), and UC (Rg = 0.419, P = 0.005). Additionally, significant positive overall correlations were observed between acute/subacute IC and IBD (Rg = 0.523, P = 0.0187) as well as UC (Rg = 0.523, P = 0.0309). In contrast, the correlation with CD (Rg = 0.3617, P = 0.0603) was not significant (Table 2). Similar results were found between IBD (encompassing CD and UC) and chronic IC (Table 2). IBD (Rg = 0.377, P = 0.0154) and UC (Rg = 0.5277, P = 0.0042) showed positive associations with chronic IC, while the correlation with CD (Rg = 0.246, P = 0.1607) did not reach statistical significance.
Regarding genetic susceptibility (Table 2), the genetic susceptibilities of IBD, CD, and UC are 0.3127, 0.4457, and 0.2387, respectively. The genetic susceptibilities of IC, Acute and subacute IC, and Chronic IC in this instance are 0.003, 0.0018, and 0.002, respectively.
Localized genetic correlations
In the local genetic correlation analysis plot, chromosome 1 exhibited enrichment for both IC and IBD, while chromosomes 1 and 6 showed enrichment for CD and UC, respectively (Fig. 2). Visually, it was observed that IBD (encompassing CD and UC) and acute/subacute IC shared local genetic overlap, primarily concentrated on chromosomes 1 and 6, with chromosome 1 being common and particularly noteworthy (Fig. 3). However, the local genetic correlation between IBD (encompassing CD and UC) and chronic IC did not reach statistical significance (Supplementary Fig. 1A-1 C).
HESS analysis of IC and IBD, CD and UC. The top and middle sections of each subgraph represent local genetic correlations and covariances, respectively, and the colored bars represent loci with significant local genetic correlations and covariances. The bottom portion represents the local snp heritability of an individual trait, and the colored bars represent loci with significant local snp heritability. (A) Local genetic correlation between IBD and IC. (B) Local genetic correlation between CD and IC. (C) Local genetic correlation between UC and IC. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
HESS analysis of acute/subacute IC and IBD, CD and UC. The top and middle sections of each subgraph represent local genetic correlations and covariances, respectively, and the colored bars represent loci with significant local genetic correlations and covariances. The bottom portion represents the local snp heritability of an individual trait, and the colored bars represent loci with significant local snp heritability. (A) Local genetic correlation between IBD and acute/subacute IC. (B) Local genetic correlation between CD and acute/subacute IC. (C) Local genetic correlation between UC and acute/subacute IC. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
ConjFDR analysis identifies shared genomic loci between two traits
From the Q-Q plot (Figs. 4, 5 and 6), it was observed that as the P-value of one disease decreases, the curve of the other disease consistently shifts to the left. This suggested a strong correlation between IBD (encompassing CD and UC) and IC (encompassing acute/subacute and chronic IC), indicating shared overlapping genetic risk loci and a common genetic background, indicative of a relationship of genetic enrichment.
Conditional quantile-quantile plot. The dashed line indicates the expected line under the null hypothesis, and the deflection to the left indicates the degree of pleiotropic enrichment. (A) IBD-IC. (B) IC-IBD. (C) CD-IC. (D) IC-CD. (E) UC-IC. (F) IC-UC. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
Conditional quantile-quantile plot. The dashed line indicates the expected line under the null hypothesis, and the deflection to the left indicates the degree of pleiotropic enrichment. (A) IBD-acute/subacute IC. (B) acute/subacute IC-IBD. (C) CD-acute/subacute IC. (D) acute/subacute IC-CD. (E) UC-acute/subacute IC. (F) acute/subacute IC-UC. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
Conditional quantile-quantile plot. The dashed line indicates the expected line under the null hypothesis, and the deflection to the left indicates the degree of pleiotropic enrichment. (A) IBD-acute/subacute IC. (B) chronic IC-IBD. (C) CD-chronic IC. (D) chronic IC-CD. (E) UC-chronic IC. (F) chronic IC-UC. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
Upon setting the screening condition to “conjFDR < 0.05”, high-confidence shared loci between the two traits were obtained. Twelve shared genetic loci and eight genes (IL23R, USP34, GPR35, ERAP1, IKZF1, EPPK1, LSP1, and GVQW3) were identified between IBD and IC (Fig. 7A, Supplementary Table 1). For CD and UC, 17 and 10 genetic risk loci corresponding to IC were identified, involving 12 and 6 genes, respectively (Fig. 7B-C, Supplementary Tables 2–3). Across all three analyses, two SNPs (rs4676410 and rs2014857) and three common genes (IL23R, GPR35, and ERAP1) were consistently identified. To further investigate the overall biological functions and mechanisms shared between IBD and IC, we conducted a functional enrichment analysis of the eight genes using GeneMania [20]. The results revealed significant enrichment in several processes, including positive T cell selection, T cell selection, positive regulation of T cell proliferation, positive regulation of T cell activation, T-helper cell lineage commitment, positive regulation of leukocyte proliferation, and positive regulation of leukocyte activation (Fig. 8).
(A) ConjFDR Manhattan plot of IBD and IC. (B) ConjFDR Manhattan plot of CD and IC. (C) ConjFDR Manhattan plot of UC and IC. The shared risk loci between IC and IBD, CD and UC were marked. The statistically significant causality is defined to be conjFDR < 0.05. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
Regarding acute/subacute IC, it was found to share 15, 19, and 9 genetic risk loci with IBD, CD, and UC, respectively (Fig. 9, Supplementary Tables 4–6). Among them, the overlapping loci included four SNPs (rs4676410, rs3124994, rs2014857, and rs9977672) and six genes (GPR35, RPL23AP12, IL23R, SNAPC4, ERAP1, and INAVA). For chronic IC, the numbers of shared genetic risk loci with IBD, CD, and UC were 5, 3, and 4, respectively (Fig. 10, Supplementary Tables 7–9). The intersecting SNPs were rs374827 and rs34920518, and the common gene was INAVA.
(A) ConjFDR Manhattan plot of IBD and acute/subacute IC. (B) ConjFDR Manhattan plot of CD and acute/subacute IC. (C) ConjFDR Manhattan plot of UC and acute/subacute IC. The shared risk loci between acute/subacute IC and IBD, CD and UC were marked. The statistically significant causality is defined to be conjFDR < 0.05. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
(A) ConjFDR Manhattan plot of IBD and chronic IC. (B) ConjFDR Manhattan plot of CD and chronic IC. (C) ConjFDR Manhattan plot of UC and chronic IC. The shared risk loci between acute/subacute chronic IC and IBD, CD and UC were marked. The statistically significant causality is defined to be conjFDR < 0.05. IC, iridocyclitis; IBD, inflammatory bowel disease; CD, Crohn’s disease; UC, ulcerative colitis
Discussion
In this study, the overall genetic correlation analysis revealed significant positive correlations in all cases except for CD and acute/subacute and chronic IC. The local genetic correlation analysis demonstrated that IC and acute/subacute IC were predominantly associated with IBD, clustering on chromosome 1, while the findings for chronic IC were less pronounced. Moreover, both aspects of the results confirmed that UC has a closer relationship with IC than CD, and acute/subacute IC displays a more intimate association with IBD compared to chronic IC. Q-Q plots were constructed at the SNP level, and the conjFDR statistical method was employed, further verifying the strong correlation between IBD and IC and identifying corresponding genetic risk loci. Consequently, our research findings deepen the understanding of the genetic structure between these two diseases and reveal their overlapping mechanisms.
Previous research has also demonstrated a certain degree of correlation between IBD and IC. A two-sample Mendelian randomization study on the two diseases mentioned that IBD (encompassing CD and UC) has a positive causal link with IC (encompassing acute/subacute and chronic IC) [21]. A nationwide cohort study in Denmark found that the incidence rate ratio of IC in individuals with CD was 8.24 (95% CI: 3.42 to 19.89), while for UC, it was 3.29 (95% CI: 1.71 to 6.29) [8]. In another cross-sectional study involving 47,325 individuals with IBD, the odds ratio for IC in CD patients was 3.6 (95% CI: 2.7 to 4.7), and the incidence rate in UC patients was 2.4 (95% CI: 2.0 to 2.9) [9]. According to a previous analysis involving more than 75,000 patients and controls, the overall genetic susceptibility of IBD was estimated to range from 0.3 to 0.5 [22]. The genetic susceptibility of CD was estimated to be approximately 0.5–0.6 [22]. The genetic susceptibility of UC was estimated to be roughly 0.2–0.3 [22]. A study in Sweden was also approximately at this level [11]. At present, there are relatively few studies conducted on the genetic susceptibility of IC. In patients with positive HLA-B27, the genetic susceptibility of IC is higher (0.4–0.7). Nevertheless, for patients with HLA-B27 negativity, the genetic susceptibility of IC is relatively low [23]. On the whole, our research is largely consistent with the previous data and possesses certain reference value. Moreover, the findings of this research, conducted at the genetic level, offer novel insights into the similarity of pathogenic mechanisms between the two diseases and confirm the existence of genetic overlap between them.
This study focused on the identified risk genes. IL23R, a protein comprising IL-12β1 and IL-23R chains, is a pro-inflammatory cytokine member. Previous GWAS studies have recognized it as an IBD gene. Additionally, a cohort study validated the role of the IL23R gene in increasing susceptibility to CD in North American and European populations [24]. IL23R has been identified as a gene that can influence the genetic progression of IC [25]. INAVA has the capability to induce MAPK and NF-κB pathways, clear cytokines, and intracellular bacteria, potentially contributing to intestinal immune homeostasis and reducing the risk of IBD [26]. Our discovered loci vary to some extent from certain previous studies. For instance, the research on GPR35 in the context of IC was not significant, but corresponding studies have been conducted in IBD. It is a susceptibility gene for the intestine and has been utilized as a target for ameliorating intestinal inflammation and restoring intestinal mucosa [27,28,29]. ERAP1 is a set of enzyme proteins that are of crucial significance for protein digestion and plays a role in autoimmune and autoinflammatory diseases [30]. It is a high-risk gene in both diseases (z > 0). Although ERAP1 was identified as a novel locus gene in a GWAS study on acute IC [31], and another extraction analysis confirmed it as an IC-related gene, it demonstrated a protective effect [32], contrary to our results. ERAP1 can serve as a risk gene for inflammation in the occurrence and development of IBD [33, 34], which is in accordance with our findings. However, research on RPL23AP12 and SNAPC4 in IBD and IC has not yet been conducted and requires further investigation.
The overall biological functions between IBD and IC are enriched in pathways such as T cell activation and leukocyte proliferation, which is an intriguing finding. T cell activation plays a vital role in the pathogenesis of IBD, promoting the occurrence and continuous development of inflammation [35]. Current studies demonstrate that antigen-presenting cells (APCs), such as dendritic cells (DCs) and macrophages, present microbial or environmental antigens in the intestine to T cells, resulting in their excessive activation and thereby driving intestinal inflammation [36]. T cell activation typically involves effector T cell subsets like Th1 and Th17, which, respectively, secrete pro-inflammatory cytokines such as IFN-γ and IL-17 to further recruit and activate neutrophils and macrophages, intensifying the local inflammatory response [37]. FN-γ and IL-17 can also stimulate the proliferation and recruitment of white blood cells (including neutrophils, monocytes, and lymphocytes), which are recruited to the intestinal inflammatory site. They further release reactive oxygen species, enzymes, and pro-inflammatory cytokines, resulting in damage to the intestinal epithelial barrier and the perpetuation of chronic inflammatory responses, aggravating the disease condition of IBD [35]. T cells, particularly CD4 + T helper cells, play a crucial role in the development and progression of IC [38]. In patients with chronic IC, the persistent activation of T cells triggers sustained inflammation, preventing the ocular tissues from returning to normal immune homeostasis [39]. Activated T cells promote the proliferation of leukocytes and their infiltration into the anterior chamber, further exacerbating local inflammation and tissue damage [39]. Additionally, this process is accompanied by the release of large amounts of reactive oxygen species (ROS) and enzymes, which intensify the damage to uveal tissues, ultimately aggravating the severity of IC [40].
The “gut-retina” axis is currently a highly popular research topic concerning the relationship between intestinal diseases and ophthalmic diseases [41]. This axis is influenced by both the ocular and intestinal immune systems, playing a crucial role in the onset and progression of various eye diseases, such as uveitis (encompassing IC), conjunctivitis, and diabetic retinopathy [42]. Moreover, the “gut-retina” axis is recognized as a new frontier in ophthalmic basic and clinical research. Some studies have proposed hypothetical pathways contributing to the disruption of the gut microbiota-ocular surface-lacrimal gland axis [43]. Both IBD and IC, being immune-mediated diseases, interact with each other, thus contributing to the formation of this “gut-retina” axis. The genetic structural overlap between these conditions may serve as a potential basis for elucidating the underlying mechanisms of action.
This study faces several unavoidable limitations. Firstly, it is challenging to rule out the possibility of LD entirely. Despite employing methods such as LDSC, HESS, and conjFDR, which mitigate sample overlap, there remains a risk of exaggerated enrichment results due to potential overlap. Additionally, uncontrollable factors such as behavioral, social, and environmental influences may impact the findings. Moreover, the GWAS conducted in this study focused solely on individuals of European ancestry, potentially limiting the generalizability of our findings to other populations. The incidence of IC does indeed present gender differences. However, the GWAS data involved in this study did not distinguish between males and females, thus preventing the conduct of gender-specific analyses. We anticipate the availability of appropriate GWAS data in the future to enable corresponding research. Lastly, since this study is based on computational simulation methods utilizing publicly available GWAS data, it lacks validation from an independent cohort.
Conclusion
In conclusion, this study demonstrates an overlap in the genetic structure between IBD and IC. SNPs and genes linked to shared risk for both conditions across different scenarios have been identified, potentially offering targets for novel immunotherapies. Future research utilizing larger and more ethnically diverse GWAS samples is expected to further clarify the underlying genetic mechanism, ultimately benefiting the treatment of comorbidities associated with these two diseases.
Data availability
All the GWAS data and statistical software used in this study were publicly available (which can be accessed through the following URLs), and all the generated results in this study were provided in the main text and supplemental data. IEU database: https://gwas.mrcieu.ac.ukFinnGen FinnGen database (https://r10.finngen.fi/) LDSC: https://github.com/bulik/ldsc conjFDR: https://github.com/precimed/pleiofdr FUMA: https://fuma.ctglab.nl.
References
McGovern DPB, Kugathasan S, Cho JH. Genetics of inflammatory bowel diseases. Gastroenterology. 2015;149:1163–e11762. https://doiorg.publicaciones.saludcastillayleon.es/10.1053/j.gastro.2015.08.001
Peloquin JM, Goel G, Villablanca EJ, Xavier RJ. Mechanisms of pediatric inflammatory bowel disease. Annu Rev Immunol. 2016;34:31–64. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev-immunol-032414-112151
Kaplan GG. The global burden of IBD: from 2015 to 2025. Nat Rev Gastroenterol Hepatol. 2015;12:720–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nrgastro.2015.150
Shah J, Shah A, Hassman L, Gutierrez A. Ocular manifestations of inflammatory bowel disease. Inflamm Bowel Dis. 2021;27:1832–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/ibd/izaa359
Harbord M, Annese V, Vavricka SR, Allez M, Barreiro-de Acosta M, Boberg KM, Burisch J, De Vos M, De Vries A-M, Dick AD, et al. The first European evidence-based consensus on extra-intestinal manifestations in inflammatory bowel disease. J Crohn’s Colitis. 2016;10:239–54. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/ecco-jcc/jjv213
Gueudry J, Muraine M. Anterior uveitis. J Français d’Ophtalmologie. 2018;41:e11–21. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jfo.2017.11.003
Greuter T, Bertoldo F, Rechner R, Straumann A, Biedermann L, Zeitz J, Misselwitz B, Scharl M, Rogler G, Safroneeva E, et al. Extraintestinal manifestations of pediatric inflammatory bowel disease: prevalence, presentation, and anti-TNF treatment. J Pediatr Gastroenterol Nutr. 2017;65:200–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1097/MPG.0000000000001455
Burisch J, Jess T, Egeberg A. Incidence of immune-mediated inflammatory diseases among patients with inflammatory bowel diseases in Denmark. Clin Gastroenterol Hepatol. 2019;17:2704–e27123. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.cgh.2019.03.040
Halling ML, Kjeldsen J, Knudsen T, Nielsen J, Hansen LK. Patients with inflammatory bowel disease have increased risk of autoimmune and inflammatory diseases. World J Gastroenterol. 2017;23:6137–46. https://doiorg.publicaciones.saludcastillayleon.es/10.3748/wjg.v23.i33.6137
Orholm M, Munkholm P, Langholz E, Nielsen OH, Sørensen TI, Binder V. Familial occurrence of inflammatory bowel disease. N Engl J Med. 1991;324:84–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1056/NEJM199101103240203
Tysk C, Lindberg E, Järnerot G, Flodérus-Myrhed B. Ulcerative colitis and Crohn’s disease in an unselected population of monozygotic and dizygotic twins. A study of heritability and the influence of smoking. Gut. 1988;29:990–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1136/gut.29.7.990
Derhaag PJ, Linssen A, Broekema N, de Waal LP, Feltkamp TE. A familial study of the inheritance of HLA-B27-positive acute anterior uveitis. Am J Ophthalmol. 1988;105:603–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/0002-9394(88)90051-7
Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh P-R, ReproGen Consortium PG, Consortium, Genetic Consortium for Anorexia Nervosa of the Wellcome Trust Case Control Consortium 3, Duncan L, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47:1236–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ng.3406
Shi H, Mancuso N, Spendlove S, Pasaniuc B. Local genetic correlation gives insights into the shared genetic architecture of complex traits. Am J Hum Genet. 2017;101:737–51. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ajhg.2017.09.022
Smeland OB, Frei O, Shadrin A, O’Connell K, Fan C-C, Bahrami S, Holland D, Djurovic S, Thompson WK, Dale AM, et al. Discovery of shared genomic loci using the conditional false discovery rate approach. Hum Genet. 2020;139:85–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00439-019-02060-2
Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, Reeve MP, Laivuori H, Aavikko M, Kaunisto MA, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. 2023;613:508–18. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41586-022-05473-8
Bulik-Sullivan BK, Loh P-R, Finucane HK, Ripke S, Yang J, Schizophrenia Working Group of the Psychiatric Genomics Consortium, Patterson N, Daly MJ, Price AL, Neale BM. LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet. 2015;47:291–5. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ng.3211
Frei O, Holland D, Smeland OB, Shadrin AA, Fan CC, Maeland S, O’Connell KS, Wang Y, Djurovic S, Thompson WK, et al. Bivariate causal mixture model quantifies polygenic overlap between complex traits beyond genetic correlation. Nat Commun. 2019;10:2417. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41467-019-10310-0
Liley J, Wallace C. A pleiotropy-informed Bayesian false discovery rate adapted to a shared control design finds new disease associations from GWAS summary statistics. PLoS Genet. 2015;11:e1004926. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pgen.1004926
Franz M, Rodriguez H, Lopes C, Zuberi K, Montojo J, Bader GD, Morris Q. GeneMANIA update 2018. Nucleic Acids Res. 2018;46:W60–4. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gky311
Meng Y, Tan Z, Liu C, Dong W, Chen C. Association between inflammatory bowel disease and iridocyclitis: a Mendelian randomization study. J Clin Med. 2023;12:1282. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/jcm12041282
Jostins L, Ripke S, Weersma RK, Duerr RH, McGovern DP, Hui KY, Lee JC, Schumm LP, Sharma Y, Anderson CA, et al. Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature. 2012;491:119–24. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nature11582
Wakefield D, Clarke D, McCluskey P. Recent developments in HLA B27 anterior uveitis. Front Immunol. 2021;11. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fimmu.2020.608134
Libioulle C, Louis E, Hansoul S, Sandor C, Farnir F, Franchimont D, Vermeire S, Dewit O, de Vos M, Dixon A, et al. Novel Crohn disease locus identified by genome-wide association maps to a gene desert on 5p13.1 and modulates expression of PTGER4. PLoS Genet. 2007;3:e58. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pgen.0030058
Huang X-F, Brown MA. Progress in the genetics of uveitis. Genes Immun. 2022;23:57–65. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41435-022-00168-6
Yan J, Hedl M, Abraham C. An inflammatory bowel disease–risk variant in INAVA decreases pattern recognition receptor–induced outcomes. J Clin Invest. 2017;127:2192–205. https://doiorg.publicaciones.saludcastillayleon.es/10.1172/JCI86282
Yansen Z, Lingang Z, Dali L, Mingyao L. Inflammatory bowel disease susceptible gene GPR35 promotes bowel inflammation in mice. Yi Chuan. 2021;43:169–81. https://doiorg.publicaciones.saludcastillayleon.es/10.16288/j.yczz.20-392
Kaya B, Melhem H, Niess JH. GPR35 in intestinal diseases: from risk gene to function. Front Immunol. 2021;12:717392. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fimmu.2021.717392
Otkur W, Wang J, Hou T, Liu F, Yang R, Li Y, Xiang K, Pei S, Qi H, Lin H, et al. Aminosalicylates target GPR35, partly contributing to the prevention of DSS-induced colitis. Eur J Pharmacol. 2023;949:175719. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ejphar.2023.175719
Țiburcă L, Zaha DC, Jurca MC, Severin E, Jurca A, Jurca AD. The role of aminopeptidase erap1 in human pathology-a review. Curr Issues Mol Biol. 2024;46:1651–67. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/cimb46030107
Huang X-F, Li Z, De Guzman E, Robinson P, Gensler L, Ward MM, Rahbar MH, Lee M, Weisman MH, Macfarlane GJ, et al. Genomewide association study of acute anterior uveitis identifies new susceptibility loci. Invest Ophthalmol Vis Sci. 2020;61:3. https://doiorg.publicaciones.saludcastillayleon.es/10.1167/iovs.61.6.3
Gelfman S, Moscati A, Huergo SM, Wang R, Rajagopal V, Parikshak N, Pounraja VK, Chen E, Leblanc M, Hazlewood R, et al. A large meta-analysis identifies genes associated with anterior uveitis. Nat Commun. 2023;14:7300. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41467-023-43036-1
Ondrejčáková L, Gregová M, Bubová K, Šenolt L, Pavelka K. Serum biomarkers and their relationship to axial spondyloarthritis associated with inflammatory bowel diseases. Autoimmun Rev. 2023;23:103512. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.autrev.2023.103512
Castro-Santos P, Moro-García MA, Marcos-Fernández R, Alonso-Arias R, Díaz-Peña R. ERAP1 and HLA-C interaction in inflammatory bowel disease in the Spanish population. Innate Immun. 2017;23:476–81. https://doiorg.publicaciones.saludcastillayleon.es/10.1177/1753425917716527
Mf N. Targeting immune cell circuits and trafficking in inflammatory bowel disease. Nat Immunol. 2019;20. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41590-019-0415-0
Kaser A, Zeissig S, Blumberg RS. Inflammatory bowel disease. Annu Rev Immunol. 2010;28:573–621. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev-immunol-030409-101225
Zhang Y-Z, Li Y-Y. Inflammatory bowel disease: pathogenesis. World J Gastroenterol. 2014;20:91–9. https://doiorg.publicaciones.saludcastillayleon.es/10.3748/wjg.v20.i1.91
Caspi RR. A look at autoimmunity and inflammation in the eye. J Clin Invest. 2010;120:3073–83. https://doiorg.publicaciones.saludcastillayleon.es/10.1172/JCI42440
Oh H-M, Yu C-R, Lee Y, Chan C-C, Maminishkis A, Egwuagu CE. Autoreactive memory CD4 + T lymphocytes that mediate chronic uveitis reside in the bone marrow through STAT3-dependent mechanisms. J Immunol. 2011;187:3338–46. https://doiorg.publicaciones.saludcastillayleon.es/10.4049/jimmunol.1004019
Rr RH. Microbiome and autoimmune uveitis. Front Immunol. 2019;10. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fimmu.2019.00232
Qin X, Zou H, Niu C. The STING pathway: an uncharacterized angle beneath the gut-retina axis. Exp Eye Res. 2022;217:108970. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.exer.2022.108970
Rowan S, Jiang S, Korem T, Szymanski J, Chang M-L, Szelog J, Cassalman C, Dasuri K, McGuire C, Nagai R, et al. Involvement of a gut-retina axis in protection against dietary glycemia-induced age-related macular degeneration. Proc Natl Acad Sci U S A. 2017;114:E4472–81. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.1702302114
Moon J, Yoon CH, Choi SH, Kim MK. Can gut microbiota affect dry eye syndrome? Int J Mol Sci. 2020;21:8443. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms21228443
Acknowledgements
The authors thank Bullet Edits Limited for the linguistic editing and proofreading of the manuscript.
Funding
This work was supported by the National Natural Science Foundation Project of China (Grant No. 82060836,82260938,82374454,82160870).
Author information
Authors and Affiliations
Contributions
Wu Liao: Conceptualization, methodology, formal analysis, data curation, writing-original draft preparation; Qinghua Luo: writing-original draft preparation, visualization; Leichang Zhang and Haiyan Wang: Data curation, writing-original draft preparation, visualization; Wei Ge and Jiawen Wang: writing-review and editing, Data curation, visualization; Zhengyun Zuo: Supervision, writing-review and editing, Funding acquisition. All authors contributed to the article and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable. The data used for analysis were obtained from published studies and public databases. The GWAS database is a database of publicly available datasets, where each study has been approved by local institutional review boards and ethics committees.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Liao, W., Luo, Q., Zhang, L. et al. Genetic overlap between inflammatory bowel disease and iridocyclitis: insights from a genome-wide association study in a European population. BMC Genom Data 25, 92 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12863-024-01274-2
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12863-024-01274-2