Associations between IL-23R gene polymorphisms and the susceptibility of rheumatoid arthritis: a meta-analysis

Abstract

Objective

This meta-analysis was aimed to verify the influences of interleukin 23 receptor (IL-23R) rs11209026 and rs2201841 polymorphisms on the susceptibility of rheumatoid arthritis (RA).

Methods

We searched potentially relevant studies on the relationships of IL-23R gene rs11209026 and rs2201841 polymorphisms with RA susceptibility from PubMed (Medline), CNKI and EMBASE web databases. Crude odds ratios (ORs) along with 95% confidence intervals (95% CIs) were calculated to assess the strength of association through the fixed- or random-effects model.

Results

Eight published articles containing 5544 RA patients and 5532 health individuals were included in our meta-analysis. In total analysis, we found IL-23R gene rs11209026 polymorphism was not related to RA susceptibility in all genetic models, but there was a significant association between rs2201841 polymorphism and RA susceptibility under G vs. A model. In stratified analysis by ethnicity, an increased RA susceptibility was detected for rs2201841 polymorphism in Caucasian group.

Conclusion

The results of meta-analysis indicated that IL-23R gene rs2201841 polymorphism might be a susceptible factor for RA under G vs. A model, especially in Caucasian population.

Keywords:

• IL-23R
• polymorphism
• RA
• meta-analysis

Introduction

Rheumatoid arthritis (RA) is a chronic autoimmune disease (AID) impacting human health severely [1]. It is characterized by synoviocyte hyperplasia, neovascularization and partial infiltration of numerous T cells [2], inducing synovial inflammation, articular cartilage damages, bone and joint capsules at different levels, and finally leads to joint deformity and functional incapacitation [3]. As one of the primary disabling diseases, the pathogenesis of RA involves in genetic factors and environmental factors [4,5]. Moreover, the genetic factors account for 50–60% of all RA cases [6,7]. However, the precise aetiology of this disease is still beyond totally understood. So far, many RA-related susceptibility genes have been identified, including HLA, PADI4 and PTPN22 [8–10]. Furthermore, other genes are still emerging to be implicated in the incidence of RA.

Interleukin 23 receptor (IL-23R) is produced by antigen-presenting cells (like macrophages and dendritic cells) with the actions of certain infections [11]. It plays a key role not only in inflammatory diseases but also in innate immunity and adaptive immunity [12,13]. In addition, IL-23R also dramatically affects the process of cell transformation from CD4⁺T to Th17, which increases the number of Th17 cells and maintains their activity, thus inducing the generation of proinflammatory cytokine IL-17 [14]. In accordance with relevant experimental models, Th17 lymphocyte subsets which secrete IL-17 may cause the damages in various positions, involving bones, cartilages, joints, brains, lungs, hearts and intestinal tracts [15]. In addition, researches also demonstrate that the variants in IL-23R gene increase the susceptibilities of AIDs in patients with inflammatory bowel disease, psoriasis or multiple sclerosis [16,17]. Sato et al. suggested that IL-23R/IL-17 axis played a vital role not only in the growth stage of autoimmune arthritis but also in osteoclasia stage, indicating that IL-23R may correlate with the susceptibility of RA [18].

Some researchers have reported the relationship between IL-23R gene polymorphisms and RA risk but got no consistent result. Therefore, we searched all eligible articles on the role of IL-23R gene rs11209026 and rs2201841 polymorphisms in individual susceptibility to RA so as to statistically discuss such relationship for a more reliable conclusion.

Materials and methods

Search strategy and selection criteria

Relevant publications were identified through searching PubMed (Medline), China National Knowledge Infrastructure (CNKI) and EMBASE web databases using the following search terms: “interleukin 23 receptor or IL-23R” and “polymorphism or mutation or variant” and “rheumatoid arthritis or RA”. What is more, we also checked the references of relevant articles to find any other potentially relevant papers. Studies for this meta-analysis were selected based on the following criteria: (1) with a case-control design; (2) involving the associations of IL-23R gene rs11209026 and rs2201841 polymorphisms with the susceptibility of RA; (3) providing adequate genotype data to calculate odds ratios (ORs) with 95% confidence intervals (95% CIs). The articles which did not meet the above-mentioned criteria were excluded.

Data extraction

The information extracted by two independent investigators from each study included first author’s name, year of publication, original country, ethnicity, source of control, genotyping method and genotype and/or allele distributions in cases and controls. Another investigator would be consulted whenever no consensus could not be reached over conflicting results of abstracted data after discussing between the two investigators.

Statistical analysis

Summary ORs and 95% CIs were calculated to evaluate the relationships of IL-23R gene rs11209026 and rs2201841 polymorphisms with the susceptibility of RA under five genetic contrasts: homozygote, dominant, recessive, allele and heterozygote models. Stratified analysis was performed based on ethnicity. Heterogeneity was measured using Q-statistic and I² test, and if p < .05 or I² > 50%, indicating significant heterogeneity, the random-effects model was applied to calculate pooled ORs; otherwise, the fixed-effects model was used. χ² test was employed to measure the deviation of genotype distribution from Hardy–Weinberg equilibrium (HWE) in controls. Publication bias was inspected by funnel plots and Egger’s test. We also performed sensitivity analysis to detect the stability of our results by sequentially omitting each single study. The STATA version 12.0 (STATA Corp, College Station, TX, USA) software was used for all data syntheses in this meta-analysis. Significant threshold was fixed at p < .05 for all tests.

Results

Study characteristics

As shown in Figure 1, 204 articles were initially searched from the databases, and 132 irrelevant publications were excluded through scanning abstracts. Then, additional 57 papers were removed in further screening without controls (17), meta-analysis (6) and not concerning IL-23R gene rs11209026 and rs2201841 polymorphisms (35). Additionally, six articles were deleted due to the lack of genotype data. Therefore, this meta-analysis included eight case-control studies [19–26], containing 5544 RA patients and 5532 controls. Table 1 demonstrates the main characteristics of included studies. The genotype distributions of the controls did not conform to HWE (p < .05) in 4 studies [20,22,23,26].

Figure 1. Flow diagram of article selection process.

Table 1. Principal characteristics of the studies included in the meta-analysis.

First author	Year	Country	Ethnicity	Control source	Genotyping method	Sample size	Case		Control		pHWE
rs11209026							GG	GA + AA	GG	GA + AA
Bogunia-Kubik	2015	Poland	Caucasian	Population	real-time PCR	89/125	83	6	111	14	.440
Chang	2008	North America	Caucasian	Population	multiplexed conventional PCR	471/475	414	57	415	60	.243
Chang	2008	North America	Caucasian	Population	allele-specific PCR	658/1318	563	95	1132	186	.305
Chang	2008	Netherlands	Caucasian	Population	allele-specific PCR	590/700	505	85	621	79	.812
Hamdy	2015	Egypt	African	Population	real-time PCR	120/120	0	120	0	120	.009
Hollis-Moffatt	2008	New Zealand	Caucasian	Population	TaqMan	850/557	749	101	490	67	.920
Orozco	2007	Spain	Caucasian	Population	TaqMan	322/342	283	39	302	40	<.001
Varade	2009	Spain	Caucasian	Population	real-time PCR	550/546	464	86	485	61	.544
rs2201841							AA	AG + GG	AA	AG + GG
Farago	2007	Hungary	Caucasian	Population	PCR-RFLP	412/220	178	234	101	119	.029
Hollis-Moffatt	2008	New Zealand	Caucasian	Population	TaqMan	644/485	294	350	228	257	<.001
Orozco	2007	Spain	Caucasian	Population	TaqMan	322/342	143	179	167	175	.084
Szabo	2013	Hungary	Caucasian	Population	PCR-RFLP	396/182	167	229	82	100	.038
Hamdy	2015	Egypt	African	Population	PCR-RFLP	120/120	104	16	104	16	.434

Notes: PCR: polymerase chain reaction; RFLP: restriction fragment length polymorphism; TaqMan: TaqManSNP; p_HWE: p values for Hardy–Weinberg equilibrium test.

Meta-analysis results

The main results of the meta-analysis are displayed in Table 2. From pooled results, we found no significant association between IL-23R gene rs11209026 and RA susceptibility under all comparison models in total analysis. However, a significant association was observed for IL-23R gene rs2201841 polymorphism in the model G vs. A (OR = 1.18, 95% CI = 1.05–1.32) (Figure 2).

Figure 2. Forest plot of RA susceptibility associated with IL-23R gene rs2201841 polymorphism under the G vs. A model. The squares and horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the weight (inverse of the variance). The diamond represents the summary OR and 95% CI.

Table 2. IL-23R Gene polymorphisms and RA susceptibility.

SNP	Comparison	Group	OR (95% CI)	ph	I^2
rs11209026
	AA vs. GG	Caucasian	1.87 (0.99, 3.56)	.299	17.80%
		Total	1.87 (0.99, 3.56)	.299	17.80%
	AA + GA vs. GG	Caucasian	1.10 (0.96, 1.27)	.343	11.40%
		Total	1.10 (0.96, 1.27)	.343	11.40%
	A vs. G	Caucasian	1.12 (0.99, 1.28)	.441	0.00%
		Total	1.25 (0.95, 1.66)	.000	75.00%
rs2201841
	GG vs. AA	Caucasian	1.84 (1.00, 3.41)	.006	76.00%
		Total	1.84 (1.00, 3.41)	.006	76.00%
	GG + AG vs. AA	Caucasian	1.11 (0.96, 1.29)	.941	0.00%
		Total	1.11 (0.96, 1.28)	.976	0.00%
	G vs. A	Caucasian	1.18 (1.05, 1.32)	.590	0.00%
		Total	1.18 (1.05, 1.32)	.590	0.00%

Note: p_h, p values of heterogeneity test.

In subgroup analysis based on ethnicity, only one article involved African populations, the combined analysis result was not accurate, thus the results were not presented in this study. Increased susceptibility of RA was found for the rs2201841 polymorphism in Caucasian group (allele G vs. allele A: OR = 1.18, 95% CI = 1.05–1.32) (Figure 2).

Heterogeneity test

We observed significant heterogeneity for IL-23R gene rs11209026 (allele A vs. allele G: p = .000) and rs2201841 (GG vs. AA: p = .006) polymorphisms, so the random-effects model was utilized to calculate combined results. The fixed-effects model was used for the rest genetic models in current meta-analysis.

Sensitivity analysis and publication bias

The results of sensitivity analysis showed that no single study impacted the overall ORs qualitatively, confirming the statistical stability of our findings (Figure 3).

Figure 3. Sensitivity analysis for IL-23R gene rs2201841 polymorphism under the allele A vs. allele G model.

None of the funnel plots showed distinct asymmetry (Figure 4), implying the absence of significant publication bias. Besides, all these results were also certified by statistical evidence from Egger’s test (p = .298).

Figure 4. Begg’s funnel plot for publication bias. Each point represents a separate study for the indicated association. Log(OR): natural logarithm of OR; Horizontal line: mean effect size.

Discussion

The pathogenesis of RA is still not completely understood, with an average prevalence rate of approximately 1% worldwide, and there are around 4 million RA patients in China. No specific medicine has been developed for the therapy of RA so far. Since conventional therapeutic approaches can only relieve patients’ pains and inflammations, retard or terminate the damages of joints, improve the health degree and action competence, they cannot radically cure such disease fundamentally. In recent years, the unceasing development of molecular biology has revealed that the mediation of cytokines may have pivotal influences on the generation and growth of joint diseases. The functions of cytokines such as interleukin and tumour necrosis factors are the critical mechanism of inflammatory responses.

IL-23, a new-found heterodimer cell-molecule, is mainly created by activated mononuclear macrophages and dendritic cells [27] and also secreted by microglias in brain [28]. It is composed of specific subunit IL-23p19 and IL-12p40 by disulfide bone, with extremely similar molecular structures as IL-12 [27,29]. IL-23, the cytokine connecting innate and adaptive immunity, can only exert its biological effects through binding to specific IL-23 receptor which has high affinity on the surface of target cells, so the latter is the key to regulate the biological activity of the former. IL-23 receptor consists of two subunits IL-23R and IL-12RB1 and is expressed on the surface of memory T cell principally [30]. IL-23R gene, encoding for IL-23R, is polymorphic and is located on human chromosome 1p31.3. Focusing on the polymorphisms in this gene, scholars have reported their roles in individual susceptibility to RA.

Varade et al. found that the A allele of IL-23R rs11209026 polymorphism was associated with increased predisposition to RA in their studied Spanish [19], and Hamdy et al. suggested a similar effect of the AA genotype of the polymorphism based on Egyptians, though they found no significant relationship for the other polymorphism rs2201841 [25]. Besides, Farago et al. detected a risk-increasing influence of the CC genotype of IL-23R rs2201841 polymorphism on RA in their study [22]. However, Hollis–Moffatt et al. insisted there was no significant relationship of RA risk either with the rs11209026 polymorphism or with the rs2201841 polymorphism [20]. Moreover, Orozco et al. reported that there was no significant difference of allele or genotype frequencies of these two polymorphisms between RA patients and controls, and concluded that they seemed not to be related to RA predisposition in a Spanish population [23].

Considering these controversial findings on the correlation between IL-23R polymorphisms and RA risk, we performed this study to combine them for a clearer perspective on this issue. Consequently, our results demonstrated that there was no significant association between IL-23R gene rs11209026 and RA susceptibility in total analysis, but the rs2201841 polymorphism increased the risk of RA under G vs. A contrast. Additionally, we also observed a similar influence of the two polymorphisms rs11209026 and rs2201841 on RA susceptibility after subgroup analysis according to ethnicity in African and Caucasian groups, respectively.

This meta-analysis had certain comprehensiveness, but there still were some limitations which should be taken into consideration. First, the number of included studies was relatively small, affecting the accuracy of results. Secondly, only published articles were included in our meta-analysis, so other relevant studies might be omitted. Third, gene–gene and loci–loci interactions were not discussed in the present study because of insufficient data. Fourth, only one article involved African populations, thus the subgroup analysis result based on ethnicity was not accurate. Finally, the majority of included studies were based on Caucasians, thus possibly introducing some selective bias.

In summary, our meta-analysis suggested that IL-23R gene rs2201841 polymorphism was related to the increased susceptibility of RA, especially in Caucasians. In future, the correlations of IL-23R gene polymorphisms with RA susceptibility require to be validated in larger studies.

Disclosure statement

No potential conflict of interest was reported by the authors.