Rheumatology Advance Access originally published online on May 11, 2006
Rheumatology 2006 45(11):1345-1348; doi:10.1093/rheumatology/kel169
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Haplotype analysis revealed no association between the PTPN22 gene and RA in a Japanese population
Institute of Rheumatology, Tokyo Women's Medical University, Tokyo, Japan.
Correspondence to: Katsunori Ikari, MD, PhD, Institute of Rheumatology, Tokyo Women's Medical University, 10-22 Kawada, Shinjuku, Tokyo 162-0054, Japan. E-mail: kikari{at}ior.twmu.ac.jp
| Abstract |
|---|
|
|
|---|
Objective. The protein tyrosine phosphatase non-receptor type 22 (PTPN22) gene is a member of the PTPs that negatively regulate T-cell activation. A missense single nucleotide polymorphism (SNP) in the PTPN22 gene known as R620W was recently reported to be associated with several autoimmune diseases including rheumatoid arthritis (RA). The association was confirmed repeatedly in the populations of North European ancestry. However, the SNP was reported to be non-polymorphic in the Asian populations. Because the gene confers an impact on autoimmune diseases, we attempt to explore an association between PTPN22 gene and RA in a Japanese population without restricting to the SNP, R620W.
Methods. We studied 1128 RA patients and 455 controls. In addition to the SNP, R620W, we selected eight testing SNPs spanning 45 kb over the PTPN22 gene using the International HapMap Project. Genotyping was performed using the TaqMan fluorogenic 5' nuclease assay. Associations between RA and each of the SNPs were estimated by the Fisher's exact test. Haplotype was constructed using the expectation-maximization algorithm.
Results. R620W was not polymorphic enough in both the patients and the controls, and was therefore excluded from further analysis. Each allele frequency for the eight other SNPs in both groups was compared and no association was detected. Haplotype analysis also revealed that PTPN22 gene was not associated with RA in a Japanese population.
Conclusion. We found no association between PTPN22 and RA in a Japanese population. The result suggests that the PTPN22 gene is associated with RA only in a specific ethnic group.
KEY WORDS: Rheumatoid arthritis, PTPN22, Association, Haplotype, IORRA
| Introduction |
|---|
|
|
|---|
Autoimmune disease (MIM 109100 [OMIM] ) is one of the common diseases affecting up to 5% of the population. It is characterized by an abnormal immune response to self-antigens. Among all the systemic autoimmune diseases, rheumatoid arthritis (RA, MIM 180300 [OMIM] ) is the most common, afflicting up to 1% of the adult population worldwide. The disease susceptibility has been estimated to have a genetic component of 60% [1].
Bottini et al. [2] first described that a missense single nucleotide polymorphism (SNP) known as R620W (rs2476601) in the protein tyrosine phosphatase non-receptor type 22 (PTPN22) gene was associated with susceptibility to type I diabetes. Begovich et al. [3] also reported that this gene was associated with RA in a North American population. Since then, the association was confirmed repeatedly in the populations of the European ancestry [47]. Furthermore, the SNP was found to be associated with other autoimmune diseases, such as systemic lupus erythematosus, Graves disease and juvenile idiopathic arthritis [57].
T-cells play a central role in the immunopathogenesis of autoimmune diseases and are also the key regulators of the destructive joint lesions [8]. The PTPN22 gene is a member of the PTPs that are involved in the negative regulation of T-cell signalling through its interaction with C-terminal Src tyrosine kinase (Csk) [9]. The amino change in PTPN22 caused by the functional SNP disrupts the binding to Csk that leads to the overactivity of the immune system [2, 3]. Thus, this susceptibility gene with the functional polymorphism is thought as one of the main players to autoimmune diseases outside the human leucocyte antigen locus [10].
Although the association of the R620W SNP and autoimmune diseases was validated repeatedly in the North European descents including the pathogenetic role of PTPN22 in the diseases, the functional SNP was non-polymorphic and the association could not be confirmed in the Asian populations [3, 11]. Recently, Carlton et al. [12] described that an SNP in the 3' untranslated region of the gene (rs3789604) is also associated with RA, independent of R620W. This finding suggests that in addition to R620W, the association between the diseases and the PTPN22 gene in the Asian populations needs to be studied more extensively.
In this study, we attempt to find an association between PTPN22 gene and RA patients in a large Japanese population using the HapMap-tagged SNPs that covers the PTPN22 gene.
| Materials and methods |
|---|
|
|
|---|
Subjects and disease criteria
Approval for this study was granted by Tokyo Women's Medical University Genome Ethics Committee. The study was part of an observational cohort project that included over 4000 Japanese RA patients, established in the year 2000 by the Institute of Rheumatology, Tokyo Women's Medical University (IORRA: Institute of Rheumatology RA cohort). The diagnosis of RA was established using the classification criteria of the American College of Rheumatology (formerly, the American Rheumatism Association) 1987 revised criteria [13]. Out of the registered RA patients, DNA samples were obtained from 1284 patients. Informed written consent was obtained from every subject. Of these, 1128 samples were randomly selected for this study. About 88% of the patients were RF positive and they were mostly females (82.6%). A total of 455 population-based control DNA samples were obtained from the Pharma SNP Consortium (http://www.jpma.or.jp/psc/index.html). All control subjects were matched for sex, ethnic origin and geographical area.
Selection of SNPs
In addition to the functional R620W SNP, we selected eight SNPs that allowed us to describe the haplotypes detected in the international HapMap project (release 19 October 2005) (Table 1) [14]. We used SNP 3 (rs3765598) instead of rs3789604, a risk polymorphism for RA independent of R620W, for the study because the supplier (Applied Biosystems, Tokyo, Japan) was unable to manufacture the assay of rs3789604 [12]. SNP3 was almost in absolute linkage disequilibrium (LD) with rs3789604 and minor alleles of these SNPs were carried by a single haplotype according to the HapMap in Japanese samples (D' = 1.0, R2 = 0.94) [14]. We finally selected nine SNPs spanning 45 kb over the PTPN22 gene for genotyping.
|
Statistical power
Our study was designed to have >90% power at a 5% significance level to detect the odds ratio (OR) of 1.40 conferred by the risk allele of SNP 3 (rs3765598, 15.6% frequency in the controls) [12]. Statistical power was calculated using a web power calculator (http://calculators.stat.ucla.edu/powercalc/).
Genotyping
Genotyping was carried out using the TaqMan fluorogenic 5' nuclease assay (Applied Biosystems). The final volume of polymerase chain reaction (PCR) was 5 µl, containing 2 ng of genomic DNA and 2.5 µl TaqMan Universal PCR Master Mix, with 0.125 µl of 40x Assay Mix or 0.25 µl of 20x Assay Mix. Thermal cycle conditions were as follows: 50°C for 2 min to activate the uracil N-glycosylase and to prevent carry-over contamination, 95°C for 10 min to activate the DNA polymerase, followed by 40 cycles of 92°C for 15 s and 60°C for 1 min. All PCR were performed using 384-well plates by a Dual 384-Well GeneAmp PCR System 9700 (Applied Biosystems) and the endpoint fluorescent readings were performed on an ABI PRISM 7900 HT Sequence Detection System (Applied Biosystems). Duplicate samples and negative controls were included to ensure accuracy of genotyping.
Statistical analysis
Allele frequencies were estimated by the gene counting method. The exact test of HardyWeinberg equilibrium was used to compare the observed numbers of each genotype with those expected for a population in the HardyWeinberg equilibrium (http://cran.r-project.org/src/contrib/Descriptions/genetics.html). Associations between RA and each of the SNPs or haplotypes were estimated by the Fisher's exact test. These tests were implemented in the R software package version 2.0.1 (http://www.r-project.org/). The expectation-maximization algorithm implemented in the LDSUPPORT program was used for the LD analysis and the haplotype estimation [15].
| Results |
|---|
|
|
|---|
From the genotyping result, the functional R620W SNP was not polymorphic enough in both the patients and the controls, and therefore it was excluded from further analysis (minor allele frequency was 0.0022 and 0.0022, respectively). This result was consistent with the results reported previously [3, 11]. Allele frequencies for the eight other SNPs that cover the gene were in HardyWeinberg equilibrium in both the patients and the controls. Each allelic frequency of the SNP in both the groups was nearly equal and no association was detected when compared independently (OR = 0.971.06) (Table 1). Allelic frequency of SNP3 that was substitute for novel risk polymorphism of RA, rs3789604, was completely equal in patients and controls (P = 1.00, OR = 1.0) [12]. Stratifying the patients by the presence of RF or sex also revealed no evidence of association with RA (data not shown).
We calculated D' values for all SNP pairs to assess the LD across the PTPN22 gene (Fig. 1). The pairwise D' values in the gene were nearly 1 among almost all SNP pairs, indicating that the SNPs were highly associated with each other and the entire PTPN22 was contained within a single LD block. Haplotype analysis predicted five common (frequency >1%) haplotypes and revealed that the PTPN22 gene was not associated with RA in this Japanese cohort (Table 2).
|
|
| Discussion |
|---|
|
|
|---|
In our study, the association between the PTPN22 gene and RA was investigated using a large Japanese RA patient cohort. Our data revealed no association between RA and the PTPN22 gene not only in the functional R620W SNP but also in the eight HapMap tagged SNPs or haplotypes that cover the gene.
Lack of association might be due to a lack of statistical power. However, based on a recent study that detected novel risk allele for RA susceptibility, a statistical power of >90% could be achieved with 2256 chromosomes of patients and 910 chromosomes of controls (P = 0.05, OR = 1.40 and a 15.6% frequency of risk allele) [12]. The statistical power is sufficient to detect a risk greater than OR 1.40 and if the OR is lower the negative result might be due to a type II error. It might be noticed that the recent findings on peptidyl arginine deiminase, type IV gene (PADI4) as a risk of RA show a different magnitude of genetic risk according to the ethnic differences [16].
Ethnic differences may play a role in the conflicting results among the genetic association studies [11]. A replication study using a population with a different ethnicity allows us to know whether a reported association is caused by a common variance through different ethnicities. In the case of PTPN22 gene and RA, the disease responsible polymorphism R620W found in the populations of North European ancestry is extremely rare in the Asian populations, and so the R620W alone cannot explain the disease susceptibility in the Asian populations [3, 11]. Since the causal variant of a disease might differ between populations, the failure to replicate an association may reflect an incomplete analysis of the variation within the gene of interest. For this reason, we conducted an expanded analysis using the HapMap to clarify the association in a Japanese population. We still did not find evidence of the association between RA and each of the eight HapMap-tagged SNPs that covered the gene nor the haplotypes.
In conclusion, our study suggests that PTPN22 gene may be associated with RA only in a specific ethnic group.
| Acknowledgements |
|---|
|
|
|---|
We thank all the DNA donors for making this study possible. We are grateful to A. C. Tang for her assistance in preparing the manuscript and K. Arai for her technical efforts. We appreciate A. Taniguchi and other members of our Institute for their efforts on the cohort project. This work was supported by a grant provided by the Japan Rheumatism Foundation (to K.I.) and a Grants-in-Aid for Scientific Research from the Japan Society for the Promotion of Science (to S.M.).
The authors have declared no conflicts of interest.
| References |
|---|
|
|
|---|
- MacGregor AJ, Snieder H, Rigby AS, et al. (2000) Characterizing the quantitative genetic contribution to rheumatoid arthritis using data from twins. Arthritis Rheum 43:307.[CrossRef][Web of Science][Medline]
- Bottini N, Musumeci L, Alonso A, et al. (2004) A functional variant of lymphoid tyrosine phosphatase is associated with type I diabetes. Nat Genet 36:3378.[CrossRef][Web of Science][Medline]
- Begovich AB, Carlton VE, Honigberg LA, et al. (2004) A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am J Hum Genet 75:3307.[CrossRef][Web of Science][Medline]
- Ladner MB, Bottini N, Valdes AM, Noble JA. (2005) Association of the single nucleotide polymorphism C1858T of the PTPN22 gene with type 1 diabetes. Hum Immunol 66:604.[Web of Science][Medline]
- Hinks A, Barton A, John S, et al. (2005) Association between the PTPN22 gene and rheumatoid arthritis and juvenile idiopathic arthritis in a UK population: further support that PTPN22 is an autoimmunity gene. Arthritis Rheum 52:16949.[CrossRef][Web of Science][Medline]
- Orozco G, Sanchez E, Gonzalez-Gay MA, et al. (2005) Association of a functional single-nucleotide polymorphism of PTPN22, encoding lymphoid protein phosphatase, with rheumatoid arthritis and systemic lupus erythematosus. Arthritis Rheum 52:21924.[CrossRef][Web of Science][Medline]
- Smyth D, Cooper JD, Collins JE, et al. (2004) Replication of an association between the lymphoid tyrosine phosphatase locus (LYP/PTPN22) with type 1 diabetes, and evidence for its role as a general autoimmunity locus. Diabetes 53:30203.
[Abstract/Free Full Text] - Klimiuk PA, Yang H, Goronzy JJ, Weyand CM. (1999) Production of cytokines and metalloproteinases in rheumatoid synovitis is T cell dependent. Clin Immunol 90:6578.[CrossRef][Web of Science][Medline]
- Cloutier JF and Veillette A. (1996) Association of inhibitory tyrosine protein kinase p50csk with protein tyrosine phosphatase PEP in T cells and other hemopoietic cells. EMBO J 15:490918.[Web of Science][Medline]
- Deighton CM, Walker DJ, Griffiths ID, Roberts DF. (1989) The contribution of HLA to rheumatoid arthritis. Clin Genet 36:17882.[Web of Science][Medline]
- Mori M, Yamada R, Kobayashi K, Kawaida R, Yamamoto K. (2005) Ethnic differences in allele frequency of autoimmune-disease-associated SNPs. J Hum Genet 50:2646.[CrossRef][Web of Science][Medline]
- Carlton VE, Hu X, Chokkalingam AP, et al. (2005) PTPN22 genetic variation: evidence for multiple variants associated with rheumatoid arthritis. Am J Hum Genet 77:56781.[CrossRef][Web of Science][Medline]
- Arnett FC, Edworthy SM, Bloch DA, et al. (1988) The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum 31:31524.[Web of Science][Medline]
- Altshuler D, Brooks LD, Chakravarti A, Collins FS, Daly MJ, Donnelly P. (2005) A haplotype map of the human genome. Nature 437:1299320.[CrossRef][Medline]
- Kitamura Y, Moriguchi M, Kaneko H, et al. (2002) Determination of probability distribution of diplotype configuration (diplotype distribution) for each subject from genotypic data using the EM algorithm. Ann Hum Genet 66:18393.[CrossRef][Web of Science][Medline]
- Iwamoto T, Ikari K, Nakamura T, et al. (2006) Association between PADI4 and rheumatoid arthritis: a meta-analysis. Rheumatology [Epub ahead of print].
This article has been cited by other articles:
![]() |
J. Bowes and A. Barton Recent advances in the genetics of RA susceptibility Rheumatology, April 1, 2008; 47(4): 399 - 402. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Okamoto, H. Kaneko, C. Terai, and N. Kamatani Protective effect of A at position 168 in the type III promoter of the MHCIITA gene in systemic lupus erythematosus Ann Rheum Dis, September 1, 2007; 66(9): 1263 - 1264. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


