Polymorphisms rs693421 and rs2499601 at locus 1q43 and their haplotypes are not associated with primary open-angle glaucoma: a case–control study
BMC Research Notes volume 12, Article number: 453 (2019)
The genetic spectrum of primary open-angle glaucoma (POAG) in middle-eastern Saudi’s is still elusive. To this end, we investigated an association between rs693421, rs2499601 and their haplotypes at chromosome 1q43 locus with POAG and its related clinical phenotypes. Genotyping was performed with TaqMan® assays. Haplotypes and their interaction analysis were carried out by SHEsis and SNPStats online tools.
The minor “T” allele frequency of rs693421 was 0.48 in controls and 0.52 in cases (odds ratio (OR) = 1.15, 95% confidence interval (CI) 0.85–1.54, p = 0.368). Similarly, for rs2499601, the minor “C” allele frequency was 0.49 in controls as compared to 0.53 in cases (OR = 1.19, 95% CI 0.89–1.60, p = 0.236). Besides, genotype distribution for both these polymorphisms was also not significant in additive, dominant and recessive models. rs693421 and rs2499601, showed significant linkage disequilibrium (D’ statistics = 0.69, p < 0.001) but haplotype association was non-significant (p = 0.698). The significance did not vary after adjustment to age and sex. No significant genotype association was observed with intraocular pressure, cup/disc ratio and number of anti-glaucoma medication in POAG group. Furthermore, age, sex and genotypes did not contribute any significant risk of POAG in regression analysis. We report no association between rs693421, rs2499601 and their haplotypes with POAG and related phenotypes.
Epidemiological studies in the past suggest that primary open-angle glaucoma (POAG), at least in part, can be caused by heritable factors [1,2,3], with five to tenfold increased risk in the first-degree relatives and an estimated heritability of 0.81 . However, given the polygenic nature of POAG, its exact genetic etiology still remains elusive. Genetic studies are a critical tool to identify genes and molecular pathways in complex human diseases. Likewise, using a similar approach, population-based genome-wide and candidate gene studies have reported a number of polymorphic variations or loci that have been associated with POAG and/or its quantitative endophenotyes [5, 6].
In an initial study, using a 2-stage genome-wide approach in a group of Japanese POAG patients, Nakano et al.  reported 3 genetic loci, harboring six single nucleotide polymorphisms (SNPs), of having significant association with the disease. Of these, 4 intergenic polymorphisms such as rs540782, rs540784, rs693421 and rs2499601 are located on chromosome 1q43. The genes flanking the locus 1q43 include Zona Pellucida glycoprotein 4 (ZP4) and Ryanodine Receptor 2 (RYR2). ZP4 plays a role in fertilization and preimplantation development; and genetic variations in ZP4 are reported to be associated with ovarian disease  and glaucoma . Ryanodine receptors (RyRs) are ubiquitous intracellular calcium (Ca2+) release channels localized in the plasma membrane of endoplasmic or sarcoplasmic reticulum of many organs, including retina . Given the important physiological role of RyR channel, defect(s) in RyR function, either due to oxidative stress or genetic mutations can render these channels leaky to Ca2+ influx. These may thereby produce disease-causing signals and contribute to severe pathologies as those observed in heart failure, muscular dystrophy, diabetes, or neurodegenerative diseases [10, 11]. Also, considering the potential role of calcium channel blockers in treatment of glaucoma , RYR2 may be a potential candidate in POAG.
Replication of case–control association studies in various ethnicities has genetic epidemiological importance, so as to enable future utilization of these genetic biomarkers in evaluation of the disease risk. Association studies using candidate gene approach are a useful tool to identify the genetic contribution of SNPs/loci, each of which may exhibit a relatively small risk but have a significant impact in a population. The genetic spectrum of POAG, especially in middle-eastern Saudi POAG patients, is still unclear. We have previously reported a lack of association between SNPs rs540782 and rs540784 with POAG in a small group of Saudi patients [13, 14]. The present study investigated an association between rs693421, rs2499601 and their haplotypes with POAG or its clinical indices in a much larger and different set of sample cohort of Saudi origin. Besides, there are no other published reports of association of these SNPs (rs693421 and rs2499601) in this population.
Study design and participants
We conducted a retrospective case–control genetic association study. The study followed the principles of the Declaration of Helsinki and had an institutional review board ethics committee approval (# 08-657). Patients (n = 185) and controls (n = 171) were recruited at King Abdul-Aziz University Hospital in Riyadh, Saudi Arabia. The inclusion–exclusion criteria for patients and controls have been detailed elsewhere . Briefly, POAG patients showed (i) optic disk or retinal nerve fiber layer changes; (ii) visual field abnormalities; (iii) bilateral open anterior chamber angles; and (iv) adult onset. Secondary forms of glaucoma and history of steroid usage or ocular trauma were excluded. Controls were > 20 years of age with normal IOP, normal optic disk, open anterior chamber angles and no history of ocular disease(s) or eye surgeries. Subjects refusing to participate were also excluded.
Genotyping of rs693421 and rs2499601 at Chr.1q43
DNA samples were genotyped for rs693421 and rs2499601using the TaqMan® genotyping assays C_80796_20 and C__26180887_10, respectively (Applied Biosystems Inc., Foster City, CA, USA; Cat#: 4351379) as described previously .
SPSS version 22 (IBM Inc. Chicago, IIinois, USA) was used to perform statistical tests as indicated. Linkage disequilibrium analysis was conducted using SHEsis online software (http://analysis.bio-x.cn/myAnalysis.php). SNPStats online software (http://bioinfo.iconcologia.net/snpstats/start.htm) was used to construct haplotypes and analyze the interactions with related factors. p < 0.05 (2-sided) was considered statistically significant.
Demographic and clinical characteristics
Except for family history of glaucoma (p = 0.019), there was no significant difference in age, gender, systemic co-morbidities and smoking habit between POAG cases and controls (Additional file 1: Table S1).
Genotype and allele frequency
rs693421 and rs2499601 did not deviate significantly from the HWE (p > 0.05). For rs693421, the minor “T” allele frequency was 0.48 in controls and 0.52 in cases, but non-significant (OR = 1.15, 95% CI 0.85–1.54, p = 0.368). Besides, genotype distribution was also not significant in additive (Chi2 = 1.74, df = 2, p = 0.417), dominant (p = 0.860) and recessive (p = 0.190) models (Table 1). Similarly, for rs2499601, the minor “C” allele frequency was 0.49 in controls as compared to 0.53 in cases (OR = 1.19, 95% CI 0.89–1.60, p = 0.236) with no significant genotype effect in additive (Chi2 = 1.69, df = 2, p = 0.429), dominant (p = 0.512) and recessive (p = 0.195) models. In addition, for both the SNPs, the genotype/allele distribution between POAG cases and controls was also non-significant after adjustment for age and sex. There was no significant interaction with sex for rs693421 (p = 0.83) and rs2499601 (p = 0.83); and the significance did not vary after adjustment for age.
Linkage disequilibrium and haplotype association
rs693421 and rs2499601 showed significant linkage disequilibrium as reflected by D’ value (D’ statistics = 0.69, p < 0.001) (Additional file 1: Figure S1). No significant haplotype association was observed (p = 0.698) (Table 2). The significance did not vary after adjustment to age and sex (p = 0.750).
Genotype effect on demographic and clinical phenotypes in POAG
Age and gender distribution was not significant for both rs693241 and rs2499601 genotypes. In addition, phenotypes such as intraocular pressure (IOP), cup/disc ratio and number of antiglaucoma medication showed no significant genotype effect (Additional file 1: Table S2). Besides, logistic regression showed that age, sex, rs693241 and rs2496601 could not significantly explain the likelihood of having POAG (Additional file 1: Table S3).
POAG is the second most common type of glaucoma in Saudi Arabia . Despite its high prevalence, the role of genetic factors contributing to the progression and/or development of POAG are largely unknown and thus warrants genetic investigation(s) to identify these factors that may plausibly be associated or contribute to POAG pathogenesis. To this end, evaluation of an association between rs693421, rs2496601 and their haplotypes in a Saudi POAG cohort in this study did not yield any significant relationship with the disease.
Nakano et al.  have previously reported a moderate association of polymorphisms rs693241 (OR = 1.48, 95% CI 1.20–1.83, p = 0.00029) and rs2499601 (OR = 1.45, 95% CI 1.17–1.79, p = 0.000589) at locus 1q43 in 1575 Japanese POAG patients. The findings were successfully replicated in stage-2 GWA study. However, except for rs693421 (OR = 1.4 (1.1–1.7), p = 0.0082) in Korean subjects, subsequent studies in a South Indian , Afro-Caribbean , Japanese , the Hong Kong cohort , Korean (for SNP rs2499601)  and very recently in the Han Chinese population  have failed to replicate this association between rs693241, rs2499601 and POAG. Besides, the functional significance of these SNPs or locus to POAG development remains unknown. In contrast to the findings of Nakano et al.  and consistent with the findings reported in other population [17, 19,20,21], no significant genotype/allelic effect of these SNPs in POAG was observed in our study. These conflicting results can be ascribed to factors such as differences in sample size, lack of adequate statistical power or differences in the phenotypes of the sample population included in the study. A comparison of minor allele frequencies of these SNPs in different ethnic groups is shown in Table 3. The minor allele frequencies observed for these variants in the Saudi population (controls) was comparable to the Oriental population’s [7, 19,20,21], but higher than that reported among the Indian  and Afro-Caribbean ethnic groups .
The genotypes of adjacent polymorphisms are often highly correlated, referred to as being in linkage disequilibrium, and could be inherited as haplotypes. Haplotypes display an interaction of these polymorphisms and can serve as a maker of disease susceptibility . In complex multifactorial diseases, such as POAG, multiple polymorphisms/loci in the form of haplotypes can be more informative as compared to single allele polymorphism analysis. Hence, to gain additional information we performed association between 1q43 locus haplotypes (for polymorphisms rs693421 and rs2499601) and POAG. The four SNPs (rs540782, rs540784, rs693421 and rs2499601) at chromosome 1q43 have been reported to exhibit strong linkage disequilibrium (D’ = 0.90–1.0) [7, 17], but data on haplotypes with these did not indicate any association with POAG among Indians  as opposed to significant association (p < 0.05) among Han Chinese POAG subjects . In our study, significant linkage disequilibrium between rs693421 and rs2499601 was found (D’ = 0.69) but analysis of distribution differences of the four haplotypes in the POAG cases and controls did not show any significant association with the disease. Unfortunately, genotype data for two other SNPs rs540782 and rs540784 at 1q43 were not available for all the subjects included in this study to assess their linkage and/or haplotype association.
Association analysis for any genotype effect of rs693421 and rs2499601 on clinical phenotypes did not provide any significant link. Age and gender showed no significant distribution and clinical phenotypes for glaucoma such as IOP, cup/disc ratio and number of anti-glaucoma medication that serves as markers of disease severity also showed no significant genotype effects in the POAG group.
Similar to Nakano et al. study in 2009, many other recent GWA-studies have reported numerous SNPs in multiple loci/genes in association with POAG and its endophenotypes among European, African and Asian ethnic groups (as reviewed elsewhere [5, 6]). In another recent large-scale study by Shiga et al.  in the Japanese population, variants near cyclin-dependent kinase inhibitor 2B-antisense 1 (CDKN2B-AS1, rs2157719), SIX Homeobox 6 (SIX6, rs33912345) and growth arrest-specific 7 (GAS7, rs9913911) were found to be associated with POAG. In contrast, with the exception of variant rs10483727 in (SIX1/SIX6) gene locus , our previous studies to identify genetic polymorphisms associated with Saudi POAG patients in genes including caveolin 1/caveolin 2 (CAV1/CAV2) , atonal homolog 7 (ATOH7) , (CDKN2B) , transmembrane and coiled-coil domain 1(TMCO1) , (GAS7)  and ZP4 [13, 14], have yielded negative results indicating that the plausible causal genetic defect(s) in POAG cases of Saudi origin may be different than those from European, African and Asian descents, including the Japanese population. There can be several potential reasons for differences in association findings. The differences could reflect an ethnic-specific genetic etiology of POAG. Furthermore, the differences could be attributed to variations in environmental exposures, clinical variability observed in patients with POAG, variations in the endophenotypic traits related to POAG or other epistatic interactions. Nonetheless, based on our findings and that of reported in other ethnic groups, it appears that polymorphisms rs693421 and rs2499601 at chromosome 1q43 (near ZP4/RYR2 genes), may not have a significant role in POAG pathogenesis and hence, cannot serve as genetic markers to assess the disease risk.
Our study could not replicate the previously observed association of rs693421 and rs2499601 in POAG. The results suggest that these SNPs and their haplotypes at locus 1q43, independently or in relation to other clinical phenotypes, are not associated with POAG and hence may not pose a significant risk of POAG in individuals of Saudi origin. However, the lack of association needs cautious interpretation since the study is limited by its sample size, performed in specific ethnicity and their interaction with other genetic or non-genetic factors cannot be ruled out. Nonetheless, based on the minor allele frequency from 1000 Genomes database or that observed in our current study, the sample size used, showed an estimated power of 90% with a detectable odds per allele of 1.60 (2-sided test) and 1.53 (one-sided test). However, assuming an OR of 1.45 as observed in the Japanese POAG patients , our study had only 70% and 71% powers (2-sided test) to detect any significant associations between POAG and rs693421 and rs2499601 polymorphisms, respectively. As in most genetic association studies, in order to detect an odds of ≤ 1.5-fold a much larger sample size would be needed.
Availability of data
The data supporting the conclusions of this article are all presented within the article.
- ATOH7 :
atonal homolog 7
- CAV1/CAV2 :
caveolin 1/caveolin 2
- CDKN2B :
cyclin-dependent kinase inhibitor 2B
genome-wide association study
- GAS7 :
growth arrest-specific 7
primary open angle glaucoma
- RYR2 :
ryanodine receptor 2
single nucleotide polymorphism
- SIX1/SIX6 :
SIX Homeobox 1/SIX Homeobox 6
- TMCO1 :
transmembrane and coiled-coil domain 1
- ZP4 :
zona pellucida glycoprotein 4
Weih LM, Nanjan M, McCarty CA, Taylor HR. Prevalence and predictors of open-angle glaucoma: results from the visual impairment project. Ophthalmology. 2001;108(11):1966–72.
Carbonaro F, Andrew T, Mackey DA, Spector TD, Hammond CJ. Heritability of intraocular pressure: a classical twin study. Br J Ophthalmol. 2008;92(8):1125–8.
Healey P, Carbonaro F, Taylor B, Spector TD, Mitchell P, Hammond CJ. The heritability of optic disc parameters: a classic twin study. Invest Ophthalmol Vis Sci. 2008;49(1):77–80.
Charlesworth J, Kramer PL, Dyer T, Diego V, Samples JR, Craig JE, Mackey DA, Hewitt AW, Blangero J, Wirtz MK. The path to open-angle glaucoma gene discovery: endophenotypic status of intraocular pressure, cup-to-disc ratio, and central corneal thickness. Invest Ophthalmol Vis Sci. 2010;51(7):3509–14.
Abu-Amero K, Kondkar AA, Chalam KV. An updated review on the genetics of primary open angle glaucoma. Int J Mol Sci. 2015;16(12):28886–911.
Wiggs JL, Pasquale LR. Genetics of glaucoma. Hum Mol Genet. 2017;26(R1):R21–7.
Nakano M, Ikeda Y, Taniguchi T, Yagi T, Fuwa M, Omi N, Tokuda Y, Tanaka M, Yoshii K, Kageyama M, et al. Three susceptible loci associated with primary open-angle glaucoma identified by genome-wide association study in a Japanese population. Proc Natl Acad Sci USA. 2009;106(31):12838–42.
Meczekalski B, Nawrot R, Nowak W, Czyzyk A, Kedzia H, Gozdzicka-Jozefiak A. Study on the zona pellucida 4 (ZP4) gene sequence and its expression in the ovaries of patients with polycystic ovary syndrome. J Endocrinol Invest. 2015;38(7):791–7.
Shoshan-Barmatz V, Zakar M, Shmuelivich F, Nahon E, Vardi N. Retina expresses a novel variant of the ryanodine receptor. Eur J Neurosci. 2007;26(11):3113–25.
Santulli G, Lewis D, des Georges A, Marks AR, Frank J. Ryanodine receptor structure and function in health and disease. Subcell Biochem. 2018;87:329–52.
Lanner JT. Ryanodine receptor physiology and its role in disease. Adv Exp Med Biol. 2012;740:217–34.
Ganekal S, Dorairaj S, Jhanji V, Kudlu K. Effect of topical calcium channel blockers on intraocular pressure in steroid-induced glaucoma. J Curr Glaucoma Pract. 2014;8(1):15–9.
Kondkar AA, Edward NB, Kalantan H, Al-Kharashi AS, Altuwaijri S, Mohamed G, Sultan T, Azad TA, Abu-Amero KK. Lack of association between polymorphism rs540782 and primary open angle glaucoma in Saudi patients. J Negat Results Biomed. 2017;16(1):3.
Azad TA, Edward NB, Kondkar AA, Kalantan H, Altuwaijri S, Sultan T, Al-Mobarak FA, Al-Obeidan SA, Abu-Amero KK. Polymorphism rs547984 on human chromosome 1q43 is not associated with primary open angle glaucoma in a Saudi cohort. J Negat Results Biomed. 2017;16(1):12.
Kondkar AA, Azad TA, Almobarak FA, Bahabri IM, Kalantan H, Abu-Amero KK, Al-Obeidan SA. Lack of association between variant rs7916697 in ATOH7 and primary open angle glaucoma in a Saudi cohort. Genet Res Int. 2018;2018:2148056.
Al Obeidan SA, Dewedar A, Osman EA, Mousa A. The profile of glaucoma in a tertiary ophthalmic University Center in Riyadh, Saudi Arabia. Saudi J Ophthalmol. 2011;25(4):373–9.
Rao KN, Kaur I, Chakrabarti S. Lack of association of three primary open-angle glaucoma-susceptible loci with primary glaucomas in an Indian population. Proc Natl Acad Sci USA. 2009;106(44):E125–6 (author reply E127).
Cao D, Jiao X, Liu X, Hennis A, Leske MC, Nemesure B, Hejtmancik JF. CDKN2B polymorphism is associated with primary open-angle glaucoma (POAG) in the Afro-Caribbean population of Barbados, West Indies. PLoS ONE. 2012;7(6):e39278.
Takamoto M, Kaburaki T, Mabuchi A, Araie M, Amano S, Aihara M, Tomidokoro A, Iwase A, Mabuchi F, Kashiwagi K, et al. Common variants on chromosome 9p21 are associated with normal tension glaucoma. PLoS ONE. 2012;7(7):e40107.
Chen LJ, Tam PO, Leung DY, Fan AH, Zhang M, Tham CC, Chiang SW, Fan BJ, Wang N, Pang CP. SNP rs1533428 at 2p16.3 as a marker for late-onset primary open-angle glaucoma. Mol Vis. 2012;18:1629–39.
Kim K, Heo DW, Kim S, Kim JS, Kim CS, Kang C. Expansive marker analysis replicating the association of glaucoma susceptibility with human chromosome loci 1q43 and 10p12.31. Eur J Hum Genet. 2014;22(3):409–13.
Tan Y, Zhang D, Gong B, Zhang X, Jiang L, Lin Y. Association of ZP4 gene polymorphism with primary open-angle glaucoma in an ethnic Chinese population from Sichuan province. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. 2018;35(5):737–40.
Zintzaras E, Lau J. Synthesis of genetic association studies for pertinent gene-disease associations requires appropriate methodological and statistical approaches. J Clin Epidemiol. 2008;61(7):634–45.
Shiga Y, Nishiguchi KM, Kawai Y, Kojima K, Sato K, Fujita K, Takahashi M, Omodaka K, Araie M, Kashiwagi K, et al. Genetic analysis of Japanese primary open-angle glaucoma patients and clinical characterization of risk alleles near CDKN2B-AS1, SIX6 and GAS7. PLoS ONE. 2017;12(12):e0186678.
Kondkar AA, Azad TA, Almobarak FA, Kalantan H, Sultan T, Alsabaani NA, Al-Obeidan SA, Abu-Amero KK. Polymorphism rs10483727 in the SIX1/SIX6 gene locus is a risk factor for primary open angle glaucoma in a Saudi cohort. Genet Test Mol Biomarkers. 2018;22(1):74–8.
Abu-Amero KK, Kondkar AA, Mousa A, Osman EA, Al-Obeidan SA. Lack of association of SNP rs4236601 near CAV1 and CAV2 with POAG in a Saudi cohort. Mol Vis. 2012;18:1960–5.
Kondkar AA, Mousa A, Azad TA, Sultan T, Almobarak FA, Alawad A, Altuwaijri S, Al-Obeidan SA, Abu-Amero KK. Analysis of polymorphism rs1900004 in atonal bHLH transcription factor 7 in Saudi patients with primary open angle glaucoma. Genet Test Mol Biomarkers. 2016;20(11):715–8.
Abu-Amero KK, Kondkar AA, Mousa A, Almobarak FA, Alawad A, Altuwaijri S, Sultan T, Azad TA, Al-Obeidan SA. Analysis of cyclin-dependent kinase inhibitor-2B rs1063192 polymorphism in Saudi patients with primary open-angle glaucoma. Genet Test Mol Biomarkers. 2016;20(10):637–41.
Kondkar AA, Mousa A, Azad TA, Sultan T, Alawad A, Altuwaijri S, Al-Obeidan SA, Abu-Amero KK. Polymorphism rs7555523 in transmembrane and coiled-coil domain 1 (TMCO1) is not a risk factor for primary open angle glaucoma in a Saudi cohort. J Negat Results Biomed. 2016;15(1):17.
Kondkar AA, Azad TA, Almobarak FA, Kalantan H, Sultan T, Al-Obeidan SA, Abu-Amero KK. Polymorphism rs11656696 in GAS7 is not associated with primary open angle glaucoma in a Saudi cohort. Genet Test Mol Biomarkers. 2017;21(12):754–8.
The authors would like to thank the Vice Deanship of Scientific Research Chair, Glaucoma Research Chair in Ophthalmology at the King Saud University for supporting this study.
This work was supported by King Saud University through Vice Deanship of Scientific Research Chair, Glaucoma Research Chair in Ophthalmology but had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The study adhered to the tenets of the Declaration of Helsinki and had received approval from the Institutional Review Board and Research Ethics Committee of College of Medicine, King Saud University, Riyadh, Saudi Arabia (approval number # 08-657). Written, informed consent was obtained from all participants prior to their inclusion in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Demographic and clinical characteristics of POAG cases and controls genotyped for polymorphism rs693421 and rs2499601 included in this study. Figure S1. Linkage disequilibrium plot between the two loci. Table S2. Genotype effect of rs693421 and rs2499601 on glaucoma specific clinical indices in PAOG cases. Table S3. Logistic regression analysis to assess the effect of age, sex, rs693421 and rs2499601 genotypes on POAG outcome.
About this article
Cite this article
Kondkar, A.A., Azad, T.A., Sultan, T. et al. Polymorphisms rs693421 and rs2499601 at locus 1q43 and their haplotypes are not associated with primary open-angle glaucoma: a case–control study. BMC Res Notes 12, 453 (2019). https://doi.org/10.1186/s13104-019-4491-x
- Chromosome 1q43
- Genetic association