<i>TMEM187-IRAK1</i> Polymorphisms Associated with Rheumatoid Arthritis Susceptibility in Tunisian and French Female Populations: Influence of Geographic Origin

Khalifa, Olfa; Balandraud, Nathalie; Lambert, Nathalie; Auger, Isabelle; Roudier, Jean; Sénéchal, Audrey; Geneviève, David; Picard, Christophe; Lefranc, Gérard; Touitou, Isabelle; Mrenda, Bakridine M’Madi; Benedito, Cécilia; Pardoux, Etienne; Gagez, Anne-Laure; Pers, Yves-Marie; Jorgensen, Christian; Mahjoub, Touhami; Apparailly, Florence

doi:https://doi.org/10.1155/2017/4915950

Journal of Immunology Research

On this page

Abstract Introduction Methods Results Discussion Conclusion Abbreviations Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Immune System and Chronic Diseases

View this Special Issue

Research Article | Open Access

Volume 2017 | Article ID 4915950 | https://doi.org/10.1155/2017/4915950

TMEM187-IRAK1 Polymorphisms Associated with Rheumatoid Arthritis Susceptibility in Tunisian and French Female Populations: Influence of Geographic Origin

Olfa Khalifa,^1,2Nathalie Balandraud,^3,4Nathalie Lambert,³Isabelle Auger,³Jean Roudier,^3,4Audrey Sénéchal,⁵David Geneviève,^1,2,6Christophe Picard,⁷Gérard Lefranc,^2,8Isabelle Touitou,^1,2,9Bakridine M’Madi Mrenda,¹⁰and Cécilia Benedito¹⁰ et al.

Academic Editor: Andréia M. Cardoso

Received16 Sept 2016

Revised21 Dec 2016

Accepted28 Dec 2016

Published08 Feb 2017

Abstract

Polymorphisms have been identified in the Xq28 locus as risk loci for rheumatoid arthritis (RA). Here, we investigated the association between three polymorphisms in the Xq28 region containing TMEM187 and IRAK1 (rs13397, rs1059703, and rs1059702) in two unstudied populations: Tunisian and French. The rs13397 G and rs1059703 T major alleles were significantly increased in RA patients () compared with age-matched controls () in both Tunisian and French women. These results were confirmed by a meta-analysis replication study including two independent Greek and Korean cohorts. The rs1059702 C major allele was significantly associated with RA, only with French women. In the French population, the GTC haplotype displayed a protective effect against RA, while the ATC, GCC, and GTT haplotypes conferred significant risk for RA. No association for these haplotypes was found in the Tunisian population. Our results replicated for the first time the association of the three Xq28 polymorphisms with RA risk in Tunisian and French populations and suggested that RA susceptibility is associated with TMEM187-IRAK1 polymorphisms in women. Our data further support the involvement of X chromosome in RA susceptibility and evidence ethnicities differences that might be explained by differences in the frequencies of SE HLA-DRB1 alleles between both populations.

1. Introduction

Rheumatoid arthritis (RA) is a systemic autoimmune disease. It is characterized by chronic destructive inflammation in synovial joints. The prevalence of RA is about 1% in the adult European population and is three times more common in women than in men [1–4]. The aetiology of this complex disease is still poorly understood but is considered as a result of interaction between susceptibility genes and environmental factors [5]. The genetic contribution to RA has been estimated to be about 50–60% [6], with the HLA (Human leukocyte antigen) classes II molecules remaining the most powerful known genetic factor [7]. There is extensive evidence for the association between certain HLA-DRB1 alleles that contain a conserved sequence of five amino acids (Q/RK/RRAA) in the third hypervariable region of the DRβ1 chain, the so-called shared epitope (SE), and RA susceptibility and severity [8–10].

Genome-wide association studies (GWAS) and high-density array Immunochip studies for single nucleotide polymorphisms (SNPs) genotyping of populations have identified over 101 RA risk loci involved among individuals of European and Asian ancestry [11–15]. Among these, IRAK1 (interleukin 1 receptor associated kinase) was the first X chromosome locus reported as associated with RA susceptibility and is thus of importance given the female predominance of the disease. IRAK1 is a serine-threonine protein kinase and an essential component of the toll/interleukin 1 receptor (TIR) signaling pathway involved in the pathogen-mediated inflammation [16]. Interestingly, the IRAK1 gene is located on the Xq28 region that harbours several SNPs that have also been associated with susceptibility to autoimmune diseases. A recent case-control study investigated numerous SNPs located on the Xq28 locus and identified rs1059703 and rs1059702, encoding for pSer532Leu and pPhe196Ser, as two IRAK1 SNPs most significantly associated with RA susceptibility in Korean families [17]. These authors also showed that the major haplotype (rs1059702 T and rs1059703 C) was associated with increased IRAK1 activity. Of note, upstream IRAK1 and within the Xq28 risk locus, Eyre et al. reported the SNP rs13397 associated with RA risk among individuals of northern European ancestry [13]. This polymorphism is located within the TMEM187 gene, which encodes a transmembrane protein of unknown function.

The relationship between these three polymorphisms on the Xq28 region and RA risk in different population remains unaddressed. Interestingly, a clear separation of different ethnic and regional populations was evidenced when considering northern and southern Europe groups separately in genetic studies, suggesting the existence of European population genetic substructures [18]. In addition, there are also clear differences between northern and southern Europe for the main genetic risk factor for RA, HLA-DRB04:04, and DRB04:01 genotypes being the most frequent in northern populations, while HLA-DRB10:01 and DRB04:05 alleles are more frequent in southern Europe, as well as in Maghreb [19, 20]. Therefore, the present study aimed at investigating the association between critical polymorphisms in Xq28, from rs13397 (A/G) through rs1059703 (C/T) and rs1059702 (C/T) in the TMEM187 and IRAK1 locus, respectively, with risk for RA in two not previously studied populations displaying different genetic background: Tunisian and French. We also examined genetic differences for the HLA-DRB1 alleles in both RA populations. A meta-analysis including two other independent cohorts was then performed to test the overall effect of these TMEM187-IRAK1 polymorphisms on RA.

2. Patients and Methods

2.1. Study Subjects

For RA patients and healthy controls, informed consents were provided in accordance with national procedures. French studies were approved by local human ethical committees: sample collection and analysis (DC-2008-327) were approved by the “Cellule Bioéthique, Direction Générale pour la Recherche et l’Innovation, Ministère de l’Enseignement Supérieur et de la Recherche” (Ministry Bioethics Unit) and by Comité de Protection des Personnes Sud Méditerrannée IV (ID RCB 2008-A01087-48). Blood samples were collected from a total of 879 unrelated female participants: 408 patients with RA and 471 controls from Tunisia and France. Detailed size of Tunisian and French cohort is shown in Table 1. RA patients fulfilled the 2010 ACR/EULAR classification criteria for RA [21]. The clinical examination consisted of questionnaire and a physical examination: age, geographical origin, family history, body mass index, smoking habits, the presence of hypertension, diabetes, levels of physical activity, and alcohol consumption. The 471 controls matched with age were recruited from Voluntary Bone Marrow Donor (VBMD). Less than 1% of healthy VBMD women come from Maghreb. All patients were positive for anti-citrullinated peptide antibodies (ACPA) as determined by commercially available ELISA kits (Orgentec®, Hamburg, Germany) and Euroimmun® for Tunisian cases and CCP2 enzyme-linked immuno sorbent assay (ELISA) (Immunoscan RA, Euro-Diagnostica, Arnhem, Netherlands) for French cases.

2.2. Genotyping

Selection for Xq28 polymorphisms from previously reported SNPs was conducted using the following criteria: (a) the minor allele frequency (MAF) was >0.05 in the Caucasian population, according to the international HapMap project databank (http://www.hapmap.org) and to the NCBI SNP database (https://www.ncbi.nlm.nih.gov), (b) it was a tagging SNP, and (c) it has been reported to be associated with RA or other autoimmune diseases in different populations or other types of studies as GWAS. Genomic DNA was extracted from total blood leukocytes using standard methods including proteinase K digestion, followed by phenol–chloroform extraction and ethanol precipitation. Finally, all subjects were genotyped for three SNPs (rs1059702 and rs1059703 and rs13397) using direct PCR sequencing with the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) and an Applied Biosystems (ABI) 3130xL genetic analyzer (Applied Bio Systems, Foster City, CA, USA). PCR reactions were performed in a Bio-Rad thermal cycler using Taq Polymerase (Prim 5 by Fisher Scientific) with initial denaturing conditions at 96°C for 5 min, followed by 30 cycles of 96°C for 30 s, 64°C for 30 s, abd 72°C for 45 s and a final extension of 72°C for 10 min.

2.3. HLA-DRB1 Genotyping

Genomic DNA from whole blood was extracted as described in the subsection “Genotyping.” HLA-DRB1 typing was carried out according to the manufacturer’s specification for LAB type SSO HD (One Lambda Inc., USA) and the retrieved output was analyzed by HLA Fusion v 1.2.1 software (One lambda Inc., USA) for allele information. The following alleles or group of alleles were genotyped: HLA-DRB01, HLA-DRB01:03, HLA-DRB03, HLA-DRB04:01, HLA-DRB04:02, HLA-DRB04:03, HLA-04:04, HLA-DRB04:05, HLA-DRB04:06, HLA-DRB04:07, HLA-DRB04:08, HLA-DRB04:11, HLA-DRB07, HLA-DRB08, HLA-DRB09, HLA-DRB1001, HLA-DRB11, HLA-DRB12, HLA-DRB13, HLA-DRB14, HLA-DRB15, and HLA-DRB16 (Table 2).

2.4. Meta-Analysis Study

A literature search using PubMed, Science Direct, and Web of Science databases was conducted to identify publications that examine the association of TMEM187 rs13397, IRAK1 rs1059703, IRAK1 rs1059702, and IRAK1 rs3027898 polymorphisms with RA. Our research strategy was based on a combination of the following keywords: variation, variant, polymorphism, SNP, Rheumatoid arthritis, RA, TMEM187, IRAK1, rs13397, rs1059703, rs1059702, and rs3027898. We also reviewed the references mentioned in the identified articles for any additional relevant studies. Association studies included in our meta-analyses had to meet the following criteria: (1) evaluation of association between TMEM187 rs13397, IRAK1 rs1059703, IRAK1 rs1059702, or IRAK1 rs3027898 polymorphisms and the susceptibility to RA, (2) use of case-control design, and (3) inclusion of available genotype frequencies or sufficient information for calculation. The following information was extracted from the selected studies: first author, years of publication, characteristics of cases and controls (mean age, distribution of gender, and ethnicity), number of cases and controls, and genotype/allele frequency information. The strength of the association between the three polymorphisms and RA was measured by pooled OR with 95% CI. The test was used to determine the significance () of the pooled OR. The heterogeneity between studies was checked by test and statistics. We used a fixed effect model when value of the test was less than 0.10 or was equal to or less than 50%; otherwise, random effect model was used. Finally this meta-analysis was performed using Review Manager 5.3 software (Cochrane Collaboration, London, UK).

2.5. Statistical Methods

In this study, the Bayesian statistical method was used to evaluate the confidence intervals (called in the Bayesian framework credibility intervals) and to compute the values associated with our tests, as in [20]. Indeed, the Bayesian statistical methods are more appropriate in case of small samples and the “approximate value” of the variance of the logarithm of the Odds Ratio. In the Bayesian framework, the unknown parameter is considered as the realization of a random variable, which forces us to choose an a priori probability distribution for that random variable. We chose an “a priori” distribution for a pair of frequencies and decided that the two random variables are independent, each with the probability density equal to 1 on the interval (0.1). It is fair to say that our prior is uninformative: it is the uniform distribution on the interval (0.1). The “a posteriori” distribution of each parameter is easily shown from the Bayes formula to be a Beta distribution. More precisely, if while observing realizations of 1 s and 0 s, we get 1 s and 0 s, then the a posteriori distribution of the proportion is the Beta (, ) probability distribution. Since the two samples come from completely unrelated sources, it is fair to assume that the two proportions are conditionally independent, given the data. Consequently a posteriori distribution of the pair of parameters (corresponding to the proportions of the given allele among the patients and among the control individuals) is the product of two Beta distributions whose density is explicitly known. Our method involves simulating one million realizations of that pair of random variables, from which we deduce one million realizations of the Odds Ratio. The value of the test is then obtained by looking at which proportion of those values is smaller (resp., larger) than 1; this gives us . Finally we want to discuss the issue of multiple testing. We have used our Bayesian test method 6 times (Table 5) where the probability of an error among six successive tests is not controlled by the probability of an error in each single test. The values were corrected for multiple comparisons using the Bonferroni correction in Table 2. PLINK 1.7 version software was used for analysis of association between haplotypes and disease, and Haploview 4.2 software was used for LD analysis. All haplotypes with frequencies less than 5% were ignored in analysis. Distributions of data were tested with the Shapiro-Wilk test. Chi² or Fisher test was used for categorical data; Yates correction was performed. For numerical data, comparisons of medians were performed using Mann–Whitney or Student’s -tests. For correlation study, statistical analyses were performed at the conventional two-tailed α level of 0.05 using R software version 3.0.2.10. Odds ratios (OR) with 95% confidence interval (CI) were also calculated. A value less than 0.05 was considered statistically significant.

3. Results

3.1. Demographics and Phenotypes of the Study Subjects

Characteristics of the 408 ACPA (anti-citrullinated protein antibody) positive RA and 471 controls from Tunisia and South Eastern France (Marseille and Montpellier) are shown in Table 1. All samples were females. The average age was 52 years for the Tunisian RA and 63 years for French RA women. Age was not significantly different between patients and control females. There were more presently smokers in controls than in patients group.

3.2. Frequency of HLA-DRB1 Alleles and SE

Allelic frequencies of the HLA-DRB1 alleles in Tunisian and French RA patients are presented in Table 2. Due to the size of cohorts, we applied Bayesian statistical method to define susceptibility or protective genotypes and calculate accurate confidence intervals for the associated Odds Ratios. HLA-DRB01, HLA-DRB03, HLA-DRB04:01, HLA-DRB04:04, HLA-DRB04:05, HLA-DRB04:06, HLA-DRB10, HLA-DRB13, HLA-DRB15, and HLA-DRB16 showed statistical differences between the French and the Tunisian RA women. No statistical differences were found for HLA-DRB01:03, HLA-DRB04:02, HLA-DRB04:03, HLA-DRB04:07, HLA-DRB04:08, HLA-DRB04:11, HLA-DRB07, HLA-DRB08, HLA-DRB09, HLA-DRB11, HLA-DRB12, and HLA-DRB14 between both populations. The RA associated HLA-DRB1 alleles that were significantly more present in the Tunisian samples compared to those in the French samples were, by order of Odds Ratios, HLA-DRB04:05 (OR = 2.15), HLA-DRB13 (OR = 2.39), HLA-DRB10 (OR = 3.65), HLA-DRB03 (OR = 4.17), and HLA-DRB04:06 (OR = +∞). Similarly, RA associated HLA-DRB1 alleles that were significantly more observed in the French than in the Tunisian samples were, by order of Odds Ratios, HLA-DRB16 (OR = 0.09), HLA-DRB04:01 (OR = 0.26), HLA-DRB01 (OR = 0.31), HLA-DRB04:04 (OR = 0.32), and HLA-DRB15 (OR = 0.45). Three HLA-DRB1 alleles (DRB04:05, DRB12:XX, and DRB15:XX) did not show significant differences between Tunisian and French RA women after the Bonferroni correction.

Two (HLA-DRB04:05 and HLA-DRB10) out of seven shared SE alleles were overrepresented in the Tunisian RA population, and three SE alleles (HLA-DRB01, HLA-DRB04:01, and HLA-DRB04:04) were overrepresented in the French population. SE single dose showed similar frequency in both Tunisian and French RA women (50% and 57.80%, resp.), while French RA population displayed higher SE double dose than Tunisian RA population (20.34% and 9.09%, resp.) (Table 3).

3.3. Association between Three TMEM187-IRAK1 Polymorphisms and RA Susceptibility in Women Samples from Tunisian and French Populations

Genotype distributions of all three investigated SNPs were in Hardy-Weinberg equilibrium in both patient and control groups ( < 3.84). Among the genotype frequencies of the three SNPs analyzed using the dominant modelling, the rs13397 and rs1059703 demonstrated a significant association with RA susceptibility in both populations. The rs13397 G allele is present in 75% of RA Tunisian patients and in 69% RA French (Table 4). The frequency of the rs13397 G major allele was significantly lower in RA patients than in controls in both Tunisian (OR 0.48 [95% CI 0.12–0.73], ) and French (OR = 0.45 [95% CI 0.21–0.59], ) women populations. The rs1059703 T major allele was significantly decreased in patients compared with controls in both Tunisian (OR = 0.46 [95% CI 0.13–0.71], ) and French (OR = 0.36 [95% CI 0.15–0.47], ) women populations. The SNP rs1059702 C showed a significant association with RA in French patients as compared with control subjects (OR 1.97 [95% CI 1.22–11.52], ), but none in the Tunisian population. Of note, no correlation was found between polymorphisms and the SE (data not shown).

3.4. Linkage Disequilibrium of TMEM187-IRAK1 Locus with RA in the Samples of the Tunisian and French Populations

To determine the extent of linkage disequilibrium (LD) among the three polymorphisms enclosed in the TMEM187-IRAK1 locus in both Tunisian and French populations, standardized LD coefficient was calculated for all pairs of the three polymorphisms. Our study showed that the IRAK1 SNPs rs1059703 and rs1059702 did not belong to the same linkage disequilibrium locus, both in Tunisian and French populations (Figure 1). We also showed that rs13397 and rs1059703 polymorphisms were in moderate linkage disequilibrium ( = 0.5) in the Tunisian population (Figure 1(a)), but not in the French population (Figure 1(b), = 0.11).

(a) Tunisian population

(b) French population

3.5. Haplotype Analysis

To examine the combined effect of three variants within the TMEM187-IRAK1 locus, the haplotype frequencies defined by the three SNPs rs13397, rs1059703, and rs1059702 were analyzed. Among the potential haplotypes, four common haplotypes were identified in our populations with frequencies greater than 5% in both cases and controls (Table 6).

Results of the analysis suggested that the GTC haplotype, which carries the rs13397 G variant, displayed a protective effect against RA in our sample of the French population (OR = 0.23 (95% CI) [0.16–0.33], ) while the ATC, GCC, and GTT haplotypes conferred significant risk for RA in the same population. French women with the ATC, GCC, or GTT haplotypes had 3.97-, 2.39-, and 2.05-fold higher risk of having RA compared to controls. In contrast, association between all 4 detected haplotypes and disease was not significant in the sample of the Tunisian women. Overall, these data further support the involvement of the Xq28 locus in RA susceptibility in the French population and underscore ethnical differences.

3.6. Meta-Analysis Study

A meta-analysis was performed and included datasets from 2076 RA cases and 2112 healthy controls. To carry out the meta-analysis, we combined our findings with three independent studies that fit the inclusion criteria and included a total of 1668 RA patients and 1641 ethnically matched healthy controls from Greece and Korea [17, 22, 23]. The samples size of these studies varied from 119 to 1318 cases and from 131 to 1016 controls. The characteristics of all studies were summarized in Table 7. Meta-analysis results showed a significant association between TMEM187 rs13397 (OR = 1.56, 95% CI = 1.25–1.94, and value < 0.0001) and IRAK1 rs1059702 (OR = 1.41, 95% CI = 1.11–1.80, and value = 0.006) polymorphisms and susceptibility to RA (Figures 2(a) and 2(c)). The IRAK1 rs1059703 polymorphism showed a significant protective effect in RA (OR = 0.68, 95% CI = 0.51–0.89, and value = 0.006) (Figure 2(b)). The IRAK1 rs3027898 polymorphism that was available in the three independent studies from Korea, China, and Greece showed no significant association with RA risk (OR = 0.96, 95% CI = 0.67–1.39, and value = 0.83).

(a)

(b)

(c)

4. Discussion

In the present study, we aimed at investigating the association between three critical polymorphisms in the TMEM187-IRAK1 locus and risk for ACPA positive RA female patients in two not previously studied populations: Tunisian and French. We found that the rs13397 A and rs1059703 C risk alleles were significantly associated with RA susceptibility for both populations. The rs1059702 T allele was only found associated with risk for RA in the French sample. Four haplotypes were detected in both populations. In the French population, the major GTC haplotype displayed a protective effect against RA, whereas the GCC, GTT, and ATC haplotypes conferred significant increased risk for RA. No association for these haplotypes was however found in the Tunisian population. Finally, a strong genetic association of RA susceptibility exists for the rs13397/rs1059703 block only in women Tunisian subjects.

Among the many loci that have been identified as associated with RA risk using either candidate gene association studies or GWAS, IRAK1 is the only gene located on the X chromosome. The IRAK1 gene encodes the interleukin 1 receptor-associated kinase 1, a serine/threonine kinase that mediates TLRs and IL-1 receptor signaling. Upon receptor stimulation, IRAK1 is recruited, phosphorylated, and responsible for TLR- and IL1-induced activation of the NF-κB signaling cascade, thus playing a critical role in the early inflammatory response [24]. The IRAK1 gene is located within the Xq28 locus that contains several polymorphisms associated with increased susceptibility to autoimmunity, including systemic lupus erythematous [25], systemic sclerosis [26, 27], and RA [13, 17, 22, 23]. Altogether, these studies demonstrated that the Xq28 locus is a shared locus for several complex diseases, thereby highlighting its role in inflammatory disease susceptibility. Recently, two studies devoted attention to the TMEM187-IRAK1 locus and identified two SNPs in a Korean population (rs1059702 and rs1059703) [17] and one SNP (rs13397) in a large cohort of northern European ancestry [13], respectively, associated with increased disease susceptibility. In an attempt to investigate new, not previously studied populations and to replicate the above studies, we examined the distribution of these three SNPs and found statistical differences between RA patients and controls for both Tunisian and French populations. Although risk alleles for both rs1059702 and rs1059703 were different, they were also associated with risk for RA in both Tunisian and French populations. Indeed, the rs1059702 T and rs1059703 C polymorphisms were in strong association with RA in Korean women while we found the rs1059702 C and rs1059703 T alleles in our cohorts. In agreement with the study of Eyre et al. [13], we found the minor rs13397 A allele to be the risk allele in both the French and Tunisian cohorts. Earlier work performed on a small Greek cohort found no association between the rs1059703 and RA [22]. These differences might have several explanations, including the sample size of cohorts used, the Bayesian statistical method being better suited than the SPSS statistical package to decide about susceptibility in small samples, the classification criteria used to select RA patients (ACR 1987 versus 2010 ACR/EULAR that discriminate or not between ACPA positive RA patients), and different ethnic and regional populations.

Concerning the haplotype and LD studies, we are the first to study the combination of LD and haplotypes of these three SNPs with RA. We found that, with the exception of rs13397 and rs1059703 polymorphisms, which were in moderate linkage disequilibrium in the Tunisian population ( = 0.5), other polymorphisms were not in linkage disequilibrium in both population. This finding seems to be in contradiction with the Korean study that previously reported a linkage disequilibrium block on Xq28 enclosing both rs1059702 and rs1059703 [17]. Discrepancies between Korean, Tunisian, and French populations might be due to differences in genetic background. Interestingly, a comparative study reported that none of the susceptibility loci identified by GWAS in Caucasian patients with RA contributed significantly to disease in Koreans [28].

We found differences between French and Tunisian populations for association of SNP haplotypes in the TMEM187-IRAK1 locus with RA susceptibility. The CTG haplotype displayed a protective effect against RA, while the ATC, GCC, and GTT haplotypes conferred significant susceptibility to RA, in French women. All four haplotypes showed no association in the Tunisian women sample. In both Tunisian and French cohorts, we found the GTC common haplotype in 57 and 68% of the control populations, respectively. This GTC haplotype carries the major alleles for the three SNPs rs13397 G, rs1059703 T, and rs1059702 C.

Although alleles and haplotypes identified among Tunisian and French subjects display differences with previous studies performed with Korean and northern European cohorts, our present work confirms significant association of the three SNPs with RA susceptibility. The minor discrepancy found for the rs1059702 in the Tunisian RA women population might be due to the size of the cohort. Indeed, pooled results for the three different ethnical groups including Tunisian confirm a significant association between rs1059702 and risk for RA.

Our meta-analysis including a larger case-control data set of two independent studies from Korea and Greece in addition to our data sets confirmed that TMEM187 rs13397, IRAK1 rs1059702, and IRAK1 rs1059703 are associated with RA. Unfortunately, data sets on the TMEM187 rs13397 from populations of European ancestry [13] were not available for inclusion in our meta-analysis. The IRAK1 rs3027898 that was not investigated in our study but available in the three independent studies from Korea, China, and Greece showed no significant association in the meta-analysis (data not shown). To the best of our knowledge, this is the first meta-analysis study on the association between RA and these four polymorphisms in TMEM187-IRAK1.

Many studies have documented that one of the key genetic elements closely linked to RA susceptibility is the HLA complex locus located in the chromosome 6 [15]. Our present study evaluated the frequency of HLA-DRB1 alleles in Tunisian and French RA women, all ACPA positive, and confirmed that both populations belonging to different ethnical groups do not share the same HLA susceptibility alleles. As expected for the French population, we find the usual association between RA and HLA-DRB01, HLA-DRB04:01, and HLA-DRB04:04 [20]. We show that two SE alleles, HLA-DRB04:05 and HLA-DRB10, are overrepresented in the RA Tunisian population. This is in agreement with a previous report showing their association with RA susceptibility in Tunisia [19]. The HLA-DRB04 allele has also been reported in RA populations from Morocco [29], Saudi Arabia [30], and Turkey [31]. The HLA-DRB10 allele has been reported to be associated with RA in Saudi Arabia [30], Italy [32], Spain [33], Turkey [34], and Asian Indian [35] populations.

A meta-analysis performed in Latin America RA population (Colombia, Peru, Brazil, Argentina, Mexico, and Chili) showed that HLA-DRB04:04 allele is associated with disease risk [36]. Moreover, this whole genome scan study also shows that the second most consistently linked region to RA is the Xq28 locus. The absence of the HLA-DRB01:03, HLA-DRB04:07, and HLA-DRB12 alleles in our RA Tunisian population sample is in accordance with the protective effect suggested by previous study [19]. Considering the observed differences between ethnical groups for HLA-DRB1, our findings suggest that the association of the TMEM187-IRAK1 locus with RA risk may play a role independently of the major genetic risk locus for RA.

5. Conclusion

In conclusion, the present study illustrates the involvement of the Xq28 locus in RA susceptibility in both Tunisian and French populations with ethnical differences. Finally, the meta-analysis combining our findings with two independent studies from Greece and Korea further supports the fact that the three polymorphisms are important in the susceptibility for RA. Although a growing number of studies suggest the contribution of the X chromosome to RA, as to many other autoimmune diseases, future works shall be performed to compare female and male RA populations to identify the exact contribution of this locus, with a focus on TMEM187-IRAK1 polymorphisms. Furthermore, functional studies on the impact of these polymorphisms on both gene expressions in RA populations might help defining new therapeutic perspectives in RA.

Abbreviations

ACPA:	Anti-citrullinated protein antibody
ACR:	American College of Rheumatology
AR:	Androgen receptor
CCP2:	Cyclic-citrullinated peptide
CD40L:	Cluster of differentiation 40 ligand
CI:	Confidence interval
ELISA:	Enzyme-linked Immunosorbent assay
EULAR:	European League against Rheumatism
FOXP3:	Forkhead box P3
GWAS:	Genome-wide association studies
HLA-DRB1:	Human leukocyte antigen-DRB1
IL2RG:	Interleukin-2 receptor gamma chain
IRAK1:	Interleukin 1 receptor associated kinase 1
LD:	Linkage disequilibrium
MAF:	Minor allele frequency
NF-kB:	Nuclear factor-kappa B
OR:	Odds Ratio
PCR:	Polymerase chain reaction
RA:	Rheumatoid arthritis
SNP:	Single nucleotide polymorphism
SE:	Shared epitope
TLR:	Toll-like receptors
TMEM187:	Transmembrane protein 187
VBMD:	Voluntary bone marrow donor.

Competing Interests

The authors declaring that they have no competing interests are Olfa Khalifa, Nathalie Balandraud, Nathalie Lambert, Isabelle Auger, Jean Roudier, Audrey Sénéchal, David Geneviève, Christophe Picard, Gérard Lefranc, Hicham Charoute, Isabelle Touitou, Bakridine M’Madi Mrenda, Cécilia Benedito, Etienne Pardoux, Anne-Laure Gagez, Yves-Marie Pers, Christian Jorgensen, Touhami Mahjoub, and Florence Apparailly.

Authors’ Contributions

Olfa Khalifa carried out the molecular genetic studies, recruited the Tunisian RA patients and healthy donors, performed the sequence alignment, analyzed data, wrote the manuscript and performed the meta-analysis study. Nathalie Balandraud, Nathalie Lambert, Jean Roudier, Isabelle Auger recruited the French RA patients and healthy donors, collected materials, analyzed data, and participate in writing the manuscript. Audrey Sénéchal performed the sequence alignment and participated in writing the manuscript. David Geneviève, Gérard Lefranc, Isabelle Touitou, Yves-Marie Pers, Christian Jorgensen overviewed the project and participated in writing the manuscript. Christophe Picard provided materials of healthy donors, participated in writing the manuscript and performed the HLA-DRB1 typing. Bakridine M’Madi Mrenda, Etienne Pardoux, Cécilia Benedito, Anne-Laure Gagez developed mathematic tools, performed statistical analyses, and participated in writing the manuscript. Touhami Mahjoub recruited the Tunisian healthy donors, overviewed the project and participated in writing the manuscript. Florence Apparailly conceived of the study, overviewed the design and coordination of the project, and participated in writing the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This work is supported by the Inserm, University of Montpellier, the EU funding IMI- BTCure (no. 115142) and ERASMUS-Mundus (no. 2888). The authors thank Elodie Sanchez and Rosana Ferreira-Lopez for technical and clinical support. They are grateful to Nikoleta Vrcelj for English editing.

References

A. Linos, J. W. Worthington, M. O'fallon, and L. T. Kurland, “The epidemiology of rheumatoid arthritis in Rochester, Minnesota: a study of incidence, prevalence, and mortality,” American Journal of Epidemiology, vol. 111, no. 1, pp. 87–98, 1980.
View at: Google Scholar
P. K. Gregersen, C. I. Amos, A. T. Lee et al., “REL, encoding a member of the NF-κB family of transcription factors, is a newly defined risk locus for rheumatoid arthritis,” Nature Genetics, vol. 41, no. 7, pp. 820–823, 2009.
View at: Publisher Site | Google Scholar
R. F. van Vollenhoven, “Sex differences in rheumatoid arthritis: more than meets the eye...,” BMC Medicine, vol. 7, article 12, 2009.
View at: Publisher Site | Google Scholar
L. Carmona, M. Cross, B. Williams, M. Lassere, and L. March, “Rheumatoid arthritis,” Best Practice and Research: Clinical Rheumatology, vol. 24, no. 6, pp. 733–745, 2010.
View at: Publisher Site | Google Scholar
P. Dieudé, “Rheumatic diseases: environment and genetics,” Joint Bone Spine, vol. 76, no. 6, pp. 602–607, 2009.
View at: Publisher Site | Google Scholar
S. Viatte, D. Plant, and S. Raychaudhuri, “Genetics and epigenetics of rheumatoid arthritis,” Nature Reviews Rheumatology, vol. 9, no. 3, pp. 141–153, 2013.
View at: Publisher Site | Google Scholar
C. Perricone, F. Ceccarelli, and G. Valesini, “An overview on the genetic of rheumatoid arthritis: a never-ending story,” Autoimmunity Reviews, vol. 10, no. 10, pp. 599–608, 2011.
View at: Publisher Site | Google Scholar
S. W. Scally, J. Petersen, S. C. Law et al., “A molecular basis for the association of the HLA-DRB1 locus, citrullination, and rheumatoid arthritis,” Journal of Experimental Medicine, vol. 210, no. 12, pp. 2569–2582, 2013.
View at: Publisher Site | Google Scholar
P. K. Gregersen, J. Silver, and R. J. Winchester, “The shared epitope hypothesis. An approach to understanding the molecular genetics of susceptibility to rheumatoid arthritis,” Arthritis and Rheumatism, vol. 30, no. 11, pp. 1205–1213, 1987.
View at: Publisher Site | Google Scholar
E. Lundström, H. Källberg, L. Alfredsson, L. Klareskog, and L. Padyukov, “Gene-environment interaction between the DRB1 shared epitope and smoking in the risk of anti-citrullinated protein antibody-positive rheumatoid arthritis: all alleles are important,” Arthritis and Rheumatism, vol. 60, no. 6, pp. 1597–1603, 2009.
View at: Publisher Site | Google Scholar
E. A. Stahl, S. Raychaudhuri, E. F. Remmers et al., “Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci,” Nature Genetics, vol. 42, no. 6, pp. 508–514, 2010.
View at: Publisher Site | Google Scholar
L. Padyukov, M. Seielstad, R. T. H. Ong et al., “A genome-wide association study suggests contrasting associations in ACPA-positive versus ACPA-negative rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 70, no. 2, pp. 259–265, 2011.
View at: Publisher Site | Google Scholar
S. Eyre, J. Bowes, D. Diogo et al., “High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis,” Nature Genetics, vol. 44, no. 12, pp. 1336–1340, 2012.
View at: Publisher Site | Google Scholar
K. Kim, S.-Y. Bang, H.-S. Lee et al., “High-density genotyping of immune loci in Koreans and Europeans identifies eight new rheumatoid arthritis risk loci,” Annals of the Rheumatic Diseases, vol. 74, no. 3, article no. e13, 2015.
View at: Publisher Site | Google Scholar
Y. Okada, C. Terao, K. Ikari et al., “Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population,” Nature Genetics, vol. 44, no. 5, pp. 511–516, 2012.
View at: Publisher Site | Google Scholar
M. U. Martin and H. Wesche, “Summary and comparison of the signaling mechanisms of the Toll/interleukin-1 receptor family,” Biochimica et Biophysica Acta—Molecular Cell Research, vol. 1592, no. 3, pp. 265–280, 2002.
View at: Publisher Site | Google Scholar
T.-U. Han, S.-K. Cho, T. Kim, Y. B. Joo, S.-C. Bae, and C. Kang, “Association of an activity-enhancing variant of IRAK1 and an MECP2-IRAK1 haplotype with increased susceptibility to rheumatoid arthritis,” Arthritis and Rheumatism, vol. 65, no. 3, pp. 590–598, 2013.
View at: Publisher Site | Google Scholar
C. Tian, R. Kosoy, R. Nassir et al., “European population genetic substructure: further definition of ancestry informative markers for distinguishing among diverse European ethnic groups,” Molecular Medicine, vol. 15, no. 11-12, pp. 371–383, 2009.
View at: Publisher Site | Google Scholar
M. B. Hamad, N. Mahfoudh, S. Marzouk et al., “Association study of human leukocyte antigen-DRB1 alleles with rheumatoid arthritis in south Tunisian patients,” Clinical Rheumatology, vol. 31, no. 6, pp. 937–942, 2012.
View at: Publisher Site | Google Scholar
N. Balandraud, C. Picard, D. Reviron et al., “HLA-DRB1 genotypes and the risk of developing anti citrullinated protein antibody (ACPA) positive rheumatoid arthritis,” PLoS ONE, vol. 8, no. 5, Article ID e64108, 2013.
View at: Publisher Site | Google Scholar
D. Aletaha, T. Neogi, A. J. Silman et al., “2010 Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative,” Arthritis and Rheumatism, vol. 62, no. 9, pp. 2569–2581, 2010.
View at: Publisher Site | Google Scholar
A. Chatzikyriakidou, P. V. Voulgari, I. Georgiou, and A. A. Drosos, “A polymorphism in the 3'-UTR of interleukin-1 receptor-associated kinase (IRAK1), a target gene of miR-146a, is associated with rheumatoid arthritis susceptibility,” Joint Bone Spine, vol. 77, no. 5, pp. 411–413, 2010.
View at: Publisher Site | Google Scholar
H. Zhang, J. Pu, X. Wang et al., “IRAK1 rs3027898 C/A polymorphism is associated with risk of rheumatoid arthritis,” Rheumatology International, vol. 33, no. 2, pp. 369–375, 2013.
View at: Publisher Site | Google Scholar
X.-F. Sun and H. Zhang, “NFKB and NFKBI polymorphisms in relation to susceptibility of tumour and other diseases,” Histology and Histopathology, vol. 22, no. 12, pp. 1387–1398, 2007.
View at: Google Scholar
K. M. Kaufman, J. Zhao, J. A. Kelly et al., “Fine mapping of Xq28: both MECP2 and IRAK1 contribute to risk for systemic lupus erythematosus in multiple ancestral groups,” Annals of the Rheumatic Diseases, vol. 72, no. 3, pp. 437–444, 2013.
View at: Publisher Site | Google Scholar
F. David Carmona, M. Carmen Cénit, L.-M. Diaz-Gallo et al., “New insight on the Xq28 association with systemic sclerosis,” Annals of the Rheumatic Diseases, vol. 72, no. 12, pp. 2032–2038, 2013.
View at: Publisher Site | Google Scholar
P. Dieudé, M. Bouaziz, M. Guedj et al., “Evidence of the contribution of the X chromosome to systemic sclerosis susceptibility: association with the functional IRAK1 196Phe/532Ser haplotype,” Arthritis and Rheumatism, vol. 63, no. 12, pp. 3979–3987, 2011.
View at: Publisher Site | Google Scholar
H.-S. Lee, B. D. Korman, J. M. Le et al., “Genetic risk factors for rheumatoid arthritis differ in caucasian and Korean populations,” Arthritis and Rheumatism, vol. 60, no. 2, pp. 364–371, 2009.
View at: Publisher Site | Google Scholar
O. Atouf, K. Benbouazza, C. Brick et al., “HLA polymorphism and early rheumatoid arthritis in the Moroccan population,” Joint Bone Spine, vol. 75, no. 5, pp. 554–558, 2008.
View at: Publisher Site | Google Scholar
R. Al-Swailem, H. Al-Rayes, S. Sobki, M. Arfin, and M. Tariq, “HLA-DRB1 association in Saudi rheumatoid arthritis patients,” Rheumatology International, vol. 26, no. 11, pp. 1019–1024, 2006.
View at: Publisher Site | Google Scholar
F. Uçar, M. Karkucak, E. Alemdaroğlu et al., “HLA-DRB1 allele distribution and its relation to rheumatoid arthritis in eastern Black Sea Turkish population,” Rheumatology International, vol. 32, no. 4, pp. 1003–1007, 2012.
View at: Publisher Site | Google Scholar
S. M. Bongi, B. Porfirio, G. Rombolà, A. Palasciano, E. Beneforti, and G. Bianucci, “Shared-epitope HLA-DRB1 alleles and sex ratio in Italian patients with rheumatoid arthritis,” Joint Bone Spine, vol. 71, no. 1, pp. 24–28, 2004.
View at: Publisher Site | Google Scholar
J. Yelamos, J. R. Garcia-Lozano, I. Moreno et al., “Association of HLA-DR4-Dw15 (DRB1*0405) and DR10 with rheumatoid arthritis in a Spanish population,” Arthritis and Rheumatism, vol. 36, no. 6, pp. 811–814, 1993.
View at: Publisher Site | Google Scholar
G. Kinikli, A. Ateş, M. Turgay, G. Akay, S. Kinikli, and G. Tokgöz, “HLA-DRB1 genes and disease severity in rheumatoid arthritis in Turkey,” Scandinavian Journal of Rheumatology, vol. 32, no. 5, pp. 277–280, 2003.
View at: Publisher Site | Google Scholar
V. Taneja, M. J. Giphart, W. Verduijn, A. Naipal, A. N. Malaviya, and N. K. Mehra, “Polymorphism of HLA-DRB, -DQA1, and -DQB1 in rheumatoid arthritis in Asian Indians: association with DRB1∗0405 and DRB1∗1001,” Human Immunology, vol. 46, no. 1, pp. 35–41, 1996.
View at: Publisher Site | Google Scholar
A. M. Delgado-Vega and J.-M. Anaya, “Meta-analysis of HLA-DRB1 polymorphism in Latin American patients with rheumatoid arthritis,” Autoimmunity Reviews, vol. 6, no. 6, pp. 402–408, 2007.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2017 Olfa Khalifa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1310

Downloads

844

Citations