International Journal of Genomics

Volume 2014 (2014), Article ID 247637, 8 pages

http://dx.doi.org/10.1155/2014/247637

Research Article

A Promoter Region Polymorphism in PDCD-1 Gene Is Associated with Risk of Rheumatoid Arthritis in the Han Chinese Population of Southeastern China

CuiPing Liu,¹ JueAn Jiang,¹ Li Gao,2,3 XiaoHan Hu,² FengMing Wang,⁴ Yu Shen,¹ GeHua Yu,² ZuoTao Zhao,⁵ and XueGuang Zhang¹

¹Jiangsu Institute of Clinical Immunology, The First Affiliated Hospital of Soochow University, Suzhou 215006, China

²Medical Biotechnology Institute, Soochow University, Suzhou 215006, China

³Wuxi Red Cross Blood Center, Wuxi 214021, China

⁴Testing Center, Center for Disease Prevention and Control, Changzhou 213000, China

⁵Department of Dermatology, First Hospital, Peking University, Beijing 100034, China

Received 21 January 2014; Revised 10 March 2014; Accepted 11 March 2014; Published 3 April 2014

Academic Editor: Ji-Fu Wei

Copyright © 2014 CuiPing Liu 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.

Objective. Programmed cell death 1 (PD-1) induces negative signals to T cells during interaction with its ligands and is therefore a candidate gene in the development of autoimmune diseases such as rheumatoid arthritis (RA). Herein, we investigate the association of PDCD-1 polymorphisms with the risk of RA among Chinese patients and healthy controls. Methods. Using the PCR-direct sequencing analysis, 4 PDCD-1 SNPs (rs36084323, rs11568821, rs2227982, and rs2227981) were genotyped in 320 RA patients and 309 matched healthy controls. Expression of PD-1 was determined in peripheral blood lymphocytes by flow cytometry and quantitative real-time reverse transcriptase polymerase chain reaction. Results. We observed that the GG genotype of rs36084323 was associated with a increased risk for developing RA (OR 1.70, 95% 1.11–2.61, ). Patients carrying G/G genotype displayed an increased mRNA level of PD-1 compared with A/A genotype and healthy controls. Meanwhile, patients homozygous for rs36084323 had induced basal PD-1 expression on activated CD4+ T cells. Conclusion. The PDCD-1 polymorphism rs36084323 was significantly associated with RA risk in Han Chinese population. This SNP, which effectively influenced the expression of PD-1, may be a biomarker of early diagnosis of RA and a suitable indicator of utilizing PD-1 inhibitor for treatment of RA.

1. Introduction

Rheumatoid arthritis (RA) is a common chronic inflammatory autoimmune disease characterized by significant disability and early mortality, which affects ~1% of adult population worldwide [1, 2]. It is generally accepted that RA is a complex autoimmune disorder, characterized by a chronic T-cell response that evaded normal control mechanisms [3, 4]. Therefore, the genes involved in the regulation of T-cell responses may be primary determinants of susceptibility to RA.

Programmed death-1 (PD-1, also called CD279) is a novel costimulatory member of B7/CD28 family, which is inducibly expressed on CD4+ T cells, CD8+ T cells, natural killer T cells, B cells, and activated monocytes [5]. PD-1 receptor has two ligands: PD-L1 (also known as B7-H1 or CD274) and PD-L2 (also called B7-DC or CD273). PD-L1 is expressed on T cells, B cells, dendritic cells (DC), macrophages, and some tumor cells and is further upregulated upon activation. PD-1 engagement by PD-L1 dephosphorylates proximal signaling molecules and augments PTEN expression, inhibiting PI3K and AKT activation [6, 7]. The critical role of PD-1 in immune regulation is highlighted by gene disruption studies demonstrating strain-specific autoimmune phenotypes [8, 9]. In addition, genetic studies revealed that there is an association between PDCD-1 gene polymorphism and susceptibility to autoimmune diseases, such as systemic lupus erythematosus (SLE) [10, 11], rheumatoid arthritis [12, 13], multiple sclerosis [14], and diabetes mellitus [15, 16].

There is mounting evidence that PD-1 is linked to human autoimmunity. In view of the pivotal role of PD-1/PD-L pathway in autoimmn immunology, it is worth considering of PDCD-1 functional SNPs, PDCD-1 −606A/G (rs36084323, also known as PD1.1) in promoter, PDCD-1 +7146A/G (rs11568821, also known as PD1.3) in intron 4, and PDCD-1 +7625G/A (rs2227982, also known as PD1.9) and PDCD-1 +7786G/C (rs2227981, also known as PD1.5) in exon 5 as prospective candidates for individuals susceptibility to RA in the Han Chinese population. Therefore, in the present study, we sought to determine the association between functional genetic variations in the PDCD-1 gene and RA risk in a Chinese population in mainland.

2. Materials and Methods

2.1. Patients and Controls

This study was approved by the Ethics Committee of Soochow University and all subjects gave informed consent for the genetic analyses. A total of 320 unrelated Chinese RA patients were recruited from the Outpatient Departments of Rheumatology in the First and the Third Affiliated Hospital of Soochow University. They were composed of 72 men and 248 women, whose mean age was 55.3 years (SD = ±12.6). Individual patients with RA were diagnosed according to the diagnosis criteria established by the American College of Rheumatology and the disease severity of individual patients was evaluated using the disease activity score 28 (DAS28) [1]. A total of 20 patients with new-onset RA (<6 months of disease duration) were recruited for expression of PD-1 protein on activated T cells. Individual RA patients were excluded if she/he received treatment with DMARDs, corticosteroids, or immunosuppressive for any reason during the past 6 months or had other chronic inflammatory and autoimmune diseases, such as diabetes, multiple sclerosis, inflammatory bowel disease, metabolic syndrome, hypertension, cardiovascular diseases, cancer, or recent infection. The controls were gender, age, and ethnically matched unrelated healthy people obtained from the checkup population in the above two hospitals (Table 1).

Table 1: Characteristics of RA patients and controls.

Characteristic
RA patients
Controls
value

Total number
320
309

Female
248 (77.5)
240 (77.7)
>0.05a

Male
72 (22.5)
69 (22.3)

Age, mean SD years

>0.05b

Duration of disease, years

RF+c
261 (81.6)

RF−c
59 (19.4)

Anti-CCP+c
244 (76.4)

Anti-CCP−c
76 (23.6)

Values are numbers (%); RA: rheumatoid arthritis; RF: rheumatoid factor; anti-CCP: anticyclic citrullinated peptide.

value calculated by Pearson chi-square test (all frequency >0.05) or Fisher’s exact test (any frequency <0.05).

b value calculated by Student’s -test.

cClinical data were not available for some cases.

2.2. DNA Extraction and Polymorphism Genotyping

Peripheral venous blood samples of 2 mL were drawn from each individual by standard venepuncture and stored at −20°C. Genomic DNA was isolated from peripheral blood leucocytes by standard procedures. The reference sequence is the human PDCD-1 sequence (GeneBank accession number AF363458). The PDCD-1 −606A/G (rs36084323) in promoter, PDCD-1 +7146A/G (rs11568821) in intron 4, and PDCD-1 +7625G/A (rs2227982) and PDCD-1 +7786G/C (rs2227981) in exon 5 polymorphisms were determined by direct sequencing in an Applied Biosystems sequencer (ABI Prism, Model 3100, Avant).

2.3. RNA Isolation and Quantitative Real-Time Reverse Transcriptase Polymerase Chain Reaction

Total RNA from PBMC of 30 RA patients and 24 healthy controls were extracted using Trizol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s protocol. Total RNA (1 μg) was used for cDNA synthesis with oligodT primers (Invitrogen, Karlsruhe, Germany) and superscript II reverse transcriptase (Takara Bio, Shiga, Japan). PCR was performed using PCR Master (Roche, Mannheim, Germany) with the following primers: for PD-1 mRNA, 5′-CTCAGGGTGACAGAGAGAAG-3′ (forward) and 5′-GACACCAACCACCAGGGTTT-3′ (reverse) and for GAPDH, 5′-GTGAAGGTCGGAGTCAACG-3′ (forward) and 5′-TGAGGTCAATGAAGGGG-TC-3′ (reverse). Fold changes were normalised based on GAPDH expression, and each assay was conducted in a 96-well ABI 7900HT real-time PCR system (Applied Biosystems, Foster City, CA, USA). This procedure was performed in triplicate.

2.4. Antibodies and Flow Cytometry

All antibodies were mouse anti-human monoclonal antibodies (mAb). PE-Cy5-conjugated anti-CD4 (OKT4) were from Beckman Coulter (Miami, FL). PE-conjugated anti-PD-1 (MIH4), fluorescein isothiocyanate (FITC)-conjugated anti-CD25 (M-A251), anti-CD69 (FN50), and anti-HLA-DR (G46-6) were from BD PharMingen (Franklin Lakes, NJ). PE- or PE-Cy5-conjugated IgG1 (679.1Mc7) (Beckman Coulter), PE-conjugated IgG2a (eBM2a; eBio-science), and FITC-conjugated IgG1 (P3) (eBioscience) were used as IgG isotype controls. PBMCs (0.5 × 10⁶) were incubated in wash buffer (PBS/2.5% fetal bovine serum [FBS]) with appropriate amounts of mAb on ice for 30 minutes. Cells were washed and were immediately analyzed on an EPICS XL-MCL flow cytometer (Beckman Coulter).

2.5. Statistics Analysis

The deviation from Hardy-Weinberg equilibrium (HWE) was examined in controls by the   test. The following statistical analyses were performed using SNPstats software (Availability: http://bioinfo.iconcologia.net/SNPstats) [17]. Based on the logistic regression method, the case-control association of genotypes was tested and the odds ratios (OR) and 95% confidence intervals (95% CI) were given. and were calculated to evaluate the magnitude of LD. Haplotype frequencies were estimated using the EM algorithm coded into the haplo.stats package (http://mayoresearch.mayo.edu/mayo/research/biostat/schaid.cfm). The association analysis of haplotypes was similar to that of genotypes with logistic regression, and results were shown as OR and 95% CI. The most frequent haplotype was automatically selected as the reference category and rare haplotypes were pooled together in a group. The log-additive inheritance model was assumed by default. The significance level of all these tests was 0.05. The nonparametric Mann-Whitney test was used for comparisons among groups with small or unequal sample sizes. Results were expressed as the mean ± SEM, and 2-tailed values less than 0.05 were considered significant.

3. Results

3.1. Single Nucleotide Polymorphism Analysis

A total of four SNPs were successfully genotyped in 320 RA patients and 309 healthy controls. Table 2 shows the allele and genotype distribution of these four SNPs. values for Hardy-Weinberg proportions of the SNPs are shown in Table 2 as well. For all four SNPs, the genotypic distribution in controls conformed to HWE. Among the four SNPs, the genotype and allele distributions of rs36084323 differed significantly between RA patients and controls (). When logistic regression was used for association analysis after modeling the SNPs’ effects as additive, dominant, or recessive, rs36084323 showed significant difference in codominant (OR, 1.70; 95% CI, 1.11–2.61), recessive (OR, 1.50; 95% CI, 1.05–2.14), and log-additive (OR, 1.30; 95% CI, 1.05–1.61) models (Table 3). The log-additive model was accepted as the best inheritance model because it showed the smallest Akaike information criterion value (869.4). The rs11568821 nucleotide is not polymorphic among Chinese population. The other two SNPs showed no association with RA in all 5 inheritance models (data not shown).

Table 2: The 4 SNPs of PDCD-1 gene investigated in the cases () and controls ().

SNP
dbSNP ID and position

Frequency (%)
HWE

Allele

Genotype

1
rs36084323

A
G

A/A
A/G
G/G

promoter (−606A/G)
Cases
47.2
52.8

24.7
45.1
30.3

Controls
54.1
45.9
0.014*
31.1
46.5
22.5
0.049*
0.12

2
rs11568821

A
G

A/A
A/G
G/G

Intron 4 (+7146A/G)
Cases
—
100

—
—
100

Controls
—
100

—
—
100
—
—

3
rs2227982

C
T

C/C
C/T
T/T

Eoxn 5 (+7625G/A)
Cases
51.5
48.5

24.8
53.5
21.8

Controls
52.5
47.5
0.74
29.6
45.8
24.6
0.19
0.21

4
rs2227981

C
T

C/C
C/T
T/T

Exon 5 (+7786G/C)
Cases
73.2
26.8

54.8
36.8
8.4

Controls
68.0
32.0
0.06
47.0
42.2
10.9
0.18
0.65

HWE indicates Hardy-Weinberg equilibrium; .

Table 3: The association analysis of rs36084323 with RA (adjusted by sex) by using logistic regression.

Model
Genotype
Cases (%)
Controls (%)
OR (95% CI)
value
AIC
BIC

Codominant
A/A
24.7
31.1
1
0.049*
871.3
884.6

A/G
45.1
46.5
1.22 (0.84–1.78)

G/G
30.3
22.5
1.70 (1.11–2.61)

Dominant
A/A
24.7
31.1
1
0.073
872.1
881

A/G-G/G
75.3
68.9
1.38 (0.97–1.96)

Recessive
A/A-A/G
69.7
77.5
1
0.026*
870.3
879.2

G/G
30.3
22.5
1.50 (1.05–2.14)

Overdominant
A/A-G/G
54.9
53.5
1
0.73
875.2
884

A/G
45.1
46.5
0.95 (0.69–1.29)

Log-additive
—
—
—
1.30 (1.05–1.61)
0.015*
869.4
878.3

OR: odds ratio; 95% CI: 95% confidence intervals; AIC: Akaike’s information criterion; BIC: Bayesian information criterion; .

3.2. Linkage Disequilibrium and Haplotypes Association Analysis

Because rs11568821 was nonpolymorphic in our population, it was excluded from the haplotype construction and LD analysis. Pairwise LD between the SNPs of the PDCD-1 (rs36084323, rs2227982, and rs2227981) was calculated for the cases and controls in the Han Chinese. We found strong LD ( > 0.75) between some pairs of markers in the PDCD-1 gene including rs36084323/rs2227982 ( = 0.7729) and rs2227982/rs2227981 ( = 0.7841).

The association analysis of the haplotypes with RA was similar to that of genotypes by logistic regression (Table 4). We found that the haplotype of rs36084323/rs2227982/rs2227981 showed significant association with the disease (); ACT and ACC haplotypes were less frequently observed in cases than in controls (OR: 0.23, 95% CI: 0.10–0.53; OR: 0.32, 95% CI: 0.11–0.96, resp.), indicating that these two haplotypes act as protective phenotype in RA. The association analysis of the haplotypes was adjusted by sex too.

Table 4: Estimated haplotype frequencies and the association analysis with RA.

Haplotype
Sequence
Cases
Controls
Total
OR (95% CI)
value

1
ATC
0.4108
0.4061
0.4081
1.00
—

2
GCC
0.2583
0.1919
0.2248
1.28 (0.94–1.73)
0.11

3
GCT
0.2127
0.2066
0.2098
0.98 (0.72–1.32)
0.87

4
GTC
0.0523
0.0462
0.0497
1.04 (0.53–2.02)
0.91

5
ACT
0.019
0.0821
0.0463
0.23 (0.10–0.53)

6
ACC
0.0138
0.0420
0.0280
0.32 (0.11–0.96)
0.043*

7
ATT
0.0269
0.0139
0.0210
1.61 (0.52–5.02)
0.41

Rare
—
—
—
0.0081
1.55 (0.27–8.96)
0.68

Global haplotype association P value: 0.0015*

OR: odds ratio; 95% CI: 95% confidence intervals..

Rare haplotypes were pooled together in a group.

3.3. Association between rs36084323 Polymorphism and PD-1 mRNA Expression

To evaluate the association between rs36084323 polymorphism and RA, we examined whether or not the rs36084323 polymorphism was associated with an altered PD-1 mRNA expression. Significant difference was observed in the relative PD-1 mRNA expression level between patients with GG and AA genotypes (Figures 1(a) and 1(b)). We concluded that PDCD-1 polymorphisms were associated with RA, and the GG genotype and G allele of rs36084323 that associated with increased PD-1 mRNA expression might be involved in RA development in Han Chinese population.

(a)

(b)

Figure 1: Analysis of PD-1 mRNA expression levels in RA patients with different genotypes of rs36084323. (a) Increased levels of PD-1 mRNA levels in PBMC of RA () patients, as compared with those from healthy controls (HC) () (). (b) The patients with GG genotype exhibited higher PD-1mRNA expression levels than those with AA genotype (). Horizontal bars indicate the mean ± SD, * = .

3.4. Upregulated Expression of PD-1 on Activated CD4+T Cells in RA Patients Was Associated with rs36084323

We next examined the expression of PD-1 on PBMCs in RA patients. Flow cytometry analyses of PBMC samples from 20 RA patients demonstrated the upregulated PD-1 expression on CD4+CD25+T cells ( versus , ), CD4+CD69+T cells (22.95 ± 2.68 versus 15.10 ± 1.432, ), and CD4+HLA-DR +T cells ( versus , ) compared with healthy controls (Figures 2(a) and 2(b)). Although we found no difference between wild-type (G/G) and heterozygous (A/G) patients, the 4 patients who were homozygous for the SNP rs36084323 (A/A) had minimal PD-1 expression on freshly isolated CD4+T cells, including the activated CD25+ and CD69+ T-cell subsets (Figure 2(c)). Collectively, these results suggest that the rs36084323 polymorphism confers increased basal PD-1 expression at early-to-intermediate stages of CD4+ T-cell activation and that this is associated with increased risk of RA.

(a)

(b)

(c)

Figure 2: Increased basal programmed death 1 (PD-1) expression on activated CD4+ T cells in RA patients. (a) A representative flow cytometry analysis of PD-1 expression on CD4+ T cells in RA patients and healthy controls is shown. (b) Upregulated expression of PD-1 on CD4+CD25+, CD4+CD69+, and CD4+HLA-DR+ T cells from RA patients, as compared with those from healthy controls. (c) RA patients homozygous for the rs36084323 SNP (G/G) have significantly increased percentages of PD-1+activated CD4+CD25+, CD4+CD69+, and CD4+HLA-DR+ T cells compared with healthy controls and rs36084323 A/A RA patients (* = , ** = by Mann-Whitney test). Horizontal bars indicate the mean ± SD.

4. Discussion

Rheumatoid arthritis is a chronic inflammatory disease that may involve extra-articular organs in addition to joints. Genetic and environmental factors are related to the pathogenesis of RA [1, 2]. PD-1 has been implicated as a critical pathway for tolerance since the first discovery of spontaneous autoimmune disease in PD-1 knockout mice. Studies using the NOD mouse model for spontaneous type I diabetes have shown that PD-1 not only is critical during the early phases of T-cell activation and expansion, but also plays a critical role in regulating T-cell effector functions and T-cell tolerance at late time points in peripheral tissue sites [16–19]. Strategies to selectively target the PD-1 pathway using antagonists have been effective in CIA model. Administration of PD-L1-Ig significantly ameliorated autoimmunity as assessed by clinical arthritis score and histology in the joints [20]. Similarly, PDL-1.Fc treatment ameliorated the severity of CIA and reduced joint swelling as well as antigen-specific T-cell proliferative responses [21]. PD-1 knockout mice have increased susceptibility and develop severe arthritis following immunization with type II collagen [21]. Collectively, these models demonstrate the importance of PD-1 in the pathogenesis of RA.

Greater than 30 single nucleotide polymorphisms (SNPs) have currently been identified in the human PDCD-1 gene, and several SNPs have been linked to development of autoimmunity in various populations [10]. Understanding how these SNPs impact PD-1 functions in human cells would be important for exploiting this pathway clinically. In our research, we selected four potentially functional polymorphisms (The PDCD-1 −606A/G (rs rs36084323) in promoter, PDCD-1 +7146A/G (rs11568821) in intron 4, and PDCD-1 +7625G/A (rs2227982) and PDCD-1 +7786G/C (rs2227981) in exon 5) of PDCD-1 gene and identified the association between the four polymorphisms and the risk of RA in Southeastern Chinese population.

A region of PDCD-1 intron 4 was described as an enhancer-like structure containing binding sites for several transcription factors [22]. G/A (PD-1.3) SNP, a regulatory SNP located in intron 4, showed to be involved in susceptibility to RA in Swedish, European American, and Mexican families and in sporadic cases. However, we found that this nucleotide is not polymorphic in our Chinese people, and this was also verified among other Chinese populations [13, 23, 24].

As it is known, mutations in the promoter region may influence the engagement between sequence motifs and transcription factors and further disrupt the activation of gene and the initiation of transcription [25]. Accordingly, as a polymorphism located in the promoter of PDCD-1 gene, rs36084323 may also affect the transcription and activation of PDCD-1 gene, thereby influencing the development of RA. These promoter polymorphisms have been considered risk factors for RA. One group observed that the frequency of SNP PD-1.1 (rs36084323) A allele is a risk allele in Italian RA patients [26]. Another group observed that the AA genotype of SNP PD-1.1 was associated with a decreased risk for developing RA in Hong Kong Chinese RA patients [13]. Previous reports also have shown that PD-1.1 was associated with risk of sporadic breast cancer [27]. We observed reduced frequencies of both the AA genotype and the A allele in RA patients compared with controls and a compensatory increase of the GG genotype in RA patients, which indicated that AA genotype might be a protective factor to decrease the risk of RA. This was consistent with the results from Hong Kong Chinese population [13]. Hence, rs36084323 may be a functional polymorphism which can influence the risk of diseases.

PD-1 is thought to be important for the “fine tuning” of lymphocyte activation at the level of synovial tissue, considering the wide pattern of expression of one of its ligands, PDL-1, in activated endothelial and epithelial cells [17, 28–30]. Furthermore, CD4+PD-1+ T cells accumulate as unique anergic cells in rheumatoid arthritis synovial fluid [31]. It is also reported that the negative costimulatory PD-1/PDL-1 pathway regulates peripheral T-cell responses in both human and murine RA [21]. In this study we found upregulated PD-1 expression on CD4+CD25+T cells, CD4+CD69+T cells, and CD4+HLA-DR+T cells compared with healthy controls. We also found that patients carrying G/G genotype displayed an increased mRNA level of PD-1 and increased basal PD-1 expression on activated CD4+ T cells. We concluded that PD-1 polymorphisms were associated with RA, and the GG genotype and G allele of rs36084323 that associated with induced PD-1 mRNA expression and PD-1 expression might be risk factor in RA development in the Han Chinese population.

In conclusion, the present study provides strong evidence that rs36084323 functional polymorphisms may contribute to the risk of RA. This SNP, which effectively influenced the expression of PD-1, may be a biomarker of early diagnosis of RA and a suitable indicator of utilizing PD-1 inhibitor for treatment of RA. However, our results were obtained from a sample-sized sample, and therefore this is a preliminary conclusion. Validation by a larger prospective study from a more diverse ethnic population is needed to confirm these findings.

Disclosure

All authors declare no commercial association, such as consultancies, stock ownership, or other equity interests or patent licensing arrangements.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

CuiPing Liu, JueAn Jiang, and Li Gao contributed equally to this work.

Acknowledgments

This work was supported by Grants from the National Basic Research Development Program of China (973 program, no. 2013CB53051), the National Natural Science Foundation of China (nos. 31170834 and 30930085), and the Science and Technology Project of Suzhou (no. SYS201331).

References

F. C. Arnett, S. M. Edworthy, D. A. Bloch et al., “The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis,” Arthritis and Rheumatism, vol. 31, no. 3, pp. 315–324, 1988.

L. Klareskog, A. I. Catrina, and S. Paget, “Rheumatoid arthritis,” The Lancet, vol. 373, no. 9664, pp. 659–672, 2009.

I. B. McInnes and G. Schett, “The pathogenesis of rheumatoid arthritis,” The New England Journal of Medicine, vol. 365, no. 23, pp. 2205–2219, 2011.

M. Gol-Ara, F. Jadidi-Niaragh, R. Sadria, et al., “The role of different subsets of regulatory T cells in immunopathogenesis of rheumatoid arthritis,” Arthritis, vol. 2012, Article ID 805875, 16 pages, 2012.

M. E. Keir, M. J. Butte, G. J. Freeman, and A. H. Sharpe, “PD-1 and its ligands in tolerance and immunity,” Annual Review of Immunology, vol. 26, pp. 677–704, 2008.

L. M. Francisco, P. T. Sage, and A. H. Sharpe, “The PD-1 pathway in tolerance and autoimmunity,” Immunological Reviews, vol. 236, no. 1, pp. 219–242, 2010.

R. V. Parry, J. M. Chemnitz, K. A. Frauwirth et al., “CTLA-4 and PD-1 receptors inhibit T-cell activation by distinct mechanisms,” Molecular and Cellular Biology, vol. 25, no. 21, pp. 9543–9553, 2005.

J. Wang, T. Yoshida, F. Nakaki, H. Hiai, T. Okazaki, and T. Honjo, “Establishment of NOD-Pdcd1-/- mice as an efficient animal model of type I diabetes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 33, pp. 11823–11828, 2005.

H. Nishimura, M. Nose, H. Hiai, N. Minato, and T. Honjo, “Development of lupus-like autoimmune diseases by disruption of the PD-1 gene encoding an ITIM motif-carrying immunoreceptor,” Immunity, vol. 11, no. 2, pp. 141–151, 1999.

T. Okazaki and T. Honjo, “PD-1 and PD-1 ligands: from discovery to clinical application,” International Immunology, vol. 19, no. 7, pp. 813–824, 2007.

G. K. Bertsias, M. Nakou, C. Choulaki et al., “Genetic, immunologic, and immunohistochemical analysis of the programmed death 1/programmed death ligand 1 pathway in human systemic lupus erythematosus,” Arthritis and Rheumatism, vol. 60, no. 1, pp. 207–218, 2009.

L. Prokunina, L. Padyukov, A. Bennet et al., “Association of the PD-1.3A allele of the PDCD1 gene in patients with rheumatoid arthritis negative for rheumatoid factor and the shared epitope,” Arthritis and Rheumatism, vol. 50, no. 6, pp. 1770–1773, 2004.

E. K. Kong, L. Prokunina-Olsson, W. H. Wong et al., “A new haplotype of PDCD1 is associated with rheumatoid arthritis in Hong Kong Chinese,” Arthritis and Rheumatism, vol. 52, no. 4, pp. 1058–1062, 2005.

A. Kroner, M. Mehling, B. Hemmer et al., “A PD-1 polymorphism is associated with disease progression in multiple sclerosis,” Annals of Neurology, vol. 58, no. 1, pp. 50–57, 2005.

C. Nielsen, D. Hansen, S. Husby, B. B. Jacobsen, and S. T. Lillevang, “Association of a putative regulatory polymorphism in the PD-1 gene with susceptibility to type 1 diabetes,” Tissue Antigens, vol. 62, no. 6, pp. 492–497, 2003.

R. Ni, K. Ihara, K. Miyako et al., “PD-1 gene haplotype is associated with the development of type 1 diabetes mellitus in Japanese children,” Human Genetics, vol. 121, no. 2, pp. 223–232, 2007.

M. E. Keir, S. C. Liang, I. Guleria et al., “Tissue expression of PD-L1 mediates peripheral T cell tolerance,” Journal of Experimental Medicine, vol. 203, no. 4, pp. 883–895, 2006.

B. T. Fife, I. Guleria, M. G. Bupp et al., “Insulin-induced remission in new-onset NOD mice is maintained by the PD-1-PD-L1 pathway,” Journal of Experimental Medicine, vol. 203, no. 12, pp. 2737–2747, 2006.

B. T. Fife, K. E. Pauken, T. N. Eagar et al., “Interactions between PD-1 and PD-L1 promote tolerance by blocking the TCR-induced stop signal,” Nature Immunology, vol. 10, no. 11, pp. 1185–1192, 2009.

G. Wang, P. Hu, J. Yang, G. Shen, and X. Wu, “The effects of PDL-Ig on collagen-induced arthritis,” Rheumatology International, vol. 31, no. 4, pp. 513–519, 2011.

A. P. Raptopoulou, G. Bertsias, D. Makrygiannakis et al., “The programmed death 1/programmed death ligand 1 inhibitory pathway is up-regulated in rheumatoid synovium and regulates peripheral T cell responses in human and murine arthritis,” Arthritis and Rheumatism, vol. 62, no. 7, pp. 1870–1880, 2010.

L. Prokunina, C. Castillejo-López, F. Öberg et al., “A regulatory polymorphism in PDCD1 is associated with susceptibility to systemic lupus erythematosus in humans,” Nature Genetics, vol. 32, no. 4, pp. 666–669, 2002.

S. Wang, Y. Chen, T. Ou et al., “Programmed death-1 gene polymorphisms in patients with systemic lupus erythematosus in Taiwan,” Journal of Clinical Immunology, vol. 26, no. 6, pp. 506–511, 2006.

Q. Meng, X. Liu, P. Yang et al., “PDCD1 genes may protect against extraocular manifestations in Chinese Han patients with Vogt-Koyanagi-Harada syndrome,” Molecular Vision, vol. 15, pp. 386–392, 2009.

K. M. K. de Vooght and W. W. van Solinge, “Gene promoter analysis in molecular diagnostics: do or don't?” Expert Review of Molecular Diagnostics, vol. 9, no. 5, pp. 403–405, 2009.

M. T. Tahoori, A. A. Pourfathollah, M. Akhlaghi, S. Daneshmandi, M. H. Nicknam, and N. Soleimanifar, “Association of programmed cell death-1 (PDCD-1) gene polymorphisms with rheumatoid arthritis in Iranian patients,” Clinical and Experimental Rheumatology, vol. 29, no. 5, pp. 763–767, 2011.

Z. Hua, D. Li, F. Xu et al., “PD-1 polymorphisms are associated with sporadic breast cancer in Chinese Han population of Northeast China,” Breast Cancer Research and Treatment, vol. 129, no. 1, pp. 195–201, 2011.

S. J. Khoury and M. H. Sayegh, “The roles of the new negative T cell costimulatory pathways in regulating autoimmunity,” Immunity, vol. 20, no. 5, pp. 529–538, 2004.

T. Okazaki and T. Honjo, “The PD-1-PD-L pathway in immunological tolerance,” Trends in Immunology, vol. 27, no. 4, pp. 195–201, 2006.

N. Martin-Orozco, Y. Wang, H. Yagita, and C. Dong, “Cutting edge: programed death (PD) ligand-1/PD-1 interaction is required for CD8+ T cell tolerance to tissue antigens,” Journal of Immunology, vol. 177, no. 12, pp. 8291–8295, 2006.

S. Hatachi, Y. Iwai, S. Kawano et al., “CD4+ PD-1+ T cells accumulate as unique anergic cells in rheumatoid arthritis synovial fluid,” Journal of Rheumatology, vol. 30, no. 7, pp. 1410–1419, 2003.