Molecular Vision 2010; 16:1874-1879 <http://www.molvis.org/molvis/v16/a203>
Received 27 July 2010 | Accepted 9 September 2010 | Published 15 September 2010

Genotyping HLA-DRB1 and HLA-DQB1 alleles in Japanese patients with normal tension glaucoma

Misa Suzuki,1 Akira Meguro,1 Masao Ota,2 Eiichi Nomura,1 Tetsuo Kato,1 Naoko Nomura,1 Kenji Kashiwagi,3 Fumihiko Mabuchi,3 Hiroyuki Iijima,3 Kazuhide Kawase,4 Tetsuya Yamamoto,4 Makoto Nakamura,5 Akira Negi,5 Takeshi Sagara,6 Teruo Nishida,6 Masaru Inatani,7 Hidenobu Tanihara,7 Makoto Aihara,8 Makoto Araie,8 Takeo Fukuchi,9 Haruki Abe,9 Tomomi Higashide,10 Kazuhisa Sugiyama,10 Takashi Kanamoto,11 Yoshiaki Kiuchi,11 Aiko Iwase,12 Shigeaki Ohno,13 Hidetoshi Inoko,14 Nobuhisa Mizuki1

The first two authors contributed equally to this work

1Department of Ophthalmology and Visual Science Yokohama City University Graduate School of Medicine, Yokohama, Kanagawa, Japan; 2Department of Legal Medicine, Shinshu University School of Medicine, Matsumoto, Nagano, Japan; 3Department of Ophthalmology, University of Yamanashi, Faculty of Medicine, Yamanashi, Japan; 4Department of Ophthalmology, Gifu University Graduate School of Medicine, Gifu, Japan; 5Department of Surgery, Division of Ophthalmology, Kobe University Graduate School of Medicine, Kobe, Hyogo, Japan; 6Department of Biomolecular Recognition and Ophthalmology, Yamaguchi University School of Medicine, Ube, Yamaguchi, Japan; 7Department of Ophthalmology and Visual Science, Graduate School of Medical Sciences, Kumamoto University, Kumamoto, Japan; 8Department of Ophthalmology, University of Tokyo School of Medicine, Tokyo, Japan; 9Division of Ophthalmology and Visual Science, Graduated School of Medical and Dental Sciences, Niigata University, Niigata, Japan; 10Department of Ophthalmology and Visual Science, Kanazawa University Graduate School of Medical Science, Kanazawa, Ishikawa, Japan; 11Department of Ophthalmology and Visual Science, Graduate School of Biomedical Sciences, Hiroshima University, Hiroshima, Japan; 12Department of Ophthalmology, Tajimi Municipal Hospital, Tajimi, Gifu, Japan; 13Department of Ophthalmology, Hokkaido University Graduate School of Medicine, Sapporo, Japan; 14Department of Genetic Information, Division of Molecular Life Science, Tokai University School of Medicine, Isehara, Kanagawa, Japan

Correspondence to: Nobuhisa Mizuki, Department of Ophthalmology and Visual Science, Yokohama City University Graduate School of Medicine, 3-9 Fukuura, Kanazawa-ku, Yokohama, Kanagawa 236-0004, Japan; Phone: +81 45 787 2683; FAX: +81 45 781 9755; email: mizunobu@med.yokohama-cu.ac.jp

Abstract

Purpose: Normal tension glaucoma (NTG) is a subtype of glaucoma in which intraocular pressure is within the statistically normal range. NTG may be associated with an immune disorder. The aim of this study was to determine whether specific alleles in the human leukocyte antigen (HLA)-DRB1 and HLA-DQB1 genes correlated with NTG in Japanese patients.

Methods: We genotyped the HLA-DRB1 and HLA-DQB1 alleles in 113 Japanese patients with NTG and in 184 healthy Japanese control subjects using the polymerase chain reaction-sequence-specific oligonucleotide probes (PCR-SSOP) Luminex method. We assessed the allelic diversity in patients and controls.

Results: There were no statistically significant differences in the allele frequency of HLADRB1 and HLA-DQB1 between NTG patients and control subjects, and no HLA-DRB1-HLA-DQB1 haplotypes demonstrated any significant association with NTG.

Conclusions: Our findings suggest that HLA-DRB1 and HLA-DQB1 polymorphisms have no significant effect on the development of NTG in Japanese patients.

Introduction

Glaucoma is a progressive optic neuropathy leading to permanent visual loss that affects approximately 70 million people worldwide [1,2]. The most common type of glaucoma is primary open angle glaucoma (POAG). Glaucoma is generally a multifactorial disorder caused by environmental and/or genetic factors. High intraocular pressure (IOP) is the most important risk factor for glaucoma [3] and has previously been used as a diagnostic criterion for the disease. A subset of patients with POAG have normal IOP, a condition called normal tension glaucoma (NTG) [4-6]; thus, factors other than IOP clearly influence the pathogenesis of glaucomatous optic neuropathy. The prevalence of NTG is 3.6% in Japan, which is higher than in other ethnicities [7]. NTG often goes undiagnosed until after vision deteriorates because the symptoms are minimal and IOP is normal. Although relatively higher IOP, myopia, and older age are known to be factors associated with the development of NTG [8], they are not pathognomonic. Other pathogenic factors such as vasospasm, ischemia, and vascular insufficiencies are also indicated to be associated with the development of NTG [9-12].

It was recently suggested that the autoimmune system plays an important role in the development of glaucoma, especially NTG. Wax [13] proposed that the role of the immune system in glaucoma was twofold and that an autoimmune mechanism may elicit damage to the optic nerve, resulting in glaucomatous injury. He additionally suggested that autoimmune-mediated glaucoma injury occurs most often in NTG patients. There have been several reports of autoantibodies against ocular antigens that are upregulated in the serum of patients with glaucoma, such as heat shock proteins (HSPs) [14], gamma-enolase [15], glutathione S-transferase [16], anti-phosphatidylserine [17], and glycosaminoglycans [18]. Therefore, NTG might actually be an organ-specific autoimmune disorder.

Human leukocyte antigen (HLA) genes regulate the immune system and the predisposition to most autoimmune disorders. HLA class II alleles have been shown to be significantly associated with glaucoma. Gil-Carrasco et al. [19] reported a higher frequency of the HLA-DR3 antigen in Mexican POAG patients compared with controls; further, an increased frequency of the HLA-DRB1*0407-DQB1*0302 haplotype was significantly associated with POAG in Mexican patients [20]. Ferreri et al. [21] showed a significant association between DR11 and DQ1 alleles and POAG in the Italian population. Therefore, HLA class II gene polymorphisms may influence the development of autoimmunity in glaucoma. On the other hand, HLA class I alleles have also been shown to be associated with glaucoma [22-24]. However numerous studies have reported no significant correlation between HLA class I alleles and glaucoma [25-28], and thus the possibility of HLA class I alleles being associated with the pathogenesis of NTG is deemed low.

The aim of this study was to clarify the role of HLA class II gene polymorphisms in the development of NTG by assessing the association of polymorphisms in HLA-DRB1 and HLA-DQB1 genes with NTG in Japanese patients.

Methods

Subjects

We recruited 113 unrelated Japanese patients with NTG and 184 unrelated healthy Japanese control subjects at the following Japanese universities: Yokohama City University, Yamanashi University, Gifu University, Kobe University, Yamaguchi University, Kumamoto University, Hokkaido University, Tokyo University, Niigata University, Kanazawa University, Hiroshima University, Tajimi municipal hospital, and Tokai University. Patients were diagnosed using the following strict inclusion criteria: the presence of glaucomatous optic neuropathy with corresponding visual field loss; normal open angle with angle width of Shaffer grade 2 or higher; absence of IOP greater than 21 mmHg on repeat pressure measurement using Goldmann applanation tonometry without medication; and lack of a pathological basis for optic nerve changes upon neurologic, rhinologic, and general medical examination, including magnetic resonance imaging. Glaucomatous optic nerve change was diagnosed when the vertical cup/disc ratio of the optic nerve head was 0.7 or higher, the rim width at the superior (11 to 1 o’clock) or inferior (5 to 7 o’clock) position was less than or equal to 10% of the disc diameter, the difference in the vertical cup/disc ratio between eyes was 0.2 or greater, or a nerve fiber layer defect was found. Glaucomatous visual field loss was defined on a hemifield basis using reliable field data examined by the Humphrey® static visual field analyzer (HFA) C-30–2 program (Carl Zeiss Meditec, Oberkochen, Germany) according to Anderson and Patella’s criteria [29]. The hemifield was considered abnormal when the pattern deviation probability plot revealed a cluster of 3 or more non-edge contiguous points having sensitivity with a probability of less than 5% in the upper or lower hemifield, with one of these points having a probability of less than 1%. In addition, the following inclusion and exclusion criteria were used. We excluded individuals who were diagnosed when they were under 20 years old or over 60 years old with a –8.0D or higher myopic refractive error of spherical equivalence. The selection criteria for the HFA mean deviation were stratified depending on the subjects’ ages to minimize the effect of aging on retinal ganglion cell loss and subsequent visual field defects. Specifically, there was no exclusion if the patient was diagnosed when under 50 years old; when the patient was –10.00 dB or worse in at least one eye but was diagnosed between the ages of 50 and 55 years; or when the patients had –15.00 dB or worse in at least one eye but the patient was over age 55 when diagnosed. Patients with relatively early onset were selected because early onset suggests a stronger involvement of genetic factors. During diagnosis, patients were excluded if their refraction values had changed due to cataract surgery, refractive surgery, etc. In cases in which glaucomatous visual field loss was present in only one eye, the refraction value and glaucomatous visual field loss of the affected eye were used. In cases in which glaucomatous visual field loss was present in both eyes, the refraction value and glaucomatous visual field loss of the more severely affected eye were used.

Patients ranged in age between 24 and 59 years (mean 47.4±6.9 years), and 47.8% were male. The mean refraction value was –3.82±3.02 D, and the mean deviation observed in the Humphrey® static visual field determination (Carl Zeiss Meditec) was –10.15±8.01 dB. The control subjects were healthy volunteers from a geographic region similar to that of the NTG patients. They had no glaucoma or ophthalmological or systemic diseases that might cause glaucoma or optic nerve changes, and they had either no myopia or mild myopia with refractive errors of –3.00 D or less. Of the control subjects, 44.6% were male.

HLA-DRB1 and HLA-DQB1 genotyping

Genotyping of the HLA-DRB1 and HLA-DQB1 alleles was performed using the PCR-SSOP Luminex method using LABType SSO (One Lambda Inc., Canoga Park, CA), a reverse SSO DNA typing system, according to the manufacturer’s instructions. Briefly, target DNAs were PCR-amplified using HLA-DRB1 or HLA-DQB1 group-specific primers. The biotinylated PCR products were denatured and hybridized to locus-specific probes conjugated to fluorescently coded microspheres. A flow analyzer, LABScanTM 100, identifies the fluorescent intensity of phycoerythrin on each microsphere. The determination of the HLA allele was based on the reaction pattern compared to patterns associated with published HLA allele sequences.

Statistical analysis

Allelic frequencies were determined by direct counting. The significance of the allelic distribution between the patients with NTG and normal controls was first analyzed using Fisher’s exact probability test. The probability of association was corrected by the Bonferroni inequality method, i.e., by multiplying the p values with the number of alleles. A corrected p (pc) value of <0.05 was considered statistically significant.

Haplotype frequencies and linkage disequilibrium (LD) in the multi-locus analyses were calculated using PyPop [30]. Haplotype frequencies were determined using the iterative expectation-maximization algorithm. LD was measured using Hedrick’s multiallelic D´ statistic [31]. D´ weights the contribution of specific allele pairs to the LD using the product of their allelic frequencies. The measure is normalized to fall between zero and one, with higher values indicating a stronger LD contribution.

Results

The gene frequencies of HLA-DRB1 and HLA-DQB1 alleles in 113 Japanese patients with NTG and in184 healthy control subjects are shown in Table 1 and Table 2, respectively. A total of 27 HLA-DRB1 alleles and 14 HLA-DQB1 alleles were identified in NTG patients and healthy controls. The frequencies of HLA-DRB1*1301 and HLA-DQB1*0603 were higher in patients with NTG compared with healthy controls, although this increase did not reach statistical significance after performing the Bonferroni correction. We found no significant differences in the frequency of other HLA-DRB1 and HLA-DQB1 alleles between the patients and controls.

The magnitude of LD between HLA-DRB1 and HLA-DQB1 was extremely high, with a pair-wise D´>0.85. Haplotype analyses indicated complete LD between the HLADRB1*1301 and HLA-DQB1*0603 alleles, and the HLA-DRB1*1301-HLA-DQB1*0603 haplotype had an increased frequency in patients; this was not significant after correction (Table 3).

Discussion

Several studies have reported associations between HLA and glaucoma, and others have suggested that HLA class II genes may have an effect on the development of POAG [19-21]. The aim of this study was to investigate whether specific alleles in the HLA class II genes affected the development of NTG, which was recently suggested to be an autoimmune disorder. Toward this end, we genotyped the HLA-DRB1 and HLA-DQB1 genes in Japanese NTG patients and in healthy control subjects.

In contrast to previous POAG studies, this study focused on NTG and did not find an allelic distribution of HLA class II genes that was statistically different between NTG patients and controls. Although the frequencies of the HLA-DRB1*1301 and HLA-DQB1*0603 alleles and the HLA-DRB1*1301-HLA-DQB1*0603 haplotype was higher in NTG patients, these increases were not statistically significance after the Bonferroni correction. These HLA types were not consistent with POAG-associated HLA types. In addition, the frequency of POAG-associated HLA types was not significantly different in patients and controls in this NTG study, and these HLA types are rare in the Japanese population.

Our NTG study did not replicate the results of previous POAG studies, most likely because NTG and POAG have distinct characteristics. The genes affecting the immune systems of NTG patients may differ from those affecting POAG patients. In addition, the distribution of the alleles of HLA genes is different for each ethnic group. It is possible that there were false-positive results in the earlier POAG studies: In HLA disease association studies, careful study design is required to observe reliable positive associations. Further, the studies must be performed using adequate sample sizes, careful multiple testing correction, and LD analysis. Although Gil-Carrasco et al. [19] and Ferreri et al. [21] reported the associations of HLA class II alleles with POAG, the sample sizes in those studies were not big enough to support their positive findings, since positive results in small studies have a high risk of being false-positives. In another study that reported a significant association between the HLA-DRB1*0407-HLA-DQB1*0302 haplotype and POAG [20], the p values were corrected by the Yates method but not by multiple testing; this, too, might lead to a false positive association. In addition, the study reported that the estimated frequency for the HLA-DRB1*0407- HLA-DQB1*0302 haplotype in the POAG patients was 25.9%, much higher than the 15.4% frequency reported for the HLA-DRB1*0407 allele in the patients. Based on the allele frequency shown in the study, it is possible that actual frequency of the haplotype was less than estimated and that the haplotype may be not associated with POAG. Recently, Yano et al. [32] reported a possible association between the HLA-DRB1*1502 allele and the presence of anti-neurofilament light subunit antibodies in a patient group that included POAG and NTG patients; no correction was performed. However, if the Bonfferoni correction is applied to that study, the adjusted p value would be greater than 0.05.

It thus remains unclear whether the HLA class II genes have an effect on the pathogenesis of glaucoma. Several clinical and experimental glaucoma studies have focused on the role of the immune system in glaucoma; results strongly suggest that activated immunity contributes to glaucomatous optic neuropathy [33]. Cartwright et al. [34] reported that 30% of NTG patients had one or more immune-related diseases compared to 8% of control patients with ocular hypertension. This suggests that the immune system plays a role in the development of NTG. It is likely that multiple immune-related genes are associated with NTG. For example, the toll-like receptor 4 gene encodes a protein that recognizes heat shock proteins and lipopolysaccharides and is critical for the innate and adaptive immune systems. We recently reported that polymorphisms in this gene affected the development of NTG [35]. However, there are no reports of a significant association between other immune-related genes and NTG. In the future, carefully designed studies will be essential to determine whether other immune-related genes contribute to the pathogenesis of NTG.

In conclusion, we found no significant associations between the genes HLA-DRB1 and HLA-DQB1 and NTG in the Japanese population. These genes do not appear to be associated with the pathogenesis of NTG. However, because there are differences in the distribution of alleles of HLA genes in different ethnic groups, the association between genes and NTG should also be assessed in other ethnic groups.

Acknowledgments

The study was approved by the ethics committee of Yokohama City University, Yamanashi University, Gifu University, Kobe University, Yamaguchi University, Kumamoto University, and Tokai University, and complied with the guidelines of the Declaration of Helsinki. The details of the study were carefully explained to all of the patients and control subjects before we obtained informed written consent for genetic screening. This work was supported by Grants-in-Aid from the Ministry of Education, Culture, Sports, Science, and Technology, Japan (No.16209052) and by a grant from the Ministry of Health, Labor, and Welfare, Japan.

References

  1. Quigley HA. Number of people with glaucoma worldwide. Br J Ophthalmol. 1996; 80:389-93. [PMID: 8695555]
  2. Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006; 90:262-7. [PMID: 16488940]
  3. Quigley HA. Open-angle glaucoma. N Engl J Med. 1993; 328:1097-106. [PMID: 8455668]
  4. Hitchings RA, Anderson SA. A comparative study of visual field defects seen in patients with low-tension glaucoma and chronic simple glaucoma. Br J Ophthalmol. 1983; 67:818-21. [PMID: 6671097]
  5. Hitchings RA. Low-tension glaucoma – its place in modern glaucoma practice. Br J Ophthalmol. 1992; 76:494-6. [PMID: 1390534]
  6. Werner EB. Normal-tension glaucoma. In: Ritch R, Shields MB, Krupin T, editors. The Glaucomas. 2nd ed. St. Louis: Mosby; 1996. p. 769–97.
  7. Iwase A, Suzuki Y, Araie M, Yamamoto T, Abe H, Shirato S, Kuwayama Y, Mishima HK, Shimizu H, Tomita G, Inoue Y, Kitazawa Y, Tajimi Study Group, Japan Glaucoma Society.. The prevalence of primary openangle glaucoma in Japanese: the Tajimi Study. Ophthalmology. 2004; 111:1641-8. [PMID: 15350316]
  8. Suzuki Y, Iwase A, Araie M, Yamamoto T, Abe H, Shirato S, Kuwayama Y, Mishima HK, Shimizu H, Tomita G, Inoue Y, Kitazawa Y. Tajimi Study Group. Risk factors for open-angle glaucoma in a Japanese population. The Tajimi Study. Ophthalmology. 2006; 113:1613-7. [PMID: 16828504]
  9. Broadway DC, Drance SM. Glaucoma and vasospasm. Br J Ophthalmol. 1998; 82:862-70. [PMID: 9828767]
  10. Fontana L, Poinoosawmy D, Bunce CV, O'Brien C, Hitchings RA. Pulsatile ocular blood flow investigation in asymmetric normal tension glaucoma and normal subjects. Br J Ophthalmol. 1998; 82:731-6. [PMID: 9924361]
  11. Ishida K, Yamamoto T, Kitazawa Y. Clinical factors associated with progression of normal-tension glaucoma. J Glaucoma. 1998; 7:372-7. [PMID: 9871858]
  12. Drance S, Anderson DR, Schulzer M, Collaborative Normal-Tension Glaucoma Study Group.. Risk factors for progression of visual field abnormalities in normal-tension glaucoma. Am J Ophthalmol. 2001; 131:699-708. [PMID: 11384564]
  13. Wax MB. Is there role for the immune system in glaucomatous optic neuropathy? Curr Opin Ophthalmol. 2000; 11:145-50. [PMID: 10848222]
  14. Wax MB, Tezel G, Saito I, Gupta RS, Harley JB, Li Z, Romano C. Anti-Ro/SS-A positivity and heat shock protein antibodies in patients with normal-pressure glaucoma (see comments). Am J Ophthalmol. 1998; 125:145-57. [PMID: 9467439]
  15. Maruyama I, Ohguro H, Ikeda Y. Retinal ganglion cells recognized by serum autoantibody against gamma-enolase found in glaucoma patients. Invest Ophthalmol Vis Sci. 2000; 41:1657-65. [PMID: 10845582]
  16. Yang J, Tezel G, Patil RV, Romano C, Wax MB. Serum autoantibody against glutathione S-transferase in patients with glaucoma. Invest Ophthalmol Vis Sci. 2001; 42:1273-6. [PMID: 11328739]
  17. Kremmer S, Kreuzfelder E, Klein R, Bontke N, Henneberg-Quester KB, Steuhl KP, Grosse-Wilde H. Antiphosphatidylserine antibodies are elevated in normal tension glaucoma. Clin Exp Immunol. 2001; 125:211-5. [PMID: 11529911]
  18. Tezel G, Edward DP, Wax MB. Serum autoantibodies to optic nerve head glycosaminoglycans in patients with glaucoma. Arch Ophthalmol. 1999; 117:917-24. [PMID: 10408457]
  19. Gil-Carrasco F, Granados J, Barojas-Weber E, Gilbert-Lucido ME, Vargas-Alarcón G. Immunogenetic aspects in primary open-angle glaucoma in family members of Mexican mestizo glaucomatous patients. Am J Ophthalmol. 1994; 118:744-8. [PMID: 7977600]
  20. Gil-Carrasco F, Vargas-Alarcón G, Zúñiga J, Tinajero-Castañeda O, Hernández-Martinez B, Hernández-Pacheco G, Rodríguez-Reyna TS, Hesiquio R, Gamboa R, Granados J. HLA-DRB and HLA-DQB loci in the genetic susceptibility to develop glaucoma in Mexicans. Am J Ophthalmol. 1999; 128:297-300. [PMID: 10511023]
  21. Ferreri G, D'Andrea A, Castagna I, Rechichi C, Pettinato G, D'Andrea D. The role of class I and class II HLA antigens in primary open angle glaucoma (POAG). Acta Ophthalmol Scand Suppl. 1998; 227:17-9. [PMID: 9972328]
  22. Aviner Z, Henley WL, Fotino M, Leopold IH. Histocompatibility (HL-A) antigens and primary open-angle glaucoma. Tissue Antigens. 1976; 7:193-200. [PMID: 941127]
  23. Shin DH, Becker B. HLA-A11 and HLA-Bw35 and resistance to glaucoma in white patients with ocular hypertension. Arch Ophthalmol. 1977; 95:423-4. [PMID: 851487]
  24. Shin DH, Becker B, Waltman SR, Palmberg PF, Bell E, Jr. The prevalence of HLA-B12 and HLA-B7 antigens in primary open-angle glaucoma. Arch Ophthalmol. 1977; 95:224-5. [PMID: 836206]
  25. Ritch R, Podos SM, Henley W, Moss A, Southern AL, Fotino M. Lack of association of histocompatibility antigens with primary open-angle glaucoma. Arch Ophthalmol. 1978; 96:2204-6. [PMID: 718509]
  26. Kass MA, Palmberg P, Becker B, Miller JP. Histocompatibility antigens and primary open-angle glaucoma. A reassessment. Arch Ophthalmol. 1978; 96:2207-8. [PMID: 718510]
  27. David R, Maier G, Baumgarten I, Abrahams C. HLA antigens in glaucoma and ocular hypertension. Br J Ophthalmol. 1979; 63:293-6. [PMID: 465401]
  28. Shaw JF, Levene RZ, Sowell JG. The incidence of HLA antigens in black primary open-angle glaucoma patients. Am J Ophthalmol. 1978; 86:501-5. [PMID: 360847]
  29. Anderson DR, Patella VM. Automated static perimetry. 2rd edition. St.Louis: Mosby; 1999. p. 121–90.
  30. Lancaster A, Nelson MP, Meyer D, Thomson G, Single RM. PyPop: a software framework for population genomics: analyzing large-scale multi-locus genotype data. Pac Symp Biocomput. 2003; 8:514-25. [PMID: 12603054]
  31. Hedrick PW. Gametic disequilibrium measures: proceed with caution. Genetics. 1987; 117:331-41. [PMID: 3666445]
  32. Yano T, Yamada K, Kimura A, Takeshita T, Minohara M, Kira J, Senju S, Nishimura Y, Tanihara H. Autoimmunity against neurofilament protein and its possible with HLADRB1* 1502 allele in glaucoma. Immunol Lett. 2005; 100:164-9. [PMID: 16005081]
  33. Tezel G, Yang J, Wax MB. Heat shock proteins, immunity and glaucoma. Brain Res Bull. 2004; 62:473-80. [PMID: 15036560]
  34. Cartwright MJ, Grajewski AL, Friedberg ML, Anderson DR, Richards DW. Immunerelated disease and normal-tension glaucoma. A case-control study. Arch Ophthalmol. 1992; 110:500-2. [PMID: 1562256]
  35. Shibuya E, Meguro A, Ota M, Kashiwagi K, Mabuchi F, Iijima H, Kawase K, Yamamoto T, Nakamura M, Negi A, Sagara T, Nishida T, Inatani M, Tanihara H, Aihara M, Araie M, Fukuchi T, Abe H, Higashide T, Sugiyama K, Kanamoto T, Kiuchi Y, Iwase A, Ohno S, Inoko H, Mizuki N. Association of toll-like receptor 4 gene polymorphisms with normal tension glaucoma. Invest Ophthalmol Vis Sci. 2008; 49:4453-7. [PMID: 18586872]