Molecular Vision 2011; 17:2248-2254 <http://www.molvis.org/molvis/v17/a244>
Received 30 May 2011 | Accepted 15 August 2011 | Published 19 August 2011

The association of hepatocyte growth factor (HGF) gene with primary angle closure glaucoma in the Nepalese population

Mona S. Awadalla,1 Suman S. Thapa,2 Kathryn P. Burdon,1 Alex W. Hewitt,1,3 Jamie E. Craig1

1Department of Ophthalmology, Flinders University, Flinders Medical Centre Adelaide, South Australia, Australia; 2Nepal Glaucoma Eye Clinics, Tilganga Institute of Ophthalmology, Kathmandu, Nepal; 3Centre for Eye Research Australia, University of Melbourne, Royal Victorian Eye and Ear Hospital, Melbourne, Australia

Correspondence to: Dr. Mona S. Awadalla, Ophthalmology Department, Flinders University, GPO Box 2100, Adelaide South Australia, Australia 5000; Phone: +61 8 8204 6985; FAX: +61 8 8277 0899; email: awad0002@flinders.edu.au

Abstract

Purpose: Genetic variation in the hepatocyte growth factor (HGF) gene has recently been associated with hyperopia, which is a known risk factor for primary angle closure glaucoma (PACG). This study aimed to investigate whether genetic variation in HGF is associated with primary angle closure glaucoma in the Nepalese population.

Methods: One hundred six Nepalese patients with primary angle closure glaucoma and 204 matched controls were recruited. Twelve tag single nucleotide polymorphisms (SNPs) were selected and genotyped to cover the majority of common variation within HGF. Genotype and haplotype analyses were conducted in PLINK.

Results: Four HGF SNPs were found to be significantly associated with PACG, rs5745718, rs12536657, rs12540393 and rs17427817 (p=0.002, 0.002, 0.0006, and 0.0006, respectively). In addition, haplotype analysis showed one common haplotype to be significantly associated with PACG (p=0.001) in this population.

Conclusions: Genetic variation in HGF is associated with PACG in the Nepalese population. Additional replication studies in other populations are necessary to confirm this association and to further explore the role of HGF in the pathogenesis of this blinding disease.

Introduction

Glaucoma represents a group of diseases with the common feature of slowly progressive destruction of the optic nerve with corresponding loss of the peripheral visual field [1]. Glaucoma is second only to cataract in causing blindness worldwide [2]. Importantly for patients with glaucoma, blindness is reported to be up to 25% higher in people with primary angle closure glaucoma (PACG), than open angle glaucoma worldwide [3].

PACG patients have been found to have particular anatomic biometric features including shallow anterior chambers [4], lens thickness and position [5], narrow iridiotrabecular drainage angles, short axial lengths [6], and hyperopic refractive error [7]. The most important risk factor is shallow anterior chamber depth [8], which has been found to correlate with older age, gender (commoner in females), and race (shallower in Eskimos and Asians than Caucasians and Africans) [9]. Asian populations are at higher risk of developing PACG than other groups [10], and the majority of bilaterally blind glaucoma patients live in China [11].

Amerasinghe et al. [12], found that siblings of Chinese patients with PACG have almost a 50% probability of having narrow angles. In another Chinese study, first degree relatives were also found to have 6–9 fold increased risk of developing ACG [13]. These studies suggest a genetic component to the risk of PACG. Several candidate genes have been studied in relation to PACG. The matrix-metalloproteinase-9 gene (MMP9), which is involved in scleral extra-cellular matrix (ECM) remodelling, was shown to be associated with PACG in Taiwanese patients [14], and in australian patients [15], however several studies failed to replicate this association [16,17]. The Methylenetetrahydrofolate reductase (MHTFR) gene was also proposed to play a role in ECM remodelling of the anterior chamber in a PACG Pakistani cohort, secondary to C677T and A1298C polymorphisms [18,19]. The Membrane frizzled-related protein (MFRP) gene, which causes recessive nanophthalmos, has been studied in relation to PACG due to some similarities in the phenotypes, however no association was detected [20,21]. The inconsistencies in genetic findings between studies suggests that further evaluation of PACG genetics is warranted.

A recent study investigated the role of hepatocyte growth factor (HGF) in causing disruption of the emmetropization process within the eye [22]. The authors identified two single nucleotide polymorphisms (SNPs; rs12536657 and rs5745718) within HGF that were significantly associated with hyperopia. Since both angle closure glaucoma and hyperopia share the same feature of short axial length [23], we hypothesized that this gene may be involved in the development of PACG. The aim of our study was to investigate the association between tag SNPs of the HGF gene and primary angle closure glaucoma in the Nepalese population.

Methods

Participants were recruited from the Nepal Glaucoma Eye Clinic, Tilganga Institute of Ophthalmology, Kathmandu, Nepal. Ethics approval was approved by the Institutional Review Committee of the Tilganga Institute of Ophthalmology (TIO), and is being conducted in accordance with the Declaration of Helsinki and its subsequent revisions. Informed consent was obtained from each individual. In total, 106 PACG cases, and 204 controls were recruited. Cases and controls were matched for sex and age although controls were slightly older than cases by design for this aging disease. All participants were from Nepal [24,25].

Each participant underwent a complete eye examination including; slit lamp examination of the anterior chamber, gonioscopy, best corrected visual acuity, measurement of intraocular pressure, fundus examination with special attention to optic disc parameters, and visual field assessment. Objective refraction was performed using a streak retinoscope (Beta 200, Heine, Germany), which was followed by a subjective refraction [24]. The diagnosis of PACG was based on the presence of glaucomatous optic neuropathy with cup:disc ratio ≥0.7, intraocular pressure more than 21 mmHg, peripheral visual loss, presence of at least 180 degrees of closed angle in which the trabecular meshwork is not visible on gonioscopy, which follow the International Society of Geographical and Epidemiological Ophthalmology (ISGEO) classification as described by Foster and colleagues [26].Controls were required to have none of the above characteristics, with no family history of glaucoma or previous glaucomatous operations. Participants with pseudophakia or secondary angle closure glaucoma caused by events such as uveitis, trauma or lens subluxation were excluded.

Genomic DNA was extracted from 2 ml of venous blood using the QiaAmp Blood Midi Kit (Qiagen, Valencia, CA). The two HGF SNPs (rs12536657, and rs5745718) identified by Veerappan et al. [22], as well as 10 other tag SNPs, were selected using the tagger program implemented in Haploview 4.2. SNPs were selected from the HapMap Han Chinese in Beijing, China (CHB) sample as the most closely related population available at the time of the study. Tag SNPs were chosen using pairwise tagging, to have an r2>0.8 with SNPs displaying a minor allele frequency of 5% in this population. SNPs previously reported to be associated with hyperopia were force included in the selection of tags. The 12 tag SNPs included: rs5745752 (A/G), rs5745718 (A/C), rs12536657 (A/G), rs2286194 (A/T), rs5745692 (C/G), rs12540393 (C/T), rs17427817 (C/G), rs12707453 (G/A), rs5745616 (A/G), rs3735520 (T/C), rs6942495 (G/C), and rs17501080 (C/A). A Bonferroni corrected p-value of 0.004 (0.05/12) was considered statistically significant.

Genotyping was conducted at the Australian Genome Research Facility, Brisbane, Australia, using the iPLEX Gold chemistry (Sequenom Inc., San Diego, CA) on an Autoflex mass spectrometer (Sequenom Inc.).

Differences in age and gender between cases and controls were assessed by t-test and χ2 test, respectively. All genetic analyses were conducted in PLINK [27]. SNPs were assessed for compliance with Hardy–Weinberg equilibrium using a χ2 test. Genetic association was assessed under allelic (allele 1 versus allele 2), dominant (1/1 genotype versus 1/2+2/2) and recessive (1/1 + 1/2 genotypes versus 2/2) models. Where fewer than 5 counts for a given genotype were observed, Fisher’s exact test was used, otherwise a χ2 test was used. Haplotypes across the two observed linkage disequilibrium blocks, as visualized in Haploview using the “solid spine” block definition, was also analyzed for association in PLINK.

Results

Samples from 310 Nepalese individuals (comprising 106 cases and 204 matched controls) were genotyped. Demographic characteristics were similar between cases and controls (Table 1). Neither age nor gender differed significantly between cases and controls. The mean spherical equivalent for PACG cases was −0.15 diopters (D)±1.46, and in controls 0.09 D±0.31 (p=0.16). The ratio of females to males in both cases and controls was approximately 3:1.

All SNPs conformed to Hardy–Weinberg equilibrium in both cases and controls (p>0.05), and the call rate was >98%. The genotype counts and frequencies of all SNPs are shown in Table 2, along with the allelic association results. All 12 tag SNPs were located in the introns of HGF, with the physical location presented in Figure 1. Four SNPs reached statistically significant association with PACG; rs5745718 (p=0.002), rs12536657 (p=0.002), rs12540393 (p=0.0006), and rs17427817 (p=0.0006). After controlling for spherical equivalent via multivariate analyses, only three SNPs remained significantly associated with PACG; rs5745718 (p=0.003), rs12540393 (p=0.001), and rs17427817 (p=0.001; Table 2). Of the 12 SNPs, rs12540393 and rs17427817, were significantly associated under the dominant genetic model (p=0.001; Table 3). SNPs rs5745718 and rs12536657 also show nominal associations but did not survive correction for multiple testing.

Haplotypic associations with PACG were also investigated. Two haplotype blocks were identified under the “solid spine” block definition, as displayed in Figure 2. Block 1 was defined by all SNPs between rs5745752 and rs5745616 (overall p-value=0.004), and block 2 as SNPs rs3735520 to rs17501080 (overall p-value of 0.024). In block 1 the frequency of the GAATGCCAG haplotype was significantly greater in cases than in controls (18.7% versus 9.5%, respectively, p=0.001) and remained significant after Bonferroni correction for the 6 haplotypes observed (p=0.006). Additionally, in block 2 the frequency of the CCA haplotype was found to be higher in cases than controls (15.2% versus 8.8%, respectively, p=0.017), but this difference was not statistically significant following Bonferronni correction (p=0.068; Table 4).

Discussion

HGF protein has been found to play an important role in stimulating the growth and migration of various eye tissues including the corneal epithelium and endothelium, iris, retinal pigment epithelium, lens epithelium, and trabecular meshwork [28-31]. Hu and Ritch [32], revealed that the concentration of HGF in the aqueous humor was significantly higher in glaucomatous eyes than in cataract eyes, with no difference between open angle and angle closure glaucoma. Another study showed that HGF mRNA concentration was increased in rabbit lacrimal gland and corneal epithelium after injury of the cornea [33]. This supports the finding of Hu and Ritch [32] that increased HGF concentration in aqueous humor of glaucomatous eyes possibly reflects the functional effects of HGF on enhancement of aqueous flow and attempt to repair trabecular injury, rather than directly causing glaucoma.

The allelic association results displayed in Table 2 show four SNPs to be significantly associated with PACG (rs5745718, rs12536657, rs12540393, and rs17427817). The former two SNPs have recently been reported to be associated with hyperopia in an Australian Caucasian population [22]. Interestingly, the risk alleles for hyperopia, rs5745718(A) and rs12536657(A), in their study were the same as in our PACG study. This finding indicates a possible common pathway, or similarities between hyperopia and PACG which are known to share similar biometric features including short axial length [22,34]. However, all four SNPs in our PACG study were significantly associated with PACG independent of spherical equivalent, so the association appears to occur through a mechanism beyond an indirect association with hyperopia in this population.

The GAATGCCAG haplotype of block 1 showed significant association with PACG (p=0.001). This haplotype contains the associated risk allele of each of the significantly associated SNPs; rs5745718(A), rs12536657(A), rs12540393(C), and rs17427817(C). The frequency of this haplotype was also found to be higher in the PACG patients (18%) than in healthy controls (9%), reflecting the allele frequencies of these four SNPs. These alleles do not occur on any other common haplotype in this population. The association of the haplotype in block 2 is likely due to the linkage disequilibrium between the associated SNPs in block 1 and SNPs in block 2.

The role of HGF in PACG remains unknown. Further work is needed to determine its involvement in the pathogenesis of this blinding disease. The similar findings in hyperopia indicate that it may be involved in influencing the structure of the eye and thus predisposing those with short axial length to the risk of angle closure. It is unlikely that the tag SNPs assayed here are the functional variants. All four tag SNPs are located in introns of HGF and are likely to be in linkage disequilibrium with actual functional variants. The causative variant will likely be found on the background of the GAATGCCAG haplotype.

The main limitation of this study is the use of CHB to tag the Nepalese samples, however we found that the allele frequencies of the SNPs were remarkably similar between the two groups and thus they are likely to be quite similar.

To our knowledge, this is the first report to identify an association between HGF and PACG. Additionally the underlying genetic etiology of PACG in people of Nepalese descent has not been previously studied. In conclusion, this study revealed an association between the HGF gene and PACG. Future replication studies in different populations are necessary to confirm this association and to further explore the role of HGF in the pathogenesis of the disease.

Acknowledgments

This work was funded by a grant from the Flinders Medical Centre Foundation. K.P.B. is funded by a National Health and Medical Research Council (NHMRC) of Australia Career Development Award and J.E.C. is an NHMRC Practitioner Fellow.

References

  1. Quigley HA, Broman AT. the number of people with glaucoma world wide in 2010 and 2020. Br J Ophthalmol. 2006; 90:262-7. [PMID: 16488940]
  2. Frick KD, Foster A. The magnitude and cost of global blindness:an increasing problem that can be alleviated. Am J Ophthalmol. 2003; 135:471-6. [PMID: 12654362]
  3. Foster PJ, Oen FT, Machin D, Ng TP, Devereux JG, Johnson GJ, Khaw PT, Seah SK. The prevalence of glaucoma in Chinese residents of Singapore: a cross-sectional population survey of the Tanjong Pagar district. Arch Ophthalmol. 2000; 118:1105-11. [PMID: 10922206]
  4. Lin YW, Wang TH, Hung PT. Biometric study of acute primary angle-closure glaucoma. J Formos Med Assoc. 1997; 96:908-12. [PMID: 9409125]
  5. Foster PJ, Baasanhu J, Alsbirk PH, Munkhbayar D, Uranchimeg D, Johnson GJ. Glaucoma in Mongolia. A population-based survey in Hovsgol province, northern Mongolia. Arch Ophthalmol. 1996; 114:1235-41. [PMID: 8859083]
  6. Abu-Amero KK, Morales J, Osman MN, Bosley TM. Nuclear and mitochondrial analysis of patients with primary angle-closure glaucoma. Invest Ophthalmol Vis Sci. 2007; 48:5591-6. [PMID: 18055808]
  7. Congdon NG, Quigley HA, Hung PT, Wang TH, Ho TC. Screening techniques for angle-closure glaucoma in rural Taiwan. Acta Ophthalmol Scand. 1996; 74:113-9. [PMID: 8739673]
  8. Alsbirk PH. Primary angle-closure glaucoma. Oculometry, epidemiology, and genetics in a high risk population. Acta Ophthalmol Suppl. 1976; 127:5-31. [PMID: 179266]
  9. Salmon JF. Predisposing factors for chronic angle-closure glaucoma. Prog Retin Eye Res. 1999; 18:121-32. [PMID: 9920501]
  10. Foster PJ, Johnson GJ. Glaucoma in China: how big is the problem? Br J Ophthalmol. 2001; 85:1277-82. [PMID: 11673287]
  11. Yip JL, Foster PJ. Ethnic differences in primary angle-closure glaucoma. Curr Opin Ophthalmol. 2006; 17:175-80. [PMID: 16552253]
  12. Amerasinghe N, Zhang J, Thalamuthu A, He M, Vithana EN, Viswanathan A, Wong TY, Foster PJ, Aung T. The heritability and sibling risk of angle closure in Asians. Ophthalmology. 118:480-5. [PMID: 21035870]
  13. Wang N, Wu H, Fan Z. Primary angle closure glaucoma in Chinese and Western populations. Chin Med J (Engl). 2002; 115:1706-15. [PMID: 12609093]
  14. Wang IJ, Chiang TH, Shih YF, Lu SC, Lin LL, Shieh JW, Wang TH, Samples JR, Hung PT. The association of single nucleotide polymorphisms in the MMP-9 genes with susceptibility to acute primary angle closure glaucoma in Taiwanese patients. Mol Vis. 2006; 12:1223-32. [PMID: 17110919]
  15. Awadalla MS, Burdon KP, Kuot A, Hewitt AW, Craig JE. Matrix metalloproteinase-9 genetic variation and primary angle closure glaucoma in a Caucasian population. Mol Vis. 2011; 17:1420-4. [PMID: 21655354]
  16. Aung T, Yong VH, Lim MC, Venkataraman D, Toh JY, Chew PT, Vithana EN. Lack of association between the rs2664538 polymorphism in the MMP-9 gene and primary angle closure glaucoma in Singaporean subjects. J Glaucoma. 2008; 17:257-8. [PMID: 18552608]
  17. Cong Y, Guo X, Liu X, Cao D, Jia X, Xiao X, Li S, Fang S, Zhang Q. Association of the single nucleotide polymorphisms in the extracellular matrix metalloprotease-9 gene with PACG in southern China. Mol Vis. 2009; 15:1412-7. [PMID: 19633731]
  18. Micheal S, Qamar R, Akhtar F, Khan MI, Khan WA, Ahmed A. MTHFR gene C677T and A1298C polymorphisms and homocysteine levels in primary open angle and primary closed angle glaucoma. Mol Vis. 2009; 15:2268-78. [PMID: 19936026]
  19. Michael S, Qamar R, Akhtar F, Khan WA, Ahmed A. C677T polymorphism in the methylenetetrahydrofolate reductase gene is associated with primary closed angle glaucoma. Mol Vis. 2008; 14:661-5. [PMID: 18385801]
  20. Aung T, Lim MC, Wong TT, Thalamuthu A, Yong VH, Venkataraman D, Venkatraman A, Chew PT, Vithana EN. Molecular analysis of CHX10 and MFRP in Chinese subjects with primary angle closure glaucoma and short axial length eyes. Mol Vis. 2008; 14:1313-8. [PMID: 18648522]
  21. Wang IJ, Lin S, Chiang TH, Chen ZT, Lin LL, Hung PT, Shih YF. The association of membrane frizzled-related protein (MFRP) gene with acute angle-closure glaucoma–a pilot study. Mol Vis. 2008; 14:1673-9. [PMID: 18781223]
  22. Veerappan S, Pertile KK, Islam AF, Schache M, Chen CY, Mitchell P, Dirani M, Baird PN. Role of the hepatocyte growth factor gene in refractive error. Ophthalmology. 2010; 117:239-45. [PMID: 20005573]
  23. Sihota R, Lakshmaiah NC, Agarwal HC, Pandey RM, Titiyal JS. Ocular parameters in the subgroups of angle closure glaucoma. Clin Experiment Ophthalmol. 2000; 28:253-8. [PMID: 11021552]
  24. Thapa SS, Rana PP, Twayana SN, Shrestha MK, Paudel I, Paudyal G, Gurung R, Ruit S, Hewitt AW, Craig JE, van Rens G. Rationale, methods and baseline demographics of the Bhaktapur Glaucoma Study. Clin Experiment Ophthalmol. 2011; 39:126-34. [PMID: 21040314]
  25. Thapa SS, Paudyal I, Khanal S, Paudel N, van Rens GH. Comparison of axial lengths in occludable angle and angle-closure glaucoma–the Bhaktapur Glaucoma Study. Optom Vis Sci. 2011; 88:150-4. [PMID: 21150676]
  26. Foster PJ, Buhrmann R, Quigley H, Johnson GJ. The definition and classification of glaucoma in prevalence surveys. Br J Ophthalmol. 2002; 86:238-42. [PMID: 11815354]
  27. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81:559-75. [PMID: 17701901]
  28. Wordinger RJ, Clark AF, Agarwal R, Lambert W, McNatt L, Wilson SE, Qu Z, Fung BK. Cultured human trabecular meshwork cells express functional growth factor receptors. Invest Ophthalmol Vis Sci. 1998; 39:1575-89. [PMID: 9699547]
  29. Li Q, Weng J, Mohan RR, Bennett GL, Schwall R, Wang ZF, Tabor K, Kim J, Hargrave S, Cuevas KH, Wilson SE. Hepatocyte growth factor and hepatocyte growth factor receptor in the lacrimal gland, tears, and cornea. Invest Ophthalmol Vis Sci. 1996; 37:727-39. [PMID: 8603858]
  30. He PM, He S, Garner JA, Ryan SJ, Hinton DR. Retinal pigment epithelial cells secrete and respond to hepatocyte growth factor. Biochem Biophys Res Commun. 1998; 249:253-7. [PMID: 9705867]
  31. Weng J, Liang Q, Mohan RR, Li Q, Wilson SE. Hepatocyte growth factor, keratinocyte growth factor, and other growth factor-receptor systems in the lens. Invest Ophthalmol Vis Sci. 1997; 38:1543-54. [PMID: 9224282]
  32. Hu DN, Ritch R. Hepatocyte growth factor is increased in the aqueous humor of glaucomatous eyes. J Glaucoma. 2001; 10:152-7. [PMID: 11442175]
  33. Wilson SE, Chen L, Mohan RR, Liang Q, Liu J. Expression of HGF, KGF, EGF and receptor messenger RNAs following corneal epithelial wounding. Exp Eye Res. 1999; 68:377-97. [PMID: 10192796]
  34. Young TL, Metlapally R, Shay AE. Complex trait genetics of refractive error. Arch Ophthalmol. 2007; 125:38-48. [PMID: 17210850]