Abstract
Background:
A missence single-nucleotide polymorphism (SNP) in the protein tyrosine phosphatase nonreceptor 22 (PTPN22) gene known as R620W (rs2476601) was recently reported to be associated with several autoimmune diseases including Graves' disease (GD). The association was repeatedly confirmed in the populations of North European ancestry. However, this amino acid was reported to be nonpolymorphic in the Asian populations. Since the gene confers an impact on autoimmune diseases, we attempt to explore an association between the PTPN22 gene and autoimmune thyroid disease (AITD) in a Japanese population without restricting to rs2476601. Previous investigations have also demonstrated that two intronic SNPs (rs706778 and rs3118470) in the interleukin-2 receptor-α (IL2RA) gene were associated with type 1 diabetes in the Japanese population.
Patients and Methods:
We genotyped the five SNPs (rs12760457, rs2797415, rs1310182, rs2476599, and rs3789604) of the PTPN22 and the two SNPs (rs706778 and rs3118470 in the IL2RA gene) in 456 Japanese patients with AITD (286 with GD, 170 with Hashimoto's thyroiditis) and 221 matched Japanese control subjects. Seven SNPs were analyzed by either the SNAPshot method or the high-resolution melting and unlabeled probe methods. Case-control association studies were performed using the χ 2 and Fisher's exact tests with Yates correction. Haplotype was conducted using the expectation-maximization algorithm.
Results:
No association was found between any of the individual SNPs of the PTPN22 gene and AITD. Permutation analysis revealed that the distribution of one haplotype is significantly different between patients with AITD and controls (p = 0.0036). A novel protective effect of a haplotype containing five SNPs was observed (p < 0.0001 for AITD, p < 0.0001 for GD, and p < 0.0001 for Hashimoto's thyroiditis, respectively). The GG allele of rs3118470 in the IL2RA gene was significantly associated with GD (p = 0.03), although the association was weak.
Conclusions:
Significant difference in the distribution of the haplotype suggests that the PTPN22 gene rather than rs2476601 is involved in the development of AITD in the Japanese population.
Introduction
The pathogenesis of AITDs is thought to involve several risk factors, including genetic risk factors [reviewed in Ref. (6)] and environmental triggers such as cigarette smoking, iodine intake, and infection (7,8) [reviewed in Ref. (3)]. However, the evidence for interactions between hereditary factors and environmental influences appears to be much stronger for cigarette smoking and iodine intake than for infections (8).
The first locus shown to be associated with AITDs was the HLA-DR (HLA-DRB1) locus [reviewed in Ref. (6)]. The HLA-DR3 (DRB1*03) has been consistently shown to be associated with GD in Caucasians, with an odds ratio of 2.0–3.0 (6). Other HLA alleles have been shown to be associated with GD in nonCaucasian populations [reviewed in Ref. (9)]. Non-HLA genes have also been shown to influence the expression of GD. These genes include the genes for cytotoxic T lymphocyte antigen-4 (CTLA-4), CD40, thyroglobulin (Tg), and TSH receptor [reviewed in Ref. (10)].
Lymphoid tyrosine phosphatase (LYP) is a protein tyrosine phosphatase (PTP) that plays a negative role in T cell signaling by dephosphorylating the Src family of kinases, thus preventing initiation of T cell activation (11). Recent studies have identified association of T allele (W620) of a functional missense single-nucleotide polymorphism (SNP), +1858C>T encoding an arginine to tryptophan (R620W) substitution at codon 620 (dbSNP accession number rs2476601) of PTP nonreceptor 22 (PTPN22) encoding LYP, with many autoimmune diseases, including type 1 diabetes (T1D) (12 –14), systemic lupus erythematosus (SLE) (15,16), GD (13,17,18), rheumatoid arthritis (RA) (19 –21), HT (16), juvenile idiopathic arthritis (21), and vitiligo (22) (Table 1).
SNP, single-nucleotide polymorphism; OR, odds ratio; T1D, type 1 diabetes; GD, Graves'disease; HT, Hashimoto's thyroiditis; SLE, systemic lupus erythematosus; RA, rheumatoid arthritis; JIA, juvenile idiopathic arthritis.
However, rs2476601 was reported to be nonpolymorphic in the Asian populations (23,24). Recently, Kawasaki et al. reported that heterozygousity (molecular heterosis) of the polymorphism in the promoter region of PTPN22 gene (dbSNP accession number rs2488457) is associated with Japanese acute onset T1D (25). In Caucasians, rs2488457 is strongly linked to the rs2476601; and several studies have disclosed that the association of rs2488457 and autoimmune diseases is a result of the linkage to the rs2476601, and real susceptible allele is 1858T (620W) (26,27).
Carlton et al. (28) identified 37 SNPs in the PTPN22 gene and reported that several SNPs including rs1310182 and rs3789604 are also associated with RA independent of rs2476601 in the Caucasian population. They also studied the haplotype and found that some haplotypes which do not contain R620W are associated with RA. However, these SNPs were not associated with Japanese RA (24). If the functional change of LYP is important for the process of autoimmune diseases, there might be disease susceptible alleles other than 620W in the Japanese population.
Another molecule with a role in the control of immune responses and the maintenance of immune homeostasis that has been demonstrated to have a role in autoimmunity is the interleukin-2 receptor-α (IL2RA). The IL2RA encodes the α-chain of the IL2R complex (also known as CD25), which is central to immune regulation as an important modulator of self-tolerance and immunity (29), and the IL2RA association with T1D was originally identified by Vella and coworkers (30). The IL2RA gene has also been associated with GD (31), RA (32), and multiple sclerosis (33) in Caucasians, which implies that this locus may have a general effect on predisposition to autoimmunity. Recently, two SNPs in intron 1 of the IL2RA, rs706778 and rs3118470, were associated with T1D in the Japanese population (34).
In this study, we examined whether the SNPs associated with Caucasian RA are also associated with Japanese AITD. We also analyzed the haplotype and found that the PTPN22 gene plays some role in the susceptibility or the protection of AITD in the Japanese population.
Materials and Methods
Subjects
This project was approved by the Institutional Review Board of Showa University. Two hundred eighty-six Japanese patients with GD, 170 patients with HT, and 221 Japanese controls were studied. GD was diagnosed based on clinical symptoms and biochemical confirmation of hyperthyroidism, including diffuse goiter, elevated radioactive iodine uptake, the presence of antithyrotropin receptor antibodies, (TRAbs) and elevated thyroid hormone levels. The patients with HT had documented clinical and biochemical hypothyroidism requiring thyroid hormone replacement therapy and showed autoantibodies against thyroid peroxidase (TPOAb) with or without antibodies against Tg (TgAb). Two hundred twenty-one age-matched Japanese volunteers served as controls in our association studies. All controls had no personal or family history of any autoimmune disease or thyroid condition.
SNP genotyping
DNA was extracted from whole blood using the Puregene kit (Gentra Systems, Minneapolis, MN). We analyzed the following five SNPs in the PTPN22 gene: rs12760457 in intron 11; rs2797415 in intron 15; rs1310182 in intron 16; rs2476599 in intron 19; and rs3789604 in putative transcription factor-binding sites downstream of the PTPN22 gene. We also analyzed the following two SNPs in the IL2RA gene: rs706778 and rs3118470 in intron 1of the IL2RA gene. The SNPs rs1310182 and rs3789604 were analyzed by the SNAPshot (PE Applied Biosystems, Foster city, CA) based on the manufacturer's protocol. The SNPs rs12760457, rs2797415, rs2476599, rs706778, and rs3118470 were genotyped by the high-resolution melting and unlabeled probe methods using LightScanner® (Idaho Technology, Salt Lake City, Utah) based on the manufacturer's protocol. Details of the primers are given in Table 2.
PTPN22, protein tyrosine phosphatase nonreceptor 22; IL2RA, interleukin-2 receptor-α.
Laboratory tests
The serum concentrations of free triiodothyronine, free thyroxine, and TSH were determined by enzyme immunoassay. The TRAb was measured by a radioreceptor assay using a commercially available kit (TRAb Cosmic III; Cosmic, Tokyo, Japan). The cutoff value for TRAb was 10%. The TPOAb and TgAb were measured by a radioimmunoassay using a commercially available kit (TPOAb Cosmic II and TgAb Cosmic II; Cosmic). The cutoff value for TPOAb and TgAb was 0.3 U/mL.
Statistical analysis
Case-control analysis and Hardy–Weinberg equilibrium (HWE) test of SNP were performed using SNPAlyze ver. 7.0 (Dynacom, Yokohama, Japan) (35). The HWE tests were carried out for all loci among subjects and controls separately. Tests in subjects and controls did not show any significant deviation from HWE for any of the SNPs. Linkage disequilibrium (LD) between SNPs was evaluated by the D′ and r 2 of pair-wise LD using SNPAlyze ver. 7.0 (Dynacom). Haplotype frequencies for multiple loci were estimated by phase estimation using the expectation-maximization algorithm. Permutation p-values were calculated by comparing haplotype frequencies between cases and controls on the basis of 10,000 replications using SNPAlyze ver. 7.0.
Results
PTPN22 variants and AITD
The frequencies of the alleles and genotypes of the five PTPN22 SNPs in our Japanese population are shown in Table 3. The individual SNPs were not associated with AITD, GD, or HT in our Japanese population (Table 3). All cases and controls were in HWE.
Values given are the number of subjects, with the percentage in parentheses.
SNP numbers corresponding to those reported by Carlton et al.28
AITD, autoimmune thyroid disease.
PTPN22 haplotype analysis
Due to the strong LD between five variants, haplotype analysis was undertaken using the computer program SNPAlyze version 7.0. Five haplotypes were identified, three of which (haplotypes 1, 2, and 4) were correlated with haplotypes 1, 4, and 5 identified in the report by Carlton et al. (Table 4). Four haplotypes (haplotypes 1–4) were relatively common and 1 haplotype was rare. Distribution of the haplotype was significantly different between AITD and control by permutation procedure (p = 0.0036). A novel protective effect of a haplotype containing five SNPs was observed (p < 0.0001 for AITD, p < 0.0001 for GD, and p < 0.0001 for HT, respectively) (Table 4).
The program, SNPAlyze ver. 7.0 Standard (Dynacom, Yokohama, Japan), was used to estimate common (frequencies >0.01) haplotypes for the five SNPs genotyped.
Each haplotype was compared with the other haplotypes combined.
IL2RA variants and AITD
The frequencies of the alleles and genotypes of the two IL2RA SNPs in our Japanese population are shown in Table 5. As shown in Table 5, the GG allele was significantly associated with GD (p = 0.03), although the association was weak. All cases and controls were in HWE.
Values given are the number of subjects, with the percentage in parentheses.
Discussion
The LYP is encoded by PTPN22 on chromosome 1p13 (36) and exerts its effect on T cell activation via association with a variety of adaptor molecules, including c-src tyrosine kinase (Csk) kinase (37), c-Cbl (36), and Grb2 (38). It is known that the presence of tryptophan at position 620 (W620) can disrupt the ability of LYP to bind to at least one of these adaptor molecules, such as Csk kinase, due to the inability of W620 to fit into the binding pocket on the Csk molecule (12). This has been shown experimentally in both Escherichia coli and COS cells, where only the construct containing R620 was precipitated by Csk (12). The combination of these genetic and functional data suggests that LYP plays an important role in regulating the autoimmune process.
The SNP rs2476601 of the PTPN22 gene was reported to be nonpolymorphic in the Asian populations (23,24). The T allele of rs2476601 of the PTPN22 gene has been shown to increase susceptibility to AITD in Caucasians (16 –18). Further, a metaanalysis showed associations of the PTPN22 1858T allele with RA, SLE, GD, and T1D in Caucasians (39). It has been proposed that the disease-associated T allele encodes a protein that does not bind to the protein tyrosine kinase Csk and may, therefore, cause general hyperresponsiveness of T cells (12). Among Japanese populations, however, the PTPN22 T variant has not been found in patients with GD, patients with T1D, patients with RA, or control subjects (23 –25,40), confirming that the PTPN22 gene polymorphism shows ethnic differences.
Recently, rs2488457 in the promoter region was reported to be associated with acute onset T1D in a Japanese population (25). However, there was no association of rs2488457 with GD (40). Further, rs3789604 of the PTPN22 gene was found to be associated with RA, independently of rs2476601 (22). The SNP rs3789604 lies 1496 bases downstream of PTPN22 at the 50 end of the round spermatid basic protein 1 gene, where it encodes either a silent mutation or putative transcription factor-binding sites, depending on the transcript. Recently, the AA-genotype and A-allele frequencies of rs3789604 were significantly higher in patients with GD than in control subjects (40), suggesting that the SNP37 polymorphism or a gene with LD may be relevant to susceptibility to GD in Japanese populations. Therefore, we further analyzed five other SNPs including rs12760457, rs2797415, rs1310182, rs2476599, and rs3789604, to clarify whether a susceptibility locus for AITD exists at another location within the PTPN22 gene. However, there was no association of any of the individual SNPs with disease.
We also analyzed the haplotype and found that the PTPN22 gene plays some role in the susceptibility or the protection of AITD in the Japanese population. There was significant LD between the SNPs tested. Five haplotypes were present in the PTPN22 gene, and a novel protective effect of a haplotype containing five SNPs (haplotype 5) was observed. Combined with the results arising from previous studies, we reveal that the association between the haplotype of the PTPN22 gene and autoimmune diseases is complicated. Identification of the precise risk or protective alleles in the region of the PTPN22 gene linked to the associated haplotypes requires further investigation. Recently, a loss-of-function mutation in the PTPN22 gene that does not exist in the Asian population was shown to confer protection against human SLE (41). It has also been proposed that significant LD may extend for several hundred kilobases surrounding PTPN22 and, thus, potentially include additional genes including round spermatid basic protein 1 and DNA cross-link repair 1B (DCLRE1B) (42). Therefore, the possibility that a gene other than PTPN22 is responsible for the association with AITD cannot be fully excluded.
In the present study, we analyzed the two SNPs in the IL2RA gene, and the GG allele of the rs3118470 was significantly associated with GD (p = 0.03), although the association was weak.
In summary, our data suggest that genetic polymorphisms within PTPN22 are influencing susceptibility in a disease-specific manner, a finding that warrants further functional analysis to elucidate their true role in AITD and the autoimmune disease process, in general. We find an evidence to support a weak effect at the IL2RA locus on susceptibility to GD, perhaps due to ethnic differences.
Footnotes
Acknowledgments
This work was supported in part by a Showa University Grant-in-aid for Innovative Collaborative Research Projects (to Yoshiyuki Ban [Y.B.]), the Showa University Medical Foundation (to Y.B.), a grant from Showa University School of Medicine Alumni Association (to Y.B.), and a grant from Yamaguchi Endocrine Research Association (to Y.B.).
Disclosure Statement
The authors declare that they have no competing financial interests.
