Polymorphism rs7555523 in transmembrane and coiled-coil domain 1 (TMCO1) is not a risk factor for primary open angle glaucoma in a Saudi cohort
© The Author(s). 2016
Received: 9 May 2016
Accepted: 22 September 2016
Published: 29 September 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.
KeywordsPOAG rs7555523 Saudi Arabia Transmembrane and coiled-coil domain 1
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.
Demographic and clinical characteristics of cases and control
Cases (N = 87) No. (%)
Controls (N = 94) No. (%)
Age in years, Mean (±SD)
Coronary artery disease
Family history of glaucoma
Awareness to glaucoma
Comparison of genotype and allele distribution of SNP rs7555523 in cases and controls
Cases (N = 87) No. (%)
Controls (N = 94) No. (%)
95 % Confidence Interval
Effect of genotypes on demographic and other clinical characteristics among POAG cases
(N = 69)
(N = 18)
Age in years, Mean (±SD)
Coronary artery disease
Health Awareness / Behavior
Family history of glaucoma
Awareness to glaucoma
Glaucoma Specific Indices
Intraocular pressure in mmHg, Mean (±SD)
Cup/disc ratio, Mean (±SD)
Number of anti-glaucoma medications, Mean (±SD)
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
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kwon YH, Fingert JH, Kuehn MH, Alward WL. Primary open-angle glaucoma. N Engl J Med. 2009;360:1113–24.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- Coleman AL, Miglior S. Risk factors for glaucoma onset and progression. Surv Ophthalmol. 2008;53 Suppl1:S3–10.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar