Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic association of -1562C>T polymorphism in the MMP9 gene with primary glaucoma in a north Indian population

  • Nanamika Thakur,

    Roles Formal analysis, Methodology, Writing – original draft

    Affiliation Department of Human Genetics, Guru Nanak Dev University, Amritsar, Punjab, India

  • Manu Kupani,

    Roles Methodology, Writing – original draft

    Affiliation Department of Human Genetics, Guru Nanak Dev University, Amritsar, Punjab, India

  • Rajeev Kumar Pandey,

    Roles Conceptualization, Writing – review & editing

    Affiliation Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America

  • Rashim Mannan,

    Roles Investigation, Supervision

    Affiliation Baba Deep Singh Charitable Hospital, Amritsar, Punjab, India

  • Archna Pruthi,

    Roles Investigation, Supervision

    Affiliation Baba Deep Singh Charitable Hospital, Amritsar, Punjab, India

  • Sanjana Mehrotra

    Roles Conceptualization, Formal analysis, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editing

    sanjana.bhu@gmail.com

    Affiliation Department of Human Genetics, Guru Nanak Dev University, Amritsar, Punjab, India

Abstract

MMP (Matrix metalloproteinase) 9 is reported to affect glaucoma pathogenesis by altering intraocular pressure (IOP) through its role in remodeling the extracellular matrix (ECM) in the trabecular meshwork. A genetic variant at the promoter region in the MMP9 gene (-1562C>T) has a putative role in regulating its transcription rate and hence can affect genetic predisposition to primary glaucoma. The present study examined the association of -1562C>T promoter polymorphism in the MMP9 gene with Primary Open Angle Glaucoma (POAG) and Primary Angle Closure Glaucoma (PACG) in a north Indian population. A total of 729 subjects (POAG = 224, PACG = 138 and 367 controls) were recruited for the study. Genotyping for the promoter sequence variant was done with PCR-RFLP method. Genotypic and allelic frequency distribution of the POAG and PACG data sets were compared to that of controls by chi-square test and genetic association was tested under different genetic models as implemented under PLINK. Statistically significant difference was observed in the genotype frequencies between PACG cases and controls (p = 0.030). However, in the POAG cases, this difference was only borderline (p = 0.052). Genetic model analysis, under the dominant model revealed 1.6 and 1.4 fold increased susceptibility to PACG and POAG (p = 0.012, p = 0.032) respectively. A higher frequency of CT genotype was observed in PACG as well as POAG males as compared to female subjects. According to the dominant model, CT+TT genotype conferred 1.8 fold higher risk of developing PACG among male patients as compared to the control group (p = 0.048, OR = 1.87;1.00–3.50). Current findings suggest significant association of MMP9 -1562C>T polymorphism with primary glaucoma in the targeted north Indian population and warrant further replication of the findings in other populations.

Introduction

Glaucoma, an optic neuropathy is the second leading cause of blindness worldwide after cataract [1, 2]. Being a complex disorder, it manifests in different clinical forms, among which primary open angle and angle closure account for bulk of the cases [36]. The vision loss occurs due to gradual damage to the retinal ganglion cells (RGCs) and the optic nerve in response to elevated intraocular pressure (IOP) which is a major risk factor for glaucoma development [7]. Genetics is another important risk factor as suggested by higher risk of both open angle and angle closure glaucoma among first-degree relatives of affected individuals in twin and family based studies [810]. Several genetic association studies have identified loci which might influence the genetic predisposition to glaucoma development and progression [9, 10], yet these variants explain only less than 10% of heritability [9]. The progressive and irreversible apoptosis of RGCs, the axons of which form the optic nerve is an important pathogenic feature in glaucoma [11]. The primary site of damage in glaucoma is controversial, it could be the optic nerve head (ONH) or retina, in either case elevation in IOP is an important contributing factor [11]. The initiating molecular events leading to high IOP conditions in the eye are not completely understood, but may occur primarily due to an imbalance of aqueous humor production by the ciliary body and its outflow resistance via the trabecular meshwork (TM) [12]. Since the TM determines the outflow resistance by homeostatic turnover of its extracellular matrix (ECM), pathways/proteins affecting ECM remodeling assume importance in glaucoma pathogenesis and can be targets for therapeutic intervention [13]. Matrix metalloproteinases (MMPs), a group of zinc proteinases are involved in degradation of ECM at TM and lamina cribosa (LC) [14]. Among different MMPs, MMP9 encodes a 92-kDa multidomain zinc dependent enzyme known as gelatinase or type V collagenase and is known to extensively affect ECM deposition and turnover in the TM and LC regions in glaucoma [14]. Numerous studies have linked changes in the expression of MMP9 in the retina, optic nerve, aqueous humor, and TM with glaucomatous eyes in humans [1517] and animal models of glaucoma [1819]. The altered expression of MMPs could be a response to elevated IOP and simultaneously contribute to it by changing the outflow resistance. In addition to their defining role in affecting IOP by remodeling of ECM of the TM in the anterior segment of the eye [20], abnormal expression of MMP9 can also affect RGCs survival as shown by Guo et alwhere MMP9 levels correlated with elevation in IOP and RGC apoptosis [7]. Aberrant MMP9 activity has also been implicated in both ischemia and excitotoxicity-mediated RGC damage [21]. In ischemic conditions wherein membrane depolarization is one of the initiating events for injury, injection of depolarizing agents like KCl into vitreous humor of mice induces up-regulation of Mmp9 activity in the retina [21]. This depolarization-induced Mmp9 up-regulation is through N-methyl-D-aspartate (NMDA) and non-NMDA type glutamate receptors as intravitreal injection of glutamate receptor antagonists along with KCl, resulted in reduced Mmp9 activity [2122]. In another study, NMDA mediated excitotoxic damage to RGCs was shown to be through Mmp9 activation via neuronal nitric oxide synthase [23]. The molecular mechanisms underlying how MMP9 contributes to RGC death are not clearly understood but may involve their role in degradation of laminin, one of the major components of the basement membrane on the inner wall of Schlemm’s canal. Experiments with Mmp9 deficient mice have indicated a direct causal relationship between MMP9 and degradation of laminin [2425]. Genetic variants in the MMP9 gene, specifically cis-regulatory elements which provide a binding site for transcription factors, can influence its expression and hence may modify overall genetic risk for glaucoma onset or progression [26]. rs3918242 (-1562C>T) variant in the promoter region exerts a functional effect on the transcription of MMP9 gene and has been found to be associated with high MMP9 levels in serum samples [27]. The present study therefore aims to assess whether -1562C>T polymorphism in the MMP9 gene promoter is associated with primary glaucoma in a north Indian population since the polymorphism has not been investigated in the targeted population, more so in any Indian population.

Materials and methods

Study participants

Study participants were recruited from Baba Deep Singh Eye hospital, Amritsar (Punjab), India after obtaining written informed consent. The study was approved by the Research Ethical Committee of Guru Nanak Dev University, Amritsar, Punjab, India and the study protocols were in accordance to the principles of Declaration of Helsinki. 362 primary glaucoma patients (POAG = 224 and PACG = 138) were enrolled. For POAG patients, the inclusion was based on following criteria: IOP of greater than 21mm Hg in either of the eyes tested using Goldmann Applanation Tonometry, glaucomatous ONH damage defined as a vertical cup-disc ratio (VCDR) 0.7 or greater as adjudged clinically on slit lamp biomicroscopy using hand held +90 D. This was confirmed using contrast enhanced fundus photograph on optical coherence tomography (OCT) as well as optic disc analysis or glaucomatous visual field defect as detected on automated perimeter using Humpherys Visual Field Analyser using Swedish Interactive Thresholding Algorithm (SITA) standard protocols. PACG cases were also recruited on the basis of above described criteria along with the presence of at least 180 degrees of closed angle in which the TM is not visible on gonioscopy. The control data set consisted of 367 unrelated age and gender matched subjects without any family history of glaucoma. The controls were all examined, prior to cataract surgery, for IOP less than 21 mmHg, normal visual field, normal optic nerve heads with CDR of <0.5. Individuals with known chronic systemic inflammatory diabetes (RBS-140mg/Dl; according to American Diabetes Association (ADA) guidelines), autoimmune or immunosuppressive disease as well as a pre-existing ocular disease (diabetic retinopathy, age-related macular degeneration) were excluded from the study. Individuals having ocular hypertension (OHT) were also excluded from the study group.

Sample collection, DNA isolation and Genotyping

Venous blood was collected in EDTA vials and genomic DNA was extracted using phenol chloroform method [28]. All samples were diluted to 100ng/μl stock concentration. Quantification of extracted DNA was done using NanoDrop ND-2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). The target region (436bp) was amplified with the primer sequences and conditions as previously reported [29]. PCR was performed in a final reaction volume of 15μl containing 50ng DNA, 1X Taq Polymerase buffer reagent, 2.5mmol of each deoxynucleotide triphosphate (dNTPs) (3B BlackBio Biotech India Ltd), 20mmol MgCl2, 0.3μmol/μl of each primer, and 0.09 units Taq DNA polymerase (3B BlackBio Biotech India Ltd; 5U/μl). The cycling conditions were initial denaturation at 95°C for 10 min, followed by 40 cycles at 95°C for 45 s, 61.5°C for 45 s and 72°C for 45 s, with a final extension at 72°C for 10 min. After amplification, the PCR products were digested with 3 units of Sph-I (New England Biolabs) for 12h at 37°C. The digested products were separated by electrophoresis on 2.5% agarose gel stained with ethidium bromide (10mg/ml) (GeNeiTM) and genotypes were scored on the basis of Restriction Fragment Length Polymorphism (RFLP) pattern as given in Fig 1. As a quality control measure, genotypes in few samples were confirmed by sanger sequencing (Fig 2).

thumbnail
Fig 1. Digested PCR product of MMP9 -1562C>T.

Well1 = 100bp Ladder; S2,10,11,12 = Heterozygous (436/242/194bp); S1,3,4,5,6,7,8,9 = Homozygous wild (436bp); UD = Undigested PCR product.

https://doi.org/10.1371/journal.pone.0192636.g001

thumbnail
Fig 2. Sanger sequencing chromatograms for validation of genotypes.

A) Homozygous CC genotype B) Heterozygous CT genotype C) Homozygous TT genotype.

https://doi.org/10.1371/journal.pone.0192636.g002

Statistical analysis

Descriptive statistical analysis for demographic and clinical characteristics for study participants was performed using GraphPad Prism 5.0. The results were tabulated as mean ± standard deviation (SD). Association analyses were performed using PLINK software (v1.07). Odds ratios (ORs) for various genetic models, dominant, recessive, additive and allelic were calculated by Medcalc software. P-values of less than 0.05 were considered to be statistically significant. To investigate the association of the variant with IOP and VCDR which are known endophenotypes of glaucoma, the values for right eye was chosen arbitrarily (as the mean for both eyes was same, Table 1) for ANOVA.

thumbnail
Table 1. Clinical and demographic details of the study groups.

Values are expressed as mean ± SD.

https://doi.org/10.1371/journal.pone.0192636.t001

Results

The present case control genetic association study was conducted on 729 individuals, comprising 224 POAG cases, 138 PACG cases and 367 control subjects. The demographic and clinical parameters of these groups are presented in Table 1.

Frequency of males (62.05%) among POAG cases was higher as compared to PACG (37.68%), while the female subjects were more prevalent in PACG (62.31%) with respect to POAG (38.39%). These results highlight the gender difference in glaucoma prevalence; women are known to be at a higher risk for PACG [30]. Although the gender predilection for POAG is not very clear, in the present study prevalence of POAG was found to be higher in men [30]. Total number of males and females in cases (both PAOG/PACG) and controls did not show any significant difference (p = 0.889) as given in Fig 3.

thumbnail
Fig 3. Frequency distribution of males and females among POAG and PACG cases with respect to control subjects.

https://doi.org/10.1371/journal.pone.0192636.g003

-1562C>T MMP9 polymorphism revealed significant association with PACG

The observed genotypes for rs3918242 followed Hardy Weinberg equilibrium (HWE) in controls (observed and expected heterozygosity = 0.28, 0.26 respectively; p = 0.244). The genotype and allele frequency distribution along with genetic model analysis for subjects of PACG and controls is given in Table 2. PACG cases showed statistically significant difference in allelic and genotypic frequencies with respect to controls. 22.82% PACG cases had minor allele, ‘T’ whereas 15.94% control individuals had T allele (p = 0.011), thereby making it the risk allele in our population. The frequencies of CC, CT and TT genotypes were 57.97%, 38.40%, 3.62% in PACG patients while in controls the corresponding frequencies were 69.75%, 28.61%, 1.63% (p = 0.030). Analysis of genotypic distribution revealed ~1.6 fold increased susceptibility to PACG under dominant (CT+TT vs CC) (p = 0.012, OR = 1.65; 1.11–2.50) and co-dominant models (CT vs CC+TT) (p = 0.035, OR = 1.55; 1.03–2.34) as shown in Table 2.

thumbnail
Table 2. Distribution of allele frequency, genotype frequency and genetic model of rs3918242 among PACG cases and controls.

https://doi.org/10.1371/journal.pone.0192636.t002

Higher prevalence of CT genotype in PACG male patients

Subsequent segregation of PACG cases on the basis of gender (Table 3) revealed a different frequency distribution of genotypes among PACG males and females. A higher number of heterozygotes were observed in PACG males (40.38%) than in PACG females (34.88%) but the difference was not statistically significant. When the comparison of PACG males and females was done with healthy male and female subjects, the study failed to obtain any significant difference in the genotypic frequencies. Genetic model investigation did not reveal a statistically significant difference in genotypic distribution among affected and control female subjects. In contrast to females, the dominant model unveiled 1.8 fold higher risk towards PACG in males (p = 0.048, OR = 1.87; 1.00–3.05).

thumbnail
Table 3. Gender wise comparison of PACG cases and controls.

https://doi.org/10.1371/journal.pone.0192636.t003

-1562C>T MMP9 polymorphism revealed borderline significance with POAG

The prevalence of CC genotype was greater in controls (69.75%) than in cases (61.16%) while frequency of CT genotype was found to be higher in cases (37.94%) as compared to the control group (28.61%). Incidence of TT genotype was lesser in both POAG cases and controls (0.89%, 1.63% respectively). The difference in genotypic distribution showed only borderline significance (p = 0.052). The prevalence of minor allele ‘T’ was higher in cases (19.86%) than controls (15.94%), however the difference was not statistically significant. The genetic model analysis showed 1.5 fold increased susceptibility toward POAG, conferred by CT genotype under co-dominant model (p = 0.018, OR = 1.52; CI = 1.03–2.17). Dominant model also revealed a statistically different genotypic distribution between cases and controls (p = 0.032, OR = 1.46; CI = 1.07–2.17) (Table 4).

thumbnail
Table 4. Distribution of allele frequency, genotype frequency and genetic model of rs3918242 among POAG cases and controls.

https://doi.org/10.1371/journal.pone.0192636.t004

Higher distribution of CT Genotype among POAG males than female patients

PAOG cases were segregated into two sub-groups of males and females. The frequency of CT genotype was found to be higher in affected males (37.41%) with respect to affected females (34.11%), while CC genotype was observed to be more in number among female patients as compared to males having POAG (64.70%, 62.58%). None of the POAG males had TT genotype in contrast to control males (2.05%) as shown in Table 5. The difference in genotypic distribution between POAG males and control male subjects was close to significance threshold (p = 0.050). However, the affected females did not show any difference in genotypic frequency as compared to control females. Genetic model analysis failed to obtain any difference in genotypic distribution among male and female patients with respect to corresponding controls. No significant correlation was observed between -1562C>T polymorphism and IOP and VCDR in POAG and PACG (Table 6).

thumbnail
Table 5. Gender-wise comparison of POAG cases and controls.

https://doi.org/10.1371/journal.pone.0192636.t005

thumbnail
Table 6. The genotype distribution of -1562C>T MMP9 gene polymorphisms in patients with open-angle glaucoma (POAG) and angle closure glaucoma (PACG) according to CDR and IOP.

https://doi.org/10.1371/journal.pone.0192636.t006

Discussion

The present study investigated the association of rs3918242 (-1562C>T), a functional variant in the promoter region of MMP9gene with primary glaucoma in a north Indian population. In glaucoma pathogenesis, MMPs are important candidate genes because they can contribute to elevation in IOP by shifting the equilibrium between ECM synthesis and proteolysis within the TM, thereby modifying outflow resistance and hindering aqueous drainage [31]. Among different members of the MMP family, MMP9expression and activity has been extensively associated with glaucomatous optic neuropathies (GONs) [32]. The exact mechanistic details of how MMP9 contributes to GONs are not well understood but may involve detachment induced RGC death due to extensive ECM remodeling [32]. This has significant implication since the vision loss in glaucoma is due to irreversible damage to RGCs. That MMP9 may have a causal role in RGCs death was shown in a rat glaucoma model wherein the neuroprotection conferred by pyrrolidine dithiocarbamate (PDTC) to RGCs was mediated by downregulation of Mmp9 [33]. Studies with Mmp9 null mice also demonstrate altered TM composition, reduced aqueous humor drainage and IOP elevation [34] thereby establishing MMP9 as an important remodeler of TM. Apart from functional evidence, molecular genetics analyses in different populations have reported a link between the MMP9gene polymorphisms and POAG and PACG. In a Taiwanese population, highly significant association was observed between PACG and a non synonymous SNP in exon 6 of the MMP9 gene (rs2664538; now merged into rs17576) [35]. Association with rs17576, along with another SNP (rs3918249) was replicated in an Australian Caucasian population [36], but not in Singaporean subjects [37] and neither in a Caucasian population from southern Austria [26]. A population based case control study on individuals of Chinese ethnicity, conducted by Shi et al revealed no significant association between the two variants of MMP9 gene (rs17576 and rs3918249) and primary angle closure glaucoma [38]. Another study in Chinese Han population, investigated 4 SNPs (rs3918249, rs3918254, rs17577 and rs3787268) in MMP9 gene, out of that only rs3918254 was found to be associated with PACG [39]. More recently, results of 2 meta-analyses seem contradictory in the sense that rs3918249 was affirmed to be associated with PACG by Rong et al [40] but not by Chen et al [41]. In the latter study, the authors in addition to performing a meta-analysis also evaluated common variants in MMP9 for association with PACG which included rs3918249 in a Chinese population but did not get significant association with any of the 6 SNPs [41]. In yet another recent meta-analysis, MMP9 rs17576 G > A polymorphism was observed to be a protective factor against the development of glaucoma [42]. However none of the above mentioned studies included promoter region polymorphisms in the MMP9, which therefore are not well investigated, both in POAG and PACG. rs3918242 is a C>T variant at position -1562 and affects the transcription of the gene; the CC genotype is associated with low promoter activity while CT, TT combinations exhibit higher activity [15]. The T substitution decreases the binding capacity of a putative transcription repressor protein, thereby causing an upregulation in the gene expression. In luciferase assays, MMP9 gene expression showed a 9-fold and 12-fold increase in the expression level for -1562 T/T genotype as compared to the -1562 C/T and -1562 C/C genotypes respectively [15]. Being a functional variant, it can act as a genetic risk factor for glaucoma development and progression. There are only two reports on the effect of -1562C>T variant on genetic predisposition to primary glaucoma, one in Caucasians from a Polish population [29] and the other is in a Pakistani population [43]. In our study, which is the first one to investigate MMP9 polymorphism in any Indian glaucomatous population we got a significant association of rs3918242 with PACG under dominant as well as co-dominant model while marginal association was observed with POAG (p = 0.052). Our findings are consistent with the results of Markiewicz et al [29] in which a statistically significant association was obtained between rs3918242 and POAG. Stratification of the dataset on the basis of sex revealed higher risk towards PACG in males (p = 0.048, OR = 1.87; 1.00–3.05). Sex specific differences in genetic susceptibility to glaucoma could be accounted by the effect of estrogens on the levels of proinflammatory cytokines like IL-1beta which is known to induce MMP9 expression in primary and immortalized cells [44]. In the Pakistani cohort of Punjabi origin, genetic association was observed between rs17576 and PACG but not with rs3918242 [43]. The contrasting results might be due to genetic admixture in the residing populations of respective areas (Punjab, India and Central Pakistan). The present study has some limitations. Only a single variant in the MMP9 gene was studied. Assessing the role of other SNPs in the gene can allow for more powerful haplotype analysis to better elucidate the role of MMP9 variants as a genetic risk factor for glaucoma. Due to lack of aqueous humor samples we were unable to correlate MMP9 expression levels with the risk genotype and therefore there exists a possibility of presence of other functional variants in linkage disequilibrium with -1562C>T variant. Due to technical constraints, other endophenotypes of glaucoma viz retinal nerve fibre layer (RNFL), rim area (RA), axial length could not be correlated with -1562C>T which has been shown to influence these quantitative traits in POAG [36]. In conclusion, the study provides a first comprehensive data on a functional polymorphism (-1562C>T) in the MMP9 gene in an Indian primary glaucomatous population. In our study population, the polymorphism was found to be associated with PACG and marginally with POAGwith CT genotype conferring a genetic risk for the condition. These results give additional impetus to investigate the role of MMP9 gene polymorphisms in primary glaucoma in other populations and also to conduct functional studies to understand the role of MMP9 in the pathology of POAG and PACG.

Acknowledgments

We are sincerely thankful to all study subjects, who participated in the study by donating their precious blood samples. The work was funded by Department of Science and Technology, Science & Engineering Research Board (SERB), Grant No. SB/YS/LS-140/2014, University Grants Commission—Special Assistance Programme and Start Up BSR-grant.

References

  1. 1. Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006; 90: 262–267. pmid:16488940
  2. 2. 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–2090. pmid:24974815
  3. 3. Friedman DS, Jampel HD, Munoz B, West SK. The prevalence of open-angle glaucoma among blacks and whites 73 years and older: the Salisbury Eye Evaluation Glaucoma Study. Arch Ophthalmol. 2006; 124: 1625–1630. pmid:17102012
  4. 4. Casson RJ, Newland HS, Muecke J, McGovern S, Abraham L, Shein WK, et al. Prevalence of glaucoma in rural Myanmar: the Meiktila Eye Study. Br J Ophthalmol. 2007; 91: 710–714. pmid:17510475
  5. 5. Song W, Shan L, Cheng F, Fan P, Zhang L, Qu W, et al. Prevalence of glaucoma in a rural northern China adult population: a population-based survey in Kailu County, Inner Mongolia. Ophthalmol.2011; 118: 1982–1988.
  6. 6. Cho HK, Kee C. Population-based glaucoma prevalence studies in Asians. Surv Opthalmol. 2014; 59: 434–447.
  7. 7. Guo L, Moss SE, Alexander RA, Ali RR, Fitzke FW, Cordeiro MF. Retinal ganglion cell apoptosis in glaucoma is related to intraocular pressure and IOP-induced effects on extracellular matrix. Invest Ophthalmol Visual Sci. 2005; 46: 175–182.
  8. 8. 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
  9. 9. Abu-Amero K, Kondkar AA, Chalam KV. An updated review on the genetics of primary open angle glaucoma. Int. J. Mol. Sci. 2015; 16: 28886–28911. pmid:26690118
  10. 10. Ahram DF, Alward WL, Kuehn MH. The genetic mechanisms of primary angle closure glaucoma. Eye (Lond) 2015; 29:1251–1259.
  11. 11. Tamm ER, Ethier CR. Biological aspects of axonal damage in glaucoma: A brief review. Exp Eye Res. 2017;157: 5–12. pmid:28223179
  12. 12. Vranka JA, Kelley MJ, Acott TS, Keller KE. Extracellular matrix in the trabecular meshwork: intraocular pressure regulation and dysregulation in glaucoma. Exp Eye Res. 2015;133: 112–125. pmid:25819459
  13. 13. Weinreb RN, Aung T, Medeiros FA. The pathophysiology and treatment of glaucoma: a review. JAMA. 2014; 311: 1901–1911. pmid:24825645
  14. 14. Page-McCaw A, Ewald AJ, Werb Z. Matrix metalloproteinases and the regulation of tissue remodelling. Nat Rev Mol Cell Biol. 2007; 8: 221–233. pmid:17318226
  15. 15. Markiewicz L, Pytel D, Mucha B, Szymanek K, Szaflik J, Szaflik JP, et al. Altered expression levels of MMP1, MMP9, MMP12, TIMP1, and IL-1β as a risk factor for the elevated IOP and optic nerve head damage in the primary open-angle glaucoma patients. Biomed Res Int. 2015;2015: 812503. pmid:26120586
  16. 16. Schlotzer-Schrehardt U, Lommatzsch J, Kuchle M, Konstas AG, Naumann GO. Matrix metalloproteinases and their inhibitors in aqueous humor of patients with pseudoexfoliation syndrome/glaucoma and primary open-angle glaucoma. Invest Ophthalmol Vis Sci. 2003;44: 1117–1125. pmid:12601038
  17. 17. Agapova OA, Ricard CS, Salvador-Silva M, Hernandez. Expression of matrix metalloproteinases and tissue inhibitors of metalloproteinases in human optic nerve head astrocytes. Glia. 2001;33: 205–216. pmid:11241738
  18. 18. Robertson JV, Siwakoti A, West-Mays JA. Altered expression of transforming growth beta 1 and matrix metalloproteinase-9 results in elevated intraocular pressure in mice. Mol Vis. 2013;19: 684–695. pmid:23559862
  19. 19. Zhang X, Cheng M, Chintala SK. Optic nerve ligation leads to astrocyte-associated matrix metalloproteinase-9 induction in the mouse retina. Neurosci Lett 2004;356: 140–144. pmid:14746883
  20. 20. Wong TT, Sethi C, Daniels JT, Limb GA, Murphy G, Khaw PT. Matrix metalloproteinases in disease and repair processes in the anterior segment. Surv Ophthalmol. 2002; 47: 239–256. pmid:12052410
  21. 21. Zhang X, Cheng M, Chintala SK. Kainic acid-mediated upregulation of matrix metalloproteinase-9 promotes retinal degeneration. Invest Ophthalmol Vis Sci. 2004; 45: 2374–2383. pmid:15223820
  22. 22. Mali RS, Cheng M, Chintala SK. Intravitreous injection of a membrane depolarization agent causes retinal degeneration via matrix metalloproteinase-9. Invest Ophthalmol Vis Sci. 2005; 46: 2125–2132. pmid:15914633
  23. 23. Manabe S, Gu Z, Lipton SA. Activation of matrix metalloproteinase-9 via neuronal nitric oxide synthase contributes to NMDA-induced retinal ganglion cell death. Invest Ophthalmol Vis Sci. 2005; 46: 4747–4753. pmid:16303975
  24. 24. Chintala SK, Zhang X, Austin JS, Fini ME. Deficiency in matrix metalloproteinase gelatinase B (MMP-9) protects against retinal ganglion cell death after optic nerve ligation. J Biol Chem. 2002; 277: 47461–47468. pmid:12354772
  25. 25. Gu Z, Cui J, Brown S, Fridman R, Mobashery S, Strongin AY, et al. A highly specific inhibitor of matrix metalloproteinase-9 rescues laminin from proteolysis and neurons from apoptosis in transient focal cerebral ischemia. J Neurosci. 2005; 25: 6401–6408. pmid:16000631
  26. 26. Mossböck G, Weger M, Faschinger C, Zimmermann C, Schmut O, Renner W, et al. Role of functional single nucleotide polymorphisms of MMP1, MMP2, and MMP9 in open angle glaucomas. Mol Vis. 2010;16: 1764–1770. pmid:20808730
  27. 27. Yadav SS, Mandal RK, Singh MK, Verma A, Dwivedi P, Sethi R, Usman K, Khattri S. High Serum Level of Matrix Metalloproteinase 9 and Promoter Polymorphism– 1562 C:T as a New Risk Factor for Metabolic Syndrome. DNA Cell Bio. 2014; 33: 816–822.
  28. 28. Adeli K, Ogbonna G. Rapid purification of human DNA from whole blood for potential application in clinical chemistry laboratories. Clin Chem. 1990; 36: 261–264. pmid:1967986
  29. 29. Markiewicz L, Majsterek I, Przybylowska K, Dziki L, Waszczyk M, Gacek M, et al. Gene polymorphisms of the MMP1, MMP9, MMP12, IL‐1β and TIMP1 and the risk of primary open‐angle glaucoma. Acta Ophthalmol. 2013; 91: 516–523.
  30. 30. Vajaranant TS, Nayak S, Wilensky JT, Joslin CE. Gender and glaucoma: what we know and what we need to know. Curr Opin Ophthalmol. 2010; 21:91. pmid:20051857
  31. 31. Singh D, Srivastava SK, Chaudhari TK, Upadhyay G. Multifaceted role of matrix metalloproteinases (MMPs). Front Mol Biosci. 2015;2:9.
  32. 32. De Groef L, Andries L, Lemmens K, Van Hove I, and Moons L. Matrix metalloproteinases in the mouse retina: a comparative study of expression patterns and MMP antibodies. BMC Ophthalmol. 2015; 15: 187. pmid:26714639
  33. 33. Sun MH, Chen KJ, Tsao YP, Kao LY, Han WH, Lin KK, et al. Down-regulation of matrix metalloproteinase-9 by pyrrolidine dithiocarbamate prevented retinal ganglion cell death after transection of optic nerve in rats. Curr Eye Res. 2011;36; 1053–1063. pmid:21999230
  34. 34. De Groef L, Andries L, Siwakoti A, Geeraerts E, Bollaerts I, et al. Aberrant collagen composition of the trabecular meshwork results in reduced aqueous humor drainage and elevated IOP in MMP-9 null mice. Invest Ophthalmol Vis Sci. 2016; 57: 5984–5995. pmid:27820954
  35. 35. Wang IJ, Chiang TH, Shih YF, Lu SC, Lin LL, Shieh JW, et al. 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–1232. pmid:17110919
  36. 36. 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–1424. pmid:21655354
  37. 37. Aung T, Yong VH, Lim MC, Venkataraman D, Toh JY, Chew PT, et al. 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–258. pmid:18552608
  38. 38. Shi H, Zhu R, Hu N, Shi J, Zhang J, Jiang L, et al. Association of frizzled-related protein (MFRP) and heat shock protein 70 (HSP70) single nucleotide polymorphisms with primary angle closure in a Han Chinese population: Jiangsu Eye Study. Mol Vis. 2013; 19: 128–134. pmid:23378726
  39. 39. Gao XJ, Hou SP, Li PH. The association between matrix metalloprotease 9 gene polymorphisms and primary angle-closure glaucoma in a Chinese Han population. Int J Ophthalmol. 2014;7: 397–402. pmid:24967180
  40. 40. Rong SS, Tang FY, Chu WK, Ma L, Yam JC, Tang SM, et al. Genetic Associations of Primary Angle-Closure Disease: A Systematic Review and Meta-analysis. Ophthalmology. 2016; 123: 1211–1121. pmid:26854036
  41. 41. Chen X, Chen Y, Wiggs JL, Pasquale LR, Sun X, Fan BJ. Association of matrix metalloproteinase-9 (MMP9) variants with primary angle closure and primary angle closure glaucoma. PloS One. 2016; 11: e0157093. pmid:27272641
  42. 42. Wu MY, Wu Y, Zhang Y, Liu CY, Deng CY, Peng L, et al. Associations between matrix metalloproteinase gene polymorphisms and glaucoma susceptibility: a meta-analysis. BMC Opthalmol. 2017; 17: 41.
  43. 43. Micheal S, Yousaf S, Khan MI, Akhtar F, Islam F, Khan WA, et al. Polymorphisms in matrix metalloproteinases MMP1 and MMP9 are associated with primary open–angle and angle closure glaucoma in a Pakistani population. Mol Vis. 2013;19: 441–447. pmid:23441116
  44. 44. Suzuki T, Sullivan DA. Estrogen stimulation of proinflammatory cytokine and matrix metalloproteinase gene expression in human corneal epithelial cells. Cornea. 2005: 24; 1004–1009. pmid:16227852