- Open Access
- Open Peer Review
Polymorphism rs7555523 in transmembrane and coiled-coil domain 1 (TMCO1) is not a risk factor for primary open angle glaucoma in a Saudi cohort
Journal of Negative Results in BioMedicinevolume 15, Article number: 17 (2016)
We investigated whether polymorphism rs7555523 (A > C) in human transmembrane and coiled-coil domain 1 (TMCO1) gene is a risk factor for primary open angle glaucoma (POAG) in a Saudi cohort.
A cohort of 87 unrelated POAG cases and 94 control subjects from Saudi Arabia were genotyped using Taq-Man® assay. The association of genotypes with POAG and other glaucoma specific clinical indices was investigated.
The genotype and allele frequency of polymorphism rs7555523 at TMCO1 did not show any statistically significant association with POAG as compared to controls. The minor allele frequency was 0.103 in cases and 0.085 in controls. Except for awareness of glaucoma (p = 0.036), no significant association of genotypes were seen with glaucoma specific clinical indices such as intraocular pressure (IOP), cup/disc ratio and number of anti-glaucoma medications used. Binary logistic regression analysis (adjusted for age and gender) showed that age was a significant indicator for the development of glaucoma in this group (adjusted odds ratio = 1.2; 95 % confidence interval = 1.078–1.157; p < 0.001).
Our study was unable to replicate the findings of previously reported association for polymorphism rs7555523 in TMCO1 with POAG and related clinical indices such as IOP and cup/disc ratio indicating that this variant is not a risk factor for POAG in the Saudi cohort.
Glaucoma, a neurodegenerative disease is characterized by progressive damage of retinal ganglion cells (RGCs) resulting in characteristic cupping of the optic nerve head and loss of peripheral vision . Primary open angle glaucoma (POAG) is the second most common form of glaucoma in Saudi Arabia which is clinically characterized by an open and normal anterior iridocorneal chamber angle . Aging, gender, African ancestry, family history, elevated intraocular pressure (IOP), central corneal thickness, and myopia are some of the well-recognized risk factors associated with POAG pathogenesis . Although POAG is clinically well-defined, the biological basis of the disease is not well understood and factors contributing to its progression are not fully characterized.
Genetic studies represent an important tool to identify genes and molecular pathways involved in disease pathogenesis. POAG is genetically complex with largely polygenic and multifactorial inheritance . Using a genome-wide and candidate gene approach, population-based genetic studies have identified several genes and genetic variants associated with POAG and related quantitative endophenotype traits . A genome-wide association study (GWAS) in Australians of European descent identified a susceptibility locus at transmembrane and coiled-coil domains 1 (TMCO1) . van Koolwijk et al. performed a GWAS for IOP in POAG patients of European descent and identified single nucleotide polymorphism (SNP) 7555523, located in TMCO1 suggesting a role in IOP regulation . TMCO1 gene is located 7.6 MB upstream of the known POAG gene, myocilin C (MYOC), on chromosome 1q24.1. It encodes a transmembrane protein with a coiled-coil domain that may localize to the Golgi apparatus and endoplasmic reticulum or to the mitochondria in different cell types with a plausible role in apoptosis of RGCs [6, 8]. TMCO1 is highly expressed in the human ciliary body, trabecular meshwork and retina [7, 9, 10]. However, precise role of TMCO1 in POAG pathogenesis is still unclear.
So far, there are no published reports of association studies at the TMCO1 locus in the Middle East population. In the current study, we investigated the association of TMCO1 SNP rs7555523 with POAG in a Saudi cohort.
Study design and setting
This case–control genetic association study was conducted between November 2015 through February 2016 at King Abdulaziz University Hospital, King Saud University, Riyadh, Saudi Arabia.
We recruited 87 Saudi adult-POAG patients who satisfied strict the following clinical criteria: i) appearance of the disc or retinal nerve fiber layer e.g., thinning or notching of disc rim, progressive changes, nerve fiber layer defect; ii) presence of characteristic abnormalities in visual field (e.g., arcuate scotoma, nasal step, paracentral scotoma, generalized depression) in the absence of other causes or explanation; iii) age >40 years, and iv) open anterior chamber angles bilaterally on gonioscopy. Exclusion criteria included evidence of secondary glaucoma, e.g., pigmentary dispersion syndrome, pseudoexfoliation, history of steroid use, or ocular trauma. All cases had onset of glaucoma after age 40 (adult onset POAG). Patients were recruited from the glaucoma clinic at King Abdulaziz University Hopsital after signing an informed consent approved by the institutional review board (proposal number # 08–657). A second group healthy controls (n = 94) of Saudi origin and free from glaucoma by examination were recruited. Inclusion criteria included: >40, normal IOP, open angles on gonioscopy and normal optic nerves upon examination.
Genotyping of intronic polymorphism, rs7555523 (g.165718979A > C), of the TMCO1 gene (NC_000001.10) was performed using the TaqMan® SNP Genotyping assay ID: C_29621671_10 (Applied Biosystems Inc., Foster City, CA, USA) on ABI 7500 real-time PCR system (Applied Biosystems). Each PCR reaction was performed in a 96-well plate in a total volume of 25 μL consisting of 1X TaqMan® Genotyping Master Mix (Applied Biosystems), 1X SNP Genotyping Assay Mix, 20 ng DNA, and two no template (negative) controls under cycling conditions recommended by the manufacturer . Genotypes of TMCO1 rs7555523 SNP were identified using the automated 2-color allele discrimination software on ABI 7500 on a two-dimensional graph.
The continuous variables were presented as mean (± Standard Deviation, SD) and tested by Student’s t-test. Categorical variables were presented as frequencies and percentages. Hardy-Weinberg Equilibrium (HWE) deviation was tested by Pearson’s Chi2 test. Odds ratio (OR) was calculated and Chi2 test was used to detect any association between different characteristics and the genetic profiles (Fisher Exact test when applicable). Mann–Whitney U test was used to investigate whether there was any significant difference between the genotypes and clinical variables. The confidence interval (CI) level was set to 95 % and a p value below 0.05 was considered statistically significant. Data were analyzed using SPSS® version 22.0 (IBM Inc., Chicago, Illinois, USA).
In the current study we recruted 87 adult-POAG patinets with confirmed diagnosis as POAG and a matching group of 94 subjects that served as controls after a confirmed clinical examination that identified them as “glaucoma free”, see methods.
As shown in Table 1 cases had a mean age of 61.1 years (ranging from 43 to 74 years) where 52 (59.8 %) of them were male and 35 (40.2 %) were female. On the other hand, controls showed a mean age of 56.5 years (range 45–70 years) where 69 (73.4 %) of them were male and 25 (26.6 %) were female. Despite the fact that males were more than females in both study groups, this diffrence was not statastcalliy significant. No significant difference was observed between cases and controls in terms of all demographic, systemic co-morbidity and glaucoma specific indices except for the family history of glaucoma and awareness of having glaucoma (p = 0.006 and <0.0001, respectively).
The genotype frequencies in both the cases and control groups did not deviate significantly from the HWE (p > 0.05). The wildtype genotype (A/A) was predominant in both cases and controls (n = 69 (79.3 %) and n = 78 (83 %), respectively) with a slightly increased frequency of heterozygous (C/A) genotype in cases (18; 20.7 %) versus controls (16; 17 %). No homozygous mutant genotype (C/C) was observed in both the groups. The cases were 1.3 times more likely to encounter a variation, however, the distribution was non-significant (OR = 1.3; 95 % CI = 0.563–2.891; p = 0.571). Likewise, the same effect was detected while comparing the frequency of the wildtype “A” allele to the mutated “C” allele, with a similar OR (1.3) and a non-significant p value (0.591). The genotype and allele frequency distribution is shown in Table 2.
Furthermore, we evaluated the effect of genotype on demographic, systemic diseases and glaucoma specific indices among POAG cases as demonstrated in Table 3. Although there was a preponderance of male subjects in both the genotype groups and the subjects with A/C genotypes were found to be slighter younger but there was no statistically significant difference in terms of age (p = 0.644) and gender (p = 0.421). Similarly, except for awareness of having glaucoma variable (p = 0.036), no statistically significant difference was observed in terms of systemic diseases and health awareness/ behavior characteristics. However, more importantly, none of the glaucoma specific indices such as IOP, cup/disc ratio and number of anti-glaucoma medications showed any statistically significant difference between the two genotype groups.
In addition, to investigate the effect of harboring a mutated genotype on having glaucoma, we performed a binary logistic regression analysis (adjusted for age and gender). The analysis showed that patients with mutated genotype seems to be 1.3 times more likely to get the disease (POAG), however, the OR was not found to be statistically significant (adjusted OR = 1.3; 95 % CI = 0.534–3.261; p = 0.548). Nevertheless, adjustment of age and gender revealed that although non-significant, females were 1.6 times more likely to get glaucoma than males (adjusted OR = 1.6; 95 % CI = 0.778–3.428; p = 0.195); and not surprisingly, age was found to be a significantly strong indicator for the development of glaucoma in this group (adjusted OR = 1.2; 95 % confidence interval = 1.078–1.157; p < 0.001).
Given the complexity and genetic mutational heterogeneity of POAG recent GWASs have identified a number of polymorphisms in multiple loci/genes including caveolin (CAV1/CAV2) , atonal homolog 7 (ATOH7) , sin oculis homeobox (SIX1/SIX6) , cyclin-dependent kinase inhibitor 2B antisense RNA 1 (CDKN2B-AS1) [6, 14] and TMCO1 [6, 14] that may contribute to the development and/or progression of POAG in various ethnic groups. In this study, we investigated whether SNP rs7555523 (A > C) in TMCO1 gene is a risk factor for POAG in a Saudi cohort.
Genetic variation in TMCO1 has been associated with POAG [6, 15]. The association of the TMCO1 locus with POAG has been replicated in another GWAS . Similarly, TMCO1 loci (including rs7555523) showed significant association with POAG and high-tension glaucoma in a Han Chinese population . However, conflicting results have been observed for populations of African ancestry [17, 18]. The frequency distribution of polymorphism(s) is known to vary significantly across the different ethnic groups. SNP rs7555523 has been reported to have a minor allele frequency (MAF) of 0.12 in Caucasians , 0.016 in Han Chinese  and 0.132 in Sub-Saharan African [http://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=7555523]. The MAF observed in our POAG Saudi cohort was 0.103 which is slightly lower than the Caucasians and Africans, but much higher than the Han Chinese (Asians). However, in contrast to the Caucasian and Chinese studies, the genotype and allele frequency of rs7555523 were not found to be an independent risk factor of POAG in our cohort indicating that this SNP may not have a significant role in Saudi POAG as compared to Caucasians and Chinese.
Family history, aging, cigarette smoking, diabetes and hypertension are well recognized risk factors of POAG [19, 20]. Sharma et al. has shown that POAG patients carrying the risk allele of SNP rs4656461 near the TMCO1 gene tend to have an earlier age at diagnosis of glaucoma . However, in our study, although the subjects with mutant heterozygous (A/C) genotypes were found to be slighter younger but the difference was statistically non-significant. Likewise, SNP rs7555523 is located in a region which has been previously suggested to be linked with blood pressure . In our study, however, none of the systemic diseases including hypertension were found to be significantly associated with this genotype.
TMCO1 is a highly evolutionary conserved gene of largely unknown function [8, 22]. A homozygous frameshift mutation in TMCO1 has been associated with a rare recessive syndrome known as “TMCO1 defect syndrome” consisting of craniofacial dysmorphism, skeletal anomalies and mental retardation . It is still unclear how this gene contributes to the pathogenesis of glaucoma. Studies suggest that TMCO1 may contribute to POAG through the pathway of IOP elevation [7, 16]. However, we did not find any significant association between heterozygous (A/C) genotype and clinical indices important for glaucoma such as IOP, cup-to-disc ratio and number of anti-glaucoma medications. It is noteworthy to mention here that TMCO1 interacts with another known POAG susceptibility gene, CAV1 via the von Hippel-Lindau (VHL) tumor suppressor protein-containing E3 ubiquitin ligase [7, 23]; and we have previously shown that SNP rs4236601 in CAV1/CAV2 is not a risk factor for POAG in Saudi population , plausibly suggesting that TMCO1 may not have an important role in POAG pathogenesis in this population. However, considering the small sample size in this study and the fact that there was no homozygous mutant (C/C) genotype observed in our sample population these observations may require further validation in a large sample population.
Our study was unable to replicate the findings of previous association reported for variant rs7555523 in TMCO1 with POAG and important clinical indices such as IOP and cup/disc ratio indicating that this SNP is not a risk factor for POAG or its important clinical indices in the Saudi cohort.
Atonal homolog 7
Cyclin-dependent kinase inhibitor 2B antisense RNA 1
Genome-wide association study
Minor allele frequency
Primary open angle glaucoma
Retinal ganglion cells
- SIX1/SIX6 :
Sin oculis homeobox
Single nucleotide polymorphism
- TMCO1 :
Transmembrane and coiled-coil domain 1
Kwon YH, Fingert JH, Kuehn MH, Alward WL. Primary open-angle glaucoma. N Engl J Med. 2009;360:1113–24.
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:373–9.
Kang JH, Loomis SJ, Rosner BA, Wiggs JL, Pasquale LR. Comparison of risk factor profiles for primary open-angle glaucoma subtypes defined by pattern of visual field loss: a prospective study. Invest Ophthalmol Vis Sci. 2015;56:2439–48.
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.
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–911.
Burdon KP, Macgregor S, Hewitt AW, Sharma S, Chidlow G, Mills RA, Danoy P, Casson R, Viswanathan AC, Liu JZ, et al. Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1. Nat Genet. 2011;43:574–8.
van Koolwijk LM, Ramdas WD, Ikram MK, Jansonius NM, Pasutto F, Hysi PG, Macgregor S, Janssen SF, Hewitt AW, Viswanathan AC, et al. Common genetic determinants of intraocular pressure and primary open-angle glaucoma. PLoS Genet. 2012;8:e1002611.
Zhang Z, Mo D, Cong P, He Z, Ling F, Li A, Niu Y, Zhao X, Zhou C, Chen Y. Molecular cloning, expression patterns and subcellular localization of porcine TMCO1 gene. Mol Biol Rep. 2010;37:1611–8.
Liton PB, Luna C, Challa P, Epstein DL, Gonzalez P. Genome-wide expression profile of human trabecular meshwork cultured cells, nonglaucomatous and primary open angle glaucoma tissue. Mol Vis. 2006;12:774–90.
Sharma S, Burdon KP, Chidlow G, Klebe S, Crawford A, Dimasi DP, Dave A, Martin S, Javadiyan S, Wood JP, et al. Association of genetic variants in the TMCO1 gene with clinical parameters related to glaucoma and characterization of the protein in the eye. Invest Ophthalmol Vis Sci. 2012;53:4917–25.
Abu-Amero KK, Kondkar AA, Mousa A, Osman EA, Al-Obeidan SA. Association of Mn-SOD mutation (c.47 T > C) with various POAG clinical indices. Ophthalmic Genet. 2014;35:85–90.
Thorleifsson G, Walters GB, Hewitt AW, Masson G, Helgason A, DeWan A, Sigurdsson A, Jonasdottir A, Gudjonsson SA, Magnusson KP, et al. Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma. Nat Genet. 2010;42:906–9.
Ramdas WD, van Koolwijk LM, Lemij HG, Pasutto F, Cree AJ, Thorleifsson G, Janssen SF, Jacoline TB, Amin N, Rivadeneira F, et al. Common genetic variants associated with open-angle glaucoma. Hum Mol Genet. 2011;20:2464–71.
Wiggs JL, Yaspan BL, Hauser MA, Kang JH, Allingham RR, Olson LM, Abdrabou W, Fan BJ, Wang DY, Brodeur W, et al. Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma. PLoS Genet. 2012;8:e1002654.
Gibson J, Griffiths H, De Salvo G, Cole M, Jacob A, Macleod A, Yang Y, Menon G, Cree A, Ennis S, Lotery A. Genome-wide association study of primary open angle glaucoma risk and quantitative traits. Mol Vis. 2012;18:1083–92.
Chen Y, Hughes G, Chen X, Qian S, Cao W, Wang L, Wang M, Sun X. Genetic variants associated with different risks for high tension glaucoma and normal tension glaucoma in a Chinese population. Invest Ophthalmol Vis Sci. 2015;56:2595–600.
Williams SE, Carmichael TR, Allingham RR, Hauser M, Ramsay M. The genetics of POAG in black South Africans: a candidate gene association study. Sci Rep. 2015;5:8378.
Liu Y, Hauser MA, Akafo SK, Qin X, Miura S, Gibson JR, Wheeler J, Gaasterland DE, Challa P, Herndon LW, et al. Investigation of known genetic risk factors for primary open angle glaucoma in two populations of African ancestry. Invest Ophthalmol Vis Sci. 2013;54:6248–54.
Fan BJ, Leung YF, Wang N, Lam SC, Liu Y, Tam OS, Pang CP. Genetic and environmental risk factors for primary open-angle glaucoma. Chin Med J (Engl). 2004;117:706–10.
Coleman AL, Miglior S. Risk factors for glaucoma onset and progression. Surv Ophthalmol. 2008;53 Suppl1:S3–10.
Ehret GB, O’Connor AA, Weder A, Cooper RS, Chakravarti A. Follow-up of a major linkage peak on chromosome 1 reveals suggestive QTLs associated with essential hypertension: GenNet study. Eur J Hum Genet. 2009;17:1650–7.
Xin B, Puffenberger EG, Turben S, Tan H, Zhou A, Wang H. Homozygous frameshift mutation in TMCO1 causes a syndrome with craniofacial dysmorphism, skeletal anomalies, and mental retardation. Proc Natl Acad Sci U S A. 2010;107:258–63.
Xie L, Xue X, Taylor M, Ramakrishnan SK, Nagaoka K, Hao C, Gonzalez FJ, Shah YM. Hypoxia-inducible factor/MAZ-dependent induction of caveolin-1 regulates colon permeability through suppression of occludin, leading to hypoxia-induced inflammation. Mol Cell Biol. 2014;34:3013–23.
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.
The authors would like to thank the Glaucoma Research Chair at the Department of Ophthalmology, College of Medicine, King Saud University for funding this study.
Availability of data and materials
The data supporting the conclusions of this article are all presented within the article.
AAK: Designed experiments, results, interpretations, writing manuscript draft; AM: statistical analyses and results; TAA, TS, SA: DNA extraction and genotyping; AA, SAA: subject recruitment, clinical examination, clinical data; KAA: study design, overall supervision, results, interpretation, and preparation of final version of the manuscript. All the authors have approved the final version of the manuscript.
Authors have no conflict of 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.
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. Written, informed consent was obtained from all participants prior to their inclusion in this study.