Copy number variation of IL17RA gene and its association with the ankylosing spondylitis risk in Iranian patients: a case-control study
BMC Medical Genetics volume 21, Article number: 147 (2020)
Ankylosing spondylitis (AS) is considered as a subtype of spondyloarthritis (SpA) that mainly leads to fatigue, stiffness, spinal ankylosis, and impaired physical functions with reduced quality of life. Interleukin (IL)-17A provokes additional inflammatory mediators and recruits immune cells to the inflamed site. IL17 expression increased in various inflammatory disorders including psoriasis, rheumatoid arthritis, multiple sclerosis, crohn’s disease, and ankylosing spondylitis. The current study aimed to evaluate the association of IL17RA copy number changes with the susceptibility to AS and their correlation to IL17RA expression in Iranian population.
IL17RA copy number genotyping assessments were carried out in 455 AS patients and 450 healthy controls, using custom TaqMan CNV assays. TaqMan primers and probe were located in Chr.22:17109553 based on pre-designed IL17RA Copy Number Assay ID, Hs02339506_cn. mRNA expression of IL17RA was also measured by SYBR Green real-time polymerase chain reaction (PCR).
A IL17RA copy number loss (< 2) was associated with AS compared to 2 copies as reference (OR:2.18, 95% CI: (1.38–3.44), P-value < 0.001) and increased the risk of AS. IL17RA mRNA expression showed a significant increase in peripheral blood mononuclear cells (PBMCs) of all AS individuals than controls. The mRNA expression level of 2 copies was significantly higher in AS patients.
Our findings revealed that a low copy number of IL17RA might confer a susceptibility risk to AS. However, it is probably not directly involved in the regulation of IL17RA mRNA expression. Epigenetic mechanisms like DNA methylation, post-transcriptional, and -translational modifications that regulate the expression of the genes may contribute in upregulation of IL17RA mRNA expression in the loss of gene copy number condition.
Ankylosing spondylitis (AS) is considered as a subtype of spondyloarthritis (SpA) that mainly leads to fatigue, stiffness, spinal ankylosis, and impaired physical functions with reduced quality of life (QoL) [1, 2]. In spite of various investigations concerning AS pathogenesis, the exact cause, and underlying mechanisms are still less clear [3,4,5,6]. However, according to the numerous studies taken in the past decades, genetic background and immunological dysfunctions are regarded as pivotal bases of the disease [7,8,9,10].
Interleukin (IL)-17A, the first known subset of the IL17 family, primarily secreted from multiple cell types including, γδ T cells, natural killer T cells (NKT cells), group 3 Innate lymphoid cells (ILC3s) and particularly Th17 [11, 12]. This pro-inflammatory cytokine provokes additional inflammatory mediators and recruits immune cells to the inflamed tissue. IL-17 mostly leads to production of several chemical messengers such as proinflammatory cytokines (IL-6, IL-8, TNF-α, IL-1β), granulocyte attracting chemokines (granulocyte-colony stimulating factor (G-CSF), granulocyte monocyte-colony stimulating factor (GM-CSF), CCL2 (MCP-1), CXCL2 (MIP-2), CXCL5, CCL20 (MIP-3A)), matrix metalloproteinases (MMP1, MMP3, MMP9, and MMP13) and anti-microbial peptides [13,14,15,16].
Taken together, the aberrant IL-17 signaling pathway contribute to developing several autoimmune disorders . It is clearly stated that the excessive production of IL-17 can lead to various inflammatory diseases such as psoriasis, rheumatoid arthritis (RA), AS, multiple sclerosis (MS) and crohn’s disease (CD) [17,18,19]. Accordingly, prior investigations have found the increased concentration of IL-17 in the serum of AS patients and the frequency of Th17 cells also elevated in peripheral blood and synovial fluids of AS patients [20, 21].
The IL-17 activation is provided after binding to receptor site on target cells. As previously reported, IL-17 binds to a complex of IL-17RA and IL-17RC to induce responses via its unique signaling pathway . IL-17RA (also termed as IL-17R), is the first member of IL-17R superfamily (IL-17RA-IL-17RE) which recognizes both IL-17A and IL-17F, but have more affinity to IL-17A than IL-17F about 100 fold greater . This family is distinctive from other known cytokine receptor families with its particular structure and signaling cascade through ACT1 (NF-κB activator 1) adaptor molecule . IL17RA showed expression on the variety of tissues, and cell types, particularly in hematopoietic cells and upon stimulation with its ligand can activate NF-κB, a major inflammatory transcription factor [25, 26]. In addition to its capacity to induce unbalanced cytokine production, it can reinforce the effect of other potent cytokines including TNF-α and IL-1β, which leads to prolonged inflammatory responses and exacerbate chronic inflammations [27, 28].
To date, there are mysteries in genetic aspects of AS. Recent discoveries have emphasized the substantial role of single nucleotide polymorphisms (SNPs) in AS pathogenesis [29,30,31,32]. The newly performed study by Vidal-Castiñeira et al. assessed the association of IL17RA polymorphisms with AS risk . Besides SNPs, copy number variations (CNVs) are believed to contribute to the etiology of ankylosing spondylitis. They are known as different copies of a specific fragment of DNA created due to deletions or duplications that captured 1 Kilobase to several megabases of the genome . Given that CNVs are harboring 12% of the human genome, they can strongly affect the expression of genes encoding proteins . CNVs exert their influence on gene expression through gene dosage imbalances in dosage-sensitive genes and cause changes in mRNA and eventually in protein level . In addition to simple gene dosage effect, rearrangements in the gene regulatory elements can also alter gene expression. These rearrangements lead to changes in the interval between genes to regulatory regions. New fused genes can form as a result of rearrangements . Moreover, a growing number of reports underscore the importance of CNVs in some complex diseases such as RA, psoriasis, Crohn’s disease, and systemic lupus erythematosus (SLE), autism, schizophrenia, type 1 diabetes (T1D), human immunodeficiency virus (HIV) and cancers [38, 39].
Recently, a genome-wide copy number variation analysis has identified 9 loci as potential candidate CNV regions . Among them, IL17RA has been significantly associated with AS and showed a protective role in susceptibility to AS . To gain more insight into the genetic structure of AS, it would be instrumental in following up initial clues provided by the copy number variation analysis study.
With regards to IL17RA expression on immune cells, it seems that changes in the copy number of the gene might influence the gene by altering its expression. However, the contribution of IL17RA gene CNVs in the pathogenesis of AS and its influence on IL17RA gene expression has not been investigated yet. Here, we set out to determine whether IL17RA copy number variations would predispose individuals to the AS in the Iranian population.
A total of 455 AS patients (85.1% males) and 450 age-, race and sex-matched healthy controls (87.3% males) were involved. All AS patients were recruited from the rheumatology clinic, Shariati Hospital, Tehran, and met the modified New York diagnostic criteria for AS . All controls had no family history of AS or other autoimmune diseases. More detailed characteristics of enrolled participants are represented in Table 1. The written informed consent was signed by all members of the case and control groups before enrolling in the study. Moreover, this study was approved by ethics committee at the Tehran University of Medical Sciences (Approval ID: IR.TUMS.MEDICINE.REC.1398.696).
Genomic DNA extraction
The blood sample was taken from participants, and the isolation of peripheral blood mononuclear cells (PBMCs) were done based on Ficoll-Hypaque gradient technique (innotrain, Germany). DNA samples were isolated from PBMCs in standard collection tube with EDTA as anticoagulants using DNSol MidiPrep DNA extraction kit (Roje Technologies Co., Iran). For evaluating the purity of DNA, the absorbance of DNA samples was measured at 260 and 280 nm using Thermo Scientific/2000C Nanodrop. Pure DNA has an A260/A280 ratio ranges from 1.8–2.0. Extracted DNA from all samples were stored at − 20 °C until copy number evaluation.
Determination of the IL17RA copy number
For analysis of IL17RA gene copy number, we used the predesigned TaqMan CNV qPCR assay (Assay ID: Hs02339506_cn), which is purchased from Applied Biosystems (Foster City, CA, USA) (http://bioinfo.appliedbiosystems.com/genomedatabase/copy-number-variation.html). IL17RA TaqMan assay includes 2 target-specific primers and a FAM/MGB dual-labeled probes which are located in Chr.22:17109553 (within exon 13). The genetic structure of IL17RA gene and the position of Hs02339506_cm TaqMan assay containing primer sets and probe illustrated in Fig. 1. Additionally, TaqMan Copy Number Reference Assay (RNase P) includes two primers and VIC/TAMRA labeled probes which used as a standard reference assay recommended by the manufacturer’s protocol. Each reaction included 5 μl of 2× TaqMan Genotyping Master Mix, 0.5 μl of IL17RA TaqMan copy number assay, 0.5 μl of TaqMan Copy Number RNAse P Reference Assay, 2.5 μl of water, and 1.5 μl of 10 ng/μl genomic DNA in a 10 μl total volume. All samples were run in duplicate in a 96 well plate, and fluorescence signals emitting from reporter dyes (FAM and VIC) were normalized to ROX (passive reference dye). Real-time PCR was done under thermal cycling program comprised of 1 holding stage of 10 min at 95 °C, followed by 40 cycles in 15 s at 95 °C and 60 s at 60 °C.
Amplification data were analyzed using Applied Biosystems StepOnePlus software (version 2.3), and raw data of CT (threshold cycle) values were imported to CopyCaller® Software v2.0 (Applied Biosystems) to calculate the gene copy number in each sample. The comparative Ct method (ΔΔCt formula) was applied to estimate the relative gene copy number. To grant validity to results, confidence values were higher than 95%, and Z-Score was < 1.75. Besides, CT of the reference gene was set < 32, and the SD of ΔCT between replicates was less than ±0.2.
mRNA expression of IL17RA
First, total RNA extraction was implemented from PBMCs of 68 cases and 53 controls with the High Pure RNA isolation kit (Roche, Germany). Total RNA concentration was determined by measuring the absorbance at 260 nm in Thermo Scientific/2000C Nanodrop. Afterward, reverse transcription (RT) was performed to synthesize cDNA from RNA by the PrimeScript RT reagent Kit (Takara, Japan). Subsequently, relative quantification of mRNA expression was conducted by SYBR Green qPCR on a StepOnePlus Real-time PCR system. Gene expression analysis was carried out with specific primers for IL17RA (Forward primer: 5′-AGTTCCACCAGCGATCCAAC-3′; Reverse primer: 5′-GTCTGAGGCAGTCATTGAGGC-3′) and for β2-microglobulin (β2M) (Forward primer: 5′-CCTGAATTGCTATGTGTCTGGG-3′; Reverse primer: 5′- TGATGCTGCTTACATGTCTCGA-3′). The mRNA expression was normalized to β2M.
The Kolmogorov–Smirnov test was implemented for computing normality of tests for continuous variables. We performed Mann-Whitney U non-parametric test to compare the gene expression data between the two groups. Correlation analysis was done by Spearman’s rank correlation. Pearson Chi-square test was used to assess the association between the categorical variables.
Furthermore, for checking the effect size odds ratio (OR) with 95% confidence intervals (95% CI) was measured. All data were analyzed using IBM SPSS version 21.0 (SPSS, Chicago, IL, USA). In addition, experimental data graphs were designed by GraphPad Prism version 6 (La Jolla, CA, USA). Results with P value less than 0.05 was regarded to be statistically significant.
Characteristics of study subjects
The demographic and clinical features of studied subjects are listed in Table 1. The mean ± SD age of patients and healthy individuals were 38.10 ± 10.46 and 36.60 ± 8.25 years, respectively. There was no statistically significant difference between two groups concerning sex (P-value = 0.37) and age (P-value = 0.86). Patients with AS had higher levels of serum ESR compared to controls (P-value < 0.001).
Analysis of IL17RA CNV
Based on our results, IL17RA copy number varied from 1 to 4 per diploid genome. As shown in Table 2, CN = 2 was the predominant copy count in both groups. One copy number of IL17RA describe gene deletion, while 3 and 4 copy counts are indicative of gene duplication relative to normal 2 copy number. The frequency of two copies of IL17RA was 76.3% in AS subjects and 85.6% in the control group. 13.4% of cases and 6.9% of controls exhibited l copy (CN < 2) of IL17RA and were remarkably significant. A high copy of CN seemed to be more in AS patients comparing to healthy subjects (10.3% vs. 7.6%), but the difference was not significant. Only 2 patients possessed 4 copies. Based on the results of chi-square test, the overall distribution of the copy-number in the cases compare to control was significantly different (P-value = 0.001).
Copy numbers were categorized into 3 groups as following: CN < 1 (losses), CN = 2 (neutral), and CN > 2 (gains). As shown in Table 2, the frequency of CN loss was significantly higher in AS patients (vs. 2 copies OR:2.18, 95% CI: (1.38–3.44), P-value < 0.001) in comparison to controls. Although there was no significant difference between the two groups for CNV gain (vs. 2 copies OR:1.53, 95% CI: (0.963–2.44), P-value = 0.071). Of note, two copies considered as the reference for estimation of AS risk.
IL17RA mRNA expression level
To address the association of IL17RA CN to the expression level of this gene, we analyzed IL17RA expression in PBMCs of 68 cases and 53 controls. The expression of IL17RA was higher in all cases compared to controls (P-value = 0.032, Fold change: AS/Control = 1.17) (Fig. 2a). Interestingly, the mRNA level displayed a significant raise in men compared to women. The detailed data are shown in Table 3 and Fig. 2b. The mRNA expression was also assessed by different copy numbers status (Table 4 and Fig. 2c). IL17RA mRNA expression fold change elevated in AS patients in 2 copies group. In contrast, there were no significant mRNA expression changes in CN loss and CN gain groups.
Correlation between copy number and mRNA expression of IL17RA with clinical features of AS
As shown in Table 5, there was no correlation between each of IL17RA copy number categories and clinical characteristics in AS patients. Moreover, IL17RA mRNA expression was not correlated significantly with clinical data (Table 5). However, BASFI and BASMI have higher non-significant effects on IL17RA mRNA expression (Table 5).
In the present study, the IL17RA CNVs and mRNA expression levels were evaluated to provide more information about its role in AS disease risk and IL-17A inflammatory pathways.
The fundamental role of IL-17A pathway in AS pathogenesis was identified by the high frequency of Th17 cells and IL-17A expression level in the peripheral blood of AS patients. Several studies showed the importance of IL-17A as a major cytokine that mediates responses in inflammation, infection, and autoimmune diseases through binding to its unique receptor (IL-17RA) and stimulates ACT1 signaling . IL-17RA signaling pathway contributes to the trigger expression of several inflammatory mediators in the active phase of inflammatory diseases .
In the recent years, the genome-wide association studies (GWASs) showed the contribution of genes outside of the MHC region in the risk of AS and CNVs received more attention due to their pathophysiological role in autoimmune diseases . A GWAS carried out by Jung et al. has determined 9 candidate genes as CNVs regions related to AS risk in the Korean population. Among other genes, the only gain in IL17RA was significantly protective against AS, and its influence is suggested in AS development . Recently, Shahba et al. reported the association signal between BMP8A copy number variation and susceptibility to AS in Iranian patients .
In line with Korean study , we showed that having only less than 2 copies of IL17RA gene significantly increased susceptibility to AS. Therefore, it suggests that there is a direct association between low copy number and AS development in both populations. These results proposed that loss of IL17RA is considered as an effective CNV in AS development in Iranian population. When stratified for sex, no sex bias was detected in CNVs categories in cases and controls. Additionally, for the first time, we examined IL17RA mRNA expression level to understand how the CNVs of the IL17RA can affect expression of the gene.
Despite the fact, we have observed the loss of gene copy number associated with disease susceptibility, and there was no significant change in mRNA expression level in cases and controls with less than 2 copies. While the increased expression of IL17RA mRNA has been observed in 2 copies group.
Therefore, a lack of significant association of IL17RA CNVs with mRNA expression levels in both groups could be attributed to epigenetic factors and modifications, dysfunction of a transcription factor, and also a specific mutation that may contribute to IL17RA upregulation in AS . DNA methylation, as one of the most epigenetic modifications, has been recognized to influence gene expression . All these strategies can naturally occur in an intuitive phenomenon called genetic compensation. It is a regulatory process that fundamentally responds to the abnormal changes in RNA or protein level .
Furthermore, it is expressed that the phenotypic effects of CNVs may have not a direct relationship to gene expression. Some biological mechanisms are proposed to make different phenotypes associated with upregulation and dysregulation of immune responses such as position effect, extra protein-coding domains, and so on [49, 50]. PBMCs isolation is a less invasive alternative method for measuring expression level as well as copy number of genes compared to the synovial biopsy isolation, therefore, we evaluated the IL17RA expression in PBMCs. The results that no significant change between the different copy number category of IL17RA gene and its mRNA level might not be true for other tissues, as the expression of this gene varies from tissue to tissue. However, a recent study has studied the genes expression profile in PBMCs of AS patients with inflamed synovial tissue. The authors found a number of immune-associated genes with similar expression change in both PBMCs and synovial tissues. This study suggested that local inflammation could show a reflection in PBMCs alterations .
Because of the possible importance of this cytokine in AS and similar conditions, as we expected, the mRNA expression level of IL17RA was significantly higher in all our patients compared to controls. Previous studies indicated increased serum IL-17 level in AS patients showing that this cytokine may play a crucial role with IL-23 in the AS pathogenesis [20, 52, 53]. Following our findings, IL17RA expression was also significantly higher in synoviocytes of RA and psoriatic arthritis (PsA) patients . Reduced production of proinflammatory mediators (e.g., IL-1, IL-6, and different MMPs) along with lessened joint inflammation has been shown in mice with IL17R-deficient streptococcal cell wall-induced arthritis .
Therefore, the upregulation of IL17RA could enhance their response to IL-17A and exacerbate inflammation of the joints . In accordance with the previous studies and the current research, it is highly possible to consider an essential role for the IL-17 family with corresponding receptors in ankylosing spondylitis.
The genetic architecture, which is an essential determinant of AS, making it a polygenic complex disease with a more heritable tendency to affect males than females [56, 57]. In our study, increased IL17RA mRNA expression is restricted to male patients, after adjusting by sex. Consistent with our results, Gracey et al. found the upregulation of IL17RA expression in male AS patients . Furthermore, there are not any significant correlation between CNVs and mRNA expression level with demographic and clinical data of patients in our study.
In conclusion, regarding the frequency distribution of copy numbers, less than 2 copies of IL17RA were significantly susceptible to AS in this cohort of Iranian patients. Epigenetic mechanisms like DNA methylation, post-transcriptional, and -translational modifications that regulate the expression of the genes may contribute in upregulation of IL17RA mRNA expression in loss of the gene copy number. On the other hand, the significant upregulation mRNA transcript level of IL17RA in all AS patients indicates the critical role of this pathway in the risk of AS. Further researches, especially in the context of functional consequences, seem to be necessary to illuminate the more details of the fundamental role of CNVs with biologic effects.
Availability of data and materials
The relevant raw data of this study are available as an Additional file 1.
Copy number variation
Genome-wide association study
Granulocyte monocyte-colony stimulating factor
Granulocyte-colony stimulating factor
Human immunodeficiency virus
Innate lymphoid cell
- NKT cells:
Natural killer T cells
Quality of life
Single nucleotide polymorphism
Systemic lupus erythematosus
Type 1 diabetes
Tam LS, Gu J, Yu D. Pathogenesis of ankylosing spondylitis. Nat Rev Rheumatol. 2010;6:399–405.
Maksymowych WP. Disease modification in ankylosing spondylitis. Nat Rev Rheumatol. 2010;6:75–81.
Fallahi S, Mahmoudi M, Nicknam MH, et al. Effect of HLA-B*27 and its subtypes on clinical manifestations and severity of ankylosing spondylitis in Iranian patients. Iran J Allergy Asthma Immunol. 2013;12:321–30.
Karami J, Mahmoudi M, Amirzargar A, et al. Promoter hypermethylation of BCL11B gene correlates with downregulation of gene transcription in ankylosing spondylitis patients. Genes Immun. 2017;18:170–5.
Aslani S, Mahmoudi M, Garshasbi M, et al. Evaluation of DNMT1 gene expression profile and methylation of its promoter region in patients with ankylosing spondylitis. Clin Rheumatol. 2016;35:2723–31.
Vanaki N, Aslani S, Jamshidi A, et al. Role of innate immune system in the pathogenesis of ankylosing spondylitis. Biomed Pharmacother. 2018;105:130–43.
Sieper J, Poddubnyy D. Axial spondyloarthritis. Lancet. 2017;390:73–84.
Dashti N, Mahmoudi M, Aslani S, et al. HLA-B*27 subtypes and their implications in the pathogenesis of ankylosing spondylitis. Gene. 2018;670:15–21.
Rezaei R, Mostafaei S, Aslani S, et al. Association study between killer immunoglobulin-like receptor polymorphisms and ankylosing spondylitis disease: an updated meta-analysis. Int J Rheum Dis. 2018;21:1746–55.
Almasi S, Aslani S, Poormoghim H, et al. Gene expression profiling of toll-like receptor 4 and 5 in peripheral blood mononuclear cells in rheumatic disorders: Ankylosing spondylitis and rheumatoid arthritis. Iran J Allergy Asthma Immunol. 2016;15:87–92.
Ranganathan V, Gracey E, Brown MA, et al. Pathogenesis of ankylosing spondylitis — recent advances and future directions. Nat Rev Rheumatol. 2017;13:359.
Beringer A, Noack M, Miossec P. IL-17 in chronic inflammation: from discovery to targeting. Trends Mol Med. 2016;22:230–41.
Jiang YX, Li PA, Yang SW, et al. Increased chemokine receptor IL-17RA expression is associated with poor survival in gastric cancer patients. Int J Clin Exp Pathol. 2015;8:7002–8.
Shen F, Gaffen SL. Structure-function relationships in the IL-17 receptor: implications for signal transduction and therapy. Cytokine. 2008;41:92–104.
Gaffen SL, Kramer JM, Yu JJ, et al. The IL-17 cytokine family. Vitam Horm. 2006;74:255–82.
Chang SH, Dong C. Signaling of Interleukin-17 family cytokines in immunity and inflammation. Cell Signal. 2011;23:1069–75.
Rezaiemanesh A, Abdolmaleki M, Abdolmohammadi K, et al. Immune cells involved in the pathogenesis of ankylosing spondylitis. Biomed Pharmacother. 2018;100:198–204.
Gu C, Wu L, Li X. IL-17 family: cytokines, receptors and signaling. Cytokine. 2013;64:477–85.
Bidad K, Salehi E, Jamshidi A, et al. Effect of all-transretinoic acid on Th17 and T regulatory cell subsets in patients with ankylosing spondylitis. J Rheumatol. 2013;40:476–83.
Mei Y, Pan F, Gao J, et al. Increased serum IL-17 and IL-23 in the patient with ankylosing spondylitis. Clin Rheumatol. 2011;30:269–73.
Jandus C, Bioley G, Rivals JP, et al. Increased numbers of circulating polyfunctional Th17 memory cells in patients with seronegative spondylarthritides. Arthritis Rheum. 2008;58:2307–17.
Gaffen SL, Jain R, Garg AV, et al. The IL-23-IL-17 immune axis: from mechanisms to therapeutic testing. Nat Rev Immunol. 2014;14:585–600.
Goepfert A, Lehmann S, Wirth E, et al. The human IL-17A/F heterodimer: a two-faced cytokine with unique receptor recognition properties. Sci Rep. 2017;7:8906.
Song X, Qian Y. IL-17 family cytokines mediated signaling in the pathogenesis of inflammatory diseases. Cell Signal. 2013;25:2335–47.
Ishigame H. Differential roles of interleukin-17A and -F in host defense against mucoepithelial bacterial infection and allergic responses. Immunity. 2009;30:108–19.
Gaffen SL. Structure and signalling in the IL-17 receptor family. Nat Rev Immunol. 2009;9:556–67.
Kehlen A, Thiele K, Riemann D, et al. Expression, modulation and signalling of IL-17 receptor in fibroblast-like synoviocytes of patients with rheumatoid arthritis. Clin Exp Immunol. 2002;127:539–46.
Abdolmohammadi K, Pakdel FD, Aghaei H, et al. Ankylosing spondylitis and mesenchymal stromal/stem cell therapy: a new therapeutic approach. Biomed Pharmacother. 2019;109:1196–205.
Soleimanifar N, Amirzargar AA, Mahmoudi M, et al. Study of programmed cell death 1 (PDCD1) gene polymorphims in Iranian patients with ankylosing spondylitis. Inflammation. 2011;34:707–12.
Daryabor G, Mahmoudi M, Jamshidi A, et al. Determination of IL-23 receptor gene polymorphism in Iranian patients with ankylosing spondylitis. Eur Cytokine Netw. 2014;25:24–9.
Mahmoudi M, Jamshidi AR, Amirzargar AA, et al. Association between endoplasmic reticulum aminopeptidase-1 (ERAP-1) and susceptibility to ankylosing spondylitis in Iran. Iran J Allergy Asthma Immunol. 2012;11:294–300.
Momenzadeh P, Mahmoudi M, Beigy M, et al. Determination of IL1 R2, ANTXR2, CARD9, and SNAPC4 single nucleotide polymorphisms in Iranian patients with ankylosing spondylitis. Rheumatol Int. 2016;36:429–35.
Vidal-Castineira JR, Lopez-Vazquez A, Diaz-Pena R, et al. A single nucleotide polymorphism in the Il17ra promoter is associated with functional severity of Ankylosing spondylitis. PLoS One. 2016;11:e0158905.
Usher CL, McCarroll SA. Complex and multi-allelic copy number variation in human disease. Brief Funct Genomics. 2015;14:329–38.
Redon R, Ishikawa S, Fitch KR, et al. Global variation in copy number in the human genome. Nature. 2006;444:444–54.
Stankiewicz P, Lupski JR. Structural variation in the human genome and its role in disease. Annu Rev Med. 2010;61:437–55.
Gamazon ER, Stranger BE. The impact of human copy number variation on gene expression. Briefings Funct Genomics. 2015;14:352–7.
Bigagli E, De Filippo C, Castagnini C, et al. DNA copy number alterations, gene expression changes and disease-free survival in patients with colorectal cancer: a 10 year follow-up. Cell Oncol (Dordr). 2016;39:545–58.
Grayson BL, Smith ME, Thomas JW, et al. Genome-wide analysis of copy number variation in type 1 diabetes. PLoS One. 2010;5:e15393.
Jung SH, Yim SH, Hu HJ, et al. Genome-wide copy number variation analysis identifies deletion variants associated with ankylosing spondylitis. Arthritis Rheumatol. 2014;66:2103–12.
Zochling J. Measures of symptoms and disease status in ankylosing spondylitis: Ankylosing Spondylitis Disease Activity Score (ASDAS), Ankylosing Spondylitis Quality of Life Scale (ASQoL), Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), Bath Ankylosing Spondylitis Functional Index (BASFI), Bath Ankylosing Spondylitis Global Score (BAS-G), Bath Ankylosing Spondylitis Metrology Index (BASMI), Dougados Functional Index (DFI), and Health Assessment Questionnaire for the Spondylarthropathies (HAQ-S). Arthritis Care Res (Hoboken). 2011;63(Suppl 11):S47–58.
Qian Y, Liu C, Hartupee J, et al. The adaptor Act1 is required for interleukin 17-dependent signaling associated with autoimmune and inflammatory disease. Nat Immunol. 2007;8:247–56.
Sadik CD, Kim ND, Alekseeva E, et al. IL-17RA signaling amplifies antibody-induced arthritis. PLoS One. 2011;6:e26342.
Olsson LM, Holmdahl R. Copy number variation in autoimmunity--importance hidden in complexity? Eur J Immunol. 2012;42:1969–76.
Shahba S, Jafari Shakib R, Jamshidi A, et al. Association study of copy number variation in BMP8A gene with the risk of ankylosing spondylitis in Iranian population. J Cell Biochem. 2018.
Myhre S, Lingjaerde OC, Hennessy BT, et al. Influence of DNA copy number and mRNA levels on the expression of breast cancer related proteins. Mol Oncol. 2013;7:704–18.
Mahmoudi M, Aslani S, Nicknam MH, et al. New insights toward the pathogenesis of ankylosing spondylitis; genetic variations and epigenetic modifications. Mod Rheumatol. 2017;27:198–209.
Veitia RA, Bottani S, Birchler JA. Gene dosage effects: nonlinearities, genetic interactions, and dosage compensation. Trends Genet. 2013;29:385–93.
Hollox EJ, Hoh BP. Human gene copy number variation and infectious disease. Hum Genet. 2014;133:1217–33.
Henrichsen CN, Chaignat E, Reymond A. Copy number variants, diseases and gene expression. Hum Mol Genet. 2009;18:R1–8.
Thomas GP, Duan R, Pettit AR, et al. Expression profiling in spondyloarthropathy synovial biopsies highlights changes in expression of inflammatory genes in conjunction with tissue remodelling genes. BMC Musculoskelet Disord. 2013;14:354.
Wang X, Lin Z, Wei Q, et al. Expression of IL-23 and IL-17 and effect of IL-23 on IL-17 production in ankylosing spondylitis. Rheumatol Int. 2009;29:1343–7.
Chen W-S, Chang Y-S, Lin K-C, et al. Association of serum interleukin-17 and interleukin-23 levels with disease activity in Chinese patients with ankylosing spondylitis. J Chinese Med Assoc. 2012;75:303–8.
Raychaudhuri SP, Raychaudhuri SK, Genovese MC. IL-17 receptor and its functional significance in psoriatic arthritis. Mol Cell Biochem. 2012;359:419–29.
Koenders MI, Kolls JK, Oppers-Walgreen B, et al. Interleukin-17 receptor deficiency results in impaired synovial expression of interleukin-1 and matrix metalloproteinases 3, 9, and 13 and prevents cartilage destruction during chronic reactivated streptococcal cell wall–induced arthritis. Arthritis Rheumatism. 2005;52:3239–47.
Lee W, Reveille JD, Davis JC, et al. Are there gender differences in severity of ankylosing spondylitis? Results from the PSOAS cohort. Ann Rheum Dis. 2007;66:633–8.
Shahlaee A, Mahmoudi M, Nicknam MH, et al. Gender differences in Iranian patients with ankylosing spondylitis. Clin Rheumatol. 2015;34:285–93.
Gracey E, Yao Y, Green B, et al. Sexual dimorphism in the Th17 signature of Ankylosing spondylitis. Arthritis Rheumatol. 2016;68:679–89.
This work was financially supported by the grant from the research deputy of Tehran University of Medical Sciences (Grant No. 95–01–30-30949). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The written informed consent was signed by all members of the case and control groups before enrolling in the study. Moreover, this study was performed based on the Declaration of Helsinki guidelines and was approved by the ethics committee at the Tehran University of Medical Sciences (Approval ID: IR.TUMS.MEDICINE.REC.1398.696).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Aghaei, H., Farhadi, E., Akhtari, M. et al. Copy number variation of IL17RA gene and its association with the ankylosing spondylitis risk in Iranian patients: a case-control study. BMC Med Genet 21, 147 (2020). https://doi.org/10.1186/s12881-020-01078-y
- Ankylosing spondylitis
- Copy number variations
- Gene expression