Plexin domain containing 2 (PLXDC2) gene polymorphism rs7081455 may not influence POAG risk in a Saudi cohort
BMC Research Notes volume 11, Article number: 733 (2018)
Plexin domain containing 2 (PLXDC2), a cell surface transmembrane protein receptor for pigment epithelium derived factor, is expressed in many tissues including the eye. Polymorphism rs7081455 flanking PLXDC2 has been associated with primary open angle glaucoma (POAG) and its clinical phenotypes and may have a role in POAG. Rs7081455 was genotyped in POAG cases (n = 188) and non-glaucomatous controls (n = 164) of Saudi origin using Taq-Man® to determine any association of this variant with POAG and its endophenotypes.
The risk variant, ‘G’ allele, frequency was 0.56 and 0.52 in controls and POAG cases, respectively (p = 0.197) with was no significant deviation from Hardy–Weinberg equilibrium. Genotype analysis between cases and controls revealed no significant distribution under additive (p = 0.482), dominant (p = 0.590) and recessive models (p = 0.228). In addition, glaucoma specific phenotypic traits such as intraocular pressure (IOP) and cup/disc ratio; and number of anti-glaucoma medications, used to assess severity of the disease, were also statistically non-significant. Furthermore, regression analysis showed no significant effect of age, sex and genotype on disease outcome. Rs7081455 was not associated with POAG or its clinical phenotypes such as IOP and cup/disc ratio and hence may not be a significant risk factor for POAG patients of Saudi origin.
Primary open angle glaucoma (POAG) is a chronic optic neuropathy clinically characterized by an open anterior chamber angle. It is the second most common type of glaucoma seen in the Saudi Arabian population . The disease is largely polygenic in nature that follows a genetically complex and multifactorial inheritance pattern leading to blindness . However, the actual pathophysiological and molecular mechanism(s) leading to glaucoma and its progression still remains elusive. Given the complexity and mutational genetic heterogeneity, POAG shows an estimated heritability of 0.81 . Genetic case–control association studies are an important tool to identify genetic features, each of which may show a relatively small effect but may exhibit significant contribution to a large number of cases. Using similar approaches based on candidate gene and genome wide association studies (GWAS), several investigators have identified a number of genetic variants in multiple loci/genes in various ethnic groups, that have been associated with POAG and related inheritable quantitative traits such as intraocular pressure (IOP) and cup/disc ratio (among others) that may have a causal effect on the development and/or progression of the disease .
Likewise, Nakano et al.  performed a two-stage GWAS in the Japanese population, and demonstrated an association of POAG with genetic variants rs7081455 flanking the plexin domain containing 2 (PLXDC2) gene at chromosomal loci 10p12.31. This association was replicated by Mabuchi et al.  in another Japanese cohort and demonstrated to influence the risk of POAG by increasing IOP. There is no evidence of any association of this gene and/or variant rs7081455 in the POAG patients of Saudi origin, so we conducted a study to determine whether variant rs7081455 near PLXD2 is associated with the disease or any of its clinical phenotypes such as IOP or cup/disc ratio and could be a risk factor for POAG in this population.
Settings and study groups
A case–control association study was performed in accordance to the Declaration of Helsinki principles and was approved by the institutional research ethics committee of College of Medicine at the King Saud University (Approval Number # 08-657). All the participants had signed a written informed consent. POAG patients and non-glaucomatous healthy controls of Saudi origin were selected from the anterior segment and general ophthalmology clinic at King Abdulaziz University Hospital, Riyadh, Saudi Arabia as detailed elsewhere . Briefly, clinical diagnosis of POAG patients (n = 188) were based on the presence of: (i) the disc or retinal nerve fiber layer (RNFL) defect; (ii) visual field abnormalities; and (iii) bilateral open angles on gonioscopy. Secondary glaucoma cases such as pigmentary glaucoma, uveitic, pseudoexfoliation, and history of steroid use or ocular trauma were excluded. Ethically-matched non-glaucomatous control subjects (n = 164) were: > 20 years of age at recruitment; normal IOP [< 21 mmHg without any anti-glaucoma medication]; open angles with normal optic disc on examination; and no history of ocular disease(s) or previous ophthalmic surgeries. Any subject refusing to participate in the study was excluded.
Genotyping of rs7081455 in PLXDC2
DNA samples were genotyped for rs7081455 using the TaqMan® SNP Genotyping Assay (assay ID: C_29390339_10, Applied Biosystems Inc., Foster City, CA, USA) on ABI 7500 Real-Time PCR System (Applied Biosystems) using PCR conditions as described by the manufacturer. The details of the PCR reaction and cycling conditions have been described previously . Genotype calling was done using the in-built automated 2-color allele discrimination software on ABI 7500.
All the statistical analyses were performed using SPSS version 22 (IBM Inc.’ Chicago, IIinois, USA). Pearson’s Chi square test was used to analyze for deviation from Hardy-Weinberg equilibrium (HWE), allelic, genotypic and model-based (additive, dominant and recessive) association with POAG as indicated. 2 × 2 table was used to calculate odds ratio (OR) and 95% confidence interval (CI). Mean differences between genotypes/groups were tested by independent sample t-test, one-way ANOVA and Kruskal–Wallis tests. In addition, effect of age, sex and genotype on the disease outcome was tested by binary logistic regression. A p value < 0.05 was considered statistically significant.
A total number of 352 subjects, consisting of 188 POAG cases and 164 controls, were genotyped for rs7081455 in this study. The demographic and other characteristics of subjects included in this study are shown in the Additional file 1: Table S1. Both the POAG patient and control groups were similar for age, gender distribution, smoking habit and other systemic co-morbidities such as diabetes, hypertension, hypercholesterolemia and heart disease status. However, there was trend towards higher individuals with in family history of glaucoma in POAG group as compared to controls (p = 0.053).
The genotype frequencies showed no significant deviation from HWE (p > 0.05). POAG patients and controls exhibited no significant genotype distribution under additive (Chi2 = 1.46, df = 2, p = 0.482), dominant (OR = 1.14, 95% CI 0.69–1.87, p = 0.590) and recessive models (OR = 1.32, 95% CI 0.87–2.07, p = 0.228). Besides, the allele frequency distribution was also non-significant (OR = 1.30; 95% CI 0.86–1.97; p = 0.197). The risk allele ‘G’ frequency was 0.56 and 0.52 in controls and patients, respectively (Table 1).
The data was further analyzed for the effect of rs7081455 genotype on different demographic and clinical variables within the POAG patient group. There was no significant difference between age (p = 0.305) and gender (p = 0.894) distribution. Similarly, there was no significant genotype effect on family history of glaucoma, smoking and systemic diseases. More importantly, clinical phenotype traits such as IOP and cup/disc ratio; and number of anti-glaucoma medications, that also serve as markers to assess disease severity were also found to be statistically non-significant (Table 2). Besides, binary logistic regression analysis also showed no significant effect of age, sex and genotype on the likelihood of POAG occurrence.
Recent GWAS’s have identified numerous genetic variants in a number of genes/loci to be associated with POAG or its related clinical endophenotypes . Using a similar genome-wide approach, Nakano et al. identified an intergenic SNP, rs7081455, near PLXDC2 gene to be associated with POAG in the Japanese population . This finding was subsequently replicated by Mabuchi et al. in POAG patients of same ethnicity and also demonstrated to be associated with increasing IOP. Since POAG is the second most prevalent type of glaucoma in Saudi Arabia , we investigated an association between SNP rs7081455 near PLXDC2 gene and POAG or its endophenotype traits in a Saudi cohort.
Our study reports a negative association between rs7081455 and POAG or its glaucoma specific indices (endophenotypes) namely, IOP and cup/disc ratio. There have been conflicting reports on association between variant rs7081455 and POAG in the literature. Rs7081455 was associated with POAG in the Japanese population in two different reports [5, 6], but not by Takamoto et al.  in the same population in patients with normal tension glaucoma, substantiating an IOP-related effect of this variant on POAG as demonstrated by Mabuchi et al. . Similarly, there was no association reported in the Chinese , Afro-Caribbean  and Korean  populations. Besides, among the two reports from India, there appeared to be a regional difference, with no association among the South Indian as compared to modest association among the POAG patients from Eastern parts of India [12, 13]. The distribution of the risk ‘G’ allele frequency observed in our study and other ethnic groups discussed above are shown in Table 3. The ‘G’ allele frequency of rs7081455 was highest in our Saudi cohort and lowest among the Chinese population.
The PLXDC2 is cell surface trans-membrane protein receptor for pigment epithelium derived factor (PEDF)  that is widely expressed in many tissues including the eye and has been implicated to have a role in POAG [5, 6]. Though the exact mechanism(s) by which PLXDC2 is involved in POAG development and/or progression is not yet clear, however, the genetic variant rs7081455 near PLXDC2 seems to contribute to glaucomatous optic neuropathy in an IOP-related mechanism, because this locus was associated with high-tension glaucoma rather than normal-tension glaucoma . The authors hypothesized that the variant may affect expression and/or activity of PLXDC2 protein leading to different level of responsiveness to PEDF, thereby causing elevated IOP . Besides, Rogers et al.  has reported that PEDF decreases the conventional outflow pathway of the aqueous fluid by increasing the barrier function of the inner wall of Schlemm’s canal. Thus, this locus/variant may be associated with IOP-regulation by decreasing the outflow pathway via the PEDF-PLXDC2 receptor [6, 15]. In contrast, our study failed to observe this link as there was no significant association of this genotype with IOP or cup/disc ratio.
Validation or replication of association studies in different ethnic groups is of genetic epidemiological relevance to evaluate the future utility of genetic markers in disease per se. However, the findings observed in our study and those reported in other populations suggests that variant rs7081455 in PLXDC2 may not have a major role to play in the pathogenesis of the disease and hence, may not be an important genetic marker that may serve as a useful tool to assess the risk of POAG.
Our study reports a negative association between rs7081455 and POAG. However, the lack of significance in our study could be due to relatively small sample size and insufficient statistical power of the study. Based on the minor allele frequency observed in our Saudi cohort, the current sample size of the study has an estimated power of > 80% (with alpha risk of 5%) to detect an odds ratio of 2.0. However, it will certainly require a much larger sample size to detect a 1.5-fold relative risk, as is the case with most genetic association studies where polymorphism(s) exhibit a relatively small effect but may show a significant contribution to a large number of cases, emphasizing the need for further validation in a large population-based cohort. In addition, only a single genetic variant in PLXDC2 gene was studied, so the role of other rare variants in this gene or gene–gene interaction cannot be ruled out.
pigment epithelium derived factor
plexin domain containing 2
primary open angle glaucoma
retinal nerve fiber layer
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.
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.
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.
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.
Mabuchi F, Mabuchi N, Takamoto M, Sakurada Y, Yoneyama S, Kashiwagi K, Iijima H, Yamagata Z, Aihara M, Iwata T, et al. Genetic variant near PLXDC2 influences the risk of primary open-angle glaucoma by increasing intraocular pressure in the Japanese population. J Glaucoma. 2017;26(11):963–6.
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.
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.
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.
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.
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).
Vishal M, Sharma A, Kaurani L, Alfano G, Mookherjee S, Narta K, Agrawal J, Bhattacharya I, Roychoudhury S, Ray J, et al. Genetic association and stress mediated down-regulation in trabecular meshwork implicates MPP7 as a novel candidate gene in primary open angle glaucoma. BMC Med Genomics. 2016;9:15.
Cheng G, Zhong M, Kawaguchi R, Kassai M, Al-Ubaidi M, Deng J, Ter-Stepanian M, Sun H. Identification of PLXDC1 and PLXDC2 as the transmembrane receptors for the multifunctional factor PEDF. Elife. 2014;3:e05401.
Rogers ME, Navarro ID, Perkumas KM, Niere SM, Allingham RR, Crosson CE, Stamer WD. Pigment epithelium-derived factor decreases outflow facility. Invest Ophthalmol Vis Sci. 2013;54(10):6655–61.
AAK: Study and experimental design, data interpretation, statistical analysis, writing and critical revision of manuscript; TS: sample collection, DNA preparation, genotyping, results and data collection; FAM, HK: subject recruitment, clinical examination, clinical data results; KAA: study design, data interpretation and preparation of final version of the manuscript; SAA: study design, subject recruitment, clinical examination, clinical data results, overall supervision and preparation of final version of the manuscript. All authors read and approved the final manuscript.
The authors would like to thank the Deanship of Scientific Research and Glaucoma Research Chair at the King Saud University for supporting this study.
The authors declare that they have no competing interests and the work was not supported or funded by any drug company. The paper has not been presented in any previous conference or scientific meeting.
Availability of data
The data supporting the conclusions of this article are all presented within the article.
Consent for publication
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. Written, informed consent was obtained from all participants prior to their inclusion in this study.
This work was supported by Deanship of Scientific Research and Glaucoma Research Chair at the King Saud University, but had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kondkar, A.A., Sultan, T., Almobarak, F.A. et al. Plexin domain containing 2 (PLXDC2) gene polymorphism rs7081455 may not influence POAG risk in a Saudi cohort. BMC Res Notes 11, 733 (2018). https://doi.org/10.1186/s13104-018-3848-x