Molecular Vision 2011; 17:2171-2176 <http://www.molvis.org/molvis/v17/a235>
Received 25 July 2011 | Accepted 8 August 2011 | Published 12 August 2011

Susceptibility to primary angle closure glaucoma in Saudi Arabia: the possible role of mitochondrial DNA ancestry informative haplogroups

Khaled K. Abu-Amero,1,3 Ana M. González,2 Essam A. Osman,1 José M. Larruga,2 Vicente M. Cabrera,2 Saleh A. Al-Obeidan1

1Department of Ophthalmology, College of Medicine, King Saud University, Riyadh, Saudi Arabia; 2Department of Genetics, Faculty of Biology, University of La Laguna, Tenerife, Spain; 3Department of Ophthalmology, College of Medicine, University of Florida, Jacksonville, FL

Correspondence to: Dr. Khaled K. Abu-Amero, Department of Ophthalmology, College of Medicine, King Saud University, P.O. Box 245, Riyadh 11411, Saudi Arabia; Phone: +96614786100; FAX: +96614775731; email: abuamero@gmail.com

Abstract

Purpose: In a previous preliminary analysis we reported that mitochondrial DNA (mtDNA) haplogroup R0a was significantly more frequent in primary angle closure glaucoma (PACG) Saudi patients than in healthy Saudi controls. This result prompted us to extend our work using a significant larger Saudi PACG cohort and more healthy controls.

Methods: We sequenced the mtDNA regulatory hypervariable region-I (HVS-I) and coding regions, comprising haplogroup diagnostic polymorphisms, in 227 PACG Saudi patients and compared their haplogroup frequencies with those obtained from 186 matched healthy controls (free of PACG by examination) and from a large sample of 810 healthy Saudi Arabs representing the general Saudi population.

Results: MtDNA Haplogroups R0a and J, the most abundant lineages in Saudi Arabia, were in significant higher frequencies in the PACG patients than in controls, while the widespread western Eurasian haplogroup U was associated with reduced risk to developing PACG.

Conclusions: Haplogroups R0a and J could be ancestry informative markers for PACG in the Saudi Arabian population. In addition, the western Eurasian haplogroup U may play a mild protective effect to this illness.

Introduction

Although glaucoma embodies a heterogeneous group of optic neuropathies, all types are defined by progressive and irreversible degeneration of the optic nerve with gradual visual field loss. Within this group, primary angle-closure glaucoma (PACG) is a type of glaucoma characterized by a narrow iridocorneal angle resulting in a blockage of the aqueous outflow structures. It seems that there is an anatomic and physiologic predisposition to PACG. Thus, shallow anterior chamber depth is considered an anatomically inheritable risk factor for this illness [1]. Congruently, Eastern Asian ascendance, hyperopia and female gender significantly predispose to this disease [2]. It has been estimated that in Saudi Arabia, 40% of glaucoma patients belong to the PACG type [3]. This percentage is closer to those calculated for Asian than for European populations [4]. Although a hereditary component for PACG exists, causative genes have not yet been identified and only occasional differences in gene frequency polymorphisms or in levels of gene expression between PACG patients and healthy controls have been reported. For instance, mutations in the myocilin (MYOC) gene were observed in PACG Chinese patients but also in Chinese controls ruling out a role of this gene in the pathogenesis of PACG [5]. The C677T polymorphism in the methylenetetrahydrofolate reductase (MTHFR) gene was found associated with PACG in a Pakistani cohort [6], but this result has not been extended to other populations. The possible role of mitochondrial abnormalities (mtDNA mutations and respiratory activity) in PACG was investigated and authors found minimal evidence of such abnormalities [7]. A significant increase of SPARC (secreted protein, acidic, and rich in cysteine) in the iris of PACG patients compared to healthy controls has also been documented pointing to a possible role of this gene in the development of PACG [8]. However, these results have not yet been confirmed. PACG is also characterized by an increase of intraocular pressure (IOP). Recently it has been demonstrated experimentally that elevated hydrostatic pressure induces mitochondrial structural changes and apoptotic cell death which suggests, once more, a possible role of mitochondrial injury in glaucoma [9]. Significant associations between mitochondrial DNA (mtDNA) haplogroups and different types of glaucoma have been detected in the Saudi Arabia population [3,10]. Furthermore, in a preliminary pilot experiment, a highly significant prevalence of mtDNA haplogroup preHV1 (now named R0a) in a small cohort of 29 PACG Saudi patients was found [11]. The aim of the present study was to evaluate the possible role of mtDNA Haplogroup(s) in large cohort of PACG Saudi patients with special emphasis to confirm the susceptibility to PACG for carriers of R0a, the most representative mtDNA haplogroup in the Arabian peninsula [12].

Methods

Patients and control subjects

We recruited 227 Saudi PACG patients who satisfied strict clinical criteria for PACG. Patients were eligible for inclusion if they had evidence of glaucomatous optic nerve damage attributable to primary angle closure [13]. Inclusion criteria included at least three of the following: 1) clinical documentation of angle closure, defined as the presence of appositional or synnechial closure of the anterior chamber angle involving at least 270 degrees by gonioscopy in either eye; 2) intraocular pressure elevated to a level ≥21 mmHg measured by Goldmann applanation tonometry before or after treatment; 3) evidence of characteristic glaucomatous optic disk damage with excavation of the disc causing a cup-to-disk ratio (c/d) vertically of at least 0.70 in at least one eye; and 4) characteristic peripheral visual field loss including nerve fiber bundle defects (nasal step, arcuate scotoma, paracentral scotoma) or advanced visual field loss (central and/or temporal island of vision) as tested by Humphrey Field Analyzer in those patients with vision better than 20/200 or Goldmann Manual Perimetry in those with worse vision. Exclusion criteria included: 1) secondary angle closure glaucoma; 2) presence of pseudoexfoliation syndrome even if coexistent with angle closure; 3) another cause of optic nerve injury affecting either eye; 4) significant visual loss in both eyes not associated with glaucoma; 5) inability to visualize the fundus for optic disk assessment; or 6) refusal to participate. Patients were recruited from the glaucoma clinic at King Abdulaziz University Hospital (KAUH) after signing an informed consent approved by the institutional review board (proposal number # 08–657).

A second group (n=186) of healthy Saudi Arabs controls (HMC group) free from glaucoma by examination were recruited. Entry criteria for those subjects were age >40, normal IOP, open angles on gonioscopy, and normal optic nerves on examination.

A third large sample (n=810) of healthy Saudi Arabs (HSA), representing all five major Saudi Arabian provinces, that were recruited previously for population genetics studies, was used for statistical comparisons. All individuals of this group were Saudi Arabs who reported no symptomatic metabolic, genetic, ocular disorders or any ophthalmic problem on an extensive questionnaire regarding family history, past medical problems, and current health. For those controls, their mother’s ancestral origin was established as a Saudi. The mtDNA haplogroup assignation for 552 subjects of this sample has been already published [11]. All patients and controls were maternally unrelated Saudi Arabs, all whose known ancestors were of Saudi Arabian origin. This research followed the tenets of the Declaration of Helsinki.

DNA extraction

Five milliters of peripheral blood were collected in EDTA tubes from all participating individuals. DNA was extracted using the illustra blood genomicPrep Mini Spin Kit from GE Healthcare (Buckinhamshire, UK), and stored at −20 °C in aliquots until required.

Mitochondrial haplogroup assortment

All samples were amplified and sequenced for the mtDNA regulatory region hypervariable segment I (HVS I) I using previously described primers and conditions [14] (Table 1). Briefly, each 25 µl PCR reaction contained 2.5 µl of 10× reaction buffer with MgCl2 (Amersham Pharmacia Biotech, Piscataway, NJ), 10 ρmol of each primer, 100 ρmol/µl each of deoxynucleoside triphophates (deoxyadenosine triphosphate, deoxyguanosine triphosphate, deoxycytidine triphosphate and deoxythymidine triphosphate; Perkin-Elmer Corporation, Foster City, CA) in Tris HCl buffer, 1 unit Taq DNA polymerase (Amersham Pharmacia Biotech, Piscataway, NJ) and 50 ng genomic DNA template. The mixture was denatured at 95 °C for 5 min and the PCR reaction was carried out for 35 cycles, in a GeneAmp 9700 PCR system (Perkin-Elmer Corporation), under the following conditions: denaturation at 95 °C for 1 min, annealing at 54 °C for 45 s, extension at 72 °C for 1 min, and final extension cycle of 72 °C for 7 min. The PCR products were electrophoresed on a 1% agarose gel and detected with 0.5 µg/ml ethidium bromide to confirm the correct amplicon size. Haplotypes were tentatively assorted into haplogroups according to their HVS diagnostic positions. This assortment was further confirmed, when necessary, by the analysis of coding-region diagnostic haplogroup polymorphisms (Table 2), following the most recent mtDNA haplogroup nomenclature [15]. To detect these polymorphisms a fragment spanning the diagnostic position was amplified and sequenced using any of the 32 overlapping pairs of primers that cover the whole mtDNA genome, and the PCR and sequencing conditions previously published for each of them [14].

Statistical analysis

Differences in haplogroup frequency between patient and control groups were tested by 2×2 contingence tests or by Fisher's exact tests when appropriate. To keep statistically adequate haplogroup sub-sample sizes within cohorts, when necessary, mtDNA haplogroups were collapsed into larger haplogroup identities following phylogenetic criteria [15]. Genetic diversity, h, within haplogroups and its variance were estimated from the unbiased formulas given by Nei [16], as implemented in the Arlequin 3.5 package. Differences in h were performed using a standard two-tailed z test. Finally the geographic origin distributions of patients and matched controls and differences in haplogroup frequencies among Saudi Arabian provinces were compared by 2×n contingency tests.

Results

As women have been reported to be at higher risk for PACG than men [17], and as gender was not a selective criteria in the recruitment of patients and controls, we compared the sex-ratio between both cohorts. Among patients, the female to male ratio (1.34:1) was significantly higher (p=0.047) than in controls (0.94), pointing to a greater female susceptibility to PACG in the Saudi population. However, mtDNA haplogroup frequencies did not differ between sexes in any group. Haplogroup frequencies obtained for the PACG patients, healthy matched controls (HMC) and the general Saudi population (HSA) is listed in Table 2. Although mtDNA haplogroup R0a was abundant in PACG (22.5%) than in HMS (14.0%) and HSA (17.4%), only the PACG*HMS comparison reached statistical significance (Table 2). Additionally, Haplogroup J, which is the most representative of the whole Saudi Arabia population [12], was found in higher frequencies in the patients group (28.2%) than in HMS (25.3%) and HSA (21.1%) controls. In this case only the PACG*HSA comparison reached statistical significance (p=0.02). Finally, haplogroup U was significantly underrepresented in PACG patients (6.2%) compared to both the HMS (13.4%) and HSA (11.2%) control groups (Table 2) pointing to a consistent protective effect of this haplogroup to PACG.

Discussion

By studying the prevalence of various mtDNA hpalogroups in 227 PACG Saudi patients and compare this to a large cohort of population controls (n=810) and to a group of individuals with no PACG by examination (healthy controls, n=186), we found that Haplogroups R0a and J could be ancestry informative markers for PACG in the Saudi Arabian population. In our previous study [11], we also found an association (p=0.00002) between mtDNA haplogroup R0a and PACG in a smaller cohort of Saudi Patients (n=29). The level of significance, for the R0a, found here (p=0.03) is lower than in our previous study [11]. The difference observed in the level of statistical significance could be explained, at first sight, by an ascertainment bias for the previous PACG sampling due to its small sample size (n=29). An associated consequence we could expect of it would be a concomitant diminution in its haplotype diversity (h) within haplogroup R0a. However, h values (Table 3) for the past and present PACG samples are very similar and not significantly different from controls (p=0.38). In addition, our PACG sample seems to be as diverse as controls for the J (p=0.20) and U (p=0.76) for which significant differences between the patient and one or both control cohorts were found (Table 2). Population genetic substructure in Saudi Arabia could also be a possible confounding factor. In fact, it was found in previous population studies that the mtDNA haplogroup frequency distribution in Arabia was geographically structured [12]. Focusing into the haplogroups that showed associations with PACG, there are significant differences in their geographic frequency distributions, across the five main provinces of Saudi Arabia, for R0a (p=0.04), with a frequency range from western (8%) to center (22%), and for J (p=0.002), with a frequency range from western (37%) to northern (16%). Only haplogroup U seems to have a uniform distribution across the country (p=0.63). However, the geographic origin distribution of the recruited PACG patients and HMS matched controls (Table 4) was strikingly similar (p=0.42) discarding wide inter-regional genetic differentiation as a main cause of the PACG and mtDNA haplogroup associations found. With the information at hand, we would attribute the R0a frequency differences found between the two independent PACG samples to biased sampling of the previous one due to its small sample size. However, other demographic population parameters, not yet uncovered, as fine-scale genetic structure and/or hidden social organization should not be discarded.

In conclusion, we have detected a significant susceptibility of Saudi R0a carriers to PACG although at a milder level than the one previously reported. It also seems likely that haplogroup J, abundanant among the Saudi population, also confers a mild susceptibility to PACG. On the contrary, the western Eurasian widespread haplogroup U seems to play a mild protective effect to PACG. In a recent analysis on primary open-angle glaucoma (POAG) in Saudi Arabia [3], it was found that mtDNA lineages of sub-Saharan Africa origin seemed to confer susceptibility to POAG, pointing to a role of ancestry informative marker for macro-haplogroup L, as it is well documented that the prevalence of POAG is higher in subjects of African ascendance [18]. In a similar vein, haplogroups R0a and J could reflect a predisposition of the genuine Saudi Arabian population to suffer PACG, which is supported by the high prevalence of this illness in this Country, near of the range observed for South and East Asians, and significantly higher than the one observed for populations of European ascendance [3].

Acknowledgments

This work was mainly supported by the Glaucoma research chair at the Department of Ophthalmology, College of Medicine, King Saud University, Riyadh, Saudi Arabia. V.C.M., J.M.L., and A.M.G. are supported by grant CGL2010–16195 from the Spanish Ministerio de Ciencia y Tecnología, Spain.

References

  1. Alsbirk PH. Anterior chamber depth and primary angle-closure glacuoma. II. A genetic study. Acta Ophthalmol (Copenh). 1975; 53:436-9. [PMID: 1174403]
  2. Nongpiur ME, Ku JY, Aung T. Angle closure glaucoma: a mechanistic review. Curr Opin Ophthalmol. 2011; 22:96-101. [PMID: 21252671]
  3. Abu-Amero KK, Gonzalez AM, Osman EA, Larruga JM, Cabrera VM, Al-Obeidan SA. Mitochondrial DNA lineages of African origin confer susceptibility to primary open-angle glaucoma in Saudi patients. Mol Vis. 2011; 17:1468-72. [PMID: 21677789]
  4. Minckler PF, Hung PT. Angle closure glaucoma-Classification and Racial variation. Asian Journal of Ophthalmology. 2001; 3:3-9.
  5. Aung T, Yong VH, Chew PT, Seah SK, Gazzard G, Foster PJ, Vithana EN. Molecular analysis of the myocilin gene in Chinese subjects with chronic primary-angle closure glaucoma. Invest Ophthalmol Vis Sci. 2005; 46:1303-6. [PMID: 15790895]
  6. 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]
  7. 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]
  8. Chua J, Seet LF, Jiang Y, Su R, Htoon HM, Charlton A, Aung T, Wong TT. Increased SPARC expression in primary angle closure glaucoma iris. Mol Vis. 2008; 14:1886-92. [PMID: 18949063]
  9. Ju WK, Kim KY, Lindsey JD, Angert M, Patel A, Scott RT, Liu Q, Crowston JG, Ellisman MH, Perkins GA, Weinreb RN. Elevated hydrostatic pressure triggers release of OPA1 and cytochrome C, and induces apoptotic cell death in differentiated RGC-5 cells. Mol Vis. 2009; 15:120-34. [PMID: 19169378]
  10. Abu-Amero KK, Cabrera VM, Larruga JM, Osman EA, Gonzalez AM, Al-Obeidan SA. Eurasian and Sub-Saharan African mitochondrial DNA haplogroup influences pseudoexfoliation glaucoma development in Saudi patients. Mol Vis. 2011; 17:543-7. [PMID: 21364909]
  11. Abu-Amero KK, Morales J, Bosley TM, Mohamed GH, Cabrera VM. The role of mitochondrial haplogroups in glaucoma: a study in an Arab population. Mol Vis. 2008; 14:518-22. [PMID: 18385785]
  12. Abu-Amero KK, Larruga JM, Cabrera VM, Gonzalez AM. Mitochondrial DNA structure in the Arabian Peninsula. BMC Evol Biol. 2008; 8:45 [PMID: 18269758]
  13. Foster PJ, Buhrmann R, Quigley HA, Johnson GJ. The definition and classification of glaucoma in prevalence surveys. Br J Ophthalmol. 2002; 86:238-42. [PMID: 11815354]
  14. Maca-Meyer N, Gonzalez AM, Larruga JM, Flores C, Cabrera VM. Major genomic mitochondrial lineages delineate early human expansions. BMC Genet. 2001; 2:13 [PMID: 11553319]
  15. van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009; 30:E386-94. [PMID: 18853457]
  16. Nei N. Molecular evolutionary genetics. New York: Columbia University Press; 1987.
  17. Congdon N, Wang F, Tielsch JM. Issues in the epidemiology and population-based screening of primary angle-closure glaucoma. Surv Ophthalmol. 1992; 36:411-23. [PMID: 1589856]
  18. Libby RT, Gould DB, Anderson MG, John SW. Complex genetics of glaucoma susceptibility. Annu Rev Genomics Hum Genet. 2005; 6:15-44. [PMID: 16124852]