Email updates

Keep up to date with the latest news and content from Arthritis Research & Therapy and BioMed Central.

Open Access Research article

Polymorphisms in peptidylarginine deiminase associate with rheumatoid arthritis in diverse Asian populations: evidence from MyEIRA study and meta-analysis

Chun Lai Too12*, Shahnaz Murad2, Jasbir Singh Dhaliwal2, Per Larsson3, Xia Jiang4, Bo Ding4, Lars Alfredsson4, Lars Klareskog1 and Leonid Padyukov1

Author Affiliations

1 Rheumatology Unit, Department of Medicine, Karolinska University Hospital, CMM L8:04, S-171 76 Stockholm, Sweden

2 Institute for Medical Research, Jalan Pahang, 50588 Kuala Lumpur, Malaysia

3 Department of Rheumatology, R 92, Karolinska University Hospital Huddinge, S-171 76 Stockholm, Sweden

4 Institute of Environmental Medicine, Karolinska Institutet, Nobel väg 13, Box 210, S-171 77 Stockholm, Sweden

For all author emails, please log on.

Arthritis Research & Therapy 2012, 14:R250  doi:10.1186/ar4093

The electronic version of this article is the complete one and can be found online at: http://arthritis-research.com/content/14/6/R250


Received:5 June 2012
Revisions received:9 October 2012
Accepted:14 November 2012
Published:19 November 2012

© 2012 Too et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Introduction

The majority of our knowledge regarding disease-related mechanisms of uncontrolled citrullination and anti-citrullinated protein antibody development in rheumatoid arthritis (RA) was investigated in Caucasian populations. However, peptidylarginine deiminase (PADI) type 4 gene polymorphisms are associated with RA in East Asian populations and weak or no association was found in Caucasian populations. This study explores the association between the PADI4 polymorphisms and RA risk in a multiethnic population residing in South East Asia with the goal of elucidating generalizability of association in non-Caucasian populations.

Methods

A total of 320 SNPs from the PADI locus (including PADI1, PADI2, PADI3, PADI4 and PADI6 genes) were genotyped in 1,238 RA cases and 1,571 control subjects from the Malaysian Epidemiological Investigation of Rheumatoid Arthritis (MyEIRA) case-control study. Additionally, we conducted meta-analysis of our data together with the previously published studies of RA from East Asian populations.

Results

The overall odds ratio (ORoverall) for the PADI4 (rs2240340) allelic model was 1.11 (95% confidence interval (CI) = 1.00 to 1.23, P = 0.04) and for the genotypic model was 1.20 (95% CI = 1.01 to 1.44, P = 0.04). Haplotype analysis for four selected PADI4 SNPs revealed a significant association of one with susceptibility (P = 0.001) and of another with a protective effect (P = 0.02). The RA susceptibility was further confirmed when combined meta-analysis was performed using these data together with data from five previously published studies from Asia comprising 5,192 RA cases and 4,317 control subjects (ORoverall = 1.23 (95% CI = 1.16 to 1.31, Pheterogeneity = 0.08) and 1.31 (95% CI = 1.20 to 1.44, Pheterogeneity = 0.32) in allele and genotype-based models, respectively). In addition, we also detected a novel association of PADI2 genetic variant rs1005753 with RA (ORoverall = 0.87 (95% CI = 0.77 to 0.99)).

Conclusion

Our study demonstrates an association between PADI4 and RA in the multiethnic population from South East Asia and suggests additional association with a PADI2 gene. The study thus provides further support for the notion that polymorphisms in genes for enzymes responsible for citrullination contribute to RA development in multiple populations of Asian descent.

Introduction

Most studies of the association between rheumatoid arthritis (RA) and genetic factors have focused on a group of human leukocyte antigens in the major histocompatibility complex [1] and a detailed account of the contributions from different major histocompatibility complex genes and their structural correlates was recently published for a number of Caucasian populations [2]. Additional contributions from more than 30 different non-HLA loci have been demonstrated, mainly in populations of Caucasian origin [3]. Important differences for RA susceptibility genes have, however, been described between Caucasian and non-Caucasian populations, as seen both from which HLA alleles are associated with disease [4-7] and from associations with non-HLA genes. A particular interesting difference between Caucasian and Asian populations has been demonstrated in genes from a peptidylarginine deiminase (PADI) locus, where a polymorphism was first demonstrated in a Japanese population [8] and later confirmed in additional Japanese and Korean populations [9-12]. These polymorphisms are of particular interest for the pathogenesis of RA since PADI4 and other PADI enzymes catalyze change from peptidylarginine to peptidylcitrulline, a target of anti-citrullinated protein antibody (ACPA), through a post-translational modification process referred to as citrullination [13,14].

The associations between PADI4 polymorphisms and RA have so far focused on the Japanese and Korean populations [9,11,15]. The effect of PADI4 polymorphisms on RA risk, however, remains unclear in the Han Chinese population [16,17]. An association between PADI4 and RA has also been observed in German and North American populations [18,19], while such an association has not been replicated in other Caucasian populations (for example, British, Spanish, Swedish and Hungarian) [19-23], despite a comparable allele frequency between these Asian and Caucasian populations. The largest study performed in the UK population with over 19,000 subjects found no evidence for association between the PADI4_94 SNP (rs2240340) and RA. In a meta-analysis on previously published European studies together with this UK study, the association between the PADI4_94 genotype and RA was weak and statistically not significant (odds ratio (OR) 1.06, 95% confidence interval (CI) = 0.99 to 1.13, P = 0.12) [23]. There is thus a need to further investigate the impact of genetic variations in PADI-associated genes in additional populations. Since ethnic differences are likely to go hand in hand with different environmental exposures, it will be helpful to study it in more divergent populations with some similarities in genetic background.

The PADI gene region is located at chromosome 1p36. This locus contains the cluster of all the PADI genes (PADI1 to PADI4 and PADI6). Of the five isotypes of PAD protein, PAD2 and PAD4 have been reported to be expressed as active enzymes in RA synovium, where citrullination of matrix proteins could potentially create antigenic peptides [24,25] and where increased citrullination has been demonstrated in inflamed synovial tissues [26,27].

In the present study, which was performed in the South East Asia region, the aim was to determine whether the association between the PADI4 polymorphisms and RA risk could be generalized to the Malaysian populations with Malay, Chinese and Indian ethnicity. We also investigated multiple SNPs from the locus (PADI1 to PADI4 and PADI6) in this multiethnic case-control study involving early RA.

Materials and methods

Study population

The source of data for our investigation was the multicenter Malaysian Epidemiological Investigation of Rheumatoid Arthritis (MyEIRA) case-control study, comprising 1,238 cases of RA and 1,571 control subjects. The demographic characteristics of RA cases and controls are shown in Table 1. Of the 1,238 cases, 516 (41.7%) were Malays, 255 (20.6%) were Chinese, 379 (30.6%) were Indians and 88 (7.1%) were of other or mixed ethnicities from South East Asia. The details of the MyEIRA study have been described elsewhere [6]. Briefly, all RA cases were diagnosed by rheumatologists according to the 1987 revised American College of Rheumatology criteria [28]. The disease duration for the RA cases was on average 1 year (interquartile range = 2 years). For each potential case, a control subject was randomly selected from the population, taking into consideration the subject's age, sex and residential area. Of the 1,571 controls, 986 (62.8%) were Malays, 206 (13.1%) were Chinese, 285 (18.1%) were Indians and 94 (6.0%) were of other sub-ethnicities. We performed case-control association analyses with regard to the influence from different genetic factors for each ethnic group separately as well as for the meta-analysis where all four ethnic groups were included. The study was approved by the Medical Research and Ethics Committee, Ministry of Health, Malaysia, and written informed consent was obtained from all participants.

Table 1. Demographic characteristics of rheumatoid arthritis patients and controls in the MyEIRA study population

DNA extraction, selection of markers and genotyping

The genomic DNA was extracted using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). All DNA samples were stored at -20°C until testing. We investigated 320 SNPs selected from the PADI locus on Immunochip and from other studies [8,19-21,29]. The list of the genotyped PADI SNPs is presented in Table S1 in Additional File 1. The SNPs were genotyped either using the TaqMan SNP genotyping assay (Applied Biosystems, Foster City, CA, USA) or by the Illumina iSelect HD custom genotyping array (Immunochip, Illumina, Inc, San Diego, CA, USA). SNPs with call rates <95%, monomorphic SNPs, SNPs with minor allele frequency <0.005, and SNPs with Hardy-Weinberg equilibrium P <1.0 × 10-4 in RA cases and controls were excluded from statistical evaluation. This resulted in 227, 238, 263 and 248 SNPs passing quality-control filters for Malay, Chinese, Indian and other or mixed ethnicities, respectively. Genetic outliers for each ethnic group (that is, Malay, four cases and two controls; Chinese, five cases and four controls; Indian, 10 cases) were removed after principal component analysis to correct for possible population stratification.

Additional file 1. Table S1 presenting a list of PADI SNPs investigated in the MyEIRA study population. A complete set of 320 SNPs selected from the PADI locus on Immunochip and from other studies. The SNPs were genotyped either using the TaqMan SNP genotyping assay (Applied Biosystems, USA) or by Illumina iSELECT HD custom genotyping array (Immunochip). Table S2 presenting the haplotype frequencies and meta-analysis of PADI4 polymorphisms in the MyEIRA study by ACPA status. The haplotype analysis and meta-analysis of PADI4 polymorphisms in the MyEIRA study were performed in different subsets of RA defined by ACPA status. Bold results indicate significant association between the PADI haplotypes and subsets of RA. Figure S1 showing the regional association plots with recombination rate on the PADI genes for the three major ethnic groups from the MyEIRA study. Regional association plots on the PADI genes including PADI1, PADI2, PADI3, PADI4 and PADI6 for the three major ethnic groups from MyEIRA study showing the peak association in each ethnic group. Graphs centered on the most significant SNP in each ethnic group. The r2 values (linkage disequilibrium between the most significant SNP and the rest of SNPs in the region) are calculated on the MyEIRA data and the recombination rates are based on the International HapMap CHB+JPT data.

Format: DOC Size: 1.1MB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

To assess genotyping robustness, comparisons were made between the TaqMan assay and the Illumina iSELECT HD custom genotyping array (Immunochip) assay regarding PADI4_94 (rs2240340), which resulted in a 100% match of the genotyping calls for both RA cases and controls.

Genotyping for HLA-DRB1 shared epitope (SE) alleles - defined by DRB1*01:01, DRB1*01:02, DRB1*01:07, DRB1*04:01, DRB1*04:04, DRB1*04:05, DRB1*04:08, DRB1*04:10, DRB1*10:01 and DRB1*10:03 alleles - was performed by HLA-DRB1 high-resolution sequence-specific oligonucleotide PCR for all the cases and controls as described previously [6]. Individuals carrying one or two SE alleles were classified as SE-positive [6].

Autoantibody measurements

ACPA and rheumatoid factors were identified and measured with ELISA kits (Immunoscan RA Mark 2; Euro-Diagnostica, Malmö, Sweden) as described elsewhere [6]. Antibody levels ≥25 AU/ml and >15 IU/ml were regarded as ACPA-positive and rheumatoid factor-positive, respectively.

Statistical analysis

Genotype, allele and haplotype frequencies were assessed with Yate's chi-square and/or Fisher's exact test when appropriate by means of IBM SPSS Statistics 20.0 software (SPSS Inc., Chicago, USA). The frequencies of the alleles and genotypes of PADI SNPs were compared between RA cases and control subjects and ORs with 95% CIs were calculated. Haplotype analysis was carried out by Haploview [30]. Power calculation was performed for a one-tail test at a significance level of 0.05. For meta-analysis, the Mantel-Haenszel method was employed with a fixed-effects model and 95% CI for cumulative or overall odds ratio [31]. The significance of the cumulative OR was determined by the Z-test. The between-study heterogeneity was assessed using the Cochran Q-statistic (P <0.10 considered significant). In addition, the I2 metric [I2 = (Q - df)/Q] was used to describe the percentage of variation across the studies due to heterogeneity. I2 values of 25%, 50% and 75% were assigned as low, moderate and high estimates, respectively.

Results

Characteristics of RA cases and controls

Among the RA cases, the overall mean ± standard deviation age was 48 ± 11.6 years, 86.3% were female, 64.5% were ACPA-positive, 40.1% were HLA-DRB1 SE-positive, and 51.5% were rheumatoid factor-positive. The mean ± standard deviation age of control subjects was 47 ± 11.4 years and 87.5% were female. The distribution of ethnic groups (Malay, Chinese, Indian and other or mixed ethnicities) is shown in Table 1.

Evaluation of the PADI gene polymorphisms and RA association in the MyEIRA study

Collectively, our study has 88% statistical power to detect genotype/minor allele frequency differences of the magnitude reported in the initial positive study in a Japanese population [8].

Overall, the single-point analyses of variations in the PADI1, PADI2, PADI3, PADI4 and PADI6 genes revealed a modest genetic effect size in the RA population and showed that the peak association varies between different ethnic groups. For example, the rs2526839 variant showed the peak association signal in Malay RA patients (OR = 1.24, 95% CI = 1.04 to 1.43, P = 0.0126), the rs3003444 variant in Chinese RA patients (OR = 1.60, 95% CI = 1.20 to 2.13, P = 0.0013) and the rs113475583 variant in Indian RA patients (OR = 1.91, 95% CI = 1.19 to 3.07, P = 0.0064) (see Figure S1 in Additional File 1).

PADI2 gene polymorphism as a risk factor in RA development

The PADI2 gene encodes PAD2 enzyme, which is the most widely expressed family member of PAD. An association between RA and the PADI2 variant was found in a Korean population [32]. We undertook this study for further investigation of whether PADI2 polymorphisms are also at risk for RA development in three independent Asian populations. Interestingly, the PADI2 variant meta-analysis demonstrated a possible novel association of the PADI2 rs1005753 variant with RA, both in the allelic model (ORoverall = 0.88, 95% CI = 0.77 to 0.99, P = 0.04) and in the dominant/recessive genotype model (ORoverall = 0.84, 95% CI = 0.72 to 0.99, P = 0.04) (Figure 1). In this multiethnic case-control study, no evidence of heterogeneity was found for the PADI2 rs1005753 variant (Pheterogeneity = 0.43, I2 = 0%).

thumbnailFigure 1. Allelic association between PADI2 rs1005753 polymorphism and risk of rheumatoid arthritis in the MyEIRA study. Meta-analysis of the minor allele (rs1005753_G) versus the common allele (rs1005753_T). CI, confidence interval; MH, Mantel-Haenszel; phet, P value for heterogeneity; I2, degree of heterogeneity (%); MyEIRA, Malaysian Epidemiological Investigation of Rheumatoid Arthritis; PADI, peptidylarginine deiminase.

Generalizability of PADI4 gene polymorphisms as a risk factor in RA development

A summary of the MyEIRA meta-analysis for the PADI4 polymorphism with RA is presented in Figure 2. The MyEIRA meta-analysis suggests an association of the PADI4 rs2240340 variant with RA also found in the Malaysian population (ORoverall = 1.11, 95% CI = 1.00 to 1.23, P = 0.04). The meta-analysis also revealed a significantly increased OR (ORoverall = 1.20, 95% CI = 1.01 to 1.44, P = 0.04) when applying a dominant/recessive genotype model for the PADI4 rs2240340 variant, suggesting that the PADI4 is associated with RA in the Malaysian population of Asian descent.

thumbnailFigure 2. Allelic association between PADI4 rs2240340 polymorphism and risk of rheumatoid arthritis in the MyEIRA study. Meta-analysis of the minor allele (rs2240340_T) versus the common allele (rs2240340_C). CI, confidence interval; MH, Mantel-Haenszel; phet, P value for heterogeneity; I2, degree of heterogeneity (%); MyEIRA = Malaysian Epidemiological Investigation of Rheumatoid Arthritis; PADI, peptidylarginine deiminase.

To clarify the role of the PADI4 gene as a possible susceptibility factor for RA development in different Asian ethnic populations, we performed a combined meta-analysis using our current data together with the previously published studies from Asia. We searched information from the PubMed database, ISI Web of Knowledge and Google. Six published original case-control studies related to PADI4 polymorphisms in RA in three different Asian populations were identified: three Japanese studies including the first positive report, one Korean study and two Chinese studies [8-11,16,17]. All studies had used the defined American College of Rheumatology criteria [28] for the diagnosis of RA as compared with our study. The strongest evidence of association reported in the first study was given by PADI4_94 (rs2240340); we therefore restricted the current combined meta-analysis to this genetic variant. Of the six published papers, the Korean study did not genotype for PADI4 rs2240340. Instead, we used genotype data of PADI4_89 (rs11203366). This choice was made because these SNPs are in strong linkage disequilibrium (r2 = 1.00) and the minor allele frequency of these SNPs are nearly equal according to the International Human Genome Project Chinese and Japanese populations (International HapMap CHB + JPT) (Asian population). Five studies were included in the combined meta-analysis together with our new data for PADI4 rs2240340 [8-11,16]. The PADI4 SNPs investigated in Fan and colleagues' study were not analyzed because they exhibited significant deviation from Hardy-Weinberg equilibrium (P <0.001) [17], which questions the validity of the genotyping data.

In total, the combined meta-analysis of the PADI4 rs2240340 included 5,192 RA cases and 4,317 control subjects from Asian populations. The Mantel-Haenszel method with fixed effects exhibited significant ORs in both the allelic model (OR = 1.23, 95% CI = 1.16 to 1.31, P <10-4) and the dominant/recessive genotype model (OR = 1.31, 95% CI = 1.20 to 1.44, P <10-4) (Figure 3a, b). A moderate level of between-study heterogeneity was observed in the allele model (Pheterogeneity = 0.08, I2 = 44%), but not in the genotype model (Pheterogeneity = 0.32, I2 = 14%). With the exclusion of the initial positive report, which tends to suggest a stronger genetic effect (winner's curse) than those in subsequent studies [9-11,16], the cumulative OR decreased slightly but remained highly significant (P <10-4) with no significant heterogeneity (Pheterogeneity = 0.47, I2 = 0%) (see Figure 3 caption).

thumbnailFigure 3. Allelic and genotypic association between PADI4 rs2240340 polymorphism and rheumatoid arthritis in Asian populations. Combined meta-analysis of (a) the allelic model and (b) the genotypic model between PADI4 rs2240340 polymorphism and risk of developing RA in Asian populations. (a) *If the first Japanese study was excluded, overall odds ratio (ORoverall) = 1.21 (1.13 to 1.29), phet = 0.15, I2 = 35%; **if the Korean study was excluded as the PADI4 rs2240340 SNP was not investigated, ORoverall = 1.21 (1.13 to 1.29), phet = 0.23, I2 = 25%; if both the first Japanese study and the Korean study were excluded, ORoverall = 1.17 (1.09 to 1.25), phet = 0.58, I2 = 0%. (b) *If the first Japanese study was excluded, ORoverall = 1.27 (1.15 to 1.40), phet = 0.47, I2 = 0%; **if the Korean study was excluded as the PADI4 rs2240340 SNP was not investigated, ORoverall = 1.28 (1.16 to 1.40), phet = 0.58, I2 = 0%; if both the first Japanese study and the Korean study were excluded, ORoverall = 1.22 (1.10 to 1.35), phet = 0.93, I2 = 0%. Minor allele (rs2240340_T) versus common allele (rs2240340_C). CI, confidence interval; MH, Mantel-Haenszel; phet, P value for heterogeneity; I2, degree of heterogeneity expressed as percentage; RA, rheumatoid arthritis.

Stratified analysis

Stratified analyses were performed to detect association between the PADI gene variants and RA risk within specific population subgroups as well as in relation to more homogeneous subsets of RA cases. Stratifying by sex, carriage of SE alleles, presence of ACPA or presence of rheumatoid factors among the cases revealed no evidence of different associations between the compared subgroups (data not shown). Nevertheless, it is noteworthy that a set of different SNPs variants were associated with risk of developing ACPA-positive and ACPA-negative RA within specific subpopulations, although the results were not significant after correction for multiple comparisons (Table 2). It would be interesting to investigate the association between RA and these SNP variants in other populations with more data.

Table 2. Significant SNP variants associated with ACPA-positive and ACPA-negative RA in the MyEIRA study population

Haplotype analysis and association test

We first performed the univariate single-population analyses for the association of PADI SNPs with RA, and subsequently those SNPs with statistical significant effects were selected for further haplotype analysis. However, we excluded the SNPs with inconsistent effects within the four ethnic groups, and also those SNPs with a significant level not low enough (P <0.05) in one or more ethnic groups. This analysis strategy led us to four PADI4 SNPs for haplotype analysis (rs79907974, rs2240340, rs1748021 and rs2240337). When haplotypes were constructed using these four PADI4 SNPs, seven haplotypes with frequency >1% out of 16 expected haplotypes were found in the study populations (Table 3). Haplotype GCGG was omitted from the association test because it was only found in the Indian ethnic group with a frequency of 5%. When meta-analyses for different ethnic groups in our study were conducted, ATAA and ATAG showed a significant difference between RA cases and controls (P = 0.02 and P = 0.001, respectively; Table 3) with no significant heterogeneity (Pheterogeneity = 0.78, I2 = 0%). These haplotypes represent either protective (ATAA) or susceptible (ATAG) variants. We further tested these haplotype associations with different RA subsets defined by ACPA status. The results showed consistent effects in the RA subsets, and the ORoverall was comparable between ACPA-positive and ACPA-negative RA subsets (see Table S2 in Additional File 1).

Table 3. Haplotype frequencies and meta-analysis of PADI4 polymorphisms in relation to RA in the MyEIRA study

Discussion

In this study, we used a multiethnic population of Asian descent in Malaysia to determine the association between PADI4 polymorphisms and risk of RA, and convincingly validated this association. Our current data, together with the previously published data on Asian populations, strongly support PADI4 as a RA susceptibility gene in different ethnic populations of Asian descent. Collectively, our data extend previous results on PADI4 and RA based on Asian populations, mainly in Japanese and Korean populations [8-11,15,33], to be observed in another population of Asian origin (that is, Malay, Chinese and Indian ethnicity). Notably, a possible novel association between the PADI2 genetic variant and RA risk was also found in the multiethnic Malaysian population.

The combined meta-analysis was performed regardless of the ACPA status, since such data were not available in the previous published papers. However, our present study showed a comparable minor allele frequency between the RA subsets defined by ACPA status, thus suggesting that the risk from PADI locus is likely to be common for the two main subgroups of RA.

Previously, the association studies of the PADI4 polymorphism and RA found in Asian populations were predominantly focused on Japanese and Korean populations [8-11,15], which are geographically and historically closely related and genetically quite similar to each other [34]. The results, however, were inconsistent across two independent studies in Han Chinese populations [16,17]. Malaysia is a multiethnic country in South East Asia representing genetic diversity across multiple large ethnic populations of Asian origin (that is, Malays, Chinese and Indians). In our study, we were able to address the question of whether PADI4 polymorphisms confer a risk of RA in this diverse population. We found that the PADI4_94 polymorphism is associated with an increased risk of developing RA in this population. Interestingly, as can be seen from Figure 3, the effect for the Han Chinese population is probably lower than for other Asian populations, which may be an explanation for previous inconsistency in published data.

Haplotype analysis did not reveal a higher effect size in comparison with univariate SNP analyses in our study. Nevertheless, meta-analysis of PADI4 haplotypes with the current data implied that association with RA was probably driven by the rs2240337 variant. This genetic variant, however, has low minor allele frequency in our materials (that is, minor allele frequency in Malay = 0.059, Chinese = 0.066, Indian = 0.032 and others = 0.042). Noteworthy is that this SNP variant, on the contrary, is in complete linkage disequilibrium (D' = 1.00) with the PADI4 SNP rs766499 variant, which was reported recently to be associated with the Japanese RA population at a genome-wide level of significance in meta-analysis [12].

An additional novel finding in our study is the discovery of a possible association between the PADI2 genetic variant (rs1005753) and RA, which was observed both in allele and genotype models. The genetic effect observed is possibly due to an indirect association, where the identified genetic variants (that is, PADI2) by themselves are not functional but are in linkage disequilibrium with a causal variant polymorphism such as the PADI4 gene. However, this would be unlikely as the PADI2 genetic variant (rs1005753) did not show any linkage disequilibrium with the numerous investigated SNPs spanning between PADI1 and PADI6 in the present study. We observed r2 <0.01 for all ethnic groups studied. A previous study by Freudenberg and colleagues reported an association between the PADI2 genetic variant (rs2075696) and RA in a Korean study [32]. Interestingly, the rs1005753 variant associated in our study and the rs2075696 variant are located in two different linkage disequilibrium blocks in the PADI2 gene locus. The International HapMap CHB+JPT project showed no linkage disequilibrium relationship between these two SNPs (D' = 0.0050, r2 = 0.0). Together, our data may suggest a possible risk effect between the rs1005753 variant and RA in the Malaysian population.

Worthy of mention is that PADI2 and PADI4 are the only two genes that are highly expressed in hemapoetic cells [35]. In RA, the expression levels of PAD2 and PAD4 were correlated with the intensity of inflammation, and both enzymes were demonstrable within or in the vicinity of citrullinated fibrins deposits [36]. The PADI2 gene encodes the PAD2 enzyme, is perhaps playing a role in RA pathogenesis on its own, independently from PADI4. Nevertheless, it is important to perform further extension of the study to gain better statistical power to investigate this association in a single-population analysis and further to replicate the association in independent cohorts of Asian origin. Since the difference between Caucasian and Asian populations may represent both genetic and environmental heterogeneity, it is logical to propose a gene-environmental interaction study for PADI genes in the development of RA.

Conclusion

This study demonstrates an association between PADI4 polymorphisms and RA in the Malaysian population. The currently updated combined meta-analysis of Asian populations further supports the hypothesis that the PADI locus contributes to the development of RA in different Asian populations and that this genetic effect is generalized to multiple ethnic populations of Asian descent.

Abbreviations

ACPA: anti-citrullinated protein antibody; CI: confidence interval; ELISA: enzyme-linked immunosorbent assay; HLA: human leukocyte antigen; I2: degree of heterogeneity; International HapMap CHB+JPT: International Human Genome Project Chinese and Japanese; MyEIRA: Malaysian Epidemiological Investigation of Rheumatoid Arthritis; OR: odds ratio; PADI: peptidylarginine deiminase; PCR: polymerase chain reaction; RA: rheumatoid arthritis; SE: shared epitope; SNP: single nucleotide polymorphism.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

TCL and LP had full access to all data in this study and take responsibility for the integrity of the data and the accuracy of the data analysis. TCL and LP drafted the manuscript, performed statistical analysis and manuscript preparation. TCL, SM, JSD, PL, XJ, BD, LA, LK and LP conceived the study and participated in the design of the study and manuscript editing. TCL and LP performed serological and molecular genetics assay. TCL and LP take responsibility for the acquisition of data, analysis and interpretation of data. All authors were involved in revising the paper critically for important intellectual content, and all authors approved the final version to be published.

Acknowledgements

The authors would like to thank the Director General of Health, Ministry of Health, Malaysia for supporting the work described in this paper and permission to publish this paper. The authors are also indebted to patients and controls for their kind participation. They also thank Nor Asiah Muhamad, Nor Aini Abdullah, Amal Nasir Mustafa, Abqariyah Yahya, Hanjeet Kaur, Kee Chee Cheong, Kasthoori Supramaniam, Masita Arip, Fatimah Kamila Abu Bakar, Ahneez Abdul Hameed, Norhazlin Mustafa and Ainun Yusniza, members of the MyEIRA study - Heselynn Hussien, Nor Suhaila Shahril, M Eashwary (Putrajaya Hospital, Wilayah Persekutuan Kuala Lumpur), Wahinuddin Sulaiman, Muhaini Othman (Tengku Bainun Hospital, Perak), Azmillah Rosman (Selayang Hospital, Selangor), Suk Chyn Gun (Tunku Jaafar Hospital, Negeri Sembilan), Mohd Shahrir M (Hospital Universiti Kebangsaan Malaysia), Ing Soo Lau (Sultanah Aminah Hospital, Johor), and Ainon Mokhtar (Tengku Ampuan Afzan Hospital, Pahang) - for the collection of data, and Maria Seddighzadeh for technical assistance. Help with Immunochip genotyping from Dr JK Sandling and Prof AC Syvänen is highly appreciated. The authors would also like to thank Dr RM Plenge for helpful discussion of the study results. In addition, they would like to thank the two anonymous reviewers for their insightful comments that helped to improve the quality of this paper. This study was supported by the Swedish Medical Research Council, the Swedish COMBINE (Vinnova) program and by the Ministry of Health, Malaysia: MRG 7/2005 and JPP-IMR 08-012 (NMRR-08-820-1975).

References

  1. Gregersen PK, Silver J, Winchester RJ: The shared epitope hypothesis. An approach to understanding the molecular genetics of susceptibility to rheumatoid arthritis.

    Arthritis Rheum 1987, 30:1205-1213. PubMed Abstract | Publisher Full Text OpenURL

  2. Raychaudhuri S, Sandor C, Stahl EA, Freudenberg J, Lee HS, Jia X, Alfredsson L, Padyukov L, Klareskog L, Worthington J, Siminovitch KA, Bae SC, Plenge RM, Gregersen PK, de Bakker PI: Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis.

    Nat Genet 2012, 44:291-296. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, Thomson BP, Li Y, Kurreeman FA, Zhernakova A, Hinks A, Guiducci C, Chen R, Alfredsson L, Amos CI, Ardlie KG, Barton A, Bowes J, Brouwer E, Burtt NP, Catanese JJ, Coblyn J, Coenen MJ, Costenbader KH, Criswell LA, Crusius JB, Cui J, de Bakker PI, De Jager PL, Ding B, Emery P, et al.: Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci.

    Nat Genet 2010, 42:508-514. PubMed Abstract | Publisher Full Text OpenURL

  4. Wakitani S, Imoto K, Murata N, Oonishi H, Ochi T, Yoneda M: An association between the natural course of shoulder joint destruction in rheumatoid arthritis and HLA-DRB1*0405 in Japanese patients.

    Scand J Rheumatol 1998, 27:146-148. PubMed Abstract | Publisher Full Text OpenURL

  5. Lee HS, Lee KW, Song GG, Kim HA, Kim SY, Bae SC: Increased susceptibility to rheumatoid arthritis in Koreans heterozygous for HLA-DRB1*0405 and *0901.

    Arthritis Rheum 2004, 50:3468-3475. PubMed Abstract | Publisher Full Text OpenURL

  6. Chun-Lai T, Padyukov L, Dhaliwal JS, Lundstrom E, Yahya A, Muhamad NA, Klareskog L, Alfredsson L, Larsson PT, Murad S: Shared epitope alleles remain a risk factor for anti-citrullinated proteins antibody (ACPA)-positive rheumatoid arthritis in three asian ethnic groups.

    PLoS One 2011, 6:e21069. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Lundstrom E, Kallberg H, Smolnikova M, Ding B, Ronnelid J, Alfredsson L, Klareskog L, Padyukov L: Opposing effects of HLA-DRB1*13 alleles on the risk of developing anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis.

    Arthritis Rheum 2009, 60:924-930. PubMed Abstract | Publisher Full Text OpenURL

  8. Suzuki A, Yamada R, Chang X, Tokuhiro S, Sawada T, Suzuki M, Nagasaki M, Nakayama-Hamada M, Kawaida R, Ono M, Ohtsuki M, Furukawa H, Yoshino S, Yukioka M, Tohma S, Matsubara T, Wakitani S, Teshima R, Nishioka Y, Sekine A, Iida A, Takahashi A, Tsunoda T, Nakamura Y, Yamamoto K: Functional haplotypes of PADI4, encoding citrullinating enzyme peptidylarginine deiminase 4, are associated with rheumatoid arthritis.

    Nat Genet 2003, 34:395-402. PubMed Abstract | Publisher Full Text OpenURL

  9. Ikari K, Kuwahara M, Nakamura T, Momohara S, Hara M, Yamanaka H, Tomatsu T, Kamatani N: Association between PADI4 and rheumatoid arthritis: a replication study.

    Arthritis Rheum 2005, 52:3054-3057. PubMed Abstract | Publisher Full Text OpenURL

  10. Kang CP, Lee HS, Ju H, Cho H, Kang C, Bae SC: A functional haplotype of the PADI4 gene associated with increased rheumatoid arthritis susceptibility in Koreans.

    Arthritis Rheum 2006, 54:90-96. PubMed Abstract | Publisher Full Text OpenURL

  11. Takata Y, Inoue H, Sato A, Tsugawa K, Miyatake K, Hamada D, Shinomiya F, Nakano S, Yasui N, Tanahashi T, Itakura M: Replication of reported genetic associations of PADI4, FCRL3, SLC22A4 and RUNX1 genes with rheumatoid arthritis: results of an independent Japanese population and evidence from meta-analysis of East Asian studies.

    J Hum Genet 2008, 53:163-173. PubMed Abstract | Publisher Full Text OpenURL

  12. Okada Y, Terao C, Ikari K, Kochi Y, Ohmura K, Suzuki A, Kawaguchi T, Stahl EA, Kurreeman FA, Nishida N, Ohmiya H, Myouzen K, Takahashi M, Sawada T, Nishioka Y, Yukioka M, Matsubara T, Wakitani S, Teshima R, Tohma S, Takasugi K, Shimada K, Murasawa A, Honjo S, Matsuo K, Tanaka H, Tajima K, Suzuki T, Iwamoto T, Kawamura Y, et al.: Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population.

    Nat Genet 2012, 44:511-516. PubMed Abstract | Publisher Full Text OpenURL

  13. Klareskog L, Ronnelid J, Lundberg K, Padyukov L, Alfredsson L: Immunity to citrullinated proteins in rheumatoid arthritis.

    Annu Rev Immunol 2008, 26:651-675. PubMed Abstract | Publisher Full Text OpenURL

  14. Schellekens GA, de Jong BA, van den Hoogen FH, van de Putte LB, van Venrooij WJ: Citrulline is an essential constituent of antigenic determinants recognized by rheumatoid arthritis-specific autoantibodies.

    J Clin Invest 1998, 101:273-281. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Bang SY, Han TU, Choi CB, Sung YK, Bae SC, Kang C: Peptidyl arginine deiminase type IV (PADI4) haplotypes interact with shared epitope regardless of anti-cyclic citrullinated peptide antibody or erosive joint status in rheumatoid arthritis: a case control study.

    Arthritis Res Ther 2010, 12:R115. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  16. Chen R, Wei Y, Cai Q, Duan S, Ren D, Shen J, He D, Fang M, Lv K, Cheng N, Sun N: The PADI4 gene does not contribute to genetic susceptibility to rheumatoid arthritis in Chinese Han population.

    Rheumatol Int 2010, 31:1631-1634. PubMed Abstract | Publisher Full Text OpenURL

  17. Fan LY, Zong M, Lu TB, Yang L, Ding YY, Ma JW: [Association of the PADI4 gene polymorphism and HLA-DRB1 shared epitope alleles with rheumatoid arthritis].

    Zhonghua Yi Xue Yi Chuan Xue Za Zhi 2009, 26:57-61. PubMed Abstract OpenURL

  18. Hoppe B, Haupl T, Gruber R, Kiesewetter H, Burmester GR, Salama A, Dorner T: Detailed analysis of the variability of peptidylarginine deiminase type 4 in German patients with rheumatoid arthritis: a case-control study.

    Arthritis Res Ther 2006, 8:R34. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  19. Plenge RM, Padyukov L, Remmers EF, Purcell S, Lee AT, Karlson EW, Wolfe F, Kastner DL, Alfredsson L, Altshuler D, Gregersen PK, Klareskog L, Rioux JD: Replication of putative candidate-gene associations with rheumatoid arthritis in >4,000 samples from North America and Sweden: association of susceptibility with PTPN22, CTLA4, and PADI4.

    Am J Hum Genet 2005, 77:1044-1060. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Harney SM, Meisel C, Sims AM, Woon PY, Wordsworth BP, Brown MA: Genetic and genomic studies of PADI4 in rheumatoid arthritis.

    Rheumatology (Oxford) 2005, 44:869-872. Publisher Full Text OpenURL

  21. Martinez A, Valdivia A, Pascual-Salcedo D, Lamas JR, Fernandez-Arquero M, Balsa A, Fernandez-Gutierrez B, de la Concha EG, Urcelay E: PADI4 polymorphisms are not associated with rheumatoid arthritis in the Spanish population.

    Rheumatology (Oxford) 2005, 44:1263-1266. Publisher Full Text OpenURL

  22. Farago B, Talian GC, Maasz A, Magyari L, Horvatovich K, Kovacs B, Cserep V, Kisfali P, Kiss CG, Czirjak L, Melegh B: Prevalence of functional haplotypes of the peptidylarginine deiminase citrullinating enzyme gene in patients with rheumatoid arthritis: no influence of the presence of anti-citrullinated peptide antibodies.

    Clin Exp Rheumatol 2007, 25:523-528. PubMed Abstract | Publisher Full Text OpenURL

  23. Burr ML, Naseem H, Hinks A, Eyre S, Gibbons LJ, Bowes J, Wilson AG, Maxwell J, Morgan AW, Emery P, Steer S, Hocking L, Reid DM, Wordsworth P, Harrison P, Thomson W, Worthington J, BIRAC Consortium; YEAR Consortium, Barton A: PADI4 genotype is not associated with rheumatoid arthritis in a large UK Caucasian population.

    Ann Rheum Dis 2010, 69:666-670. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Vossenaar ER, Radstake TR, van der Heijden A, van Mansum MA, Dieteren C, de Rooij DJ, Barrera P, Zendman AJ, van Venrooij WJ: Expression and activity of citrullinating peptidylarginine deiminase enzymes in monocytes and macrophages.

    Ann Rheum Dis 2004, 63:373-381. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Kinloch A, Lundberg K, Wait R, Wegner N, Lim NH, Zendman AJ, Saxne T, Malmstrom V, Venables PJ: Synovial fluid is a site of citrullination of autoantigens in inflammatory arthritis.

    Arthritis Rheum 2008, 58:2287-2295. PubMed Abstract | Publisher Full Text OpenURL

  26. Baeten D, Peene I, Union A, Meheus L, Sebbag M, Serre G, Veys EM, De Keyser F: Specific presence of intracellular citrullinated proteins in rheumatoid arthritis synovium: relevance to antifilaggrin autoantibodies.

    Arthritis Rheum 2001, 44:2255-2262. PubMed Abstract | Publisher Full Text OpenURL

  27. Makrygiannakis D, af Klint E, Lundberg IE, Lofberg R, Ulfgren AK, Klareskog L, Catrina AI: Citrullination is an inflammation-dependent process.

    Ann Rheum Dis 2006, 65:1219-1222. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA Jr, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis.

    Arthritis Rheum 1988, 31:315-324. PubMed Abstract | Publisher Full Text OpenURL

  29. Hoppe B, Haupl T, Egerer K, Gruber R, Kiesewetter H, Salama A, Burmester GR, Dorner T: Influence of peptidylarginine deiminase type 4 genotype and shared epitope on clinical characteristics and autoantibody profile of rheumatoid arthritis.

    Ann Rheum Dis 2009, 68:898-903. PubMed Abstract | Publisher Full Text OpenURL

  30. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps.

    Bioinformatics 2005, 21:263-265. PubMed Abstract | Publisher Full Text OpenURL

  31. Review Manager (RevMan) [Computer program]: Version 5.1. Copenhagen: The Nordic Cochrane Centre, The Cochrane Collaboration; 2011.

  32. Freudenberg J, Lee HS, Han BG, Shin HD, Kang YM, Sung YK, Shim SC, Choi CB, Lee AT, Gregersen PK, Bae SC: Genome-wide association study of rheumatoid arthritis in Koreans: population-specific loci as well as overlap with European susceptibility loci.

    Arthritis Rheum 2011, 63:884-893. PubMed Abstract | Publisher Full Text OpenURL

  33. Cha S, Choi CB, Han TU, Kang CP, Kang C, Bae SC: Association of anti-cyclic citrullinated peptide antibody levels with PADI4 haplotypes in early rheumatoid arthritis and with shared epitope alleles in very late rheumatoid arthritis.

    Arthritis Rheum 2007, 56:1454-1463. PubMed Abstract | Publisher Full Text OpenURL

  34. Kim JJ, Verdu P, Pakstis AJ, Speed WC, Kidd JR, Kidd KK: Use of autosomal loci for clustering individuals and populations of East Asian origin.

    Hum Genet 2005, 117:511-519. PubMed Abstract | Publisher Full Text OpenURL

  35. Vossenaar ER, Zendman AJ, van Venrooij WJ, Pruijn GJ: PAD, a growing family of citrullinating enzymes: genes, features and involvement in disease.

    Bioessays 2003, 25:1106-1118. PubMed Abstract | Publisher Full Text OpenURL

  36. Foulquier C, Sebbag M, Clavel C, Chapuy-Regaud S, Al Badine R, Mechin MC, Vincent C, Nachat R, Yamada M, Takahara H, Simon M, Guerrin M, Serre G: Peptidyl arginine deiminase type 2 (PAD-2) and PAD-4 but not PAD-1, PAD-3, and PAD-6 are expressed in rheumatoid arthritis synovium in close association with tissue inflammation.

    Arthritis Rheum 2007, 56:3541-3553. PubMed Abstract | Publisher Full Text OpenURL