Single Nucleotide Polymorphism of TYK2 Gene and Susceptibility to Rheumatoid Arthritis in Iranian Population

- - - -

PMID: 31057722 (PubMed) - PMCID: PMC6490405 - View online: PubReader
Volume 11, Issue 2, April-June 2019, Page 187 to 191
Wednesday, April 11, 2018 :Received , Monday, June 11, 2018 :Accepted

• Mohamadhosseini, Azadeh
  - Department of Immunology, Faculty of Medicine, Shahid Sadoughi University of Medical Sciences, Yazd, Iran

• Mansouri, Reza Immunology Department, Faculty of Medicine, Shahid Sadoughi University of Medical Sciences, Yazd, Iran, Tel/Fax: +98 2188220067, +98 3538203410, E-mail: rmm542003@yahoo.com
  - Department of Immunology, Faculty of Medicine, Shahid Sadoughi University of Medical Sciences, Yazd, Iran

• Javinani, Ali
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

• Ashraf-Ganjouei, Amir
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

• Akhlaghi, Massoumeh
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

• Aslani, Saeed
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

• Hamzeh, Elham
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

• Jamshidi, Ahmadreza
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

• Ahmadzadeh, Nooshin
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

• Mahmoudi, Mahdi Rheumatology Research Center, Shariati Hospital, Tehran, Iran, Tel/Fax: +98 2188220067, +98 3538203410, E-mail: mahmoudim@tums.ac.ir
  - Rheumatology Research Center, Tehran University of Medical Sciences, Tehran, Iran

Rheumatoid Arthritis (RA) is an auto-immune disease characterized by inflammatory cascades primarily happening in joints but the other organs, albeit less so, can be affected ¹. According to the World Health Organization (WHO) regions classification, the prevalence of RA was reported in a range from 0.37% in Eastern Mediterranean region to 1.25% in American region ². RA exhibits a strong gender predilection with the female to male ratio of 4.7:1 in low and middle-income countries ². RA has the estimated disability-adjusted life years of 4.8 million in 2010 which was ranked as the 42^nd contributor to the global disability ³. Thus, according to its relatively high prevalence and morbidity, RA is one of the major fields of interests in Rheumatology and Immunology.
In RA, arthritis-associated antigens are presented by dendritic cells to naïve CD4+ T cells ⁴. Following their activation and differentiation toward the T Helper (TH) 17 cells, the immune system orchestrates an inflammatory cascade with a high spectrum of cytokines production. These cytokines affect a wide range of receptors, which are linked to the signaling mediators from Janus Kinase (JAK) family ⁵. JAK family has four major members; JAK1, JAK2, JAK3,and tyrosine kinase 2 (TYK2) ⁶. They are non-receptor tyrosine kinases, which are activated after the binding of cytokines to their relative receptors. Following their activation, they make a docking site for signaling molecules, mostly members of the Signal Transducer and Activator of Transcription (STAT) family ⁶. STAT molecules make a dimer after phosphorylation which then enter the nucleus and induce the expression of further molecules of the inflammatory cascade ⁵.
TYK2 is linked to the cytoplasmic domain of type I and II cytokine receptors, which are involved in the signaling pathway of interleukin (IL)-6, IL-12 and IL-23 ⁵. Its gene is located in 19p13.2 position andits polymorphismshave been linked to a variety of auto-immune diseases ⁷. It has been reported previously that TYK2 gene polymorphisms influenced the susceptibility to RA ⁸, Systemic Lupus Erythematosus (SLE) ⁹, systemic sclerosis ¹⁰, Inflammatory Bowel Disease (IBD) ¹¹ and Multiple Sclerosis (MS) ¹². The present study aimedto investigate the rs34536443 Single Nucleotide Polymorphism (SNP) of TYK2 gene in Iranian RA patients. Also, its association with demographic and serological features of the patients was described as well.

Study population: This study was conducted on 700 RA patients from the outpatient’s rheumatology clinic of Shariati   Hospital who were diagnosed based on the American College of Rheumatology (ACR)/European League Against Rheumatism (EULAR) 2010 RA classification criteria 13. Control group of this study composed of 700 age, sex, and ethnicity matched healthy individuals. Written informed consent was obtained from all the participants. The present study was approved by the ethics committee of Yazd University of Medical Sciences.
DNA preparation and polymorphism analysis: For all participants, peripheral blood was collected into tubes with EDTA as an anticoagulant and Erythrocyte Sedimentation Rate (ESR) was measured. Genomic DNA was extracted using standard proteinase K digestion- phenol/chloroform extraction method ¹⁴. The extracted DNA samples were stored at -20°C. Approximately, 20 ng of the genomic DNA of each sample was used for genotyping. Genotyping of TYK2 gene rs34536443 was performed using the MGB-TaqMan Allelic Discrimination method (Applied Biosystems, Foster City, CA, USA). Amplification was performed in 10 μl reaction volumes, containing 5 μl of the TaqMan Genotyping master mix, 0.25 μl of TaqMan Genotyping assay mix, 0.25 μl of distilled water, and 4.5 μl of genomic DNA.
Serologic analysis: The plasma samples were obtained after centrifuging the patients׳ whole blood samples, which were drained previously into EDTA-containing test tubes. The Rheumatoid Factor (RF) and anti-Cyclic CitrullinatedPeptide antibody (anti-CCP), and C-Reactive Protein (CRP) were measured by commercially available kits. 
Statistical analysis: Control group was tested to be in the Hardy-Weinberg Equilibrium. The frequency of alleles and genotypes were compared between RA patients and control group using chi-square test. The mean comparison of ESR level, RF and anti-CCP titers of different genotypes in case group was done by ANOVA. p-values less than 0.05 were considered statistically significant. All analyses were performed using SPSS package for Windows (Version 19.0, IBM SPSS Inc., USA).

The demographic and clinical characteristics of our study population are shown in table 1. The mean age of the RA patients was 50.32±11.27 years whereas the mean age of control group was 48.44±12.37 years. The female to male ratio was 3.48:1 in patients group and 6.52:1 in healthy individuals. The mean ESR level was higher in patients (20.64±17.14 vs. 10.27±7.95), but the healthy individuals had a higher percentage of positive CRP (61.1 vs. 46%). The means of RF and anti-CCP titer were 85.5±7.63 and 138.7±18.96 international unit (IU)/ml, respectively among RA patients.
The allele and genotype frequencies of rs34536443 in RA patients and control group are shown in table 2. The frequency of minor allele (C allele) was not different among patients and healthy individuals (1.7 vs. 2.61%, OR (95% CI)=1.35 (0.78-2.33), p-value= 0.27). The frequencies of CC and GC genotypes were negligible; however, according to our analyses, there was not a statistically significant difference between the genotypic distribution of the case and control groups. Moreover, the dominant model of CC+CG had no significant differences between patients and controls. Finally, the association of genotype frequencies with clinical characteristics of RA patients is shown in table 3. There was not any correlation between genotype frequencies and patients' age and gender, ESR level, and the titer of RF and anti-CCP.

In the present study, an attempt was made to investigate the association of TYK2 gene rs34536443 polymorphism with RA susceptibility. According to our results in Iranian population, rs34536443 polymorphisms do not have any influence on RA development. Furthermore, neither demographic features nor RF and the anti-CCP titer was correlated with rs34536443 polymorphisms.
The importance of TYK2 gene polymorphisms in autoimmune diseases was supported by a vast majority of the studies. The previously studied works on TYK2 gene polymorphisms in RA are shown briefly in table 4. According to a meta-analysis published in 2011, the rs34536443 and rs2304256 SNPs were associated with auto-immune inflammatory disorders including RA, SLE, MS and IBD 15. But, at that time, there was not any published evidence on rs34536443 polymorphism and RA susceptibility. In 2012, the Genome-Wide Association Study (GWAS) on European RA patients, have reported the association of rs34536443 with RA for the first time ⁸. According to results from the Immunochip array, G allele of rs34536443 had a protective effect on RA development. After regrouping the patients based on anti-CCP status, the association remained significant only in seropositive individuals. These findings were also replicated for combined analysis of Immunochip results and imputed GWAS data from non-overlapping samples. In line with mentioned discoveries, Diogo et al showed that the three protein-coding variants in TYK2 gene, including P1104A (rs34536443) have a protective effect against RA in European patients ¹⁶. Additionally, they showed that the same three variants also protect against SLE and IBD. Moreover, based on the phenome-wide association study, they reported that the P1104A was not associated with non-autoimmune phenotypes. Their findings can support further the role of TYK2 targeting as the possible promising treatment for these auto-immune diseases.
According to the analyses of the Immunochip array on 2816 European patients with Juvenile Idiopathic Arthritis (JIA), it was showen that G allele of rs34536443 had a protective effect against oligoarticular and RF-negative polyarticular JIA ¹⁷. In contrast to this study, no association was found between TYK2 gene rs34536443 SNPs and JIA development either at an allelic or genotypic level among Greek patients ¹⁸. Unexpectedly, in a study performed on Greece adult RA patients, Myrthianou et al revealed that the GC genotype of rs34536443 SNP had a protective role against RA with an OR of 0.33 (0.123-0.90) 19.
In a functional point of view, it has been shown previously that the rs34536443 SNP, located in protein-coding region of TYK2 gene (Assembly: GRCh37.p13, annotation release: 105, Chromosome: 19, Location: NC_000019.9), can influence the three-dimensional structure of TYK2, which can consequently affect its function ¹⁸. Moreover, studies on MS patients revealed that the protective C allele of rs34536443 SNP mediates the divergence of T lymphocytes differentiation toward a TH2 phenotype ²⁰. They showed that the downstream pathways of TYK2, which were activated by interferon (IFN)-β, IL-6 and IL-10were less pronounced in carriers of GC genotype compared to GG carriers. In addition, they reported that the GATA3 expression was higher in GC genotype carriers, which can favor the TH2 modulation pathway. Besides, the secretion level of IL-5 and IL-13 was also higher in GC carrier T cells upon activation by CD3/CD28 signaling. Nonetheless, no correlation between rs34536443 polymorphisms and M1/M2 balance of macrophages phenotype was identified in their work. Although these data are governed by MS, according to the crucial role of TH1 and TH17 in RA pathophysiology 21, this scenario might also be implicated in RA. Confusingly, a certain piece of evidence described that despite the catalytic impairment of TYK2 P1104A, it has a comparable IFN-α signaling potency relative to the wild-type protein ²². 
According to the intimate contribution of JAK family proteins in inflammatory pathways, there was a great interest to target them for the treatment of inflammatory disorders. Tofacitinib, a monoclonal antibody, is a JAK inhibitor approved for RA treatment with long-term safety and efficacy ²³. Nevertheless, Tofacitinib targets more specifically the JAK3 and JAK1 with less affinity for TYK2 ²⁴. Moreover, it was described that the Tofacitinib could bind to the adenosine triphosphate-binding cavity of TYK2, which can further diminish its potential to activate the STAT molecules ²⁵. Recent studies have indicated promising outcomes about the efficacy of TYK2 inhibitors in SLE and IBD treatment via blocking the IL-12, IL-23, and IFN-I pathways in animal models with short-term safety profile in human ^26,27. 

In general, JAK family has a key role in mediating the inflammatory responses in auto-immune diseases. TYK2 is one of the family members that its gene polymorphisms were associated with RA susceptibility in different populations (Table 4). However, in the present work, no association was found between the rs34536443 polymorphism and RA. These findings can further support the determinative role of ethnicity in the etiopathogenesis of such diseases.

The authors are grateful of the patients for their participation in the study.

None.

Larger Image Table 1. Demographic and clinical data of RA patients and healthy controls

Larger Image Table 2. Allele and genotype distribution of TYK2 gene rs34536443 SNP in RA patients and healthy controls RA: Rheumatoid Arthritis; HWE: Hardy–Weinberg Equilibrium

Larger Image Table 3. Association of genotypes of TYK2 gene rs34536443 SNP with characteristics of the RA patients * The data was available for only one patient with CC genotype for TYK2 gene rs34536443 SNP.

Larger Image Table 4. Association of TYK2 gene polymorphism with RA susceptibility in different populations JIA: Juvenile Inflammatory Arthritis