Molecular Vision 2009; 15:386-392 <http://www.molvis.org/molvis/v15/a40>
Received 20 December 2008 | Accepted 13 February 2009 | Published 20 February 2009

PDCD1 genes may protect against extraocular manifestations in Chinese Han patients with Vogt-Koyanagi-Harada syndrome

Qianli Meng,1 Xiaoli Liu,1,2 Peizeng Yang,2 Shengping Hou,1 Liping Du,1 Hongyan Zhou,1 Aize Kijlstra3,4

The first two authors contributed equally to this work.

1Zhongshan Ophthalmic Center, Sun Yat-sen University, Uveitis Study Center of Sun Yat-sen University and International Uveitis Study Laboratory of Guangdong Province, State Key Laboratory of Ophthalmology of Sun Yat-sen University, Guangzhou, P.R. China; 2The First Affiliated Hospital, Chongqing Medical University, Chongqing key laboratory of ophthalmology, Chongqing, P. R. China; 3Eye Research Institute Maastricht, Department of Ophthalmology, University Hospital Maastricht, Maastricht, the Netherlands; 4Animal Sciences Group, Wageningen UR, Lelystad, the Netherlands

Correspondence to: Dr. Peizeng Yang, MD, PhD, The First Affiliated Hospital, Chongqing Medical University, Chongqing Key Laboratory of Ophthalmology, Youyi Road 1, Chongqing, 400016, P. R. China; Phone: 0086-23-89012851; FAX: 0086-23-89012851; email: peizengycmu@126.com or peizengy@126.com

Abstract

Purpose: To analyze the potential association of programmed cell death 1 (PDCD1) with Vogt-Koyanagi-Harada (VKH) syndrome in a Chinese Han population.

Methods: Three single nucleotide polymorphism (SNPs), PD-1.3G/A, PD-1.5C/T, and PD-1.6G/A, were genotyped in 247 VKH patients and 289 age-, sex-, and ethnically-matched healthy controls using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay. The associations of genotypes and alleles with VKH syndrome were analyzed.

Results: All genotype distributions in healthy controls were in Hardy–Weinberg equilibrium. The genotype and allele frequencies of PD-1.3, PD-1.5, and PD-1.6 were not different between patients with VKH syndrome and healthy controls. No significant difference was observed according to the status of human leukocyte antigen (HLA)-DR4 and HLA-DRw53. Compared to the controls, lower frequencies of the PD-1.5C genotype and allele frequencies were observed in VKH patients with extraocular findings.

Conclusions: PD-1.3 and PD-1.6 polymorphisms are not associated with the susceptibility to VKH syndrome in the Chinese Han population. However, PD-1.5 may be negatively associated with the occurrence of extraocular manifestations of VKH syndrome.

Introduction

Vogt-Koyanagi-Harada (VKH) syndrome is an autoimmune disease characterized by a bilateral granulomatous panuveitis and systemic disorders including poliosis, vitiligo, alopecia, and central nervous system and auditory signs [1]. Although the etiology of VKH syndrome remains unclear, several studies have suggested that an autoimmune response against melanocytes or tyrosinase family proteins may play a key role in this disease [2,3]. Genetic predisposition is an important element as evidenced by familial cases [3] and a strong association with the human leukocyte antigen (HLA) system [4,5], especially HLA-DR4 and HLA-DRw53. However, the association between VKH syndrome and the HLA system does not fully explain the genetic risk for this disease. Moreover, little is known about the association of non-HLA genetic factors with VKH syndrome. Therefore, studies on the disease association with certain genes involved in the immune response may highlight the genetic background of VKH syndrome.

Programmed cell death 1 (PD-1), originally identified as a molecule linked to in vitro induction of apoptotic cell death in murine lymphoid cell lines [6], is a member of the cluster of differentiation (CD)28/B7 family. PD-1 contains an immunoreceptor tyrosine-based inhibiting motif (ITIM) and an immunoreceptor tyrosine-based switch motif (ITSM) [7] and is transcriptionally induced in activated T cells, B cells, and myeloid cells [8,9]. As an immunoinhibitory receptor, PD-1 has been shown to inhibit lymphocyte activation and cytokine production after interacting with its ligands, PD-L1 (B7-H1) and PD-L2 (B7-DC) [10-12]. The human gene encoding PD-1, i.e., PDCD1, is located on chromosome 2q37. Studies on PDCD1 polymorphisms have demonstrated its association with several autoimmune diseases including systemic lupus erythematosus (SLE) [13-17], rheumatoid arthritis (RA) [18], type I diabetes (TID) [19], multiple sclerosis (MS) [20], ankylosing spondylitis (AS) [21], and Graves' disease [22], although there are also some conflicting results [23-27]. Our present study was designed to investigate the association of PDCD1 single nucleotide polymorphisms (SNPs) including PD-1.3, PD-1.5, and PD-1.6 with VKH syndrome in a Chinese Han population.

Methods

Patients and healthy controls

Two hundred and forty-seven VKH patients with an average age at onset of 33.6 years (109 female and 138 male) and 289 age-, sex-, ethnically-matched healthy controls with an average age of 35.4 years (129 female and 160 male) were enrolled in this study. All patients were recruited from the Uveitis Study Center of the Sun Yat-sen University (Guangzhou, China) and the First Affiliated Hospital, Chongqing Medical University (Chongqing, China) and fulfilled the First International Workshop criteria for VKH disease [28]. These new criteria mainly include bilateral diffuse choroiditis at early phage and sunset glow funds and nummular chorioretinal depigmented scars at late phase in association with neurologic, auditory, or integumentary findings. Additionally but importantly, there is no history of penetrating ocular trauma or surgery preceding the initial onset of uveitis and there is no clinical or laboratory evidence suggestive of other ocular disease entities. The clinical characteristics of the tested VKH patients were summarized in Table 1. Controls were mainly the accompanying person or spouses of the patients. The local institutional ethics committee approved the study, and written informed consent was obtained from all the study subjects. Blood samples were collected in EDTA tubes and kept at −70 °C until use.

Genomic DNA extraction and genotyping

Genomic DNA was isolated from peripheral blood of patients and controls using standard proteinase K digestion and phenol-chloroform extraction. Amplification of the target DNA in PDCD1 was performed by polymerase chain reaction (PCR). We used those primers previously determined by Ferreuris-Vidal, Velazquez-Cruz, and Kong (Table 2) [14,17,18]. A 15 μl reaction mixture, which consisted of 7.5 μl Premix Taq (Ex Taq Version; TaKaRa Biotechnology Co. Ltd., Dalian, China), 50 pmoles primers, and 0.2 μg of genomic DNA, was amplified by PCR. PCR conditions were as follows: initial denaturation at 95 °C for 5 min followed by 35 cycles of denaturation at 94 °C for 30 s, annealing at different temperatures (67 °C for PD-1.3, 60 °C for PD-1.5, and 61 °C for PD-1.6) for 30 s, extension at 72 °C for 30 s, and a final extension at 72 °C for 5 min.

The SNPs were genotyped by PCR-restriction fragment length polymorphism (RFLP) analysis. PCR products of PD-1.3, PD-1.5, and PD-1.6 polymorphisms were respectively digested with 2 U of PstI (Fermentas international Inc., Ontario, Canada), Pvu II (Fermentas international Inc.), and Nla III (New England Biolabs, Inc., Ontario, Canada) restriction enzymes (Table 2) in a 10 µl reaction volume overnight. Digestion products were visualized on agarose gels of appropriate concentration and stained with GoldView™ (SBS Genetech, Beijing, China). To confirm the accuracy of the method employed, randomly selected subjects (20% of all samples) were analyzed by direct sequencing (Invitrogen Biotechnology Co., Guangzhou, China). Appropriate controls (no template and known genotype) were included in each typing run. HLA-DR4 and HLA-DRw53 typing was performed using a PCR sequence specific primers (SSP) method as previously described [29].

Statistical analysis

Hardy–Weinberg equilibrium (HWE) was tested using the χ2 test. Allele frequencies were estimated by direct counting. Allele and genotype frequencies were compared between patients and controls by the χ2 test using SPSS (version 10.0; SPSS Inc., Chicago, IL). The p values were corrected (pc) with the Bonferroni correction by multiplying with the number of analyses performed. pc<0.05 was considered significant.

Results

Association analysis of single nucleotide polymorphism and haplotype

Three SNPs were successfully genotyped in 247 patients with VKH syndrome and 289 healthy controls. Table 3 shows the genotype and allele distribution of these SNPs. The distribution of genotype frequencies of each SNP in the healthy controls did not show any significant deviation from the Hardy–Weinberg equilibrium. The frequency of the C allele of SNP PD-1.5 tended to be lower in VKH patients than in healthy controls (71.5% versus 78.0%, p=0.013). However, no difference was found when the Bonferroni correction was applied (pc=0.078, n=6). Similarly, although a lower frequency of SNP PD-1.5CC genotype was observed in VKH patients as compared with that in healthy controls, the difference was not statistically significant following the Bonferroni correction (p=0.034, pc=0.306, n=9). The genotype and allele frequencies of PD-1.3 and PD-1.6 in patients with VKH syndrome were not different from those in healthy controls. In particular, only the GG genotype and the G allele of PD-1.3 were identified in all patients and controls.

Haplotype analysis using Haploview software showed no linkage disequilibrium in the tested three SNPs in both patients with VKH syndrome and healthy controls.

Stratification analysis according to HLA status and clinical characteristics

The frequencies of HLA-DR4 and HLA-DRw53 were shown to be significantly increased in 231 VKH patients as compared with those in 289 healthy controls (Table 4). The allele and genotype frequencies of PD-1.3, PD-1.5, and PD-1.6 were not different between VKH patients and healthy controls when a stratification analysis was performed according to the status of HLA-DR4 and HLA-DRw53.

When the genotype frequencies were analyzed according to the clinical features, our results showed significantly lower frequencies of PD-1.5CC in VKH patients either with poliosis or with dysacusis as compared to that observed in healthy controls. Concerning the allele frequencies analyzed according to the clinical features within the VKH cohort, it was found that significantly lower frequencies of the PD-1.5C allele were associated with VKH patients with extraocular findings including alopecia, poliosis, dysacusis, and tinnitus (Table 5). Neither PD-1.3 nor PD-1.6 was found to be associated with any of the extraocular findings including alopecia, poliosis, dysacusis, tinnitus, and neck stiffness.

Discussion

Our results showed that a lower frequency of PD-1.5 genotype or allele was associated with certain extraocular findings within the VKH cohort. However, PD-1.3 and PD-1.6 polymorphism differences were identified between VKH patients and healthy controls.

VKH syndrome is a multifactorial disease that may result from interactions among susceptibility genes, environmental factors, and immunological responses. The association of HLA-DR4 and HLADRw53 with VKH syndrome has been reported in the Chinese [4] and Japanese [5] populations. Our present study confirmed these observations. The non-HLA genetic background of the disease has not been completely understood. It is reported that there are no genetic susceptibilities for the tyrosinase gene family and interferon-γ in Japanese VKH patients [30,31]. Recently, we reported an association of cytotoxic T lymphocyte-associated antigen-4 (CTLA-4) genetic polymorphisms with the susceptibility to VKH syndrome in a Chinese population [32]. PD-1 is an immunoglobulin (Ig) superfamily member related to CD28 and CTLA-4. One hundred and thirty-five SNPs (found in the NCBI-Entrez SNP database) have been identified in the human PDCD1 region. Several studies showed that polymorphisms of PD-1.3, PD-1.5, and PD-1.6 were associated with autoimmune diseases [13-22], although inconsistent results were also reported [23-27]. In this study, we tested whether the polymorphisms of these three SNPs were associated with VKH syndrome in Han Chinese patients.

To ensure the analysis results, the following attempts were made. First, all genotype distributions of three SNPs in healthy controls were tested and found to be in Hardy–Weinberg equilibrium. Second, the controls and patients were strictly matched according to the places where they were born to exclude the possible influence of stratification of the population. VKH patients of Chinese Han descendents were also strictly selected in this study to avoid the influence of gene background. Third, a total of 247 VKH patients and 289 age-, sex-, ethnically-matched healthy controls were used in this study, and the number of tested samples was large enough to avoid a bias of the results. Finally, direct sequencing was performed to validate the genotype findings.

As genetic interaction between HLA genes and non-HLA genes may influence the susceptibility to autoimmune disease, a stratification analysis was performed according to the status of HLA-DR4 and HLA-DRw53. Unexpectedly, no association with PDCD1 was observed in HLA-DR4 and HLA-DRw53 positive or negative VKH patients. As VKH syndrome affects not only the eye but also skin and hair, we further analyzed the association of polymorphisms of three SNPs with extraocular findings in VKH patients. Lower frequencies of genotype and allele of PD-1.5 were observed in VKH patients with any of the extraocular findings including alopecia, poliosis, dysacusis, or tinnitus. These results seem to suggest that the genotype and allele of PD-1.5 may possibly protect VKH patients from extraocular findings.

Several previous studies have shown that there is a large variation in the frequencies of PDCD1 polymorphisms among different ethnic groups. In the present study, we found a higher frequency of the PD-1.5C allele and PD-1.6A allele in Chinese (78%, 67%, respectively) than in Swedish (60%, 8%, respectively), European American (57%, 11%, respectively), Mexican (54%, 39%. respectively) and African American individuals (42%, 48%, respectively) [13]. A higher frequency of the PD-1.3A allele was reported in the European population (5%–12%) compared to that observed in the Mexican and African American populations (2%–7%) [13,16] Interestingly, we did not detect the PD-1.3A allele in either the healthy controls or the VKH patients in this study, which is consistent with earlier results in healthy Hong Kong Chinese [18] and Japanese populations [27]. With regards to the association of PD-1.5 polymorphism with autoimmune disease, a study by Lee et al. [21] showed that a PD-1.5 polymorphism was positively associated with ankylosing spondylitis in Korean patients. Kong and Iwamoto found that a polymorphism of PD-1.5 was not associated with rheumatoid arthritis in Hong Kong Chinese [18] and Japanese [27] patients.

In conclusion, this study showed that all genotype distributions of PD-1.3, PD-1.5, and PD-1.6 in Chinese Han healthy controls were in Hardy–Weinberg equilibrium. Polymorphisms of PD-1.3 and PD-1.6 were not associated with the susceptibility to VKH syndrome in a Chinese Han population. However, a genotype or allele of PD-1.5 was negatively associated with VKH patients having accompanying extraocular findings. It is not clear whether this association is present in other ethnic VKH patients and whether other PDCD1 polymorphisms are associated with the susceptibility to VKH syndrome. More studies are needed to clarify these issues within the VKH cohort.

Acknowledgments

This work was supported by the Key Project of Natural Science Foundation (30630064), National Supporting Project of P.R. China, Natural Science Foundation of Guangdong Province (7301447), Key Project of Chongqing Health Bureau, Clinical Key Project of Ministry of Heath, and Natural Science Foundation for Research Groups of Guangdong Province (05200176). We would like thank all subjects participated in this study.

References

  1. Yang P, Ren Y, Li B, Fang W, Meng Q, Kijlstra A. Clinical characteristics of Vogt-Koyanagi-Harada syndrome in Chinese patients. Ophthalmology. 2007; 114:606-14. [PMID: 17123618]
  2. Norose K, Yano A. Melanoma specific Th1 cytotoxic T lymphocyte lines in Vogt-Koyanagi-Harada disease. Br J Ophthalmol. 1996; 80:1002-8. [PMID: 8976730]
  3. Yamaki K, Gocho K, Hayakawa K, Kondo I, Sakuragi S. Tyrosinase family proteins are antigens specific to Vogt-Koyanagi-Harada disease. J Immunol. 2000; 165:7323-9. [PMID: 11120868]
  4. Zhao M, Jiang Y, Abrahams IW. Association of HLA antigens with Vogt-Koyanagi-Harada syndrome in a Han Chinese population. Arch Ophthalmol. 1991; 109:368-70. [PMID: 2003797]
  5. Islam SM, Numaga J, Fujino Y, Hirata R, Matsuki K, Maeda H, Masuda K. HLA class II genes in Vogt-Koyanagi-Harada disease. Invest Ophthalmol Vis Sci. 1994; 35:3890-6. [PMID: 7928186]
  6. Ishida Y, Agata Y, Shibahara K, Honjo T. Induced expression of PD-1, a novel member of the immunoglobulin gene superfamily, upon programmed cell death. EMBO J. 1992; 11:3887-95. [PMID: 1396582]
  7. Shlapatska LM, Mikhalap SV, Berdova AG, Zelensky OM, Yun TJ, Nichols KE, Clark EA, Sidorenko SP. CD150 association with either the SH2-containing inositol phosphatase or the SH2-containing protein tyrosine phosphatase is regulated by the adaptor protein SH2D1A. J Immunol. 2001; 166:5480-7. [PMID: 11313386]
  8. Agata Y, Kawasaki A, Nishimura H, Ishida Y, Tsubata T, Yagita H, Honjo T. Expression of the PD-1 antigen on the surface of stimulated mouse T and B lymphocytes. Int Immunol. 1996; 8:765-72. [PMID: 8671665]
  9. Greenwald RJ, Freeman GJ, Sharpe AH. The B7 family revisited. Annu Rev Immunol. 2005; 23:515-48. [PMID: 15771580]
  10. Salama AD, Chitnis T, Imitola J, Ansari MJ, Akiba H, Tushima F, Azuma M, Yagita H, Sayegh MH, Khoury SJ. Critical role of the programmed death-1 (PD-1) pathway in regulation of experimental autoimmune encephalomyelitis. J Exp Med. 2003; 198:71-8. [PMID: 12847138]
  11. Kitazawa Y, Fujino M, Wang Q, Kimura H, Azuma M, Kubo M, Abe R, Li XK. Involvement of the programmed death-1/programmed death-1 ligand pathway in CD4+CD25+ regulatory T-cell activity to suppress alloimmune responses. Transplantation. 2007; 83:774-82. [PMID: 17414712]
  12. Sharpe AH, Wherry EJ, Ahmed R, Freeman GJ. The function of programmed cell death 1 and its ligands in regulating autoimmunity and infection. Nat Immunol. 2007; 8:239-45. [PMID: 17304234]
  13. Prokunina L, Castillejo-Lopez C, Oberg F, Gunnarsson I, Berg L, Magnusson V, Brookes AJ, Tentler D, Kristjansdottir H, Grondal G, Bolstad AI, Svenungsson E, Lundberg I, Sturfelt G, Jonssen A, Truedsson L, Lima G, Alcocer-Varela J, Jonsson R, Gyllensten UB, Harley JB, Alarcon-Segovia D, Steinsson K, Alarcon-Riquelme ME. A regulatory polymorphism in PDCD1 is associated with susceptibility to systemic lupus erythematosus in humans. Nat Genet. 2002; 32:666-9. [PMID: 12402038]
  14. Ferreiros-Vidal I, Gomez-Reino JJ, Barros F, Carracedo A, Carreira P, Gonzalez-Escribano F, Liz M, Martin J, Ordi J, Vicario JL, Gonzalez A. Association of PDCD1 with susceptibility to systemic lupus erythematosus: evidence of population-specific effects. Arthritis Rheum. 2004; 50:2590-7. [PMID: 15334473]
  15. Sanghera DK, Manzi S, Bontempo F, Nestlerode C, Kamboh MI. Role of an intronic polymorphism in the PDCD1 gene with the risk of sporadic systemic lupus erythematosus and the occurrence of antiphospholipid antibodies. Hum Genet. 2004; 115:393-8. [PMID: 15322919]
  16. Thorburn CM, Prokunina-Olsson L, Sterba KA, Lum RF, Seldin MF, Alarcon-Riquelme ME, Criswell LA. Association of PDCD1 genetic variation with risk and clinical manifestations of systemic lupus erythematosus in a multiethnic cohort. Genes Immun. 2007; 8:279-87. [PMID: 17344889]
  17. Velazquez-Cruz R, Orozco L, Espinosa-Rosales F, Carreno-Manjarrez R, Solis-Vallejo E, Lopez-Lara ND, Ruiz-Lopez IK, Rodriguez-Lozano AL, Estrada-Gil JK, Jimenez-Sanchez G, Baca V. Association of PDCD1 polymorphisms with childhood-onset systemic lupus erythematosus. Eur J Hum Genet. 2007; 15:336-41. [PMID: 17228327]
  18. Kong EK, Prokunina-Olsson L, Wong WH, Lau CS, Chan TM, Alarcon-Riquelme M, Lau YL. A new haplotype of PDCD1 is associated with rheumatoid arthritis in Hong Kong Chinese. Arthritis Rheum. 2005; 52:1058-62. [PMID: 15818672]
  19. Nielsen C, Hansen D, Husby S, Jacobsen BB, Lillevang ST. Association of a putative regulatory polymorphism in the PD-1 gene with susceptibility to type 1 diabetes. Tissue Antigens. 2003; 62:492-7. [PMID: 14617032]
  20. Kroner A, Mehling M, Hemmer B, Rieckmann P, Toyka KV, Maurer M, Wiendl HA. PD-1 polymorphism is associated with disease progression in multiple sclerosis. Ann Neurol. 2005; 58:50-7. [PMID: 15912506]
  21. Lee SH, Lee YA, Woo DH, Song R, Park EK, Ryu MH, Kim YH, Kim KS, Hong SJ, Yoo MC, Yang HI. Association of the programmed cell death 1 (PDCD1) gene polymorphism with ankylosing spondylitis in the Korean population. Arthritis Res Ther. 2006; 8:R163 [PMID: 17064404]
  22. Newby PR, Roberts-Davies EL, Brand OJ, Heward JM, Franklyn JA, Gough SC, Simmonds MJ. Tag SNP screening of the PDCD1 gene for association with Graves' disease. Clin Endocrinol (Oxf). 2007; 67:125-8. [PMID: 17490403]
  23. Lin SC, Yen JH, Tsai JJ, Tsai WC, Ou TT, Liu HW, Chen CJ. Association of a programmed death 1 gene polymorphism with the development of rheumatoid arthritis, but not systemic lupus erythematosus. Arthritis Rheum. 2004; 50:770-5. [PMID: 15022318]
  24. Abelson AK, Johansson CM, Kozyrev SV, Kristjansdottir H, Gunnarsson I, Svenungsson E, Jonsen A, Lima G, Scherbarth HR, Gamron S, Allievi A, Palatnik SA, Alvarellos A, Paira S, Graf C, Guilleron C, Catoggio LJ, Prigione C, Battagliotti CG, Berbotto GA, Garcia MA, Perandones CE, Truedsson L, Steinsson K, Sturfelt G, Pons-Estel B. Argentinean Collaborative Group, Alarcon-Riquelme ME. No evidence of association between genetic variants of the PDCD1 ligands and SLE. Genes Immun. 2007; 8:69-74. [PMID: 17136123]
  25. Asad S, Nikamo P, Torn C, Landin-Olsson M, Lernmark A, Alarcon-Riquelme M, Kockum I, Diabetes Incidence in Sweden Study Group.. No evidence of association of the PDCD1 gene with Type 1 diabetes. Diabet Med. 2007; 24:1473-7. [PMID: 18042085]
  26. Ferreiros-Vidal I, D'Alfonso S, Papasteriades C, Skopouli FN, Marchini M, Scorza R, Migliaresi S, Sebastiani GD, Endreffy E, Mavromati M, Kappou-Rigatou I, Ruzickova S, Dostal C, Schmidt RE, Witte T, Gomez-Reino JJ, Gonzalez A. Bias in association studies of systemic lupus erythematosus susceptibility due to geographical variation in the frequency of a programmed cell death 1 polymorphism across Europe. Genes Immun. 2007; 8:138-46. [PMID: 17230193]
  27. Iwamoto T, Ikari K, Inoue E, Toyama Y, Hara M, Yamanaka H, Tomatsu T, Momohara S, Kamatani N. Failure to confirm association between PDCD1 polymorphisms and rheumatoid arthritis in a Japanese population. J Hum Genet. 2007; 52:557-60. [PMID: 17468813]
  28. Read RW, Holland GN, Rao NA, Tabbara KF, Ohno S, Arellanes-Garcia L, Pivetti-Pezzi P, Tessler HH, Usui M. Revised diagnostic criteria for Vogt-Koyanagi-Harada disease: report of an international committee on nomenclature. Am J Ophthalmol. 2001; 131:647-52. [PMID: 11336942]
  29. Bunce M, O'Neill CM, Barnardo MC, Krausa P, Browning MJ, Morris PJ, Welsh KI. Phototyping: comprehensive DNA typing for HLA-A, B, C, DRB1, DRB3, DRB4, DRB5 & DQB1 by PCR with 144 primer mixes utilizing sequence-specific primers (PCR-SSP). Tissue Antigens. 1995; 46:355-67. [PMID: 8838344]
  30. Horie Y, Takemoto Y, Miyazaki A, Namba K, Kase S, Yoshida K, Ota M, Hasumi Y, Inoko H, Mizuki N, Ohno S. Tyrosinase gene family and Vogt-Koyanagi-Harada disease in Japanese patients. Mol Vis. 2006; 12:1601-5. [PMID: 17200659]
  31. Horie Y, Kitaichi N, Takemoto Y, Namba K, Yoshida K, Hirose S, Hasumi Y, Ota M, Inoko H, Mizuki N, Ohno S. Polymorphism of IFN-gamma gene and Vogt-Koyanagi-Harada disease. Mol Vis. 2007; 13:2334-8. [PMID: 18199975]
  32. Du L, Yang P, Hou S, Lin X, Zhou H, Huang X, Wang L, Kijlstra A. Association of the CTLA-4 gene with Vogt-Koyanagi-Harada syndrome. Clin Immunol. 2008; 127:43-8. [PMID: 18282809]