1. Introduction
Rheumatoid arthritis (RA) is a systemic, chronic, and progressive inflammatory disease. This disease mainly affects multiple joints as well as exhibits extra-articular manifestations such as rheumatoid nodules, systemic comorbidities, and pulmonary involvement of vacuities. The therapeutic revolution in the past few decades that includes the advancement of the new therapeutic options, the development of new criteria for classification, the introduction of early therapy, and the application of new effective strategies has transformed the articular and systemic outcomes in RA [Reference Aletaha1, Reference Smolen and Aletaha2]. The modern genetic technologies combined with the large and well-characterized clinical cohorts have advanced our understanding of the genetics of this disease. Genome-wide association studies (GWAS) have reported more than a hundred loci associated with the RA risk [Reference Roberson and Bowcock3–Reference Zhu5]. The presence of genetic polymorphisms in HLA (human leukocyte antigen) and non-HLA genes contributes about 60% to the pathogenesis of RA. Out of HLA genes, the HLA-DRB1 and associated polymorphisms are crucial in genetic susceptibility of RA, whereas other important non-HLA genes are PTPN22, IL23R, TRAF1, CTLA4, IRF5, AFF3, STAT4, CCR6, and PADI4, etc. [Reference Goëb, Dieudé and Vittecoq6–Reference Plenge and Rioux9].
The AFF3 (lymphoid nuclear protein related to the AF4 gene) is located at 2q11.2 position and encodes a protein of 1227 amino acids. This gene is highly expressed in the lymphoid tissue and has been suggested to be involved in its development, while a lower level of expression has been reported in other tissues like the brain and lungs [Reference Aletaha1, Reference Hiwatari, Taki and Taketani10–Reference Viatte and Barton12]. This gene forms a nuclear factor that binds to DNA through its transcriptional activation domain. This property of the AFF3 gene has made it a strong candidate for autoimmunity in the human cell [Reference Barton, Eyre and Ke13–Reference Orozco and Barton15]. A large number of studies, including the genome-wide association studies and their meta-analysis, conducted in different populations, have investigated the association of AFF3 rs10865035, rs1160542, rs17023158, and rs1437377 variants in the pathogenesis of RA [Reference Barton, Eyre and Ke13, Reference Hinks and Eyre16–Reference Plant, Flynn and Mbarek20].
The AFF3 gene is not only important in the susceptibility to RA but also is a good candidate in the assessment of the therapeutic response in individuals [Reference Viatte and Barton12]. The SNPs in any population are helpful in understanding its genetic diversity and to design a population-specific therapy for RA. Previously, European RA risk loci including AFF3 rs10865035 were replicated in a small sample set of Pakistani individuals [Reference Jalil, Bhatti and Demirci7]. However, in this study, the association of AFF3 rs10865035 with RA was investigated in another sample set of Pakistani origin and was merged with our previous data for statistical significance.
2. Methodology
2.1. Study Subjects
A cohort of 703 individuals of Pakistani origin (409 RA cases and 294 healthy controls) were included in this study after getting their informed consent. The study cohort was divided into group A (n = 203) and group B (n = 500). The dataset of group A was taken from our previous study [Reference Jalil, Bhatti and Demirci7]. The details of group B are given in Table 1. The cases were clinically diagnosed by certified rheumatologists following American College of Rheumatology (ACR) criteria [Reference Aletaha1], and the data including their clinical and demographic information were recorded on a specially designed questionnaire. A comparison of the RA cases and controls showed that the level of inflammatory markers like rheumatoid factor (RF) and erythrocyte sedimentation rate (ESR) was significantly higher in cases. The participants having no immunological disease/symptoms were counted as controls in this study. The study protocol was approved by the Ethical Review Board of LRH and AWKUM (Abdul Wali khan University Mardan).
∗ Data obtained from our previous study [Reference Jalil, Bhatti and Demirci7].
2.2. DNA Extraction and Genotyping
Whole blood samples were collected and processed for group B (n = 500) as previously done by Jalil et al., 2013, for group A (n = 203). The genomic DNA was extracted using an organic phenol-chloroform method [Reference Sambrook, Fritsch and Maniatis21]. The genotyping of the AFF3 gene variant rs10865035 was carried out using TaqMan genotype SNP assay (Applied Biosystems) and ARMS-PCR techniques. The individuals in group A (n = 203) were analyzed through TaqMan assay where the PCR amplification was performed in 384-well plates on dual block GeneAmp PCR system 9700 (Applied Biosystems), following the manufacture’s protocol [Reference Jalil, Bhatti and Demirci7]. The remaining 500 samples were analyzed through ARMS-PCR using a set of three primers, including two forward and a common reverse primer (F1 (for A allele): TTTAAAACCTCTATCTGGGGAAAAA, F2 (for G allele): TAAAACCTCTATGGGGAAAAG and R: CCCCTCTAATAGTCAATCAATCAAAATA).
The PCR amplification was performed in 96-well plates on a thermocycler (T100, Bio-Rad). The amplification conditions were set as the initial denaturation at 94°C for 5 minutes, followed by 35 cycles of denaturation at 94°C for 30 seconds, annealing for 30 seconds at 58°C, and extension at 72°C for 1 minute. The final elongation was performed at 72°C for 5 minutes. The PCR products were resolved in 2% agarose gel, and the genotype calls for each subject were recorded using the visual inspection method of the gel.
2.3. Statistical Analysis
The association of AFF3 rs10865035 with RA was tested through different statistical models such as codominant, homozygous dominant, homozygous recessive, heterozygous, and additive. The data were further checked for deviation from the Hardy–Weinberg Equilibrium (HWE), and the error in genotyping was measured by repeating 10% of the samples. The association of AFF3 rs10865035 with RA was measured by calculating the odds ratio (OR) at 95% confidence interval (CI) using chi-square (χ 2) and Fisher’s exact test. The P value < 0.05 was considered statistically significant.
3. Results
The data of the current study (n = 500) were combined with our previous findings (n = 203) to make the sample size more effective. In the codominant model, the frequency of the genotypes was found almost similar in both the cases and controls (cases: A/A 133 (32.52%), A/G 210 (51.35%), and G/G 66 (16.14%); controls: A/A 105 (35.72%), A/G 159 (54.09%), and G/G 30 (10.21%), χ 2 = 5.169; P = 0.07). Similarly, in the homozygous dominant model, no significant difference was observed in the distribution of A/A vs. G/G + A/G genotypes in the study subjects (OR = 0.867 (0.636–1.187); P = 0.41). However, a significant distribution was determined in the homozygous recessive model by comparing G/G vs. A/G + A/A (OR = 1.693 (1.06–2.648); Pe = 0.02). Furthermore, an insignificant distribution of A/G vs. A/A + GG genotypes was observed in the heterozygous model (OR = 0.491 (0.667–1.215); P = 0.49) whereas in the additive model, the alleles (A vs. G) were distributed insignificantly among the cases and controls (OR = 0.826 (0.665–1.027); P = 0.08). The association results, P values along with ORs, are given in Table 2.
4. Discussion
Rheumatoid arthritis (RA) is an autoimmune disease that disrupts the normal physical activities of the patients [Reference Firestein22, Reference McInnes and Schett23]. It causes inflammation in joints which results in damaging the articular cartilage along with the synovial hyperplasia that ultimately leads to consistent pain and permanent disability. In RA, the genetic risk factors (50–60%) are triggered by various environmental factors (40%) [Reference McInnes and Schett23–Reference Klareskog, Stolt and Lundberg25]. Among the genetic factors, the major histocompatibility complex (MHC) and non-MHC loci are accountable for ∼23% of the genetic risk which indicates that there are still undiscovered risk alleles [Reference Stahl, Raychaudhuri and Remmers18]. Furthermore, the association studies have also reported >150 SNPs located at more than 70 gene loci [Reference Stahl, Raychaudhuri and Remmers18, Reference Kim, Bang, Lee and Bae26–Reference Mohammadi, Aslani, Mostafaei, Jamshidi, Riahi and Mahmoudi28]. Another gene-based association study that was conducted on European (14,361 cases and 43,923 controls) and Asian populations (4,873 cases and 17,642 controls) identified 221 newly RA-associated genes including 76 European-specific, 74 Asian-specific, and 71 genes which were found overlapped among both the populations [Reference Zhu5]. Other studies that followed different study protocols have either validated the role of some of the reported genes or have determined the association of novel genes and polymorphisms in both European and Asian individuals [Reference Aletaha1, Reference Jalil, Bhatti and Demirci7, Reference Prasad17, Reference Stahl, Raychaudhuri and Remmers18, Reference Arya29, Reference Saad, Mabrouk, Eldeib and Shaker30].
The aim of this study was to assess the association of the AFF3 rs10865035 with the genetic background of RA in the Pakistani population. The replication of known genetic polymorphism will enhance our understanding of the ethnogenetic heterogeneity and homogeneity of the Pakistani population with others. The understanding of replication patterns among different populations would help us in the development of a common therapy for disease management.
Previously, we determined an insignificant association (P = 0.117) of AFF3 gene polymorphism rs10865035 in a small sample set (n = 213) of Pakistani origin [Reference Jalil, Bhatti and Demirci7]. Therefore, the current study was aimed to replicate these findings in a larger sample (n = 500) set of the same ethnicity. To enhance the statistical significance of this study, our previous data of 203 individuals were merged with the current data of 500 subjects. These statistical analyses showed that AFF3 gene polymorphism rs10865035 has no association with RA in the Pakistani population (Table 2).
The AFF3 gene encodes a DNA interacting nuclear factor having a transcriptional activation domain that makes it a strong candidate gene in autoimmunity. Previously, more than 50 genetic polymorphisms have been mapped in the same gene [Reference Barton, Eyre and Ke13]. Also, an extended study was conducted on a sample set of 6819 RA patients and 12650 healthy controls, provided a convincing proof for the association of AFF3 gene variants rs10865035 (OR: 1.12 (1.07–1.17); P = 2.8 × 10−7) and rs1160542 (OR: 1.12 (1.05–1.20); P = 0.001) as novel susceptibility targets [Reference Barton, Eyre and Ke13]. Using these findings, another group of researchers investigated the role of AFF3 gene polymorphisms (rs10865035 and rs1160542) with respect to anti-TNF treatment in RA (coefficient −0.14 (95% CI −0.25 to −0.03), P = 0.01) [Reference Tan, Gibbons and Potter31]. Furthermore, a GWAS meta-analysis by Stahl et al. (2010) also confirmed the association of AFF3 gene polymorphism rs10865035 (OR: 1.12 (1.07–1.17); P = 2.0 × 10−6) with RA in a sample set of 41,282 individuals (12,307 RA cases and 28,975 controls) of European descent [Reference Stahl, Raychaudhuri and Remmers18].
In order to explore population-specific as well as European-overlap susceptibility loci/genes, a GWAS was carried out in the Korean population. In this study, 1519 RA cases and 1476 healthy controls were genotyped for 441,398 single-nucleotide polymorphisms, and along with many other susceptibility signals, AFF3 rs10865035 was successfully associated with the RA at genome-wide significance level (P < 5 × 10−08) [Reference Freudenberg, Lee and Han19]. Another replication study, conducted by Prasad et al. [Reference Prasad17] on 983 RA cases and 1007 controls of North Indian descent, established the association of seven genes with the RA. They selected 42 candidate genes/loci (3 Asian and 39 European) and tested 603 SNPs, which were either index SNPs or were surrogate SNPs, using Infinium Human 660w-quad microarray methods for genotyping. Out of 12 SNPs, tested in AFF3 gene, 3 SNPs were found to be significantly associated with RA (rs17023158, OR: 1.45 (1.12–1.88), P = 0.005; rs6706188, OR: 0.81 (0.7–0.94), P = 0.005; rs1437377, OR: 0.77 (0.65–0.92), P = 0.003) [Reference Prasad17]. However, rs10865035 was not included in their SNPs array. Similarly, a sequenom mass array-based meta-analysis conducted on European individuals (3311 RA cases and 3709 controls) investigated another SNP rs1160542 and obtained a significant association with RA (OR: 1.08 (1.01–1.16), P = 0.029) [Reference Plant, Flynn and Mbarek20]. In addition, the AFF3 gene has been detected to be significantly associated (OR: 1.25 (1.13–1.39); P = 2.05 × 10−5) with juvenile idiopathic arthritis (JIA) as well [Reference Hinks and Eyre16]. It was further reported that, besides, its role in RA, the AFF3 gene shares the same genetic basis with the systemic lupus erythematosus (SLE). A Chinese case-control study (868 SLE patients and 975 controls) observed significant association of the AFF3 rs10865035 with SLE (OR: 1.26 (1.11–1.44); P = 4.81 × 10−4) [Reference Cen, Leng and Wang32]. Furthermore, studies have shown that the AFF3 gene was initially linked to type 1 diabetes (T1D), and then, later, its role was determined in 16 different autoimmune diseases including RA and SLE [Reference Ramos33, Reference Todd and Walker34].
In conclusion, our previous study [Reference Jalil, Bhatti and Demirci7] regarding the association of AFF3 rs10865035 with RA was further validated using different statistical models in a larger sample set (n = 703). This reaffirmation with different statistical tools suggests the nonsignificance of AFF3 rs10865035 in a Pakistani cohort. However, there is a possibility of the existence of ethnic variability and we urge the need for multiethnic large population-based study in order to understand the exact mechanism of pathogenicity as well as the evolutionary background of the genetic factors involved.
Data Availability
All the data used to support the findings of this study are cited in the text.
Disclosure
Yasir Ali and Suleman Khan are shared first authors of this article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
YA analyzed the data and participated in manuscript writing. YA, SK, NF, and AA collected blood samples and clinical data. YA, MA, and AA performed the laboratory experiments. FJ being a supervisor at AWKUM, Pakistan, helped in writing and editing the manuscript. YC, ZI, AAS, and MJ edited the manuscript. All authors read and approved the final manuscript.
Acknowledgments
The authors are highly thankful to the study subjects who participated in this study. The study was sponsored by Higher Education Commission, Pakistan, and School of Biomedical Sciences, The Chinese University of Hong Kong.