Molecular Vision 2015; 21:1340-1344 <http://www.molvis.org/molvis/v21/1340>
Received 12 March 2015 | Accepted 28 December 2015 | Published 30 December 2015

Analysis of functional polymorphisms in apoptosis-related genes in primary open angle glaucoma

Ewald Lindner, Wilfried Glatz, Christoph Schwab, Yosuf El-Shabrawi, Georg Mossböck

Department of Ophthalmology, Medical University of Graz, Austria

Correspondence to: Ewald Lindner, Department of Ophthalmology, Medical University of Graz, Auenbruggerplatz 4, 8036 Graz, Austria; Phone: +43 316 385 80810; FAX: +43 316 385 13261; email: ewald.lindner@medunigraz.at

Abstract

Purpose: Glaucoma is a disease with high heritability in which the degradation of retinal ganglion cells occurs via apoptosis. Therefore, we investigated the role of four functional apoptosis-related gene variants (Akt1 rs1130233, Bax rs4645878, Fas rs223476, and FasL rs763110) in patients with primary open angle glaucoma.

Methods: 334 patients with primary open angle glaucoma and 334 controls were recruited for this case–control study. The main outcome measures were genotype distribution and allelic frequencies determined with PCR.

Results: After adjustment for multiple testing, no significant difference in either the genotype distribution or the allelic frequencies of any investigated gene variant was found.

Conclusions: Our findings indicate that the investigated gene polymorphisms are unlikely to be major risk factors for primary open angle glaucoma in Caucasian patients.

Introduction

Glaucoma is defined as a group of progressive optic neuropathies characterized by a steady loss of retinal ganglion cells and consequently visual field defects. The number of people with glaucoma is about 60 million worldwide and will exceed 110 million by 2040 [1]. The pathomechanisms that lead to retinal ganglion cell (RGC) death are not fully understood. High intraocular pressure (IOP) is the main risk factor, but in the case of normal tension glaucoma, glaucomatous nerve degeneration occurs although the IOP is within normal limits.

Thus, other factors, such as ocular blood flow [2], inflammation [3], and genetics, must be considered. A positive family history for glaucoma is known to increase individual risk substantially [4]. There is no clear Mendelian inheritance pattern, but some portion of glaucoma is caused by defects in distinct genes, for example, myocilin [5]. The genetics of the larger portion of glaucoma is more complex, and a vast number of gene variants have been implicated [6].

Apoptosis or programmed cell death is a complex process, which ultimately leads to cellular degradation [7]. In glaucoma, several studies emphasized the role of apoptotic factors, such as apoptosis stimulating fragment (FAS) [8], BCL2-associated X protein (BAX) [9], and v-akt murine thymoma viral oncogene homolog 1 (Akt1) [10]. The binding of FAS ligand (FASL) to FAS forms the death-inducing signaling complex, which triggers apoptosis. FAS ligand (FASL) deficiency protects RGCs from cell death [10]. AKT is a serine-threonine protein kinase, which affects the components of the intrinsic cell death machinery [11,12]. The BAX protein is a member of the BCL2 family and has proapoptotic properties by antagonizing BCL2 [13]. Recently, C*HSDGIC*, a novel cyclopeptide from the cyclization of pituitary adenylate cyclase-activating polypeptide, was shown to inhibit apoptosis in retinal ganglion cells by decreasing BAX and increasing BCL2 [14].

For this candidate gene analysis, we selected polymorphisms that had already been demonstrated to alter gene expression or the activity of apoptosis-related genes [15-18]. Identifying the associations of gene polymorphisms with primary open angle glaucoma (POAG) may enhance our understanding of the pathogenesis and may pave the way for future neuroprotective therapies. In addition, a more precise individual risk profile can provide assistance for screening and foreseeing the progression rate of the condition. Therefore, we investigated a hypothesized association between apoptosis-related functional gene polymorphisms and POAG.

Methods

In the present institutional, retrospective case-control study we investigated a total of 668 unrelated Caucasian subjects comprising 334 patients with POAG, and 334 control subjects. In each group 187 females and 147 males were involved. The mean age of patients with POAG was 74.1 ± 8.1 years and 73.2 ± 7.8 years in control subjects, respectively. All participants were of Caucasian origin, living in the same geographical area and were seen at the local Department of Ophthalmology, Medical University of Graz. The study was approved by the Institutional Review Board of the Medical University of Graz and followed the principles of the Declaration of Helsinki. The study adhered to the ARVO statement on human subjects. Prior to enrollment, written informed consent was obtained from all participants.

All patients underwent slit-lamp biomicroscopy, testing for best-corrected visual acuity, Goldmann applanation tonometry, gonioscopy, and standard automated perimetry (Interzeag Octopus 500, program G2) or, in cases of profoundly decreased visual acuity, Goldmann perimetry. In all patients, photographs of the optic discs were taken.

POAG was defined by the IOP before the initiation of a pressure-lowering therapy of at least 21 mmHg, an open anterior chamber angle, optic disc changes characteristic of glaucoma (notching, thinning of the neuroretinal rim, increased cup/disc ratio in relation to the optic disc size), visual field defects characteristic of glaucoma (inferior or superior arcuate scotoma, nasal step, paracentral scotoma), and the absence of conditions that lead to secondary glaucoma. A narrow angle was defined as Schaffer’s grade 2 or less and was excluded with gonioscopy.

The control group consisted of 334 unrelated patients with no morphological or functional damage indicative of primary or secondary open angle or angle closure glaucoma. The control subjects had been admitted to the department for cataract surgery. All participants were Caucasians from the same geographic area (southern Austria).

Genotype determination

Venous blood was collected in 5 ml EDTA tubes. DNA was isolated form peripheral lymphocytes using a QIAamp DNA blood mini-kit (Qiagen, Venlo, Netherlands) following the manufacturer’s protocol and stored at −20 °C. Genotype determination was performed using high-resolution melting curve analysis on the LightCycler® 480 PCR system (Hoffmann-La Roche, Basel, Switzerland). The samples were amplified in duplicate 20 µl reactions using the Light Cycler 480 High Resolution Melting Master kit (Roche Diagnostics, Wien, Austria) and analyzed on an LC480 instrument I (Roche Diagnostics GmbH, Mannheim, Germany). The final reaction mix contained 1X Master Mix, 3 mM MgCl2, 4 µM forward and reverse primer, and 50 ng of genomic DNA. Primer sequences are shown in Table 1. For PCR, the following cycling conditions were chosen: one cycle of 95 °C for 10 min followed by 45 cycles of 95 °C for 10 s, 60 °C for 15 s, and 72 °C for 20 s. The amplicons were then denaturated at 95 °C for 1 min, cooled down to 40 °C for 1 min, and then melted from 65 °C to 95 °C with 25 signal acquisitions per degree. To detect sequence variations, Gene Scanning Software v 1.5 (Roche Diagnostics) was used. Using the Auto Group, mode samples were automatically grouped because of their melting curves.

Statistical analysis

Descriptive statistics were used to calculate the frequencies and percentages of discrete variables. Continuous data are given as mean ± standard deviation (SD). Means were compared using the Mann–Whitney test. Proportions of groups were compared with a chi-square test. The criterion for statistical significance was p value less than 0.05. The Bonferroni correction was used to adjust for multiple testing. Hardy–Weinberg equilibrium was calculated using HW Diagnostics-Version 1.beta (Fox Chase Cancer Center). Statistical analysis was performed using the SPSS statistical package (SPSS, version 17.0, Chicago, IL). Power calculation was done using PS Power and Sample Size Calculation software version 2.1.30.

Results

Our study included 334 patients with POAG and an equal number of control subjects. In each group, 187 women and 147 men were included. The mean age of the patients with POAG and the control subjects was 74.1±8.1 years and 73.2±7.8 years, respectively.

The observed genotype distributions did not deviate from those predicted by Hardy–Weinberg equilibrium. Table 2 shows the genotype and minor allele frequencies of the investigated polymorphisms in the patients with POAG and the control subjects. Genotype AA of AKT1 rs1130233 was significantly more common in POAG, but this significance did not remain after the Bonferroni adjustment. Therefore, no significant difference in either the genotype or allelic frequencies of any investigated gene variant was found. Quantitative trait analysis with age at presentation, mean defect in visual fields, vertical cup/disc ratio, or intraocular pressure did not reveal any association (p>0.05).

The present study had a statistical power of 0.80 to detect an odds ratio of 1.86 at a significance level of 0.05.

Discussion

To the best of our knowledge, the present study is the first to investigate a hypothesized association between POAG and functional apoptosis-related gene polymorphisms. Genotypes of four functional apoptosis-related polymorphisms were analyzed in 334 patients with POAG and 334 controls matched for age and sex. After correction for multiple testing, the allelic frequencies and genotype distribution did not differ significantly between the groups.

AKT inhibits apoptosis by phosphorylation and, thus, inhibition of transcription factors of the forkhead family, which trigger programmed cell death by inducing FASL [19]. Furthermore, AKT1 influences nerve cell plasticity and development [20]. The influence of rs1130233 on the expression of AKT1 has been demonstrated in lymphoblasts [21]. Although gene variants of AKT1 have been shown to affect apoptosis in a cell model [18], we could not find an association with POAG.

Gene polymorphisms have been shown to affect nuclear protein binding affinity and thus the transcription of BAX [22]. Mice deficient in BAX were protected from retinal ganglion cell death suggesting a major role of BAX in glaucoma development [23]. However, we did not find an association of the BAX gene polymorphism rs4645878 with POAG, which is in line with the findings of a recent meta-analysis in chronic lymphatic leukemia [24]. Despite the functionality of this polymorphism, it was not shown to alter the predisposition of POAG or chronic lymphatic leukemia, two diseases in which apoptosis plays a pivotal role.

Binding of FAS to FASL leads to caspase 8 activation and is a key mechanism of normal tissue homeostasis [25]. Polymorphisms in FAS and FASL have been found to be associated with various inflammatory diseases [26-28] and cancer [29-31]. Moreover, the expression of FAS and FASL is altered in neurodegenerative disorders. The FAS/FASL pathway has already been demonstrated to play a critical role in mediating beta-amyloid-induced death of neurons in Alzheimer’s disease [32,33]. Interestingly, in a study that used enzyme-linked immunosorbent assay (ELISA), FASL levels were not significantly altered in patients with POAG [34], and with multiplex beads analysis, the FASL levels were below the limits of detection in aqueous humor [35]. However, other studies identified the upregulation of FAS/FASL in glaucoma [10,36,37]. Due to these conflicting results, the role of FAS/FASL in glaucoma remains elusive. Still, our findings do not support the pivotal influence of FAS/FASL on glaucoma pathogenesis.

The following potential limitations should be kept in mind when interpreting our results. First, only a small number of single nucleotide polymorphisms (SNPs) were investigated in the present study. We cannot rule out associations with other polymorphisms of the investigated genes. A tagged SNP approach would be a possible strategy here, but it was beyond the scope of this survey. All investigated SNPs had already been demonstrated to be functional. Second, as genetic polymorphisms have been shown to vary between populations, our findings do not necessarily apply to populations other than Caucasian. In conclusion, our findings suggest that the investigated gene polymorphisms are unlikely to be major risk factors for POAG in Caucasian patients.

Acknowledgments

The authors thank Mrs. Fischl, Mrs. Trummer, Mrs. Wachswender and Mrs. Kirchengast for their skilful technical assistance. No funding was received.

References

  1. Tham YC, Li X, Wong TY, Quigley HA, Aung T, Cheng CY. Global Prevalence of Glaucoma and Projections of Glaucoma Burden through 2040: A Systematic Review and Meta-Analysis. Ophthalmology. 2014; 121:2081-90. [PMID: 24974815]
  2. Flammer J, Orgul S, Costa VP, Orzalesi N, Krieglstein GK, Serra LM, Renard JP, Stefansson E. The impact of ocular blood flow in glaucoma. Prog Retin Eye Res. 2002; 21:359-93. [PMID: 12150988]
  3. Tezel G, Wax MB. The immune system and glaucoma. Curr Opin Ophthalmol. 2004; 15:80-4. [PMID: 15021215]
  4. Tielsch JM, Katz J, Sommer A, Quigley HA, Javitt JC. Family history and risk of primary open angle glaucoma. The Baltimore Eye Survey. Arch Ophthalmol. 1994; 112:69-73. [PMID: 8285897]
  5. Mimivati Z, Nurliza K, Marini M, Liza-Sharmini A. Identification of MYOC gene mutation and polymorphism in a large Malay family with juvenile-onset open angle glaucoma. Mol Vis. 2014; 20:714-23. [PMID: 24883016]
  6. Janssen SF, Gorgels TG, Ramdas WD, Klaver CC, van Duijn CM, Jansonius NM, Bergen AA. The vast complexity of primary open angle glaucoma: disease genes, risks, molecular mechanisms and pathobiology. Prog Retin Eye Res. 2013; 37:31-67. [PMID: 24055863]
  7. Hengartner MO. The biochemistry of apoptosis. Nature. 2000; 407:770-6. [PMID: 11048727]
  8. Deckwerth TL, Elliott JL, Knudson CM, Johnson EM, , Jr Snider WD, Korsmeyer SJ. BAX is required for neuronal death after trophic factor deprivation and during development. Neuron. 1996; 17:401-11. [PMID: 8816704]
  9. Dudek H, Datta SR, Franke TF, Birnbaum MJ, Yao R, Cooper GM, Segal RA, Kaplan DR, Greenberg ME. Regulation of neuronal survival by the serine-threonine protein kinase Akt. Science. 1997; 275:661-5. [PMID: 9005851]
  10. Gregory MS, Hackett CG, Abernathy EF, Lee KS, Saff RR, Hohlbaum AM, Moody KS, Hobson MW, Jones A, Kolovou P, Karray S, Giani A, John SW, Chen DF, Marshak-Rothstein A, Ksander BR. Opposing roles for membrane bound and soluble Fas ligand in glaucoma-associated retinal ganglion cell death. PLoS ONE. 2011; 6:e17659 [PMID: 21479271]
  11. Yao R, Cooper GM. Requirement for phosphatidylinositol-3 kinase in the prevention of apoptosis by nerve growth factor. Science. 1995; 267:2003-6. [PMID: 7701324]
  12. Datta SR, Brunet A, Greenberg ME. Cellular survival: a play in three Akts. Genes Dev. 1999; 13:2905-27. [PMID: 10579998]
  13. Oltvai ZN, Milliman CL, Korsmeyer SJ. Bcl-2 heterodimerizes in vivo with a conserved homolog, Bax, that accelerates programmed cell death. Cell. 1993; 74:609-19. [PMID: 8358790]
  14. Cheng H, Ding Y, Yu R, Chen J, Wu C. Neuroprotection of a Novel Cyclopeptide C*HSDGIC* from the Cyclization of PACAP (1–5) in Cellular and Rodent Models of Retinal Ganglion Cell Apoptosis. PLoS ONE. 2014; 9:e108090 [PMID: 25286089]
  15. Starczynski J, Pepper C, Pratt G, Hooper L, Thomas A, Milligan D, Bentley P, Fegan C. Common polymorphism G(−248)A in the promoter region of the bax gene results in significantly shorter survival in patients with chronic lymphocytic Leukemia once treatment is initiated. J Clin Oncol. 2005; 23:1514-21. [PMID: 15735127]
  16. Sibley K, Rollinson S, Allan JM, Smith AG, Law GR, Roddam PL, Skibola CF, Smith MT, Morgan GJ. Functional FAS promoter polymorphisms are associated with increased risk of acute myeloid leukemia. Cancer Res. 2003; 63:4327-30. [PMID: 12907599]
  17. Wu J, Metz C, Xu X, Abe R, Gibson AW, Edberg JC, Cooke J, Xie F, Cooper GS, Kimberly RP. A novel polymorphic CAAT/enhancer-binding protein beta element in the FasL gene promoter alters Fas ligand expression: a candidate background gene in African American systemic lupus erythematosus patients. J Immunol. 2003; 170:132-8. [PMID: 12496392]
  18. Harris SL, Gil G, Robins H, Hu W, Hirshfield K, Bond E, Bond G, Levine AJ. Detection of functional single-nucleotide polymorphisms that affect apoptosis. Proc Natl Acad Sci USA. 2005; 102:16297-302. [PMID: 16260726]
  19. Brunet A, Bonni A, Zigmond MJ, Lin MZ, Juo P, Hu LS, Anderson MJ, Arden KC, Blenis J, Greenberg ME. Akt promotes cell survival by phosphorylating and inhibiting a Forkhead transcription factor. Cell. 1999; 96:857-68. [PMID: 10102273]
  20. Tan HY, Chen AG, Chen Q, Browne LB, Verchinski B, Kolachana B, Zhang F, Apud J, Callicott JH, Mattay VS, Weinberger DR. Epistatic interactions of AKT1 on human medial temporal lobe biology and pharmacogenetic implications. Mol Psychiatry. 2012; 17:1007-16. [PMID: 21788944]
  21. Tan HY, Nicodemus KK, Chen Q, Li Z, Brooke JK, Honea R, Kolachana BS, Straub RE, Meyer-Lindenberg A, Sei Y, Mattay VS, Callicott JH, Weinberger DR. Genetic variation in AKT1 is linked to dopamine-associated prefrontal cortical structure and function in humans. J Clin Invest. 2008; 118:2200-8. [PMID: 18497887]
  22. Semaan SJ, Li Y, Nickells RW. A single nucleotide polymorphism in the Bax gene promoter affects transcription and influences retinal ganglion cell death. ASN Neuro. 2010; 2:e00032 [PMID: 20360947]
  23. Libby RT, Li Y, Savinova OV, Barter J, Smith RS, Nickells RW, John SW. Susceptibility to neurodegeneration in a glaucoma is modified by Bax gene dosage. PLoS Genet. 2005; 1:17-26. [PMID: 16103918]
  24. Zintzaras E, Kitsios GD. Synopsis and synthesis of candidate-gene association studies in chronic lymphocytic leukemia: the CUMAGAS-CLL information system. Am J Epidemiol. 2009; 170:671-8. [PMID: 19700502]
  25. Curtin JF, Cotter TG. Live and let die: regulatory mechanisms in Fas-mediated apoptosis. Cell Signal. 2003; 15:983-92. [PMID: 14499341]
  26. Xiang N, Li XM, Wang GS, Tao JH, Li XP. Association of Fas gene polymorphisms with systemic lupus erythematosus: a meta-analysis. Mol Biol Rep. 2013; 40:407-15. [PMID: 23065220]
  27. Sezgin M, Barlas IO, Yildir S, Turkoz G, Ankarali HC, Sahin G, Erdal ME. Apoptosis-related Fas and FasL gene polymorphisms' associations with knee osteoarthritis. Rheumatol Int. 2013; 33:2039-43. [PMID: 23392773]
  28. Yildir S, Sezgin M, Barlas IO, Turkoz G, Ankarali HC, Sahin G, Erdal ME. Relation of the Fas and FasL gene polymorphisms with susceptibility to and severity of rheumatoid arthritis. Rheumatol Int. 2013; 33:2637-45. [PMID: 23749041]
  29. Wang X, Xing GH, Fan CC. Association between the FAS rs2234767G/A polymorphism and cancer risk: a systematic review and meta-analysis. DNA Cell Biol. 2014; 33:320-7. [PMID: 24568648]
  30. Sunter NJ, Scott K, Hills R, Grimwade D, Taylor S, Worrillow LJ, Fordham SE, Forster VJ, Jackson G, Bomken S, Jones G, Allan JM. A functional variant in the core promoter of the CD95 cell death receptor gene predicts prognosis in acute promyelocytic leukemia. Blood. 2012; 119:196-205. [PMID: 22084312]
  31. Zhang Z, Qiu L, Wang M, Tong N, Li J. The FAS ligand promoter polymorphism, rs763110 (−844C>T), contributes to cancer susceptibility: evidence from 19 case-control studies. Eur J Hum Genet. 2009; 17:1294-303. [PMID: 19337311]
  32. Morishima Y, Gotoh Y, Zieg J, Barrett T, Takano H, Flavell R, Davis RJ, Shirasaki Y, Greenberg ME. Beta-amyloid induces neuronal apoptosis via a mechanism that involves the c-Jun N-terminal kinase pathway and the induction of Fas ligand. J Neurosci. 2001; 21:7551-60. [PMID: 11567045]
  33. Ethell DW, Kinloch R, Green DR. Metalloproteinase shedding of Fas ligand regulates beta-amyloid neurotoxicity. Curr Biol. 2002; 12:1595-600. [PMID: 12372252]
  34. Razeghinejad MR, Kamali-Sarvestani E. Aqueous humor levels of soluble Fas and Fas-ligand in patients with primary open angle and pseudoexfoliation glaucoma. Iran J Immunol. 2007; 4:215-9. [PMID: 18057579]
  35. Borkenstein A, Faschinger C, Maier R, Weger M, Theisl A, Demel U, Graninger W, Irene H, Mossböck G. Measurement of tumor necrosis factor-alpha, interleukin-6, Fas ligand, interleukin-1alpha, and interleukin-1beta in the aqueous humor of patients with open angle glaucoma using multiplex bead analysis. Mol Vis. 2013; 19:2306-11. [PMID: 24265545]
  36. Zalewska R, Zalewski B, Reszec J, Mariak Z, Zimnoch L, Proniewska-Skretek E. The expressions of Fas and caspase-3 in human glaucomatous optic nerve axons. Med Sci Monit. 2008; 14:BR274-8. [PMID: 19043361]
  37. Ju KR, Kim HS, Kim JH, Lee NY, Park CK. Retinal glial cell responses and Fas/FasL activation in rats with chronic ocular hypertension. Brain Res. 2006; 1122:209-21. [PMID: 17045251]