The Association of Interleukin 6 Receptor Gene Variant (Rs4129267) with Susceptibility to Ankylosing Spondylitis and Soluble IL-6 Receptor Level in Iranian Patients

Abstract

Objectives:

Ankylosing spondylitis (AS) is an axial skeleton-related auto-inflammatory arthritis. Interleukin 6 (IL-6) and its receptor (IL-6R) play a role in the etiopathology of AS. The recent genome-wide association studies (GWASs) revealed that rs4129267; an intrinsic variant of IL6-R is associated with the risk of AS. We aimed to study the rs4129267 IL6R variant in Iranian AS patients and the correlation of this variant with the soluble interleukin-6 receptor (sIL-6R) serum level.

Materials and Methods:

In this case-control study, we genotyped rs4129267 IL6-R variant in extracted DNA from 498 AS patients and 495 matched healthy controls by Amplification-Refractory Mutation System Polymerase Chain Reaction (ARMS-PCR). The serum level of sIL-6R and IL-6 was also analysed in 39 controls and 42 AS patients by the ELISA method. The correlation between rs4129267 IL-6R variants and patients’ clinical manifestations and sIL-6R level was also investigated.

Results:

In comparison to the AS individuals, C allele were significantly less prevalent in the healthy control group (OR = 1.23 [1.01 –1.48], P value = .034). Patients had significantly elevated levels of IL-6 serum levels. There was a significant correlation between the presence of the rs4129267 T allele and increased sIL-6R level.

Conclusions:

Our study confirmed the significant association between the rs4129267 T variant and the protection against AS disease susceptibility. Increased IL-6 level and reduced frequency of sIL-6R rs4129267 T allele in patients which contributes to reduced sIL-6R serum level, highlights the significant role of IL-6 signalling pathway in AS pathogenesis.

Keywords

Ankylosing spondylitis interleukin 6 interleukin 6 receptor single nucleotide poly

Introduction

Ankylosing spondylitis (AS) is an inflammatory rheumatic condition that often affects the sacroiliac joints, spine, and hip joints. Typically, AS is characterised by spondylitis and sacroiliitis accompanied by severe joint pain, stiffness, and functional disability in damaged joints. In the severe form of the disease, persistent inflammation of these areas causes osteogenesis in intervertebral discs and fusion of the vertebrae together. The result of these changes will be a loss of movement and flexibility in these parts.¹ AS is more common in young people between the ages of 16 and 36 and usually lasts a lifetime. Although it’s less prevalent and the symptoms are milder among women, the involvement of joints other than the spine is more likely in women than in men. In 5% of patients with AS, the disease begins in childhood.² A study of the Iranian population shows that the disease is 3 to 1 more common in men. Besides, the average age of onset of symptoms is 24.7 years and the incidence of the disease in the south of Iran was reported approximately 8.03 per 100,000 people.³

While the cause of AS is unknown, studies have shown the contribution of both genetic and environmental factors in the development of AS. There is a strong genetic relationship between the incidence of the disease and the presence of human leukocyte antigen (HLA)-B27 polymorphisms. About 90% of patients are HLA-B27 positive, nonetheless only about 5% of respective patients develop the disease.^4,5 Genome-wide association studies (GWAS) suggested that genetic variants in both HLA and non-HLA genes are associated with susceptibility to AS and its progression.⁶

Genetic studies reveal that several variants in key genes of the pro-inflammatory IL-23/IL-17 pathway are associated with AS disease and its pathogenesis.^6,7 Interleukin (IL)-6, an important cytokine related to the IL-23 pathway, is involved in the pathogenesis of AS. IL-6 pro-inflammatory properties are mediated by binding to the IL-6 receptor (IL-6R). IL-6R has two forms by alternative splicing⁸: (a) The membrane-bound form (mIL-6R) that is expressed in neutrophils, monocytes, CD4⁺ T cells, B cells, and hepatocytes; (b) The soluble form (sIL-6R) circulating in body fluids which are mostly generated by proteolysis from the mIL-6R. The IL-6/sIL-6R complex could be a mediator for the IL-6 signalling pathway in non-expressing mIL-6R cells by binding to gp130 transmembrane glycoprotein.^9,10

Many studies have reported elevated serum levels of IL-6 in AS patients and its association with the severity of patients’ clinical manifestations.^11,12 Although targeting IL-6 signalling in AS patients by tocilizumab (an IL-6R antagonist) has been repeatedly investigated, the results were inconsistent and some studies demonstrated tocilizumab as a non-efficient treatment.^13,14 The differences in clinical responses to tocilizumab could be due to variations in the IL-6R gene.

The rs4129267 single nucleotide polymorphism (SNP) of the IL-6R gene is a candidate risk variant for AS and consistently, GWAS studies revealed that the T allele of this SNP is significantly associated with increased susceptibility to AS progression,^15-17 which is a non-synonymous intronic variant that has a gene-regulatory effect.¹⁸ This non-coding variant had an impact on sIL-6R availability and serum levels.¹⁹ In previous studies, rs4129267 showed a strong association with AS risk in a European cohort including European, South Asian, and Latinx populations¹⁹ and a Chinese cohort.²⁰ Moreover, a significant association with asthma is reported, while a negative association with AS.²¹

In this study, we aimed to examine the association of the rs4129267 variant of the IL-6R gene, as a susceptible loci, with the risk of AS in the Iranian population, to check the reproducibility of the results in a population from the Middle East and investigate the correlation between rs4129267 variants and sIL-6R serum level.

Materials and Methods

Participants

A total of 498 AS patients (388 men and 110 women, mean age: 37.5 years) were referred to the rheumatology clinic of Shariati Hospital of Tehran University of Medical Sciences between November 2013 and June 2015 based on modified New York criteria²² plus 495 ethnicity-, age- and sex-matched healthy controls (mean age: 36.8 years) without any familial history of autoimmune diseases participated in this study. Clinical and demographical characteristics of AS patients and healthy individuals were gathered at the time of sample collection in a questionnaire (Table 1). Informed consent was taken from all participants before enrollment.

Table 1.

Baseline Characteristics of AS Patients and Healthy Controls.

Variables	AS patient (N = 498)	Control (N = 495)	P Value
Age (years)	37.5 ± 10.1	36.7 ± 10.6	.111
Sex (female)	106 (21%)	113 (23%)	.702
BASFI score (0–10)	3.7 ± 2.6	–	–
BASMI score (0–10)	4.0 ± 1.8	–	–
BASDAI score	4.7 ± 2.5	–	–
ASQOL score	7.7 ± 5.3	–	–
NRS	7.7 ± 5.4	–	–
CRP* (µg/ml)	17.4 ± 23.2	1.20 ± 3.4	<.001
ESR* (mm/h)	19.97 ± 19.8	5.33 ± 4.7	<.001

Notes: BASFI: Bath Ankylosing Spondylitis Functional Index; BASMI: Bath Ankylosing Spondylitis Metrologic Index; BASDAI: Bath Ankylosing Spondylitis Disease Activity Index; ASQOL: Ankylosing Spondylitis Quality of Life; NRS: Numerical Rating Scale of pain.

Data are presented as mean (SD) or counts (percentages). *P value <.05 was considered statistically significant.

DNA Extraction and Genotyping

Whole DNA was extracted from 4 mL venous blood samples by phenol-chloroform method.²³ The optical density values were used to evaluate the concentration and purity of the extracted DNA (NanoDrop 2000C). Subsequently, it was stored at -20°C. The rs4129267 is located on position 154453788 (intron variant) of chromosome 1 and we employed the ARMS-PCR technique for its genotyping. To design the assay, three primers were used: two forward primers (one for the mutant and one for the wild-type) and a single reverse primer common to both, which were created using web-based allele-specific PCR (WASP). The ARMS-PCR was carried out in a reaction mixture with a total volume of 15 µL, which included 1.5 µL of genomic DNA (150 ng), 7.1 µL of distilled water, 0.25 µL of 10 mM dNTP mix, 1.2 µL of each specific primer at a concentration of 10 pmol/µL, 2.5 µL of 10X PCR buffer, 0.75 µL of 1.5 mM MgCl2, and 0.3 µL of Taq DNA polymerase. Additionally, two internal control primers for growth hormone (GH) were added at a volume of 0.1 µL each (with a concentration of 10 pmol/µL) to every reaction. The PCR cycling conditions were consisting of the initial denaturation (94°C, 3 min) followed by 40 cycles of 94°C for 30 s, 62°C for 30 s, 72°C for 30 s, final extension at 72°C for 10 min and eventually 10 min of storage at 4°C. Amplified DNA products were separated on a 2% agarose gel, stained with a DNA-safe dye, and visualised using a UV transilluminator. To verify the accuracy of the amplification, Sanger sequencing was performed on the PCR products. The sequence of primers for ARMS-PCR genotyping and Sanger sequencing is listed in Table 2.

Table 2.

Primers for ARMS-PCR Genotyping.

Primer Name	Primer Sequence (IL6R/rs4129267)	Length	Tm	GC%	Product Size
IL6R-F Seq	5’-GTTTAGATGTGCCCTCTGTCTAGT-3’	24	59.84	45.83	302
IL6R-NF	5’- AGCTTGGAGTGGGGTCAATTCGC -3’	23	62.46	52.17	296
IL6R-MF	5’- AGCTTGGAGTGGGGTCAATTCGT -3’	23	62.46	52.17
IL6R-R	5’- ATGTCAGAGAGGATTGGATGCTTTA -3’	25	59.87	40.00
IL6R-R Seq	5’-CTCCTGTCTGTTTTCTCTACCTTTG-3’	25	59.30	44.00	371
GH (F)	5’-TGCCTTCCCAACCATTCCCTTA-3’	22	61.65	50.00
GH(R)	5’-CCACTCACGGATTTCTGTTGTGTTTC-3’	26	62.85	46.15	434

Note: IL6R-F: Forward primer for IL-6R gene; IL6R-NF: Forward primer for non-mutant IL-6R gene; IL6R-MF: Forward primer for mutant IL-6R gene; IL6R-R: Reverse primer for IL-6R gene; IL6R-Seq: Reverse primer for Sanger sequencing of IL-6R gene; GH (F): Forward primer of growth hormone; GH(R): Reverse primer of growth hormone.

Serum Protein Measurements

Serum was isolated from 10 ml blood samples of 42 AS patients and 39 healthy controls by centrifuging at 4,500 rpm for 7 minutes. The ELISA kits (R&D Systems, Minneapolis, MN, USA) were utilised to measure sIL-6R and IL-6 levels based on the manufacturer’s procedures.

Statistical Analysis

Data were described as mean ± SD or n (%). We used the Shapiro-Wilks test to assess the normality distribution for continuous variables. Association between categorical variables was assessed by Chi-square test, odds ratio (OR), and 95% confidence interval (CI). The genotype distributions of SNP were examined for statistical deviation from Hardy-Weinberg equilibrium. We used the Kruskal–Wallis test and one-way ANOVA to compare the means of continuous variables among groups; Also, the Bonferroni post hoc test was used for pairwise comparison if needed. We performed this statistical analysis using SPSS software (version 27.0, IBM SPSS Inc., USA), and a P value less than .05 was considered statistically significant.

Results

Clinical Features and Inflammatory Markers

Clinical features of AS patients and healthy individuals are presented in Table 1. There is no notable difference between patients and healthy individuals regarding their age and sex (P = .111 and P = .702, respectively). The AS group had significantly elevated levels of erythrocyte sedimentation rate (ESR) and C-reactive protein (CRP) (P <.001 for all) as expected.

The rs4129267 Genotyping

Table 3 displays the associations between rs4129267 alleles and genotype frequency with AS susceptibility. In the healthy group, there was no significant deviation from the Hardy-Weinberg Equilibrium for rs4129267 genotypes distribution (P = .714). According to our results, the CC genotype is significantly more common in the AS group (CC vs. TT+TC P = .038), and the C allele is significantly increased in AS patients [P = .034, OR (95% CI) = 1.23 (01.01–1.48)] (Data not shown).

Table 3.

Alleles and Genotypes Distribution of rs4129267 Variant in AS and Healthy Groups.

SNP	Allele/Genotype	AS (N = 498) N (%)	Control (N = 495) N (%)	OR (95% CI)	P Value
rs4129267	T	292 (29.3)	334 (33.7)
	C*	704 (70.7)	656 (66.3)	1.23(1.01, 1.48)	.034
	TT	45 (9.0)	56 (11.3)
	TC	202 (40.6)	222 (44.8)	1.13(0.73, 1.75)	.576
	CC+TC	453 (91%)	439 (88.7%)	1.28(0.85, 1.94)	.235
	CC	251 (50.4)	217 (43.8)	1.44(0.93, 2.22)	.098
HWE					.714

Note: AS: Ankylosing Spondylitis; OR: Odd Ratio; CI: Confidence Interval; HWE: Hardy-Weinberg equilibrium;^*P value <0.05 was considered statistically significant.

IL-6 and sIL-6R Serum Level

We also assessed the level of IL-6 and sIL-6R in serum samples of 39 healthy and 42 randomly selected patients. Although IL-6 level was significantly upregulated in AS patients (P <.001) there was no statistically significant difference in the level of sIL-6R between the two groups (P = .182) (Table 4).

Table 4.

IL-6 Cytokine and sIL-6R Levels in AS Patients and Healthy Controls.

Variables	AS Patient (N = 42)	Control (N = 39)	P Value
IL-6 (pg/ml)	5.35 ± 2.59	1.99 ± 2.59	<.001*
sIL-6R (ng/ml)	30.6 ± 6.9	28.63 ± 5.5	.182**

Note: *Calculated based on Mann Whitney test; the Median (Q1, Q3) of IL-6 is 5.34 (3.21, 6.99) and 1.61 (0.97, 2.35) in AS patients and controls, respectively. **Calculated based on independent t-test.

Association of rs4129267 Genotypes and Clinical Features

The association between rs4129267 genotypes with HLA-B27 positivity and the score of patients’ clinical features by Bath Ankylosing Spondylitis Metrologic Index (BASMI), Bath Ankylosing Spondylitis Functional Index (BASFI), Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), Ankylosing Spondylitis Quality of Life (ASQOL), Numerical Rating Scale of pain (NRS), ESR and CRP were analysed and are shown in Table 5. No significant association between the mentioned clinical indices with different genotypes of rs4129267 was observed.

Table 5.

Association of rs4129267 Variant with Disease Severity Indicators and Inflammatory Factors.

Clinical Manifestations	Mean ± SD			P Value*
Clinical Manifestations	CC (n = 251)	TC (n = 202)	TT (n = 45)	P Value*
BASFI score	3.59 ± 2.6	3.91 ± 2.6	3.43 ± 2.6	.343
BASMI score	4.01 ± 1.8	4.13 ± 1.9	3.86 ± 1.6	.671
BASDAI score	4.69 ± 2.5	4.94 ± 2.4	4.13 ± 2.3	.167
ASQOL score	7.33 ± 5.5	8.30 ± 4.8	6.85 ± 5.6	.110
ESR	19.38 ± 18.85	21.02 ± 20.85	18.26 ± 19.76	.058
NRS	7.95 ± 5.4	7.58 ± 5.3	6.27 ± 5.2	.140
CRP	19.94 ± 27.11	13.79 ± 15.42	17.47 ± 24.81	.239
HLA-B27 (positive)	173 (68.9%)	128(63.4%)	27(60.0%)	.698

Note: *P value(s) calculated based on Kruskal–Wallis test.

Association of CRP Levels with IL-6 and sIL-6R in Different rs4129267 Genotypes

Figure 1 shows that, unlike TT, levels of CRP and IL-6 are significantly correlated (P <.05) in patients with CC and TC genotypes. However, no significant association was observed between CRP and sIL-6R levels in different genotypes (Figure 2).

Figure 1.

The Association of CRP and IL-6 Levels Divided by Different Genotypes.

Figure 2.

The Association of CRP and sIL-6R Levels Divided by Different Genotypes.

Association of rs4129267 Genotypes with sIL-6R and IL-6 Level

Later on, we investigated the association between rs4129267 variants and sIL-6R level in participants. The findings indicated a significant correlation between the rs4129267 variants and the levels of sIL-6R in the serum. In particular, those possessing the TT genotype were observed to have the highest levels of sIL-6R, whereas individuals with the CC genotype exhibited the lowest levels. Furthermore, there was a noteworthy decrease in the serum levels of sIL-6R among AS patients and control subjects who have the TC genotype compared with TT.

Conversely, while the genotype did not display a significant link to the serum levels of IL-6, AS patients, regardless of their genotype, had a significantly higher concentration of IL-6 in their serum compared to the control group (Table 6).

Table 6.

Comparison of IL-6 and sIL-6R Levels Among CC, TC, and TT Genotypes Separately in AS Patients and Healthy Controls.

Variables		Mean ± SD			P Value
Variables		CC (n = 31)	TC (n = 38)	TT (n = 11)	P Value
IL-6 (pg/ml)	AS patient	5.78 ± 2.86	4.87 ± 2.49	5.54 ± 1.51	.557
	Control	2.33 ± 1.66	1.99 ± 1.48	1.37 ± 0.53	.303
	P value	<.001	<.001	<.001	–
sIL-6R (ng/ml)	AS patient	25.37 ± 4.91	34.26 ± 5.01*	38.15 ± 6.20*	<.001
	Control	24.54 ± 3.55	29.14 ± 4.15*^Ä	34.52 ± 7.00*	<.001
	P value	.615	.002	.413	–

Note: (*: TT or TC vs. CC; Δ: TC vs. TT) = Statistically significant (P value <.05) based on multiple comparisons by Bonferroni post hoc test.

Discussion

Several GWASs reveal that IL-23 pathway-related genes are associated with susceptibility to AS disease.²⁴ IL-6R is an effector signalling molecule in the IL-23 pathway^25,26 and implicates in the pathogenesis and auto-inflammatory responses of AS.⁵ Association of IL-6R rs4129267 intronic variant with AS was previously reported.¹⁷ Here, we studied the association of rs4129267 with the susceptibility to AS in the Iranian population for the first time and our findings highlighted the role of the IL-6R gene in AS disease.

The IL-6R variant rs4129267 is described as a biomarker for inflammation in a GWAS study in the Italian population.²⁷ Ferreira et al. also reported the association of rs4129267 with the susceptibility of asthma in an Australian population and found that T allele increases the risk of asthma.¹⁸ Cortes et al. studied 10470 AS patients and 12338 controls, of European, South Asian, and Latin American origin. They found that the rs4129267 C allele is significantly associated with an elevated risk of AS in individuals of European ancestry individuals, while the T allele plays a protective role. Besides they observed that this variant is associated with significantly higher concentrations of CRP in the TT genotype than CC, which highlights its impact on the severity and clinical features of AS.¹⁷ Furthermore, Ruan et al. found the frequency of TT genotype rising significantly in the AS group of the Chinese‎ population from Wuhan. Additionally, in line with clinical characteristics, significantly elevated BASFI scores were observed in the TT genotype compared to CC. Apart from the genotype, they found out that the rs4129267 T allele could also serve as a prognostic factor for AS pathogenesis because of its influence on the worsening of treatment response during the follow-up.¹⁵ In our study, the frequency of the T allele is significantly reduced in the AS group, and consistent with Cortes et al., our results showed that the T allele has a protective role against AS incidence in the Iranian population.

Previously it was found that rs4129267 is a regulatory variant and has an impact on the level of sIL-6R.¹⁷ The study by Melzer et al. showed sIL-6R serum concentration to be higher in people with rs4129267 TT genotype than other genotypes and each T allele increases 1.4-fold the serum concentration of sIL-6R.²⁸ Consistent with previous reports, we showed a strong correlation between increased serum concentration of sIL-6R and rs4129267 TT genotype and found that each T allele increases the receptor serum concentration. Although the small sample size for this measurement was one of our limitations in this study; considerable significant results were achieved that could be increased with enhanced power of analysis by enrolling more serum samples of participants. In a study by Aparicio-Siegmund et al., it was found that sIL-6R and sgp130 protein act as IL-6-neutralising buffers and are reduced in diabetic patients.²⁹ They found that the higher level of sIL-6R can bind to IL-6 and neutralise low amounts of it and consequently block IL-6 signalling through the mIL-6R.²⁹ These findings confirmed the regulatory role of the rs4129267 variant on the sIL-6R level.

In our study, we found that the T allele, which is reduced in AS patients, is strongly and significantly associated with an elevated level of sIL-6R serum level. It seems that the rs4129267 T allele by increasing the serum level of sIL-6R, acts as an IL-6 neutralising buffer and by blocking the IL-6 signalling pathway may prevent AS disease incidence. Blocking of IL-6 signalling using its anti-IL-6R antibody (tocilizumab) is one of the treatments for this rheumatic disease; however, the response to tocilizumab for AS is still obscure. A Pharmacogenetics study of the correlation between IL-6R genotypes and response to tocilizumab could help to resolve ambiguities.

However, we observed an increase in sIL-6R levels with increasing T allele; the effect of this increase is not propagated in disease severity indicators or inflammatory factors. This rather indicates that rs4129267 does not have an impact on AS pathogenesis.

Conclusion

In conclusion, we identified the association between the rs4129267 T variant and the protection against susceptibility to AS disease in the Iranian population for the first time. We found that AS patients had significantly elevated IL-6 levels compared to healthy individuals. Increased IL-6 level and reduced frequency of IL-6R rs4129267 T allele in patients which contributes to reduced sIL-6R serum level, highlights the regulatory role of this variant and the importance of IL-6 signalling pathway in AS pathogenesis. However, more investigations are necessary to determine the exact function of IL-6R variants in AS pathogenesis.

Footnotes

Authors’ Contributions

SE, MA, FD, NA, SM, YV, and EM: Acquisition of data, drafting the article, analysis and interpretation of data, and final approval of the article. EF, MA, AJ, and MM: The conception and design of the study, revising the article critically, interpretation of data, final approval of the article.

Availability of Data and Material

All data generated or analyzed during this study are available upon request.

Consent to Participate

The written informed consent was signed by all participants before enrolling in the study.

Consent for Publication

Consent was not applicable, as this is a review article compiled from various research articles and guidelines and not from patients directly.

Declaration of Conflicting Interests

The authors declared no potential conflicts of interest with respect to the research, authorship and/or publication of this article.

Ethical Approval

The procedures were approved by the ethics committee of the National Institute for Medical Research Development (Approval No: IR.TUMS.MEDICIN.REC.1397.774).

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by a grant from the Deputy of Research, Tehran University of Medical Sciences (Grant No. 93-11-28-8001).

References

Khan

, van der Linden

. Ankylosing spondylitis and other spondyloarthropathies. Rheum Dis Clin North Am. 1990;16:551–579.

Dean

, Jones

, MacDonald

, Downham

, Sturrock

, Macfarlane

. Global prevalence of ankylosing spondylitis. Rheumatology (Oxford). 2014;53:650–657.

Shahlaee

, Mahmoudi

, Nicknam

, Farhadi

, Fallahi

, Jamshidi

. Gender differences in Iranian patients with ankylosing spondylitis. Clin Rheumatol. 2015;34:285–293.

Brown

, Kenna

, Wordsworth

. Genetics of ankylosing spondylitis–insights into pathogenesis. Nat Rev Rheumatol. 2016;12:81–91.

Mahmoudi

, Aslani

, Nicknam

, Karami

, Jamshidi

. New insights toward the pathogenesis of ankylosing spondylitis; genetic variations and epigenetic modifications. Mod Rheumatol. 2017;27:198–209.

Tsui

, Tsui

, Akram

, Haroon

, Inman

. The genetic basis of ankylosing spondylitis: new insights into disease pathogenesis. Appl Clin Genet. 2014;7:105–115.

Jethwa

, Bowness

. The interleukin (IL)-23/IL-17 axis in ankylosing spondylitis: new advances and potentials for treatment. Clin Exp Immunol. 2016;183:30–36.

Lust

, Donovan

, Kline

, Greipp

, Kyle

, Maihle

. Isolation of an mRNA encoding a soluble form of the human interleukin-6 receptor. Cytokine. 1992;4:96–100.

Müllberg

, Schooltink

, Stoyan

, . The soluble interleukin-6 receptor is generated by shedding. Eur J Immunol. 1993;23:473–480.

10.

Chalaris

, Rabe

, Paliga

, . Apoptosis is a natural stimulus of IL6R shedding and contributes to the proinflammatory trans-signaling function of neutrophils. Blood. 2007;110:1748–1755.

11.

Przepiera-Będzak

, Fischer

, Brzosko

. Serum IL-6 and IL-23 levels and their correlation with angiogenic cytokines and disease activity in ankylosing spondylitis, psoriatic arthritis, and SAPHO syndrome. Mediators Inflamm. 2015;2015: 785705.

12.

Liu

, Wu

, Zhang

, . Elevated serum levels of IL-6 and IL-17 may associate with the development of ankylosing spondylitis. Int J Clin Exp Med. 2015;8:17362–17376.

13.

Sieper

, Porter-Brown

, Thompson

, Harari

, Dougados

. Assessment of short-term symptomatic efficacy of tocilizumab in ankylosing spondylitis: results of randomised, placebo-controlled trials. Ann Rheum Dis. 2014;73:95–100.

14.

Henes

, Horger

, Guenaydin

, Kanz

, Koetter

. Mixed response to tocilizumab for ankylosing spondylitis. Ann Rheum Dis. 2010;69:2217–2218.

15.

Ruan

, Xie

, Jin

, Wang

, Ping

. The diagnostic and prognostic role of interleukin 12B and interleukin 6R gene polymorphism in patients with ankylosing spondylitis. J Clin Rheumatol. 2018;24:18–24.

16.

Australo-Anglo-American

Spondyloarthritis C

, Reveille

, Sims

, . Genome-wide association study of ankylosing spondylitis identifies non-MHC susceptibility loci. Nat Genet. 2010;42:123–127.

17.

Cortes

, Hadler

, Pointon

, . Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat Genet. 2013;45:730–738.

18.

Ferreira

, Matheson

, Duffy

, . Identification of IL6R and chromosome 11q13.5 as risk loci for asthma. Lancet. 2011;378:1006–1014.

19.

Cortes

, Hadler

, Pointon

, . Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat Genet. 2013;45:730.

20.

Ruan

W-F

, Xie

J-T

, Jin

, Wang

W-D

, Ping

A-S

. The diagnostic and prognostic role of interleukin 12B and interleukin 6R gene polymorphism in patients with ankylosing spondylitis. J Clin Rheumatol. 2018;24:18–24.

21.

Ferreira

, Matheson

, Duffy

, . Identification of IL6R and chromosome 11q13. 5 as risk loci for asthma. Lancet. 2011;378:1006–1014.

22.

Deodhar

. Axial spondyloarthritis criteria and modified NY criteria: issues and controversies. Clin Rheumatol. 2014;33:741–747.

23.

Köchl

, Niederstätter

, Parson

. DNA extraction and quantitation of forensic samples using the phenol-chloroform method and real-time PCR. Methods Mol Biol. 2005;297:13–30.

24.

Adamopoulos

, Tessmer

, Chao

, . IL-23 is critical for induction of arthritis, osteoclast formation, and maintenance of bone mass. J Immunol. 2011;187:951–959.

25.

Yen

, Cheung

, Scheerens

, Poulet

, . IL-23 is essential for T cell-mediated colitis and promotes inflammation via IL-17 and IL-6. J Clin Invest. 2006;116:1310–1316.

26.

Bianchi

, Rogge

. The IL-23/IL-17 pathway in human chronic inflammatory diseases-new insight from genetics and targeted therapies. Genes Immun. 2019;20:415–425.

27.

Naitza

, Porcu

, Steri

, . A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation. PLoS Genet. 2012;8:e1002480.

28.

Melzer

, Perry

, Hernandez

, . A genome-wide association study identifies protein quantitative trait loci (pQTLs). PLoS Genet. 2008;4:e1000072.

29.

Aparicio-Siegmund

, Garbers

, Flynn

, . The IL-6-neutralizing sIL-6R-sgp130 buffer system is disturbed in patients with type 2 diabetes. Am J Physiol Endocrinol Metab. 2019;317:E411–e20.