International Journal of Genomics
Volume 2017 (2017), Article ID 8753498, 5 pages
https://doi.org/10.1155/2017/8753498

Research Article

Association of PTPN22 Haplotypes (−1123G>C/+1858C>T) with Rheumatoid Arthritis in Western Mexican Population

Yeniley Ruiz-Noa,¹ Jorge Ramón Padilla-Gutiérrez,¹ Jorge Hernández-Bello,¹ Claudia Azucena Palafox-Sánchez,¹ Yeminia Valle,¹ Edith Oregón-Romero,¹ Ana Laura Pereira-Suárez,¹ Ana Guilaisne Bernard-Medina,² and José Francisco Muñoz-Valle¹

¹Instituto de Investigación en Ciencias Biomédicas, Centro Universitario de Ciencias de la Salud, Universidad de Guadalajara, Guadalajara, JAL, Mexico
²Servicio de Reumatología, OPD Hospital Civil de Guadalajara “Fray Antonio Alcalde”, Guadalajara, JAL, Mexico

Correspondence should be addressed to José Francisco Muñoz-Valle; biologiamolecular@hotmail.com

Received 30 March 2016; Revised 6 November 2016; Accepted 18 December 2016; Published 22 January 2017

Academic Editor: Mohamed Salem

Copyright © 2017 Yeniley Ruiz-Noa 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.

Abstract.

Rheumatoid arthritis (RA) is an autoimmune disease characterized by the presence of antibodies against cyclic citrullinated peptide (anti-CCP), a consequence of the breakdown of immune tolerance. The lymphoid tyrosine phosphatase (Lyp) protein has significant effects on maintenance of peripheral immune tolerance. Two polymorphic variants (−1123G>C and +1858C>T) at PTPN22 gene that encodes this protein have been associated with autoimmune disorders and found in strong linkage disequilibrium in Caucasian population. We evaluated whether PTPN22 haplotypes (−1123G>C/+1858C>T) are associated with anti-CCP antibodies, as well as susceptibility to RA in a Western Mexican population. A total of 315 RA patients and 315 control subjects (CS) were included. The polymorphisms were genotyped by PCR-RFLP and the anti-CCP antibodies were determined by ELISA. The PTPN22 polymorphisms were in strong linkage disequilibrium (D′ = 1.00 in CS). The susceptibility haplotype CT was significantly more frequent in RA patients than in CS (OR 2.18, 95% CI 1.15–4.16, ). No association between haplotypes and anti-CCP antibodies levels was observed. In conclusion, this study confirmed that −1123G>C and +1858C>T PTPN22 polymorphisms are in strong linkage disequilibrium and the CT haplotype is a susceptibility marker to RA in Western Mexico. However, the PTPN22 haplotypes are not associated with anti-CCP antibodies.

1. Introduction

Rheumatoid arthritis (RA) is a complex autoimmune disease, characterized by immune cell-mediated destruction of the joints. Genetic factors appear as major contributors to its pathogenesis [1]. Single nucleotide polymorphisms (SNP) in the protein tyrosine phosphatase nonreceptor type 22 gene (PTPN22) encoding the lymphoid protein tyrosine phosphatase LYP [2] were found to be associated with several autoimmune disorders [3–10]. LYP is a potent negative regulator of the lymphocytic activation; it acts at the initial steps of T-cell receptor (TCR) signaling [11].

PTPN22 +1858C>T SNP (rs2476601) has been consistently associated with RA [12–14]. This SNP consists of a single-base change resulting in a substitution of Arg to Trp in the 620 position at P1 domain of the phosphatase, which plays an important role in the interaction to Csk kinase in order to downregulate the activation signaling pathway, causing by this an hyperresponsive T and B cell phenotype [15].

Nevertheless, in Asian population, the polymorphic variant of this SNP was not found, suggesting that the presence of another variant is associated with the disease. In that study, the authors identified a novel promoter polymorphism, the −1123G>C (rs2488457), which was associated with acute-onset but not slow-onset type 1 diabetes (T1DM) in the Japanese and Korean population [16]. Similar studies in Caucasian and Asian patients revealed a strong association between −1123G>C SNP with RA [17, 18].

The haplotype constructions of both SNPs have been tested in several populations and strong linkage disequilibrium was found between them [16, 19–21]. The CT haplotype that includes the risk alleles from both PTPN22 polymorphisms (−1123G>C and +1858C>T) has shown a significant association with T1DM [16, 19, 20], RA [21], and juvenile idiopathic arthritis (JIA) [19].

The most common feature of RA is the presence of antibodies to citrullinated proteins/peptides (ACPAs). ACPAs comprise a group of antibodies that are highly specific for the disease; among them are anti-CCP antibodies. The presence of high levels of anti-CCP at early diagnosis predicts more pronounced radiographic progression, as has been demonstrated by several studies [22–24].

Based on this knowledge, the aim of our study was to investigate the association of PTPN22 haplotypes (−1123G>C/+1858C>T) with RA and their relationship with anti-CCP antibodies in patients from Western Mexico.

2. Material and Methods

2.1. RA Patients and Control Subjects

Our study enrolled three hundred fifteen RA patients and the same number of subjects with no medical record of autoimmune or chronic diseases, named as control subjects, all of them with Western Mexican origin. Patients were recruited from the Rheumatology Department of Hospital Civil “Fray Antonio Alcalde,” Guadalajara, Jalisco, Mexico, and they were classified according to the 2010 ACR/EULAR criteria [25]. Clinical and demographic characteristics of this group are shown in Table 1. We observed a female predominance in our RA population (93%); their mean age was years, and most patients had disease duration of ten years and mean disease activity score DAS28 of 4.8. All individuals signed the written informed consent and this study was approved by the ethical committee of Universidad de Guadalajara.

Table 1: Demographic and clinical characteristics of RA patients.

Variables
RA
()
CS
()

Demographic and clinic characteristics

Age
48 (18–82)
36 (18–74)

Sex%  

Males
7 (23)
40 (127)

Females
93 (292)
60 (188)

29 (71)
9 (29)

Disease status

Duration of disease
10 (0.1–52)
—

—

Spanish HAQ-DI
0.76 (0–2.28)
—

Autoantibodies

Rheumatoid factor
178 (15–717)
15 (0–62)

Negative (<20 
16 (50)
95 (300)

Low positive (20–59.9 
21 (65)
5 (15)

High positive (≥60 
63 (200)
0 (0)

Anti-CCP
119 (2–460)
1.5 (0.9–2.5)

Negative (<5 
17 (52)
100 (315)

Low positive (5.1–14.9 
9 (27)
0 (0)

High positive (≥15 
74 (236)
0 (0)

presented in median (p25–p75). provided in percentage and n. provided in mean ± SD. RA: rheumatoid arthritis; CS: control subject; Anti-CCP: anticyclic citrullinated peptide antibody; DAS28: disease activity score 28; Spanish HAQ-DI: Spanish version of the Health Assessment Questionnaire Disability Index.

2.2. PTPN22 −1123G>C and +1858C>T SNPs Genotyping

Genomic DNA was extracted from peripheral blood cells according to a modified Miller technique. Genotyping was performed using polymerase chain reaction followed by a restriction fragment length polymorphism (PCR-RFLP) technique for the PTPN22 −1123G>C and +1858C>T polymorphisms. PCR products were digested for one hour at 37°C with appropriate restriction enzymes and then resolved on a 6% polyacrylamide gel.

For the −1123G>C polymorphism a 205 bp fragment was amplified (following our previously reported method) [26] and digested with Sac I restriction enzyme (New England Biolabs, Ipswich, MA, USA). The polymorphic allele −1123C was cleaved in two fragments of 183 bp and 22 bp, respectively, while the −1123G allele remained intact. For the +1858C>T SNP we used previously reported conditions [14]. Amplified products of 412 bp were digested with Xcm I restriction enzyme (New England Biolabs, Ipswich, MA, USA); only the polymorphic +1858T allele was cleaved into two fragments of 246 bp and 166 bp. Results were confirmed by automatized sequencing of one random sample of each genotype (Applied Biosystems, USA).

2.3. Enzyme-Linked Immunoassays

Anti-CCP antibodies were measured in serum samples from RA patients and CS by enzyme-linked immunosorbent assay (ELISA) (Anti-CCP, DIASTAT™, Axis Shield Diagnostics, Dundee, UK). The cut-off level was 5 U/mL for semiquantitative analysis. The immunoassays were performed following the manufacturer’s instructions. The sensitivity limit of the assay was 1.04 U/mL.

2.4. Statistical Analysis

Hardy-Weinberg equilibrium was determined by a Chi-squared test. Haplotype frequencies were estimated based on expectation–maximization algorithm according to Slatkin and Excoffier [27]. Linkage disequilibrium (LD) was expressed as Lewontin’s coefficient (D′) using Arlequin software version 3.0. Haplotype frequency comparisons between RA patients and CS were performed using odds ratios (OR) with 95% confidence intervals (CI) by Fisher’s exact test. p values < 0.05 refer to a two-sided test which was considered statistically significant. For other comparisons we used Mann–Whitney U test, SPSS 18.0 software (SPSS Inc., Chicago, IL, USA).

3. Results

3.1. Clinical and Demographic Features

Clinical and demographic data of the RA patients and control subjects (CS) groups are shown in Table 1. The RA patients consisted predominantly of women of 48 (18–82) years with a long-term disease evolution (median ten years), while the median of year in CS was 36 years. The RA presented a condition with moderate clinical activity defined by the DAS28 index ( score) and some functional disability degree, according to Spanish HAQ-DI index ( score). Furthermore, 83% of total patients were positive for anti-CCP antibodies, and 84% were positive for RF. None of the control subjects are positive to anti-CCP antibodies, and only 5% are positive for RF.

3.2. Linkage Disequilibrium and Haplotypes of PTPN22 Gene

The genotypes and alleles distribution of the PTPN22 −1123G>C and +1858C>T polymorphisms were in Hardy-Weinberg equilibrium among control population (). Haplotype constructions of two SNPs revealed three of four possible haplotypes (with a frequency >1%); GC, CC, and CT haplotypes. The linkage disequilibrium analysis showed the strongest linkage between −1123G>C and +1858C>T PTPN22 polymorphisms (D′ = 1, ) calculated for CS. The PTPN22 susceptibility haplotype (−1123C/+1858T or CT), which included the polymorphic alleles from both SNPs, was more frequent in RA patients than in CS (OR 2.18, 95% CI 1.15–4.16, ) (Table 2).

Table 2: Association analyses of the PTPN22 haplotypes constructed for the −1123G>C and +1858C>T SNPs in RA patients and CS.

Haplotypes −1123G>C/+1858C>T
RA ()%
CS ()%
OR
(95% CI)
p

GC
66.26
63.33
NA
NA
NA

CC
28.18
34.44
0.78
(0.62–0.99)
0.05

CT
5.15
2.22
2.18
(1.15–4.16)
0.01

GT
0.40
0
6.69
(0.34–130.09)
0.25

Values represent a frequency of cases with respect to total number of haplotypes, for RA patients and control subject. SNP: single nucleotide polymorphism; RA: rheumatoid arthritis; CS: control subjects; OR: odds ratios; CI: confidence intervals; level of significance for Fisher exact test; NA: not applicable.

3.3. Association of PTPN22 Haplotypes with Anti-CCP Antibodies

Furthermore, we compared anti-CCP antibodies levels between RA patients homozygous for PTPN22 susceptibility haplotype CT and homozygous carriers of (−1123G/+1858C or GC) nonsusceptibility PTPN22 haplotype, which included the wild alleles from both SNPs. Our results did not show statistical differences among levels of anti-CCP antibodies among those patients (Figure 1).

Figure 1: Levels of anti-CCP autoantibodies according to homozygous haplotypes of the PTPN22 −1123G>C and +1858C>T SNPs. GC nonrisk haplotype and CT risk haplotype. The medians and interquartile ranges of the parameters are represented as a black line. Levels of anti-CCP are given in international units (U/mL) according to the Axis Shield Company.

Later, we selected anti-CCP positive RA patients and stratified them according to the antibodies levels. Patients who had anti-CCP antibodies levels higher than the upper limit of normal (ULN) laboratory assay, but less or equal to 3 times this value, were considered as low positive and patients with anti-CCP antibodies values that are higher than 3 times ULN were considered high positive (according to 2010 ACR/EULAR criteria) [25]. These groups were compared according to the presence of CT susceptibility or GC nonsusceptibility PTPN22 haplotypes (in individuals homozygous for both SNPs), to determine whether haplotypes influence the anti-CCP antibodies status. In our comparison, we did not find significant differences between carriers of susceptibility or nonsusceptibility haplotypes (Table 3). The relationship between anti-CCP levels and these haplotypes was also compared among RA patients with or without a history of smoking habit, but we did not find a significant correlation (data not shown).

Table 3: Frequencies of −1123G>C and +1858C>T PTPN22 haplotypes in RA seropositive patients stratified according to anti-CCP antibodies status.

GC% (n)
CT% (n)
OR
95% CI
p value

Anti-CCP low
3.9 (5)
0 (0)
—
—
—

Anti-CCP high
82.8 (106)
100 (3)
1.028
(0.996–1.061)
0.873

to antibodies values that are > than the upper limit of normal (ULN) and ≤ 3 times the ULN for laboratory assay; to antibodies values that are > than 3 times ULN for laboratory assay (according to ACR/EULAR 2010 criteria). OR: odds ratio; 95% CI: confidence interval; anti-CCP: anti-cyclic citrullinated peptide; Ab: antibody; RA: rheumatoid arthritis. The haplotypes include only homozygous individuals for each polymorphism.

4. Discussion

Haplotype frequencies and D′ values in our study indicate strong linkage disequilibrium between both PTPN22 SNPs (D′ = 1.00) which is in concordant with previous reports from Caucasian [16, 19, 21], Japanese [16], Norwegian [21], and Polish [20] populations. The extent of linkage disequilibrium found by these authors was very similar to ours (with D ′ values in CS: 0.99, 0.97, 0.99, and 0.97, resp.) [16, 19, 20].

The frequency distribution analysis revealed a significant increase in the susceptibility haplotype comprised by the risk alleles (−1123C and +1858T) of both SNPs, compared among RA patients and controls subjects. There is only one previous study in which the association between haplotypes formed by the PTPN22 −1123G>C and +1858C>T polymorphisms and RA gene was established [21]. In this report, the authors demonstrated an association between the CT haplotype consisting of the risk alleles of both polymorphisms and RA (OR 1.54, 95% CI 1.23–1.93, ). Our results agree with this work, as haplotype association between CT and RA susceptibility was established.

Altogether, these results confirm the direct involvement of PTPN22 gene independently of HLA locus in the pathogenesis of autoimmune diseases. However, considering the limited effect that has shown individually the −1123 G>C SNP in some reports, it is evident that the potential effects of the CT haplotype are mainly due to +1858 C>T PTPN22 SNP [28]. Thus, that is the reason why there was no evidence for an independent effect at +1858 C>T in some studies [19, 29, 30].

In established comparisons between the levels of anti-CCP antibodies between individuals carrying the susceptibility haplotype CT with those who had nonsusceptibility haplotype GC, we did not find any significant difference (). However, there is a marked tendency in patients with CT haplotype to have higher anti-CCP levels (median value of 240.5 U/mL) than those obtained for patients carrying GC haplotype (median value of 83.5 U/mL). These negative results can be due to a low frequency of the individuals homozygous to the risk haplotype CT in our population. Another important point to consider is the treatment and the evolution time of these patients, since they can modify the levels of anti-CCP antibodies; so it would be interesting to evaluate the outcome in patients with early RA, if possible, prior to treatment with disease modifying drugs (DMARDs).

Anti-CCP status was defined by 2010 ACR/EULAR criteria as negative, low, or high positive anti-CCP levels [25]. Later, we stratified seropositive RA patients group according to this classification and did not find statistical differences () between subjects carrying the CT risk haplotype compared to patients carrying GC haplotype, which indicates that CT risk haplotype did not predispose to high positive or any particular anti-CCP status; this observation remains constant even in smokers.

It is important to note that until now there are no studies that compare status or levels of anti-CCP antibodies with PTPN22 haplotypes in rheumatoid arthritis patients; in this respect, our study is the first of its type.

In conclusion, −1123G>C and +1858C>T PTPN22 SNPs are in strong linkage disequilibrium and (−1123C/+1858T) haplotype, which comprises polymorphic variants from both polymorphisms, is associated with susceptibility to RA but is not associated with status or higher levels of anti-CCP antibodies in a Western Mexican population. However, further studies are needed to elucidate the relationship between the risk haplotype CT and anti-CCP antibodies.

Competing Interests

The authors declare that they have no competing interests.

References

J. Ding, S. Eyre, and J. Worthington, “Genetics of RA susceptibility, what comes next?” RMD Open, vol. 1, no. 1, Article ID e000028, 2015.

S. Cohen, H. Dadi, E. Shaoul, N. Sharfe, and C. M. Roifman, “Cloning and characterization of a lymphoid-specific, inducible human protein tyrosine phosphatase, Lyp,” Blood, vol. 93, no. 6, pp. 2013–2024, 1999.

Y. Kochi, “Genetics of autoimmune diseases: perspectives from genome-wide association studies,” International Immunology, vol. 28, no. 4, pp. 155–161, 2016.

N. A. C. Tavares, M. M. S. Santos, R. Moura et al., “Association of TNF-α, CTLA4, and PTPN22 polymorphisms with type 1 diabetes and other autoimmune diseases in Brazil,” Genetics and Molecular Research, vol. 14, no. 4, pp. 18936–18944, 2015.

A. B. Begovich, V. E. H. Carlton, L. A. Honigberg et al., “A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis,” The American Journal of Human Genetics, vol. 75, no. 2, pp. 330–337, 2004.

N. Bottini, L. Musumeci, A. Alonso et al., “A functional variant of lymphoid tyrosine phosphatase is associated with type I diabetes,” Nature Genetics, vol. 36, no. 4, pp. 337–338, 2004.

C. Kyogoku, C. D. Langefeld, W. A. Ortmann et al., “Genetic association of the R620W polymorphism of protein tyrosine phosphatase PTPN22 with human SLE,” American Journal of Human Genetics, vol. 75, no. 3, pp. 504–507, 2004.

G. Orozco, D. Pascual-Salcedo, M. A. López-Nevot et al., “Auto-antibodies, HLA and PTPN22: susceptibility markers for rheumatoid arthritis,” Rheumatology, vol. 47, no. 2, pp. 138–141, 2008.

A. Skórka, T. Bednarczukt, E. Bar-Andziak, J. Nauman, and R. Ploski, “Lymphoid tyrosine phosphatase (PTPN22/LYP) variant and Graves' disease in a Polish population: association and gene dose-dependent correlation with age of onset,” Clinical Endocrinology, vol. 62, no. 6, pp. 679–682, 2005.

D. Smyth, J. D. Cooper, J. E. Collins et al., “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, vol. 53, no. 11, pp. 3020–3023, 2004.

J. A. Hamerman, J. Pottle, M. Ni, Y. He, Z.-Y. Zhang, and J. H. Buckner, “Negative regulation of TLR signaling in myeloid cells-implications for autoimmune diseases,” Immunological Reviews, vol. 269, no. 1, pp. 212–227, 2016.

G. Nabi, N. Akhter, M. Wahid et al., “Meta-analysis reveals PTPN22 1858C/T polymorphism confers susceptibility to rheumatoid arthritis in Caucasian but not in Asian population,” Autoimmunity, vol. 49, no. 3, pp. 197–210, 2016.

G.-P. Tang, L. Hu, and Q.-H. Zhang, “PTPN22 1858C/T polymorphism is associated with rheumatoid arthritis susceptibility in Caucasian population: a meta-analysis,” Journal of Zhejiang University. Medical sciences, vol. 43, no. 4, pp. 466–473, 2014.

N. M. Torres-Carrillo, Y. Ruiz-Noa, G. E. Martínez-Bonilla et al., “The +1858C/T PTPN22 gene polymorphism confers genetic susceptibility to rheumatoid arthritis in Mexican population from the Western Mexico,” Immunology Letters, vol. 147, no. 1-2, pp. 41–46, 2012.

J. Zhang, N. Zahir, Q. Jiang et al., “The autoimmune disease-associated PTPN22 variant promotes calpain-mediated Lyp/Pep degradation associated with lymphocyte and dendritic cell hyperresponsiveness,” Nature Genetics, vol. 43, no. 9, pp. 902–907, 2011.

E. Kawasaki, T. Awata, H. Ikegami et al., “Systematic search for single nucleotide polymorphisms in a lymphoid tyrosine phosphatase gene (PTPN22): association between a promoter polymorphism and type 1 diabetes in Asian populations,” American Journal of Medical Genetics, vol. 140, no. 6, pp. 586–593, 2006.

J. Huang, Y. Qiu, H. Li, D. Sun, J. Yang, and C. Yang, “A PTPN22 promoter polymorphism −1123G>C is associated with RA pathogenesis in Chinese,” Rheumatology International, vol. 32, no. 3, pp. 767–771, 2012.

X. Feng, Y.-Z. Li, Y. Zhang et al., “Association of the PTPN22 gene (−1123G > C) polymorphism with rheumatoid arthritis in Chinese patients,” Tissue Antigens, vol. 76, no. 4, pp. 297–300, 2010.

O. Cinek, O. Hradsky, G. Ahmedov et al., “No independent role of the -1123 G > C and + 2740 A > G variants in the association of PTPN22 with type 1 diabetes and juvenile idiopathic arthritis in two Caucasian populations,” Diabetes Research and Clinical Practice, vol. 76, no. 2, pp. 297–303, 2007.

M. Fichna, M. Zurawek, D. Januszkiewicz-Lewandowska, P. Fichna, and J. Nowak, “PTPN22, PDCD1 and CYP27B1 polymorphisms and susceptibility to type 1 diabetes in Polish patients,” International Journal of Immunogenetics, vol. 37, no. 5, pp. 367–372, 2010.

M. K. Viken, M. Olsson, S. T. Flåm et al., “The PTPN22 promoter polymorphism -1123G>C association cannot be distinguished from the 1858C>T association in a Norwegian rheumatoid arthritis material,” Tissue Antigens, vol. 70, no. 3, pp. 190–197, 2007.

L. I. Sakkas, D. P. Bogdanos, C. Katsiari, and C. D. Platsoucas, “Anti-citrullinated peptides as autoantigens in rheumatoid arthritis-relevance to treatment,” Autoimmunity Reviews, vol. 13, no. 11, pp. 1114–1120, 2014.

F. van Gaalen, A. Ioan-Facsinay, T. W. J. Huizinga, and R. E. M. Toes, “The devil in the details: the emerging role of anticitrulline autoimmunity in rheumatoid arthritis,” The Journal of Immunology, vol. 175, no. 9, pp. 5575–5580, 2005.

A. H. M. van der Helm-van Mil, J. Z. Wesoly, and T. W. J. Huizinga, “Understanding the genetic contribution to rheumatoid arthritis,” Current Opinion in Rheumatology, vol. 17, no. 3, pp. 299–304, 2005.

D. Aletaha, T. Neogi, A. J. Silman et al., “Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative,” Arthritis Rheum, vol. 62, pp. 2569–2581, 2010.

J. R. Padilla-Gutiérrez, Y. Valle, M. Vázquez-Del Mercado, M. Maldonado, and J. F. Muñoz-Valle, “A new PCR-RFLP assay for -1123 G>C polymorphism in the PTPN22 gene: allele and genotype frequencies in a western Mexican population,” Clinical Chemistry and Laboratory Medicine, vol. 47, no. 4, pp. 491–493, 2009.

M. Slatkin and L. Excoffier, “Testing for linkage disequilibrium in genotypic data using the Expectation-Maximization algorithm,” Heredity, vol. 76, no. 4, pp. 377–383, 1996.

J.-E. Martín, B. Z. Alizadeh, M. A. González-Gay et al., “Evidence for PTPN22 R620W polymorphism as the sole common risk variant for rheumatoid arthritis in the 1p13.2 region,” Journal of Rheumatology, vol. 38, no. 11, pp. 2290–2296, 2011.

W. R. Wan Taib, D. J. Smyth, M. E. Merriman et al., “The PTPN22 locus and rheumatoid arthritis: no evidence for an effect on risk independent of Arg620Trp,” PLoS ONE, vol. 5, no. 10, Article ID e13544, 2010.

M. Zoledziewska, C. Perra, V. Orrù et al., “Further evidence of a primary, causal association of the PTPN22 620W variant with type 1 diabetes,” Diabetes, vol. 57, no. 1, pp. 229–234, 2008.