Skip to main content
  • Research article
  • Open access
  • Published:

Genotype–phenotype associations of polymorphisms within the gene locus of NOD-like receptor pyrin domain containing 3 in Swiss inflammatory bowel disease patients

This article has been updated

Abstract

Background

Genetic variations within the regulatory region of the gene encoding NOD-like receptor pyrin domain containing 3 (NLRP3) have been associated with Crohn’s Disease (CD). NLRP3 is part of the NLRP3-inflammasome that mediates the maturation of IL-1β and IL-18. Carrying the major allele of the single nucleotide polymorphisms (SNPs) rs10733113, rs4353135 and rs55646866 is associated with an increased risk for CD. We here studied the impact of these polymorphisms on clinical characteristics in patients of the Swiss IBD Cohort Study (SIBDCS).

Methods

We included 981 Crohn’s disease (CD) patients and 690 ulcerative colitis (UC) patients of the SIBDCS. We analyzed whether three CD-associated NLRP3 polymorphisms have an impact on the clinical disease course in these patients.

Results

In CD patients presence of the major allele (G) of rs10733113 was associated with less surgeries and lower maximal CDAI and a similar trend was observed for rs55646866 and rs4353135. Presence of the major allele of all three SNPs was negatively correlated to maximal CDAI. In UC patients homozygous genotype for the major allele (CC) for rs55646866 was associated with a higher age at diagnosis and a higher MTWAI index. Homozygous genotype for the major allele of all three polymorphisms was associated with a higher number of ambulatory visits and longer hospital stays.

Conclusions

In CD patients presence of the major allele of all three polymorphisms was associated with markers of a less severe disease course, while in UC the homozygous genotype for all major alleles suggested a more severe disease activity.

Peer Review reports

Background

Genetic variations in a predicted regulatory region of the NOD-like receptor pyrin domain containing 3 (NLRP3) gene locus have been associated with an increased risk to develop Crohn’s disease (CD) [1]. CD risk has been linked to the major alleles of the NLRP3 single nucleotide polymorphisms (SNPs) rs10733113, rs4353135 and rs55646866; however, 64–85% of the healthy European or Canadian population carry the disease-associated alleles [1]. For an overview of the major allele frequency in different populations and for details on SNP location see Additional file 2: Table S1 and Additional file 1: Figure S1.

The protein NLRP3 belongs to the family of nucleotide binding oligomerization domain (NOD) and leucine rich repeat (LRR) containing receptors (NLRs). It has an important role in the innate immune response, namely sensing of pathogen associated molecular patterns (PAMPs) and danger associated molecular patterns (DAMPs).

Upon activation, NLRP3 assembles with pro-caspase-1 and apoptosis-associated speck-like protein (ASC) to the NLRP3 inflammasome, resulting in autoproteolytic cleavage of pro-caspase-1 to active caspase-1 (CASP1). CASP1 is responsible for the conversion of pro-interleukin (IL)-1β and pro-IL-18 to the active cytokines [2]. Inflammatory bowel disease (IBD) patients tend to have higher levels of these cytokines in the serum [3]. Elevated secretion of IL-1β by lamina propria mononuclear cells isolated from inflamed IBD patients has been reported in several studies [4,5,6], and IL-1β levels in colonic perfusion fluid correlated with disease activity [7] or risk of relapse within the next year [8]. Yet, experimental colitis in mice deficient in NLRP3 could not prove an unequivocal pro-inflammatory role of NLRP3. Depending on the animal model, environment and hygienic conditions, lack of NLRP3 either ameliorated or aggravated experimental colitis [9,10,11,12].

So far, little is known about the functional consequences of human NLRP3 associated genetic variants. The major alleles of rs4353135 and rs10733113 have been shown to lead to lower NLRP3 mRNA expression in human peripheral blood mononuclear cells (PBMCs) [1, 13]. The major allele (GG) of rs6672995 was associated with lower IL-1β secretion from stimulated PBMC [1], suggesting that the NLRP3 polymorphisms that are associated with CD might lead to lower NLRP3-inflammasome activity. In line with the proposed pro-inflammatory role of the NLRP3 inflammasome during established IBD, polymorphisms that affect NLRP3 activity might also affect the clinical course of IBD patients, but such a relationship has not been addressed, so far.

We here studied for the first time whether the CD associated variants of the three SNPs (rs10733113, rs55646866, rs4353135) located in a regulatory region on chromosome 1q44 downstream of NLRP3 affect the disease characteristics in IBD patients of the Swiss IBD cohort study (SIBDCS). In CD patients, the major allele of all three SNPs was negatively correlated to maximal CDAI, suggesting a less severe course of disease. In contrast, in UC patients, presence of the major allele of rs55646866 was associated with a higher MTWAI while homozygous genotype for the major allele of all three polymorphisms was associated with a higher number of ambulatory visits. Our findings might lead to a better understanding of disease pathogenesis depending on NLRP3 genotype and identify patients that might benefit from a NLRP3 targeted therapy.

Methods

Patient data

Data were obtained from the nationwide Swiss Inflammatory Bowel Disease Cohort Study (SIBDCS). The SIBDCS is a prospective multicenter observational population-based study into which patients with IBD from all regions of Switzerland have been included since 2006. The goals and the methodology of the cohort have been described in detail previously [14, 15].

Study design

We included 1671 IBD patients (981 CD and 690 UC) that were enrolled in the SIBDCS at time of data acquisition and had been genotyped for the CD-associated single nucleotide polymorphisms (SNP) rs4353134, rs55646866 and rs10733113 within a regulatory region on chromosome 1q44 downstream of the NLRP3 gene. Genotyping was performed as part of an analysis of the whole SIBDCS for selected SNPs that are currently known to be associated with IBD [16]. SNP Genotyping of SIBDCS samples was performed by MALDI-TOFF based analysis [17].

The aim of this study was to analyze whether the CD-associated risk variants for NLRP3 rs4353134, rs55646866 and rs10733113 are associated with clinical parameters defining the severity/course of disease. For rs4353134, the T-allele is the major allele and the G-allele the minor allele. The major and minor alleles for rs55646866 are C/T and for rs10733113 G/A, respectively.

For analysis of clinical characteristics, we included gender, diagnosis, age at diagnosis, disease duration, complications, Crohn’s disease activity index (CDAI) or modified Truelove and Witts activity index (MTWAI), reported flare (including possible cause of flare, flare management and hospitalization), fistula, abscess or anal fissure, stenosis, presence of extraintestinal manifestations, intestinal surgery and medication.

Statistical analysis

Clinical data were obtained from the data center of the SIBDCS at the University of Lausanne. These data and additional data obtained from a review of the patients' files were entered into a database (Access 2000; Microsoft Switzerland Ltd Liab. Co., Wallisellen, Switzerland). Stata 14 software (StataCorp, 2015, College Station, TX was used for the statistical analysis.

We analyzed for associations between clinical characteristics and genotype for rs4353135, rs55646866 and rs10733113 with Fisher’s exact test for discrete variables and with the Kruskal–Wallis test for continuous variables. To analyze for a cumulative risk of the three variants, we calculated a score based on the odds ratio for the risk allele as described by Villani et al. [1] for all patients (Score = x * log(OR rs4353135 T-allele) + y * log(OR rs55646866 C-allele) + z * log (OR rs10733113 G-allele); x, y, z = number of major allele (0, 1 or 2)) resulting in a score between 0 and 2.58 for each patient. We performed linear regression analyses with this score as predictor and clinical parameters as response. Log(maxCDAI) was used for the analysis to reduce the impact of very high values.

For calculation of linkage disequilibrium (LD) we used the square of the correlation coefficient (r2) between SNP pairs (R, genetics package, with the function LD()) https://0-doi-org.brum.beds.ac.uk/10.1086/381000 [1].

Results

In this study, we investigated a total of 1671 IBD patients of the SIBDCS. 981 (58.7%) patients were classified as CD and 690 (41.3%) as UC patients.

Distribution of NLRP3 SNP genotypes in CD and UC patients

The distribution of the alleles for rs4353135, rs55646866 and rs10733113713 and the joint distribution of alleles in CD and UC patients are presented in Tables 1 and 2. There were no significant differences in distribution of the genotypes between CD and UC patients. The major allele frequency in SIBDCS patients was 0.68 (T) for rs4353135, 0.89 (C) for rs55646866, and 0.74 (G) for rs10733113713, corresponding well with the major allele frequency in the European population. For an overview of the major allele frequencies in different populations see Additional file 2: Table S1. 422 (43.0%, CD) and 285 (41.3%, UC) patients are homozygous carriers for the major allele of all three SNPs. As CD and UC are distinct forms of IBD with diverging disease pathogenesis, we performed individual analyses for associations between NLRP3 SNPs and disease characteristics in CD and UC patients.

Table 1 Cross table of joint SNP distribution in CD patients
Table 2 Cross table of joint SNP distribution in UC patients

Presence of the major allele of all three SNPs is associated with lower maximal CDAI in CD patients

Presence of the major allele (G) of rs10733113 was associated with a significantly lower maximal CDAI value throughout follow-up (median CDAI = 63 in GG or AG vs. 76 in AA patients, p = 0.011, and we observed similar trends, yet not significant, for lower CDAI in carriers of the major alleles of rs4353135 and rs55646866 (Fig. 1A). The investigated SNPs have been reported to be linked due to close proximity of their loci [1] (they all lie within a 5321 bp region, for details see Additional file 1: Figure S1). In our cohort, however, the square of the correlation coefficient r2 (Additional file 2: Table S2) was below 0.2, suggesting no relevant linkage of the studied SNPs. To analyze whether there is an additive effect on the association between NLRP3 SNPs and clinical parameters, we calculated a score for the cumulative genetic risk for carrying the major allele of all 3 polymorphisms. To calculate this score we added the log of the odds ratio (as identified by Villani et al. [1]) for the CD associated allele of the given SNP (log (1.78) = 0.58 for rs10733113, log (1.21) = 0.19 for rs43353135, log (1.69) = 0.52 for rs55646866) resulting in a score between 0 and 2.58 for each patient. In a linear regression model with this genetic score as predictor and log (max. CDAI) as response we observed a significant negative association (Coefficient − 0.097; 95% CI − 0.182 to − 0.011; p = 0.026), indicating that carrying the CD-associated alleles for the 3 investigated NLRP3 SNPs is associated with a less severe course of disease as evaluated by disease symptoms (Fig. 1B), and that there might be an additive effect of the three SNPs.

Fig. 1
figure 1

Effect of the major alleles of rs10733113, rs4353135 and rs55646866 on clinical parameters in CD patients. The graphs show median Log (max. CDAI) values for maximal CDAI within the 25% and 75% percentile (box borders) and outliers (dots) in CD patients A homozygous for the major allele, heterozygous or homozygous for the minor allele for rs10733113, rs4353135 and rs55646866, B depending on the cumulative genetic risk score for carrying the major allele for rs1073313, rs4353135 and rs55646866. In a linear regression model, a significant negative correlation between the cumulative genetic score and log (max. CDAI) can be observed (Coefficient − 0.097; 95% CI − 0.182 to − 0.011; p = 0.026) C CD duration in months within the 25% and 75% percentile (box borders) and outliers (dots) and D percentage of patients with stenosis that required an operation for treatment of stenosis in CD patients homozygous for the major allele, heterozygous or homozygous for the minor allele for rs10733113, rs4353135 and rs55646866

Carrying the major allele for rs55646844 is associated with longer disease duration in CD patients

Duration of disease, the time from diagnosis until last patient visit (i.e., the most recent patient visit, either at the regular yearly visit for the SIBDCS or the latest disease associated visit) 

was significantly longer in patients carrying the major allele for rs55646866 (12.3 years for CC or CT vs. 6.1 years for TT, p = 0.013). For rs4353135 a similar trend, yet not significant, could be observed (12.4 years for GT or TT vs. 9.9 years for GG, p = 0.090) while there was no effect for rs10733113. Yet, age at diagnosis was not affected by the genotype for rs55646844 (Fig. 1C).

The major allele (GG) for rs10733113 is associated with fewer operations due to stenosis in CD patients

While the overall incidence of stenosis did not differ between genotypes, homozygous carriers of the major allele (GG) for rs10733113 had significantly less operations for treatment of stenosis (153 (49.4%) in GG, 0 (0%) in AG, 68 (58.1%) in AA, p = 0.045). The genotype of the SNPs rs4353235 and rs55646866 was not associated to the frequency of operations needed for treatment of stenosis (Fig. 1D).

NLRP3 SNPs are not associated with further clinical disease characteristics in CD patients.

The major alleles of rs10733113, rs4353235 and rs55646866 were associated with a lower maximal CDAI value, suggesting a less severe course of disease in carriers of the major alleles. We analyzed whether individual clinical characteristics e.g. flare frequency, extraintestinal manifestations or medication are associated with carrying the major allele of the investigated polymorphisms. No further associations between genotype for rs10733113, rs4353235, rs55646866 and disease characteristics were observed in CD patients. Tables 3, 4 and 5 show the demographic distribution and clinical characteristics investigated for the different genotypes of rs10733113, rs4353235 and rs55646866.

Table 3 Demographic and clinical characteristics for the different genotypes of rs10733113 in CD patients
Table 4 Demographic and clinical characteristics for the different genotypes of rs4353135 in CD patients
Table 5 Demographic and clinical characteristics for the different genotypes of rs55646866 in CD patients

Carrying the major alleles for rs10733113, rs4353235 and rs55646866 is associated with a higher number of ambulatory hospitalizations in UC patients

Additional ambulatory visits are strong indicators for occurrence of a moderate to severe flare in UC patients. Homozygous genotype for the major allele of all three polymorphisms was associated with a higher number of ambulatory visits for flare management in UC patients (Fig. 2A). Despite their close proximity and earlier reports of linkage disequilibrium, the investigated polymorphisms showed low linkage disequilibrium in the SIBDCS (Additional file 2: Table S2) and we observed a significant association between the number of risk-associated alleles and the number of visits. We applied the same genetic score as calculated for the cumulative effect of all three SNPs on CDAI in CD patients for analysis of ambulatory flare management in UC patients. Based on a logistic regression, we find that this score is significantly associated with the occurrence of at least one ambulatory flare measure (OR: 2.081 (95% CI 1.393–3.111; p < 0.001)) (Fig. 2B).

Fig. 2
figure 2

Individual and cumulative effect of the major alleles of rs10733113, rs4353135 and rs55646866 on ambulatory flare management in UC patients. The graphs show the percentage of UC patients with at least one ambulatory visit for flare management A Homozygous for the major allele, heterozygous or homozygous for the minor allele for rs10733113, rs4353135 and rs55646866. B Depending on the cumulative genetic risk score for carrying the major allele for rs10733113, rs4353135 and rs55646866. In a linear regression model a significant association between the cumulative genetic score and the percentage of patients needing ambulatory flare management could be observed (Odds ratio 2.081; 95% CI − 1.393 to 3.11; p < 0.001). C Age at disease onset is shown within the 25% and 75% percentile (box borders) and outliers (dots) in UC patients.

Homozygous genotype for the major allele (CC) of rs55646866 is associated with a higher age at disease onset and a higher MTWAI index

The MTWAI index characterizes the severity of disease as defined by symptoms in UC patients. Homozygous genotype (CC) for rs55646866 was associated with a higher MTWAI index in UC patients (mean MTWAI: 5 (CC), 5 (CT), 1.5 (TT), p = 0.009). We further observed that homozygous genotype for the major allele of rs55646866 was associated with a higher age at diagnosis (mean age at diagnosis: 28 (CC), 30 (CT), 18.3 (TT), p = 0.004) (Fig. 2C).

NLRP3 SNPs are not associated with further clinical disease characteristics in UC patients.

The major alleles of rs10733113, rs4353235 and rs55646866 were associated with a higher number of ambulatory hospitalizations for flare management, suggesting a more severe course of disease in carriers of the major alleles. We therefore investigated whether individual clinical characteristics, e.g. flare frequency, extraintestinal manifestations or medication, that are usually affected during a more severe disease course, are associated with carrying the major allele for rs10733113, rs4353235 and rs55646866. No further associations between homozygous genotype for the major allele of rs10733113, rs4353235, rs55646866 and disease characteristics were observed in UC patients. Tables 6, 7 and 8 show the demographic distribution and clinical characteristics investigated for the different genotypes of rs10733113, rs4353235 and rs55646866.

Table 6 Demographic and clinical characteristics for the different genotypes of rs10733113 in UC patients
Table 7 Demographic and clinical characteristics for the different genotypes of rs4353135 in UC patients
Table 8 Demographic and clinical characteristics for the different genotypes of rs55646866 in UC patients

Discussion

Polymorphisms in the regulatory region of NLRP3 have been associated with an increased risk to develop CD, yet the impact of these variants on established IBD has not been clarified, so far. We here analyzed for possible associations between the SNPs rs10733113, rs4353135 and rs55646866 and the clinical characteristics of patients from the SIBDCS. Homozygous genotype for the major allele for rs10733113 was associated with a lower CDAI in CD patients. The same trend could be observed for rs4353135 and rs55646866. The SNPs rs4353135, rs55646866 and rs10733113 are in close proximity in the downstream region of NLRP3. Possible linkage disequilibrium between SNP pairs has been reported previously [1]: while D’ values indicated linkage disequilibrium for all SNP pairs, r2 showed very low linkage for 2 of the 3 SNP pairs. However, our calculations of r2 values with SIBDCS data showed no indication of substantial linkage disequilibrium among the three studied SNPs. Furthermore, linear regression models showed log-additive allele effects on maximal CDAI, suggesting independent effects of SNPs’ alleles Similarly, additive effects and a better diagnostic value of a genetic risk score in comparison to individual SNPs have recently been reported by Cleynen et al. [18]. As the CDAI describes the disease activity in CD patients, our findings suggest that the investigated CD-associated NLRP3 variants are associated with a less severe course of disease. In addition, homozygous genotype for the major allele for rs10733113 was associated with fewer operations due to stenosis supporting our finding of an overall less severe course of disease in patients of this genotype.

The investigated polymorphisms in the regulatory region of NLRP3 had been associated to CD, only, but not to UC [1]. Still, as similarly elevated IL-1β levels have been reported for UC patients, suggesting that alterations in NLRP3 inflammasome activity and IL-1β levels might also affect clinical disease characteristics in UC patients, UC patients were also included in our analysis. In contrast to CD, UC patients homozygous for the major allele of the studied NLRP3 polymorphisms seem to have a more severe course of disease. Homozygous genotype for all three polymorphisms was significantly associated with more days of hospitalization indicating most likely an acute flare. Again, the association with all three SNPs could not be explained by linkage disequilibrium, as we observed an additive risk for the genetic risk score for the number of ambulatory visits and did not find r2 values indicative for linkage disequilibrium of these SNPs in the SIBDCS. In addition, homozygous genotype (GG) for the major allele of rs55646866 was associated with a higher MTWAI. As the MTWAI characterizes disease severity in UC patients, this finding further suggests a more severe course of disease in UC patients with GG genotype. In contrast, this genotype was associated with a higher age at diagnosis of disease. Higher age at diagnosis might suggest that the GG variant of rs55646866 is protective against development of UC, but we observed no further results that would support such a finding.

Further clinical characteristics, e.g. other parameters indicative of more severe disease like type of medication or need of treatment intensification, effects on CRP or calprotectin levels or nutrient deficiencies, were not affected by NLRP3 variants, neither in CD nor in UC patients.

Genetic determination of high interindividual variation of in vitro IL-1β secretion and association with disease course in CD and UC patients has been previously reported [19, 20]. So far, little is known about the functional consequences of the investigated SNPs in the regulatory region of NLRP3. NLRP3 mRNA expression was lower in PBMCs isolated from carriers of the major alleles of rs4353135 and rs10733113 in the CD-cohorts studied by Villani et al. and in patients with myocardial infarction and healthy controls [1, 13]. IL-1β secretion from isolated PBMCs was lower in patients with homozygous genotype for the major allele (GG) of SNP rs6672995 [1] and the same trend could be observed for homozygous genotype for the major alleles of rs55646866 and rs10733113. Homozygous genotype for the minor allele of rs4353135 was further associated with a higher risk for oligoarticular and polyarticular juvenile idiopathic arthritis and higher levels of inflammatory markers in a Taiwanese population [21], while presence of the G-allele of rs10733113 was associated with early onset disease in psoriatic arthritis [22].

These findings suggest that the NLRP3 polymorphisms that are associated with CD might lead to lower NLRP3-inflammasome activity and lower IL-1β secretion. In line with the proposed pro-inflammatory role for NLRP3 and IL-1β, such lower NLRP3 activity and IL-1β secretion, might indeed be causally related to a less severe disease course resulting in a lower CDAI in CD patients homozygous for the major allele for these SNPs. As we observed a cumulative effect of the investigated NLRP3 variants on CDAI it would be of particular interest to investigate whether being homozygous for the major allele for more than one of the investigated SNPs has a cumulative effect on IL-1β secretion, thus that patients homozygous for all three major alleles would have the lowest IL-1β levels. Furthermore, polymorphisms in the IL1B gene have been reported to affect course and severity of IBD [23] and nucleotide-binding domain and leucine-rich repeat caspase recruitment domain 4 (NLRC4) and NLRP3 can be recruited to the same macromolecular inflammasome complex [24]. Although no associations between NLCR4 or CASP1 polymorphisms and IBD have been reported so far, it would be very interesting to investigate whether IL-1β levels are affected by cumulative effects of polymorphisms in IL1B itself and further genes involved in IL-1β processing. Unfortunately, such an analysis is beyond the scope of our study.

CD and UC are distinct clinical entities with differences in disease location, histology and only partly shared genetic susceptibility [25]. Still, the diverging effect of CD-associated NLRP3 polymorphisms on disease characteristics in CD patients and UC patients is striking, as similarly elevated secretion levels of IL-1β have been observed in CD and UC patients [4,5,6, 8, 26]. The diverging role of NLRP3 polymorphisms on clinical parameters in CD and UC patients remains unexplained and the role of NLRP3 and IL-1β in IBD is still elusive, as reviewed in detail recently [27]. Although consistently elevated levels of IL-1β secretion have been reported for both CD and UC patients, IL-1β does not seem to play a decisive role in the inflammatory process of IBD patients, as only patients with IL-10R deficiency responded to IL-1R-antagonist treatment [27]. Data from murine models of experimental colitis with animals lacking NLRP3 show diverging effects depending on the model, experimental conditions and environment. Both, aggravated as well as less severe course of disease have been reported for Nlrp3−/− mice during acute DSS and TNBS colitis [9, 10, 12, 28, 29]. Itani et al. describe a more severe course of disease of Nlrp3−/− mice during oxazolone colitis [11], while Mak’Anyengo et al. observed reduced colitis in Nlrp3−/− with the T-cell transfer colitis model [30]. These findings may be explained by the different roles of IL-β during intestinal inflammation: Mak’Anyengo et al. showed NLRP3 dependent IL-1β promoted Th17 differentiation and GM-CSF production of T-cells acting pro-inflammatory and resulting in aggravated colitis. Bersudsky et al. on the other hand, demonstrated that not IL-1β but IL-1α is the decisive pro-inflammatory mediator in the DSS colitis model, whereas IL-1β supports epithelial cell proliferation and restoration of colon barrier [31].

Data from patients or animals with NLRP3 hyper-activity provide a similar controversial picture for the role of NLRP3 in colitis. Both aggravation as well as amelioration of colitis have been shown in patients and mice with various reasons of NLRP3 hyperactivation. Lack of negative regulation by IL-10 leads to NLRP3 overactivation in Il10−/− mice that develop spontaneous colitis [32, 33]. In this model the pro-inflammatory role of IL-1β seems to be non-redundant as treatment with a CASP1 or a NLRP3 inhibitor ameliorated established disease and prevented or delayed development of spontaneous colitis. Similarly, lack of negative regulation of NLRP3 by mutated CARD8 was found to be responsible for development of CD in three related patients, as disease was responsive to IL-1β inhibition [34].

On the contrary, situations with protective effects of NLRP3 hyperactivation have been observed in mice carrying an activating, autoimmune-associated mutation of Ptpn22 [29]. As NLRP3 is activated by PTPN22-dephosphorylation, the Ptpn619W variant leads to increased NLRP3 activation and IL-1β secretion. Mice carrying this mutation were protected from DSS colitis. Likewise, mice carrying the same activating mutation (Nlrp3R258W) that is responsible for increased NLRP3 activation and IL-1β levels in patients with cryoporin associated periodic syndrome (CAPS) were resistant to DSS colitis and T cell transfer colitis [35]. In these animals NLRP3 hyperactivation led to remodelling of the gut microbiota and increased production of regulatory T cells and antimicrobial peptides.

As gut microbiota alterations have been reported for animal models of both lack [10] or hyperactivation [35] of NLRP3, the question arises, whether NLRP3 variant induced modulation of the gut microbiota might be involved in the effect on IBD disease course. Yet, Yilmaz et al. [36] performed an extensive analysis of microbiota disturbances in the SIBDCS, but did not report effects of polymorphisms and we are not aware of any other study investigating the impact of NLRP3 SNPs on the gut microbiota.

So far, it is not clear whether polymorphisms that affect IBD risk also contribute to disease course and severity. While O’Donnell et al. [37] report association of 8 SNPs with time to-abdominal surgery, Lee et al. [38] could not identify any associations of disease susceptibility loci with disease prognosis, but identified four loci with significant association to the course of CD without any association to disease susceptibility. In addition, TLR SNPs that do not affect disease susceptibility, were found to be associated with extensive colonic disease in both CD and UC patients [39], and polymorphisms associated with smoking behavior could predict the number of surgeries and the time to first intestinal surgery in smoking Crohn’s disease patients [40].

Studies addressing the association of genetic loci to disease prognosis included a smaller number of patients in contrast to large genome wide association studies addressing IBD risk. We here studied 981 CD and 690 UC patients leading to significant associations, not robust to Bonferroni correction. Therefore, the identified associations between NLRP3 polymorphisms and IBD prognosis need to be confirmed in additional cohorts.

Based on the currently available (sparse) data, that homozygous genotype for the major allele of CD associated NLRP3 variants leads to lower NLRP3 mRNA levels and activity resulting in lower IL-1β secretion, our findings suggest that during established CD the pro-inflammatory role of IL-1β might dominate, thus that patients homozygous for the major allele of these SNPs might profit from a less severe course of disease. In UC patients on the other hand, our findings support a protective role of NLRP3 such that homozygous genotype for the major alleles of the investigated NLRP3 polymorphisms might favour a more severe course of disease.

The fact, that the major alleles of NLRP3 polymorphisms are the CD associated alleles implies that a considerable number of CD patients—in our study 43%—might profit from the beneficial effect on disease severity. On the other hand, in CD patients homozygous for the minor allele, potentially resulting in a more active NLRP3 inflammasome, newly available NLRP3 inhibitors, e.g. MCC950 or Cy-09 [41, 42] might present a novel treatment option. In contrast, based on our findings that the major alleles of CD-associated NLRP3 variants were associated with more days of hospitalization in UC patients, treatment with NLRP3 inhibitors might be detrimental in UC patients. Nevertheless, a detailed analysis of the effects of NLRP3 variants on NLRP3-inflammasome activity and on IL-1β cytokine levels in IBD patients is necessary, before considering NLRP3 variant data for treatment decisions.

Conclusions

In summary, we here show that SNPs in the regulatory region of the NLRP3 gene have diverging effects on the course of disease in CD and UC patients. Our findings should be taken into account, if treatment with novel NLRP3 inhibitors is considered. Yet, further confirmation of the impact of NLRP3 polymorphisms on IBD disease course in additional and/or larger cohorts and a more detailed analysis of the consequences of NLRP3 variants is necessary.

Availability of data and materials

Data may be made available upon request due to ethical restrictions imposed by the Swiss IBD cohort study group which collected this data. Interested researchers may request access to the collected data from this study in the same manner as the authors did. Data is available upon ethical approval and a request to the Head of the cohort, Dr. Gerhard Rogler, at gerhard.rogler@usz.ch. Data access requests may also be made to the Swiss IBD Cohort at sibdcs-submission@chuv.ch or http://ibdcohort.ch/index.php/informationen-fuer-forscher.html.

Change history

  • 11 January 2022

    In the original HTML version of this article, the given and family names of Diana Bakker de Jong, Ulrike Graf Bigler, Carolina Khalid-de Bakker, Christina Knellwolf(-Grieger) and Sybille Schmid(-Uebelhart) were incorrectly structured. The article has been corrected.

Abbreviations

NLRP3:

NOD-like receptor pyrin domain containing 3

SNP:

Single nucleotide polymorphism

SIBDCS:

Swiss IBD cohort study

CD:

Crohn’s disease

UC:

Ulcerative colitis

PAMP:

Pathogen associated molecular pattern

DAMP:

Danger associated molecular pattern

NOD:

Nucleotide binding oligomerization domain

LRR:

Leucine rich repeat

NLR:

NOD-like receptor

CDAI:

Crohn’s disease activity index

MTWAI:

Truelove and Witts activity index

NLRC4:

Nucleotide-binding domain and leucine-rich repeat caspase recruitment domain 4

References

  1. Villani AC, Lemire M, Fortin G, et al. Common variants in the NLRP3 region contribute to Crohn’s disease susceptibility. Nat Genet. 2009;41:71–6.

    Article  CAS  PubMed  Google Scholar 

  2. Guo H, Callaway JB, Ting JP. Inflammasomes: mechanism of action, role in disease, and therapeutics. Nat Med. 2015;21:677–87.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Sartor RB. Cytokines in intestinal inflammation: pathophysiological and clinical considerations. Gastroenterology. 1994;106:533–9.

    Article  CAS  PubMed  Google Scholar 

  4. Grottrup-Wolfers E, Moeller J, Karbach U, et al. Elevated cell-associated levels of interleukin 1beta and interleukin 6 in inflamed mucosa of inflammatory bowel disease. Eur J Clin Invest. 1996;26:115–22.

    Article  CAS  PubMed  Google Scholar 

  5. Mahida YR, Wu K, Jewell DP. Enhanced production of interleukin 1-beta by mononuclear cells isolated from mucosa with active ulcerative colitis of Crohn’s disease. Gut. 1989;30:835–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Reinecker HC, Steffen M, Witthoeft T, et al. Enhanced secretion of tumour necrosis factor-alpha, IL-6, and IL-1 beta by isolated lamina propria mononuclear cells from patients with ulcerative colitis and Crohn’s disease. Clin Exp Immunol. 1993;94:174–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Guimbaud R, Bertrand V, Chauvelot-Moachon L, et al. Network of inflammatory cytokines and correlation with disease activity in ulcerative colitis. Am J Gastroenterol. 1998;93:2397–404.

    Article  CAS  PubMed  Google Scholar 

  8. Schreiber S, Nikolaus S, Hampe J, et al. Tumour necrosis factor alpha and interleukin 1beta in relapse of Crohn’s disease. Lancet. 1999;353:459–61.

    Article  CAS  PubMed  Google Scholar 

  9. Bauer C, Duewell P, Mayer C, et al. Colitis induced in mice with dextran sulfate sodium (DSS) is mediated by the NLRP3 inflammasome. Gut. 2010;59:1192–9.

    Article  CAS  PubMed  Google Scholar 

  10. Hirota SA, Ng J, Lueng A, et al. NLRP3 inflammasome plays a key role in the regulation of intestinal homeostasis. Inflamm Bowel Dis. 2011;17:1359–72.

    Article  PubMed  Google Scholar 

  11. Itani S, Watanabe T, Nadatani Y, et al. NLRP3 inflammasome has a protective effect against oxazolone-induced colitis: a possible role in ulcerative colitis. Sci Rep. 2016;6:39075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Zaki MH, Boyd KL, Vogel P, et al. The NLRP3 inflammasome protects against loss of epithelial integrity and mortality during experimental colitis. Immunity. 2010;32:379–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Paramel Varghese G, Folkersen L, Strawbridge RJ, et al. NLRP3 inflammasome expression and activation in human atherosclerosis. J Am Heart Assoc. 2016;5:e003031.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Pittet V, Juillerat P, Mottet C, et al. Cohort profile: the Swiss Inflammatory Bowel Disease Cohort Study (SIBDCS). Int J Epidemiol. 2009;38:922–31.

    Article  PubMed  Google Scholar 

  15. Pittet V, Michetti P, Mueller C, et al. Cohort profile update: the Swiss inflammatory bowel disease cohort study (SIBDCS). Int J Epidemiol. 2019;48:385–386f.

    Article  PubMed  Google Scholar 

  16. Jostins L, Ripke S, Weersma RK, et al. Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature. 2012;491:119–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Storm N, Darnhofer-Patel B, van den Boom D, et al. MALDI-TOF mass spectrometry-based SNP genotyping. Methods Mol Biol. 2003;212:241–62.

    CAS  PubMed  Google Scholar 

  18. Cleynen I, Boucher G, Jostins L, et al. Inherited determinants of Crohn’s disease and ulcerative colitis phenotypes: a genetic association study. Lancet. 2016;387:156–67.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Endres S, Cannon JG, Ghorbani R, et al. In vitro production of IL 1 beta, IL 1 alpha, TNF and IL2 in healthy subjects: distribution, effect of cyclooxygenase inhibition and evidence of independent gene regulation. Eur J Immunol. 1989;19:2327–33.

    Article  CAS  PubMed  Google Scholar 

  20. Nemetz A, Nosti-Escanilla MP, Molnar T, et al. IL1B gene polymorphisms influence the course and severity of inflammatory bowel disease. Immunogenetics. 1999;49:527–31.

    Article  CAS  PubMed  Google Scholar 

  21. Yang CA, Huang ST, Chiang BL. Association of NLRP3 and CARD8 genetic polymorphisms with juvenile idiopathic arthritis in a Taiwanese population. Scand J Rheumatol. 2014;43:146–52.

    Article  PubMed  Google Scholar 

  22. Juneblad K, Kastbom A, Johansson L, et al. Association between inflammasome-related polymorphisms and psoriatic arthritis. Scand J Rheumatol. 2021;50:206–12.

    Article  CAS  PubMed  Google Scholar 

  23. Nemetz A, Kope A, Molnar T, et al. Significant differences in the interleukin-1beta and interleukin-1 receptor antagonist gene polymorphisms in a Hungarian population with inflammatory bowel disease. Scand J Gastroenterol. 1999;34:175–9.

    Article  CAS  PubMed  Google Scholar 

  24. Man SM, Hopkins LJ, Nugent E, et al. Inflammasome activation causes dual recruitment of NLRC4 and NLRP3 to the same macromolecular complex. Proc Natl Acad Sci U S A. 2014;111:7403–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Khor B, Gardet A, Xavier RJ. Genetics and pathogenesis of inflammatory bowel disease. Nature. 2011;474:307–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Lazaridis LD, Pistiki A, Giamarellos-Bourboulis EJ, et al. Activation of NLRP3 inflammasome in inflammatory bowel disease: differences between Crohn’s disease and ulcerative colitis. Dig Dis Sci. 2017;62:2348–56.

    Article  CAS  PubMed  Google Scholar 

  27. Mao L, Kitani A, Strober W, et al. The role of NLRP3 and IL-1beta in the pathogenesis of inflammatory bowel disease. Front Immunol. 2018;9:2566.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Bauer C, Duewell P, Lehr HA, et al. Protective and aggravating effects of Nlrp3 inflammasome activation in IBD models: influence of genetic and environmental factors. Dig Dis. 2012;30(Suppl 1):82–90.

    Article  PubMed  Google Scholar 

  29. Spalinger MR, Kasper S, Gottier C, et al. NLRP3 tyrosine phosphorylation is controlled by protein tyrosine phosphatase PTPN22. J Clin Invest. 2016;126:4388.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Mak'Anyengo R, Duewell P, Reichl C, et al. Nlrp3-dependent IL-1beta inhibits CD103+ dendritic cell differentiation in the gut. JCI Insight 2018;3.

  31. Bersudsky M, Luski L, Fishman D, et al. Non-redundant properties of IL-1alpha and IL-1beta during acute colon inflammation in mice. Gut. 2014;63:598–609.

    Article  CAS  PubMed  Google Scholar 

  32. Liu L, Dong Y, Ye M, et al. The pathogenic role of NLRP3 inflammasome activation in inflammatory bowel diseases of both mice and humans. J Crohns Colitis. 2017;11:737–50.

    PubMed  PubMed Central  Google Scholar 

  33. Zhang J, Fu S, Sun S, et al. Inflammasome activation has an important role in the development of spontaneous colitis. Mucosal Immunol. 2014;7:1139–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Mao L, Kitani A, Similuk M, et al. Loss-of-function CARD8 mutation causes NLRP3 inflammasome activation and Crohn’s disease. J Clin Invest. 2018;128:1793–806.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Yao X, Zhang C, Xing Y, et al. Remodelling of the gut microbiota by hyperactive NLRP3 induces regulatory T cells to maintain homeostasis. Nat Commun. 2017;8:1896.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Yilmaz B, Juillerat P, Oyas O, et al. Microbial network disturbances in relapsing refractory Crohn’s disease. Nat Med. 2019;25:323–36.

    Article  CAS  PubMed  Google Scholar 

  37. O’Donnell S, Borowski K, Espin-Garcia O, et al. The Unsolved link of genetic markers and Crohn’s disease progression: a North American cohort experience. Inflamm Bowel Dis. 2019;25:1541–9.

    Article  PubMed  Google Scholar 

  38. Lee JC, Biasci D, Roberts R, et al. Genome-wide association study identifies distinct genetic contributions to prognosis and susceptibility in Crohn’s disease. Nat Genet. 2017;49:262–8.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Pierik M, Joossens S, Van Steen K, et al. Toll-like receptor-1, -2, and -6 polymorphisms influence disease extension in inflammatory bowel diseases. Inflamm Bowel Dis. 2006;12:1–8.

    Article  PubMed  Google Scholar 

  40. Lang BM, Biedermann L, van Haaften WT, et al. Genetic polymorphisms associated with smoking behaviour predict the risk of surgery in patients with Crohn’s disease. Aliment Pharmacol Ther. 2018;47:55–66.

    Article  CAS  PubMed  Google Scholar 

  41. Coll RC, Robertson AA, Chae JJ, et al. A small-molecule inhibitor of the NLRP3 inflammasome for the treatment of inflammatory diseases. Nat Med. 2015;21:248–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Jiang H, He H, Chen Y, et al. Identification of a selective and direct NLRP3 inhibitor to treat inflammatory disorders. J Exp Med. 2017;214:3219–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank all the patients for their collaboration and the members of the Swiss Inflammatory Bowel Disease Cohort Study for their contribution. Members of the Swiss IBD Cohort Study Group: Karim Abdelrahman, Gentiana Ademi, Patrick Aepli, Amman Thomas, Claudia Anderegg, Anca-Teodora Antonino, Eva Archanioti, Eviano Arrigoni, Diana Bakker de Jong, Bruno Balsiger, Polat Bastürk, Peter Bauerfeind, Andrea Becocci, Dominique Belli, José M. Bengoa, Luc Biedermann, Janek Binek, Mirjam Blattmann, Stephan Boehm, Tujana Boldanova, Jan Borovicka, Christian P. Braegger, Stephan Brand, Lukas Brügger, Simon Brunner, Patrick Bühr, Bernard Burnand, Sabine Burk, Emanuel Burri, Sophie Buyse, Dahlia-Thao Cao, Ove Carstens, Dominique H. Criblez, Sophie Cunningham, Fabrizia D’Angelo, Philippe de Saussure, Lukas Degen, Joakim Delarive, Christopher Doerig, Barbara Dora, Susan Drerup, Mara Egger, Ali El-Wafa, Matthias Engelmann, Jessica Ezri, Christian Felley, Markus Fliegner, Nicolas Fournier, Montserrat Fraga, Yannick Franc, Pascal Frei, Remus Frei, Michael Fried, Florian Froehlich, Raoul Ivano Furlano, Luca Garzoni, Martin Geyer, Laurent Girard, Marc Girardin, Delphine Golay, Ignaz Good, Ulrike Graf Bigler, Beat Gysi, Johannes Haarer, Marcel Halama, Janine Haldemann, Pius Heer, Benjamin Heimgartner, Beat Helbling, Peter Hengstler, Denise Herzog, Cyrill Hess, Roxane Hessler, Klaas Heyland, Thomas Hinterleitner, Claudia Hirschi, Petr Hruz, Pascal Juillerat, Carolina Khalid-de Bakker, Stephan Kayser, Céline Keller, Christina Knellwolf(-Grieger), Christoph Knoblauch, Henrik Köhler, Rebekka Koller, Claudia Krieger(-Grübel), Patrizia Künzler, Rachel Kusche, Frank Serge Lehmann, Andrew Macpherson, Michel H. Maillard, Michael Manz, Astrid Marot, Rémy Meier, Christa Meyenberger, Pamela Meyer, Pierre Michetti, Benjamin Misselwitz, Patrick Mosler, Christian Mottet, Christoph Müller, Beat Müllhaupt, Leilla Musso, Michaela Neagu, Cristina Nichita, Jan Niess, Andreas Nydegger, Nicole Obialo, Diana Ollo, Cassandra Oropesa, Ulrich Peter, Daniel Peternac, Laetitia Marie Petit, Valérie Pittet, Daniel Pohl, Marc Porzner, Claudia Preissler, Nadia Raschle, Ronald Rentsch, Alexandre Restellini, Sophie Restellini, Jean-Pierre Richterich, Frederic Ris, Branislav Risti, Marc Alain Ritz, Gerhard Rogler, Nina Röhrich, Jean-Benoît Rossel, Vanessa Rueger, Monica Rusticeanu, Markus Sagmeister, Gaby Saner, Bernhard Sauter, Mikael Sawatzki, Michael Scharl, Martin Schelling, Susanne Schibli, Hugo Schlauri, Dominique Schluckebier, Daniela Schmid, Sybille Schmid(-Uebelhart), Jean-François Schnegg, Alain Schoepfer, Vivianne Seematter, Frank Seibold, Mariam Seirafi, Gian-Marco Semadeni, Arne Senning, Christiane Sokollik, Joachim Sommer, Johannes Spalinger, Holger Spangenberger, Philippe Stadler, Peter Staub, Dominic Staudenmann, Volker Stenz, Michael Steuerwald, Alex Straumann, Bruno Strebel, Andreas Stulz, Michael Sulz, Aurora Tatu, Michela Tempia-Caliera, Joël Thorens, Kaspar Truninger, Radu Tutuian, Patrick Urfer, Stephan Vavricka, Francesco Viani, Jürg Vögtlin, Roland Von Känel, Dominique Vouillamoz, Rachel Vulliamy, Paul Wiesel, Reiner Wiest, Stefanie Wöhrle, Samuel Zamora, Silvan Zander, Tina Wylie, Jonas Zeitz, Dorothee Zimmermann. Head of the Swiss IBD study group is Prof. Gerhard Rogler; Division of Gastroenterology and Hepatology, University Hospital Zurich, University of Zurich, Zurich, Switzerland; e-mail: Gerhard.rogler@usz.ch

Funding

This research was supported by a research grant from the Swiss National Science Foundation to GR for the Swiss IBD Cohort (Grant No. 3347CO-108792). Funding enabled data acquisition and the genetic analysis of the SIBDC. The analysis at hand was performed by the authors without external funding. The funder did not influence the design of the study, the collection, analysis, and interpretation of data, nor the writing of the manuscript.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

BM, MH, GR, MS and IFW: concept and design of the study; PY, JBR, SBUJ, YF, LB, BM, IFW: acquisition and analysis of the data; JBR, SBUJ, YF and IFW: statistical analysis and visualization, PY, JBR and IFW drafting the article; JBR, SBUJ, YF, LB, BM, MH, GR, MS and IFW interpretation of the data and critical review. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Isabelle Frey-Wagner.

Ethics declarations

Ethics approval and consent to participate

The SIBDCS has been approved by the respective local ethics committees in Switzerland (BASEC-No. 2018-02068, approved on 09.03.2020 by the Cantonal Ethics Committee of the Canton Zurich, Switzerland) conforming to the ethical guidelines of the 1975 Declaration of Helsinki. All patients signed an informed consent and confirmed their participation in the cohort study at the time of enrolment and gave informed consent for data collection and analysis for research purposes. The current substudy has been evaluated and approved by the scientific board of SIBDCS.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests relevant to this study.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yoganathan, P., Rossel, JB., Jordi, S.B.U. et al. Genotype–phenotype associations of polymorphisms within the gene locus of NOD-like receptor pyrin domain containing 3 in Swiss inflammatory bowel disease patients. BMC Gastroenterol 21, 310 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12876-021-01880-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12876-021-01880-9

Keywords