eISSN: 1644-4124
ISSN: 1426-3912
Central European Journal of Immunology
Current issue Archive Manuscripts accepted About the journal Special Issues Editorial board Abstracting and indexing Subscription Contact Instructions for authors Publication charge Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
1/2019
vol. 44
 
Share:
Share:
Clinical immunology

Ankylosing spondylitis: analysis of gene-gene interactions between IL-12β, JAK2, and STAT3 in Han Chinese and Algerian cohorts

Ahlem Saadi
1
,
Jie Dang
2
,
Shan Shan
2
,
Aicha Ladjouze-Rezig
3
,
Salima Lefkir-Tafiani
3
,
Yaoqin Gong
2
,
Qiji Liu
2
,
Traki Benhassine
4

  1. Laboratory of Cellular and Molecular Biology, Faculty of Biological Sciences, University of Sciences and Technology Houari Boumediene (USTHB), Algiers, Algeria
  2. Key Laboratory for Experimental Teratology of the Ministry of Education and Department of Medical Genetics, Shandong University School of Medicine, Jinan, China
  3. Rhumatology Service, Etablissement Hospitalier Spécialisé (EHS Mohamed Boukhroufa) Ben Aknoun, Algiers, Algeria
  4. Faculty of Biological Sciences, University of Sciences and Technology Houari Boumediene (USTHB), Algiers, Algeria
(Centr Eur J Immunol 2019; 44 (1): 65-74)
Online publish date: 2019/04/15
Article file
Get citation
 
PlumX metrics:
 

Introduction

Ankylosing spondylitis (AS) is a common inflammatory arthritis characterised by acute and chronic inflammation of the axial skeleton, causing pain, stiffness, and progressive joint ankylosis. AS has been strongly associated with human leukocyte antigen-B27 (HLA-B27), and its contribution to the genetics of AS is best characterised at present in various populations [1-4]. For instance, HLA-B27 is found in more than 80% of Chinese patients with AS [5], while only 1-5% of HLA-B27 carriers develop AS in the Chinese population [6]. However, the reported prevalence of AS ranges between 2% and 5% among major Arab populations [7]. Moreover, association studies on AS in the Maghreb population have mainly focused on the implication of the HLA-B27 gene in disease susceptibility [8-10], thus showing a weaker association between AS and HLA-B27 for Algerian patients as compared to European patients (63% vs. 80-90%, respectively) [8].

Although the contribution of HLA-B27 to the genetics of AS is best characterised at present, the pathogenesis of AS remains poorly understood, and other genetic factors encoded by non-major-histocompatibility-complex (non-MHC) genes may contribute to this disease [11, 12]. Accordingly, many previous association studies tagging a number of non-MHC candidate genes have identified significant associations with AS in the Chinese population [13-17]. Recently, a genome-wide association study in Caucasian European populations from Australia, North America, and the United Kingdom provided the first evidence that interleukin-23 (IL-23) is involved in the pathogenesis of ankylosing spondylitis, and variants in several genes involved in the IL-23 pro-inflammatory cytokine pathway have been shown to be associated with the disease [18].

Interleukin-23 and interleukin-12 are critical cytokines that bridge innate and adaptive immunities, although they have different activities. IL-23 promotes a T-cell population characterised by the production of IL-17, IL-17F, TNF, and IL-6 [19]. IL-12β is on chromosome 5q31 and encodes the IL-12 p40 subunit that heterodimerises with the IL-12 p35 subunit (IL-12A gene) to form IL-12 cytokine, or with the IL-23 p19 subunit (IL-23A gene) to form IL-23 cytokine. IL-12β is important for the differentiation of naive T cells into IFN-γ-producing T-helper 1 cells (Th1), which is essential for the antimicrobial response [20], and this could explain its role in AS and other inflammatory-mediated processes [21]. The functional roles of IL-12β polymorphisms are still unclear, but most of the case-control studies on inflammatory diseases, including ankylosing spondylitis, have analysed one important variant in the IL-12β gene: the 3’ untranslated region (3’UTR) variant rs3212227 (or A1188C) [22-24].

The aim of our work was to confirm whether there is any association between ankylosing spondylitis and four SNPs (single-nucleotide polymorphisms) within three genes (IL-12β, JAK2, and STAT3), all involved in the IL-23 signalling pathway, for two ethnically different populations: Han Chinese and Algerian. In addition, we performed for each population a gene-gene interaction analysis to determine whether these four sequence variants could have any substantial effect increasing the incidence of AS.

Material and methods

Subjects

Two different populations were included in this case-control study, Han Chinese and Algerian. Experiments were carried out separately for the two cohorts in the Institute of Medical Genetics of Shandong University. We recruited 1010 unrelated Chinese individuals (430 patients with ankylosing spondylitis and 580 healthy controls) living in the same urban area and collected from the Department of Rheumatology and Immunology of Qilu Hospital in Shandong, an eastern coastal province in the north of China. We also recruited 250 Algerian individuals (130 unrelated AS patients and 120 healthy controls) from EHS Ben Aknoun, a central hospital in Algiers, Algeria.

Ankylosing spondylitis was diagnosed according to the modified 1984 New York criteria [25]. Magnetic resonance imaging (MRI) and computerised tomography (CT) were used as imaging modalities for sacroiliitis. Diagnosis of AS was confirmed by a qualified rheumatologist. All patients were only affected by ankylosing spondylitis, while those displaying other diseases such as inflammatory bowel disease (IBD) or psoriasis were excluded.

All healthy, unrelated controls recruited for this study were matched for age. Detailed information, clinical characteristics, and original medical institution were obtained for all participants, who gave written informed consent. The study was approved by the relevant Ethics Committee.

SNP selection and genotyping

Four single-nucleotide polymorphisms (SNPs) within three genes were selected for the Han Chinese cohort from the international HapMap project database (http://www.hapmap.org/) with the Chinese Beijing Han population (CHB), namely rs3212227 (A1188C) and rs6887695 for the same gene IL-12β, rs7857730 for JAK2, and rs2293152 for STAT3. All these SNPs are involved in the IL-23 signalling pathway (Fig. 1). These four SNPs were also investigated for the first time in an Algerian cohort, since no previous association studies between AS and the three studied genes were available for the Algerian population.

Fig. 1

Schematic representation of the three candidate genes and the four tag SNPs selected for this study in the IL-23 signalling pathway (SNPs rs3212227 and rs6887695 in IL-12β, SNP rs7857730 in JAK2, and SNP rs2293152 in STAT3)

/f/fulltexts/CEJOI/36206/CEJI-44-84019-g001_min.jpg

Genomic DNA from all Chinese subjects was extracted from peripheral blood mononuclear cells (PBMCs; leukocytes) by a standard phenol-chloroform method and was diluted to 20 mg/ml. QIAamp DNA Mini Kit® was used to extract genomic DNA from peripheral blood cells for all Algerian individuals [26].

Genotyping of the two cohorts involved allelic discrimination by TaqMan® SNP Genotyping Assays with an ABI prism 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA), using four pre-made, double-dye probes (C_2084293_10 for rs3212227), (C_1994992_10 for rs6887695), (C_3140302_10 for rs2293152), and (C_29340600_20 for rs7857730). The endpoint fluorescence was read using the Roche 480 Real-Time PCR System (Hoffmann-La Roche Ltd, Basel City, Switzerland) with 96-well plates. The PCR reaction was run in a total solution volume of 5 µl containing 2.5 µl of 2 × c Premix EX Taq (Takara, Otsu, Shiga, Japan), at least 10 ng of genomic DNA, 2.5 µM of each primer, and 0.05 µl of probe. The PCR amplification conditions were 95°C for 30 sec, 45 cycles of 95°C for 5 sec, 60°C for 20 sec, and 40°C for 30 sec for cooling. Genotyping accuracy in the samples was confirmed by direct sequencing of PCR products for 5% of randomly chosen samples. Only those individuals with 100% genotype success for all markers were included for final analysis.

Statistical analysis

Single nucleotide polymorphisms were assessed for genotypic and allelic association analysis among cases and controls, separately for the Chinese and the Algerian cohorts. The genotypic data of the four SNPs were tested for Hardy-Weinberg equilibrium by comparing genotype frequencies within the two groups (AS patients and controls) by Chi-squared test. The odds ratios (ORs) with 95% confidence intervals (95% CIs) were also calculated for single site analysis to estimate the risk associated with ankylosing spondylitis. Statistical analyses were carried out to calculate p-values using SPSS statistical software (SPSS Statistics 16.0; 2008 SPSS Inc., Chicago, IL).

In order to analyse gene-gene interactions, multiplicative interactions were estimated separately for the two cohorts using SPSS 16.0 software for logistic regression analyses. The multiplicative variable model was defined as disease status (1 is unaffected, 2 is affected) and SNP genotypic values ranging from 0 to 2 indicating the number of risk alleles in an individual subject (0 is zero risk allele, 1 is one risk allele, 2 is two risk alleles). For each SNP, random combinations were constructed to estimate a binary logistic regression model based on the case-control status (dependent variable) and the other SNP variables studied (independent variables or covariates).

Because the effect of any single genetic variation was likely to be dependent on other genetic variations, multifactor dimensionality reduction (MDR) software analysis (http://sourceforge.net/projects/mdr/) was performed to test various interaction models for one-variable, two-variable, three-variable, and four-variable gene-gene epistasis on ankylosing spondylitis. We analysed all data using 10-fold cross-validation.

Results

Association results for the Han Chinese cohort

Genotypic and allelic association analysis of Chinese cases and controls is shown in Table 1. The frequency distributions of the four SNPs within patients and controls were in Hardy-Weinberg equilibrium (p > 0.05). For the Chinese Han cohort, the SNP IL-12β rs3212227 (or A1188C) showed the strongest association with AS (p genotype = 0.007, p allele = 0.005 [OR 0.77, 95% CI: 0.65-0.93]). Another association with AS was also observed for SNP rs7857730 of the JAK2 gene (genotypic p = 0.029 and allelic p = 0.014). The minor allele frequency results revealed that C and T might be protective alleles for rs3212227 and rs7857730, respectively. Thus, subjects carrying A and G alleles showed increased risk of AS for the same variants. However, no significant differences were observed between AS patients and controls for the two other variants (IL-12β rs6887695 and STAT3 rs2293152) within the Han Chinese cohort.

Table 1

Genotypic and allelic association analysis of the four single-nucleotide polymorphisms in the Han Chinese and Algerian ankylosing spondylitis (AS) cohorts

 
Chinese cohort
Algerian cohort
Cytogenetic locationGeneSNPsFunctionGenotype/AlleleNumber (%) AS patientsNumber (%) ControlsOR(95% CI)p-valueGenotype/AlleleNumber (%) AS patientsNumber (%) ControlsOR (95% CI)p-value
chr5q33IL-12βrs3212227
(A1188C)
UTR 3’CC
CA
AA
C
A
54 (0.125)
217 (0.505)
159 (0.370)
325 (0.378)
535 (0.622)
115 (0.198)
280 (0.482)
185 (0.320)
510 (0.440)
650 (0.560)
1.00
0.61 (0.42-0.88)
0.55 (0.37-0.80)
1.00
0.77 (0.65-0.93)
0.007
0.005
CC
CA
AA
C
A
6 (0.046)
54 (0.415)
70 (0.540)
66 (0.254)
194 (0.746)
17 (0.142)
48 (0.400)
55 (0.458)
82 (0.342)
158 (0.658)
1.00
0.31 (0.11-0.86)
0.28 (0.10-0.75)
1.00
0.66 (0.45-0.96)
0.029
0.031
chr5q33IL-12βrs6887695PromoterCC
CG
GG
C
G
60 (0.140)
196 (0.460)
168 (0.400)
316 (0.370)
532 (0.630)
92 (0.160)
272 (0.470)
216 (0.370)
456 (0.393)
704 (0.607)
1.00
0.91 (0.62-1.31)
0.84 (0.57-1.23)
1.00
0.92 (0.76-1.10)
0.650
0.351
CC
CG
GG
C
G
23 (0.177)
60 (0.462)
47 (0.361)
106 (0.408)
154 (0.592)
28 (0.230)
55 (0.460)
37 (0.310)
111 (0.460)
129 (0.540)
1.00
0.75 (0.39-1.46)
0.65 (0.32-1.30)
1.00
0.80 (0.56-1.14)
0.472
0.216
chr9p24JAK-2rs7857730IntronTT
TG
GG
T
G
80 (0.186)
230 (0.535)
120 (0.279)
390 (0.453)
470 (0.547)
147 (0.253)
296 (0.510)
137 (0.237)
590 (0.509)
570 (0.491)
1.00
0.70 (0.51-0.97)
0.62 (0.43-0.90)
1.00
0.99 (0.83-1.18)
0.029
0.014
GG
GT
TT
G
T
21 (0.161)
53 (0.408)
56 (0.431)
95 (0.365)
165 (0.635)
15 (0.125)
63 (0.525)
42 (0.350)
93 (0.388)
147 (0.612)
1.00
1.66 (0.78-3.55)
1.05 (0.48-2.28)
1.00
0.91 (0.63-1.31)
0.176
0.610
chr17q21STAT-3rs2293152IntronGG
GC
CC
G
C
102 (0.237)
201 (0.468)
127 (0.295)
405 (0.470)
455 (0.530)
124 (0.214)
287 (0.495)
169 (0.291)
535 (0.461)
625 (0.539)
1.00
1.17 (0.85-1.61)
1.09 (0.77-1.55)
1.00
1.04 (0.87-1.24)
0.606
0.664
GG
GC
CC
G
C
8 (0.061)
48 (0.370)
74 (0.569)
64 (0.246)
196 (0.754)
18 (0.150)
50 (0.417)
52 (0.433)
86 (0.358)
154 (0.642)
1.00
0.46 (0.18-1.16)
0.31 (0.13-0.77)
1.00
0.58 (0.40-0.86)
0.025
0.006

[i] OR – odds ratio, 95% CI: 95% confidence interval

[ii] p values < 0.05 are indicated in bold.

Association results for the Algerian cohort

Association analysis results of Algerian patients and controls are represented in Table 1. All calculated frequencies were in Hardy-Weinberg equilibrium. Two variants showed associations with AS: rs3212227 of IL-12β (allelic p = 0.031) and rs2293152 of STAT3 (allelic p = 0.006). These results might reveal that A and C represent risk alleles for Algerian AS patients within the IL-12β and STAT3 genes, respectively. For variants IL-12β rs6887695 and JAK2 rs7857730, no significant association was found between AS and Algerian individuals.

Multiplicative interaction analysis results for the two cohorts

The analysis of random combinations of two, three, or four SNPs by binary logistic regression model revealed 11 interactions for each cohort, as summarised in Table 2. Significantly different results between Chinese AS patients and controls were shown by the combinations rs3212227 * rs6887695 (p = 0.008), rs6887695 * rs7857730 * rs2293152 (p = 0.006), and rs3212227 * rs6887695 * rs7857730 * rs2293152 (p = 0.045). For Algerian AS patients and controls, only rs3212227 * rs7857730 * rs2293152 was found as a significant combination (p = 0.035).

Table 2

Multiplicative interaction analysis reported by logistic regression for the two populations investigated

PopulationChineseAlgerian
Interactionχ2OR (95% CI)p valueχ2OR (95% CI)p value
rs3212227 * rs68876956.9300.783 (0.652-0.939)0.0082.1430.722 (0.143-1.117)0.143
rs6887695 * rs78577300.6851.111 (0.865-1.427)0.4080.7440.736 (0.367-1.476)0.388
rs3212227 * rs78577301.2751.168 (0.892-1.530)0.2590.4320.899 (0.653-1.236)0.511
rs3212227 * rs78577300.1771.056 (0.819-1.363)0.6741.3221.301 (0.831-2.038)0.250
rs3212227 * rs22931521.7061.161 (0.928-1.452)0.1920.6621.155 (0.816-1.635)0.416
rs6887695 * rs22931520.0341.016 (0.855-1.208)0.8532.1921.408 (0.895-2.216)0.139
rs3212227 * rs6887695 * rs78577300.9371.119 (0.892-1.403)0.3330.6241.205 (0.758-1.916)0.430
rs6887695 * rs7857730 * rs22931527.4310.801 (0.683-0.940)0.0060.0960.949 (0.680-1.324)0.757
rs3212227 * rs6887695 * rs22931520.4901.064 (0.895-1.264)0.4840.4901.064 (0.895-1.264)0.484
rs3212227 * rs7857730 * rs22931520.9891.050 (0.954-1.156)0.3204.4681.227 (1.015-1.482)0.035
rs3212227 * rs6887695 * rs7857730 * rs22931524.0051.136 (1.003-1.287)0.0450.0230.981 (0.765-1.258)0.880

[i] OR – odds ratio, 95% CI: 95% confidence interval

[ii] p values < 0.05 are indicated in bold

MDR analysis results and gene-gene epistasis for the two cohorts

Various interaction models were tested by MDR analysis in order to detect any potential interaction effect of gene-gene epistasis represented in Table 3 for the two cohorts. All data were based on the highest testing accuracy and cross-validation consistency. For the Han Chinese cohort, the one-variable model SNP IL-12β rs3212227 showed the highest testing accuracy among the four SNPs. The two-variable gene-gene interaction model (IL-12β rs3212227, JAK2 rs7857730) and the three-variable gene-gene interaction model (IL-12β rs3212227, JAK2 rs7857730, and STAT3 rs2293152) showed the highest testing accuracy as compared to the other two-locus and three-locus interaction models (p = 0.0001 and p < 0.0001; respectively). MDR also showed a high testing accuracy for the four-variable gene-gene interaction model (p < 0.0001) for the Chinese cohort (Fig. 2).

Table 3

Gene-gene interaction model results of MDR for each population

ModelTrainingaccuracy (%)Testingaccuracy (%)Cross-validation consistency
Chinese
IL-12β (rs3212227)
IL-12β (rs3212227), JAK2 (rs7857730)
IL-12β (rs3212227), JAK2 (rs7857730), STAT3 (rs2293152)
IL-12β (rs3212227), IL-12β (rs6887695), JAK2 (rs7857730), STAT3 (rs2293152)
53.63
55.51
58.02
61.91
53.63
53.90
55.52
54.80
10/10
10/10
8/10
10/10
Algerian
STAT3 (rs2293152)
IL-12β (rs3212227), JAK2 (rs7857730)
IL-12β (rs3212227), JAK2 (rs7857730), STAT3 (rs2293152)
IL-12β (rs3212227), IL-12β (rs6887695), JAK2 (rs7857730), STAT3 (rs2293152)
57.68
62.63
66.73
61.91
53.62
59.49
62.44
54.80
8/10
10/10
10/10
10/10

[i] MDR – multifactor dimensionality reduction

Fig. 2

Graphic models determined by multifactor dimensionality reduction (MDR), representing the following different interactions in the Chinese cohort: A) IL-12β rs3212227, B) IL-12β rs3212227 and JAK2 rs7857730, C) IL-12β rs3212227, JAK2 rs7857730, and STAT3 rs2293152, D) IL-12β rs3212227, IL-12β rs6887695, JAK2 rs7857730, and STAT3 rs2293152. SNP genotypic values range from 0 to 2 (0 = zero risk allele, 1 = one risk allele, 2 = two risk alleles). Each cell shows counts of “Class 1 = unaffected” on the left and “Class 2 = affected” on the right, and shades on cells represent the risk degree of the disease (dark cell = high risk, light cell = low risk)

/f/fulltexts/CEJOI/36206/CEJI-44-84019-g002_min.jpg

For the Algerian cohort, the one-locus model SNP STAT3 rs2293152 represented the model with highest testing accuracy among the four SNPs. The two-locus (IL-12β rs3212227 and JAK2 rs7857730), the three-locus (IL-12β rs3212227, JAK2 rs7857730, and STAT3 rs2293152), and the four-locus interaction models all showed the highest testing accuracy (p < 0.0001) within the other interaction models (Fig. 3).

Fig. 3

Graphic models determined by multifactor dimensionality reduction (MDR), representing the following different interactions in the Algerian cohort: A) STAT3 rs2293152, B) IL-12β rs3212227 and JAK2 rs7857730, C) IL-12β rs3212227, JAK2 rs7857730, and STAT3 rs2293152, D) IL-12β rs3212227, IL-12β rs6w887695, JAK2 rs7857730, and STAT3 rs2293152. SNP genotypic values range from 0 to 2 (0 = zero risk allele, 1 = one risk allele, 2 = two risk alleles). Each cell shows counts of “Class 1 = unaffected” on the left and “Class 2 = affected” on the right, and shades on cells represent the risk degree of the disease (dark cell = high risk, light cell = low risk)

/f/fulltexts/CEJOI/36206/CEJI-44-84019-g003_min.jpg

Discussion

While the precise molecular aetiology of ankylosing spondylitis remains unclear, a variety of genetic factors have been involved in the development and incidence of the disease. Thus, many studies have reported a number of new non-MHC genes to be involved in the pathogenesis of AS, in particular those involved in the IL-23 pathway [27-29]. In the present study, our data indicated that SNP IL-12β rs3212227 (+1188A/C) was significantly associated with AS in the Han Chinese cohort. This polymorphism was also associated with AS in the Algerian cohort, and it is worth noting that the AA genotype may increase the risk of AS in both Chinese and Algerian subjects. Previously, linkage analysis studies in the Chinese population have considered IL-12B as a candidate gene for AS, and several polymorphisms were tested within this gene, but their results were contentious regarding whether some of these polymorphisms are associated with AS susceptibility and others are associated with disease severity [30, 31]. Nevertheless, in a recent association study of IL-12B polymorphism susceptibility with ankylosing spondylitis in the Mainland Han population, no significant difference was found for IL-12β rs3212227 variant, while only the SNP rs6871626 from the eight polymorphisms studied showed a significant difference in genotype distribution between AS and healthy controls [32]. Another meta-analysis studying the relationship between rs3212227 polymorphism of the IL-12B gene and susceptibility to multiple autoimmune diseases in East Asian populations failed to find any association with Graves’ disease and ankylosing spondylitis [23]. Additionally, the International Genetics of Ankylosing Spondylitis Consortium (IGAS) reported other SNPs within IL-12β that conferred risk of ankylosing spondylitis in Europeans, but the two SNPs analysed in our study were not genotyped [28].

In fact, IL-12β may play a key role in the pathogenesis of several inflammatory diseases because it encodes the p40 subunit of IL-23 and IL-12 cytokines and is thus involved in both the IL-12/Th1 and IL-23/Th17 pathways. Many studies were recently conducted to confirm the susceptibility of polymorphisms in IL-12β, in particular rs3212227 (+1188A/C) polymorphism, to cause ankylosing spondylitis (AS) and/or other autoimmune diseases such as rheumatoid arthritis (RA), type 1 diabetes (T1D), Behcet’s disease (BD), Graves’ disease (GD), and multiple sclerosis (MS) [23, 33]. Hence, taken together, our results may accentuate the importance of IL-12β rs3212227 (+1188A/C) as a novel risk polymorphism in the development of ankylosing spondylitis in both Chinese and Algerian subjects. Unfortunately, the second polymorphism in IL-12β rs6887695 was not associated with susceptibility to AS, neither in the Chinese nor in the Algerian individuals. The two polymorphisms rs3212227 and rs6887695 in IL-12β have been previously found to be associated with risk of psoriasis and psoriatic arthritis in the Chinese population [24, 34].

We further analysed the two other polymorphisms rs2293152 in STAT3 and rs7857730 in JAK2, but our results were only partially in agreement with the findings of Chen et al. [14], who reported no difference in the distribution of alleles and genotypes between AS groups and controls in their Chinese cohort in both JAK2 and STAT3 genes. In fact, our findings only showed a significant association between Chinese AS patients and controls for rs7857730 in the JAK2 gene. The difference in these findings is probably due to the difference in the recruitment of Chinese subjects for the two studies. Chen et al. also reported that haplotype analysis revealed an association of haplotype rs1536798/rs10119004/rs7857730-CGT in the JAK2 locus with AS, which is consistent with our results. Therefore, our findings may confirm the association of rs7857730 polymorphism in JAK2 with the susceptibility to AS and may emphasise again the importance of the IL-23/IL-17 axis as a causal pathway of AS in the Chinese population, for which earlier studies have also investigated its role in AS [27] and related diseases such as IBD [35] and psoriasis [36].

Another case-control study, carried out by Davidson et al. [15], in which a number of genes implicated in the pathogenic mechanisms of AS were tagged in the Han Chinese population, showed that SNP rs2293152 in STAT3 had a significant association with AS in their Chinese cohort, a result that corroborates the one previously found in STAT3 rs2293152 variant by the Australo-Anglo-American Spondyloarthritis Consortium (TASC) in Caucasian European populations [18]; however, this was in contrast to our results.

Unlike the Chinese cohort, we reported no significant difference between Algerian AS patients and controls for the variant JAK2 rs7857730, but a significant association for the variant STAT3 rs2293152 was found in this group. This discrepancy in our findings could be explained by differences in the genetic pools between the two populations investigated.

Our analysis also considered different interactions among one-variable, two-variable, three-variable, and four-variable gene-gene epistasis with ankylosing spondylitis, which were analysed to address concerns about DNA sequence variations that could have an interaction effect but not a statistically significant independent main effect on AS. Remarkably, logistic regression analysis showed significant interactions between IL-12β rs6887695 and the three other variables in the Han Chinese cohort, although this SNP itself showed no significant association with AS. MDR analysis also revealed a gene-gene interaction model for the Han Chinese cohort, in which IL-12β rs6887695 interacts with the other variables. This difference for IL-12β rs6887695 between individual association results and multiple interaction findings could be explained by the relatively small sample size. In addition, IL-12β rs3212227 showed the strongest association with AS in the Chinese population, so the significant interaction observed between the two variables of the IL-12β gene (rs3212227 and rs6887695) could confirm once again the important role played by this gene in the IL-23 signalling pathway and the pathogenesis of AS.

Analysis of interaction by logistic regression in the Algerian cohort revealed a significant three-variable interaction between IL-12β (rs3212227), JAK2 (rs7857730), and STAT3 (rs2293152) SNPs for AS susceptibility (p = 0.035). Because both IL-12β and STAT3 polymorphisms individually analysed showed a significant association with AS, this interaction suggests that IL-12β and STAT3 may be involved in the risk of developing AS in Algerian patients. In fact, STAT3 is predominantly activated by IL-23, which signals through IL-12Rβ1 and IL-23R to activate JAK and STAT signalling molecules, and induces IL-17A, IL-17F, and/or IL-22 and stabilises TH17 cells [19]. Thus, the additive interaction effect of the two SNPs in STAT3 and IL-12β may increase the risk of AS for Algerian patients.

In conclusion, we confirmed in our study a significant association between ankylosing spondylitis and the IL-12β gene for the rs3212227 variant in both the Han Chinese and Algerian populations. Also, JAK2 rs7857730 polymorphism showed an association in the Han Chinese population, while rs2293152 polymorphism in the STAT3 gene was associated with AS in the Algerian population. Despite the relatively small sample sizes and the number of SNPs selected for each gene, our findings remain novel. Additionally, the gene-gene interaction analysis through the effect of different SNPs was combined and may therefore provide a better understanding of AS and its pathogenic mechanisms.

Notes

[6] Conflicts of interest The authors declare no conflict of interest.

References

1 

Pile KD, Kennedy LG, Calin A, et al. , authors. HLA class I associations of ankylosing spondylitis in the white population in the United Kingdom. Ann Rheum Dis. 1996. 55:p. 268–270

2 

Wei JC, Tsai WC, Lin HS, et al. , authors. HLA-B60 and B61 are strongly associated with ankylosing spondylitis in HLA-B27-negative Taiwan Chinese patients. Rheumatology (Oxford). 2004. 43:p. 839–842

3 

Spencer CC, Pointon JJ, Su Z, et al. , authors. Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility. Nat Genet. 2011. 43:p. 761–767

4 

Yi L, Wang J, Guo X, Espitia MG, et al. , authors. Profiling of hla-B alleles for association studies with ankylosing spondylitis in the Chinese population. Open Rheumatol J. 2013. 7:p. 51–54

5 

Liu Y, Jiang L, Cai Q, et al. , authors. Predominant association of HLA-B*2704 with ankylosing spondylitis in Chinese Han patients. Tissue Antigens. 2010. 75:p. 61–64

6 

Brown MA, Kennedy LG, MacGregor AJ, et al. , authors. Susceptibility to ankylosing spondylitis in twins: the role of genes, HLA, and the environment. Arthritis Rheum. 1997. 40:p. 1823–1828

7 

Mustafa KN, Hammoudeh M, Khan MA , authors. HLA-B27 prevalence in Arab populations and among patients with ankylosing spondylitis. H J Rheumatol. 2012. 39:p. 1675–1677

8 

Amroun H, Djoudi H, Busson M, et al. , authors. Early-onset ankylosing spondylitis is associated with a functional MICA polymorphism. Hum Immunol. 2005. 66:p. 1057–1061

9 

Kchir MM, Hamdi W, Laadhar L, et al. , authors. HLA-B, DR and DQ antigens polymorphism in Tunisian patients with ankylosing spondylitis (a case-control study). Rheumatol Int. 2010. 30:p. 933–939

10 

Akassou A, Yacoubi H, Jamil A, et al. , authors. Prevalence of HLA-B27 in Moroccan healthy subjects and patients with ankylosing spondylitis and mapping construction of several factors influencing AS diagnosis by using multiple correspondence analysis. Rheumatol Int. 2015. 35:p. 1889–1894

11 

Reveille JD , author. Recent studies on the genetic basis of ankylosing spondylitis. Curr Rheumatol. 2009. 11:p. 340–348

12 

Reveille JD , author. Genetics of spondyloarthritis: beyond the MHC. Nat Rev Rheumatol. 2012. 8:p. 296–304

13 

Davidson SI, Wu X, Liu Y, et al. , authors. Association of ERAP1, but not IL23R, with ankylosing spondylitis in a Han Chinese population. Arthritis Rheum. 2009. 60:p. 3263–3268

14 

Chen C, Zhang XS, Wang Y, et al. , authors. Analysis of JAK2 and STAT3 polymorphisms in patients with ankylosing spondylitis in Chinese Han population. Clin Immunol. 2010. 136:p. 442–446

15 

Davidson SI, Liu Y, Danoy PA, et al. , authors. Association of STAT3 and TNFRSF1A with ankylosing spondylitis in Han Chinese. Ann Rheum. 2011. 70:p. 289–292

16 

Liu YC, Zhang H, Li J, et al. , authors. Association of common variants in KIF21B and ankylosing spondylitis in a Chinese Han population: a replication study. Immunogenetics. 2013. 65:p. 835–839

17 

Shan S, Dang J, Li JX, et al. , authors. ETS1 variants confer susceptibility to ankylosing spondylitis in Han Chinese. Arthritis Res Ther. 2014. 16:p. R87

18 

Reveille JD, Sims AM, Danoy P, et al. , authors. Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci. Nat Genet. 2010. 42:p. 123–127

19 

Parham C, Chirica M, Timans J, et al. , authors. A receptor for the heterodimeric cytokine IL-23 is composed of IL-12Rβ1 and a novel cytokine receptor subunit, IL-23R. J Immunol. 2002. 168:p. 5699–5708

20 

O’Garra A, Arai N , authors. The molecular basis of T helper 1 and T helper 2 cell differentiation. Trends Cell Biol. 2000. 10:p. 542–550

21 

Windsor L, Morahan G, Huang D, et al. , authors. Alleles of the IL12B 3’UTR associate with late onset of type 1 diabetes. Hum Immunol. 2004. 65:p. 1432–1436

22 

Tsunemi Y, Saeki H, Nakamura K, et al. , authors. Interleukin-12 p40 gene (IL12B) 30-untranslated region polymorphism is associated with susceptibility to atopic dermatitis and psoriasis vulgaris. J Dermatol Sci. 2002. 30:p. 161–166

23 

Zhang C, Wang JB, Chen SS, et al. , authors. Relationship between the IL12B (rs3212227) gene polymorphism and susceptibility to multiple autoimmune diseases: a meta-analysis. Mod Rheumatol. 2016. 26:p. 1–27

24 

Zhu K J, Zhu CY, Shi G, Fan YM , authors. Meta-analysis of IL12B polymorphisms (rs3212227, rs6887695) with psoriasis and psoriatic arthritis. Rheumatol Int. 2013. 33:p. 1785–1790

25 

Van der Linden S, Valkenburg HA, Cats A , authors. : Evaluation of diagnostic criteria for ankylosing spondylitis: a proposal for modification of the New York criteria. Arthritis Rheum. 1984. 27:p. 361–368

26 

QIAamp® DNA Mini and Blood Mini Handbook. 2010. Third Edition. p. 26–29. DNA purification from blood or body fluids (spin protocol).

27 

Burton PR, Clayton DG, Cardon LR, et al. , authors. Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat Genet. 2007. 39:p. 1329–1337

28 

International Genetics of Ankylosing Spondylitis Consortium (IGAS); Cortes A, Hadler J, Pointon JP, et al. , authors. Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat Genet. 2013. 45:p. 730–738

29 

Zhai J, Rong J, Li J, et al. , authors. Immunogenetic study in Chinese population with ankylosing spondylitis: are there specific genes recently disclosed? Clin Develop Immunol. 2013. 2013:p. 1–6. ID 419357.

30 

Zhou L, Zhang JY, Liu Y, et al. , authors. Novel association between IL-12B gene polymorphism and risk of ankylosing spondylitis in Ningxia population. Int J Lab Med. 2013. 34:p. 1372–1375

31 

Wong RH, Wei JC, Huang CH, et al. , authors. Association of IL-12B genetic polymorphism with the susceptibility and disease severity of ankylosing spondylitis. J Rheumatol. 2012. 39:p. 135–140

32 

Zhang L, Fan D, Liu L, et al. , authors. Association study of IL-12Bpolymorphisms susceptibility with ankylosing spondylitis in Mainland Han population. PLoS One. 2015. 10:p. 1–9

33 

Huang J, Yang Y, Zhou F, et al. , authors. Meta-analysis of the IL23R and IL12B polymorphisms in multiple sclerosis. Int J Neurosci. 2016. 126:p. 205–212

34 

Zheng HF, Zuo XB, Lu WS, et al. , authors. Variants in MHC, LCE and IL12B have epistatic effects on psoriasis risk in Chinese population. J Dermatol Sci. 2011. 61:p. 124–128

35 

Duerr RH, Taylor KD, Brant SR, et al. , authors. A genome- wide association study identifies IL23R as an inflammatory bowel disease gene. J Science. 2006. 314:p. 1461–1463

36 

Cargill M, Schrodi SJ, Chang M, et al. , authors. A large-scale genetic association study confirms IL12B and leads to the identification of IL23R as psoriasis-risk genes. Am J Hum Genet. 2007. 80:p. 273–290. 2007.

Copyright: © 2019 Polish Society of Experimental and Clinical Immunology This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0) License (http://creativecommons.org/licenses/by-nc-sa/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material, provided the original work is properly cited and states its license.
 
Quick links
© 2024 Termedia Sp. z o.o.
Developed by Bentus.