SNP located in an AluJb repeat downstream of TMCO1, rs4657473, is protective for POAG in African Americans.

AIMS: To determine the association of single nucleotide polymorphisms (SNPs) downstream from the TMCO1 gene with primary open-angle glaucoma (POAG) in African Americans (AA).
METHODS: AA subjects were recruited for the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study from the Scheie Eye Institute and its satellite sites in Philadelphia. A region containing an AluJb repeat and seven SNPs, including rs4656461 near the TMCO1 gene, were PCR-Sanger sequenced from POAAGG cases (n=1537) and controls (n=1570). Association between POAG and SNPs near TMCO1 was investigated by logistic regression analysis. Phenotypic trait associations with these SNPs were assessed by analysis of variance. Electrophoretic mobility shift assay (EMSA) was performed to assess the affinity of human T-box 5 (TBX5) protein for a predicted binding motif in the TMCO1 region. Dual Luciferase assays were performed by transfecting recombinant plasmids containing the region surrounding the above SNPs in HEK293T and trabecular meshwork cells.
RESULTS: The SNP rs4657473 (C>T) was associated with POAG; the TT genotype was protective (OR 0.20, 95% CI 0.09 to 0.42; p<0.001). No significant associations were found between the TMCO1 variants and phenotypic traits. EMSA confirmed the affinity of TBX5 for a predicted binding motif containing TMCO1 SNP rs4657475. Luciferase assays demonstrated a regulatory function for the genomic region around SNP rs4656561, located within AluJb repeat.
CONCLUSION: Our results demonstrate that a SNP downstream of TMCO1, rs4657473, is associated with POAG in an AA population. Our studies suggest a regulatory role for the previously POAG-associated locus near the TMCO1 gene that may affect gene expression. © Author(s) (or their employer(s)) 2019. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Introduction

Glaucoma is a complex neurodegenerative disease and the leading cause of irreversible blindness worldwide.1 Primary open-angle glaucoma (POAG) is a progressive degeneration of retinal ganglion cells and their axons in the optic nerve that leads to characteristic changes in the optic disc and retinal nerve fibre layer (RNFL) and corresponding visual field loss.2 Risk factors for POAG include elevated intraocular pressure (IOP), advanced age, family history and African American (AA) race; however, the exact aetiology of POAG remains unknown.3 In recent years, family linkage and genome-wide association studies (GWAS) have made increasingly clear the importance of genetic risk factors in the development and progression of glaucoma.4

The TMCO1 region has been identified as a likely contributor to the genetic risk for POAG. A single nucleotide polymorphism (SNP), rs4656461, located at the 1q24 locus in a non-coding region ~6.5 kb downstream to TMCO1 was linked (p=6.0×10–14) to severe POAG-induced visual field loss in a GWAS of a Caucasian cohort.5 The same TMCO1 locus was associated with POAG in case-control association studies in Pakistani and Han Chinese patients with high tension glaucoma (HTG).6 7 The SNP rs755523, an intronic variant of TMCO1, has also been associated with HTG in the same Han Chinese population and with high IOP in a GWAS of Caucasians (p=1.6×10-–8).7

Few studies have investigated the genetic risk factors of POAG in AA, despite this population having a higher disease prevalence than other populations.5 8 9 A candidate gene case-control study in AAs (1150 cases and 999 controls) and Ghanaian (West African) subjects (483 cases and 593 controls) was performed for eight SNPs in the TMCO1 region including rs4656461. No significant associations were found for TMCO1 SNPs or other loci that have been previously associated with POAG in populations of non-African ancestry, suggesting that these susceptibility alleles may not be relevant in AAs.10

We investigated whether SNP rs4656461 and other SNPs downstream from TMCO1 (in an Alu region) are associated with POAG in a large AA cohort. We also examined the impact of TMCO1 SNPs in this Alu region on predicted TBX5 transcription factor binding and interrogated these for a regulatory role.

Methods

Patient recruitment and data collection

Subjects for our case-control study were recruited from the Scheie Eye Institute, University of Pennsylvania and its research affiliate sites in Philadelphia for the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study. Information regarding POAAGG study patient recruitment and data collection and the baseline demographics of the 2520 subjects enrolled as of March 2014 has been previously described.11 In brief, subjects over 35 years who self-identified as Black were enrolled. We defined POAG as having an open iridocorneal angle and at least one characteristic glaucomatous optic nerve finding (notching, neuroretinal rim thinning, excavation or a nerve fibre layer defect) that corresponded to a visual field defect. All secondary causes of glaucoma were excluded. Informed consent was obtained from all participating subjects. Approval was obtained from the University of Pennsylvania institutional review board and research was conducted in accordance with the Declaration of Helsinki.

Genotyping and association analysis

Sanger sequencing was performed on genomic DNA from blood and saliva samples from 1537 AA cases and 1570 controls to genotyp e seven SNPs downstream of TMCO1 (rs11430311, rs4657473, rs4657474, rs4657475, rs116461237, rs61058666 and rs4656461) as previously described.12 All the seven candidate SNPs in our study were identified from sequencing around 1 kb of the intergenic region between TMCO1 and ALDH9A1 surrounding the previously associated POAG GWAS SNP, rs4656461 in both AA cases and controls. Briefly, BigDye Terminator chemistry (ThermoFisher, Waltham, Massachusetts, USA) was used, followed by capillary electrophoresis on an ABI 3130xl genetic analyzer (Applied Biosystems, Foster City, California, USA). Linkage disequilibrium (LD) patterns of the TMCO1 SNPs were plotted using LDlink.13

Association between disease status and the genotype frequencies and allele frequencies of the seven SNPs was assessed using logistic regression analysis to calculate ORs with 95% CI while adjusting for age and sex. Correction for multiple comparisons was performed using a Bonferroni correction, requiring a p≤0.007 to be considered as statistically significant after accounting for the testing of seven SNPs. Information regarding the quantitative traits of glaucoma, including central corneal thickness (CCT), IOP, cup-to-disc ratio (CDR), mean deviation (MD), pattern standard deviation (PSD) and RNFL thickness, was collected from the electronic medical records for POAG cases only. We did not collect information about quantitative traits from controls. Associations between each SNP and the quantitative traits among POAG cases were performed using analysis of variance (ANOVA). In these analyses, when measures from left and right eyes were available, averaged values for CCT, maximum-recorded IOP, CDR, MD, PSD and RNFL thickness were calculated and used for statistical comparisons. All statistical analyses were conducted using SAS V.9.4 (SAS Institute, Cary, North Carolina, USA).

Electrophoretic mobility shift assay (EMSA) binding studies

The seven SNPs near TMCO1 from the genotype-phenotype study were queried in RegulomeDB to determine whether any lie within predicted transcription factor binding sites.14 Three of the SNPs were predicted to lie within a TBX5 binding motif. EMSA was performed using Biotinylated oligonucleotide probes (online supplementary table 1) and purified recombinant protein of human T-box 5 (TBX5) with C-terminal MYC/DDK tag (Origene, Rockville, MD; TP316520) and the Odyssey Infrared EMSA kit (LiCor, Lincoln, Nebraska; 829–07910). The binding reactions were assembled according to the manufacturer’s protocol. TBX5 protein was incubated with the reaction mixture for 10 min at 4°C followed by incubation with the DNA probes for an additional 30 min. For the competition assay, 2000-fold excess of unlabeled oligonucleotides was added to the reaction mixture before the addition of the labelled oligonucleotide probe. The DNA-protein complexes were electrophoresed through a TBE gel on an ice bath, visualised with a 1:10 000 dilution of IR dye in PBS and 0.1% Tween 20 and analysed on an Odyssey infrared fluorescence scanner (LiCor Biosciences, Lincoln, Nebraska, USA).

Cloning of the TMCO1 SNP region and generation of luciferase reporter constructs

Luciferase reporter constructs were generated by cloning a 3337 bp region of genomic DNA surrounding the TMCO1 SNP region (including rs4656461) from four AA individuals that represented naturally occurring haplotypes (figure 1A). Genomic DNA samples were selected for cloning based on the prior sanger sequencing of POAAGG subjects. Haplotype construct 1_3337 contained wildtype alleles at all seven TMCO1 SNPs (rs11430311, rs4657473, rs4657474, rs4657475, rs116461237, rs61058666 and rs4656461). The haplotypes associated with the SNP rs4657473 variant (C>T), rs116461237 variant (G>A) and SNP rs4656461 variant (A>G) were associated with constructs 2, 3 and 4_3337 respectively. The construct 4_3337 also had variants at SNPs rs11430311, rs4657474 and rs4657475, which are in LD with rs4656461.

Figure 1

(A) Haplotypes of inserts used in luciferase assay. (B) Luciferase assay was performed with inserts of three sizes surrounding SNP rs4656461 and Alu element. Inserts were formed by serial deletion of the initial 3337 bp construct. 703 bp insert was cloned in both a forward and reverse direction. SNP, single nucleotide polymorphism.

PCR amplification was performed using Phusion High-Fidelity DNA polymerase (NEB) and primers containing Acc65I and EcoRV restriction enzyme sites (online supplementary table 2) to generate the 3337 bp sequences containing TMCO1 SNPs. Each PCR construct was ligated into a vector upstream of a minimal promoter (pGL4.23[luc2/minP]; Promega). The resulting plasmids were transfected into 10-beta Competent E. coli (High Efficiency) cells (NEB) and plasmid DNA was quantitated using Qubit Fluorometer 2.0 (ThermoFisher).

In order to identify the minimal regulatory region surrounding the SNPs, serial deletion constructs were developed by selective amplification of the 3337 bp plasmids containing the desired SNPs using primer sets mentioned in online supplementary table 2 to generate 1274 bp and 703 bp constructs, which were again cloned into pGL4.23 vectors (figure 1). The first number for each of the luciferase clones indicates the set of SNPs in the clone (figure 1A) and the second number indicates the size of the clone (figure 1B). All recombinant plasmids were verified by Sanger sequencing. Nearly 300–400 bp of insert sequence surrounding SNP rs4656461 was compared with previous genotyping.

Dual luciferase reporter assay

Transcriptional regulatory activity of all inserts was interrogated by dual luciferase reporter assay in HEK293T cells (human embryonic kidney ATCC-CRL3216) and primary trabecular meshwork (TM) cells (gift from Dr Markus Kuehn, University of Iowa). To validate the authenticity of the TM cells, we performed the MYOC induction experiment with DEX in the TM cells. We used Passage 2 (P2) to Passage 4 (P4) primary TM cultures for all our experiments. To characterise the TM cultures, we induced 65%–70% confluent P4 TM cultures with 100 mM dexamethasone (Sigma) for 3 days. The RNA was isolated from the induced and uninduced TM cultures and was converted to cDNA using Superscript III One-Step RT-PCR kit (Invitrogen, CA) and a quantitative Real-time PCR was performed using primers (5'-CATCTGGCTATCTCAGGAGTGG-3'; Rev-5'- TCTCGCATCCACACACCATAC-3') to amplify Myocilin transcript and the gene expression was compared with GAPDH. We found that there was a 7-fold increase in the Myocilin expression in TM cells induced with Dexamethasone when compared with uninduced controls (online supplementary figure S1).

Both HEK293T and TM cells were plated in a 24-well format and cultured at 37°C with 5% CO2. Cells were seeded on a 96-well white opaque bottomed luminometric plates (Falcon) at ~60% confluency 24 hours prior to transfection. Cells were transiently transfected with reporter plasmids using Lipofectamine 3000 following manufacturer’s instructions (ThermoFisher). Renilla luciferase vector (pGL4.74[hRluc/TK]; Promega) was cotransfected with experimental plasmids to control for transfection efficiency. Luciferase activity was measured following protocol from Pierce Renilla-Firefly Luciferase Dual Assay Kit (Thermo Scientific) on an Infinite 200 PRO plate reader (Tecan). Three replicates were analysed for each transfection set, and the ratio of Firefly luciferase activity to Renilla luciferase activity was used to determine the transfection efficiency. Plasmids were compared by one-way ANOVA and student’s t-test for differential luciferase activity. In HEK293T cells, three replicates of 3337 bp and 703 bp_Rev inserts were tested in triplicate transfection sets, while 1274 bp and 703 bp inserts were performed as a single transfection set. In TM cells, all inserts were tested in duplicate transfection sets.

Results

TMCO1 genotype association with glaucoma phenotype

The present case-control study included 1537 AA POAG cases and 1570 AA controls. Patients with POAG had a mean age of 70.9 years compared with 61.5 years in the controls (p<0.001). 60.4% of the POAG cases were female, while 69.4% of the controls were female (p<0.001). Demographic and clinical characteristics of the cohort are shown in table 1.

Table 1

Demographic and clinical characteristics of the cohort

		Glaucoma status
		Case (n=1537)	Control (n=1570)	P value
Female, %	F	929 (60.4%)	1090 (69.4%)	<0.001
Age, years	n	1537	1569	<0.001
	Mean (SD)	70.9 (11.3)	61.5 (11.8)
Maximum-recordedIOP (mean), mm Hg	n	1535	1515	<0.001
	Mean (SD)	25.9 (8.0)	17.1 (3.2)
CDR (mean)	n	1494	1247	<0.001
	Mean (SD)	0.71 (0.16)	0.33 (0.13)
CCT (mean)	n	1470	–	–
	Mean (SD)	531.0 (38.8)	–
RNFL (mean)	n	1096	–	–
	Mean (SD)	72.9 (13.3)	–
PSD (mean)	n	1106	–	–
	Mean (SD)	5.1 (3.1)	–
MD (mean)	n	1106	–	–
	Mean (SD)	−7.4 (7.8)	–

CCT, central corneal thickness; CDR, cup-to-disc ratio; IOP, intraocular pressure; MD, mean deviation; PSD, pattern standard deviation; RNFL, retinal nerve fibre layer.

Since our cohort enrolment was done by self-reporting as Black, we have performed a PCA of the cases and controls from our study to assess our cohorts for outliers in heredity. Boxplots results showing the comparison of the PC1 score between the cases and controls showed no statistically significant difference (p=0.06) (online supplementary figure S2). Genotyping efficiency was ≥99.8% for the seven SNPs. LD analysis revealed that SNPs rs11430311, rs4657474 and rs4657475 are in LD with SNP rs4656461 (figure 2B). Genotype frequencies and p values for the association of seven TMCO1 SNPs in POAG cases and controls are shown in table 2.

Figure 2

(A) Schematic showing POAG-associated SNP rs4656461 and six SNPs in an adjacent Alu element lie in an intergenic region between genes TMCO1 and ALDH9A1. (B) LD relationships are shown for TMCO1 SNPs. SNP rs4656461 is in LD with rs11430311, rs4657474 and rs4657475. LD, linkage disequilibrium; SNP, single nucleotide polymorphism.

Table 2

Genotype frequencies for 7 TMCO1 SNPs and association with POAG

SNP	Genotypesalleles	Cases (%)	Controls (%)	Hardy-Weinberg	Age and sexadjustedOR (95% CI)	P value
rs11430311	InsT/insT	925 (60.3)	938 (59.8)	0.99	Reference	0.76
	InsT/WT	533 (34.8)	546 (34.8)		0.99 (0.84 to 1.17)
	WT/WT	75 (4.9)	84 (5.4)		0.88 (0.62 to 1.25)
	InsT	2383 (77.7)	2422 (77.2)		Reference	0.60
	WT	683 (22.3)	714 (22.8)		0.97 (0.85 to 1.10)
rs4657473	CC	1251 (81.4)	1272 (81.0)	0.02	Reference	<0.001
	CT	275 (17.9)	265 (16.9)		1.05 (0.86 to 1.29)
	TT	10 (0.7)	33 (2.1)		0.20 (0.09 to 0.42)
	C	2777 (90.4)	2809 (89.5)		Reference	0.08
	T	295 (9.6)	331 (10.5)		0.85 (0.71 to 1.02)
rs4657474	GG	926 (60.3)	941 (59.9)	0.76	Reference	0.81
	GA	537 (35.0)	549 (35.0)		1.00 (0.85 to 1.18)
	AA	73 (4.8)	80 (5.1)		0.89 (0.62 to 1.27)
	G	2389 (77.8)	2431 (77.4)		Reference	0.70
	A	683 (22.2)	709 (22.6)		0.97 (0.86 to 1.11)
rs4657475	CC	926 (60.3)	941 (60.0)	0.65	Reference	0.81
	CT	538 (35.0)	549 (35.0)		1.00 (0.85 to 1.18)
	TT	72 (4.7)	79 (5.0)		0.89 (0.62 to 1.28)
	C	2390 (77.8)	2431 (77.5)		Reference	0.72
	T	682 (22.2)	707 (22.5)		0.98 (0.86 to 1.11)
rs116461237	GG	1492 (97.2)	1527 (97.3)	0.60	Reference	0.68
	GA	43 (2.8)	41 (2.6)		1.03 (0.64 to 1.64)
	AA	0 (0.0)	1 (0.1)		N/A
	G	3027 (98.6)	3095 (98.6)		Reference	0.98
	A	43 (1.4)	43 (1.4)		0.99 (0.63 to 1.57)
rs61058666	GG	1336 (87.0)	1355 (86.4)	0.47	Reference	0.95
	GC	192 (12.5)	209 (13.3)		1.01 (0.81 to 1.27)
	CC	7 (0.5)	5 (0.3)		1.21 (0.34 to 4.27)
	G	2864 (93.3)	2919 (93.0)		Reference	0.84
	C	206 (6.7)	219 (7.0)		1.02 (0.82 to 1.27)
rs4656461	AA	925 (60.3)	940 (59.9)	0.65	Reference	0.81
	AG	538 (35.0)	549 (35.0)		1.01 (0.86 to 1.18)
	GG	72 (4.7)	79 (5.0)		0.89 (0.62 to 1.28)
	A	2388 (77.8)	2429 (77.5)		Reference	0.75
	G	682 (22.2)	707 (22.5)		0.78 (0.86 to 1.11)

*P values for any difference between three genotypes, and linear difference by number of alleles, respectively.

POAG, primary open-angle glaucoma; SNP, single nucleotide polymorphism.

We observed a significant association of POAG with SNP rs4657473, with a protective effect shown by the TT variant at this locus (OR 0.20, 95% CI 0.09 to 0.42; p<0.001). The allele association analysis found no statistically significant association between any of the seven SNPs with POAG (table 2). The Hardy-Weinberg (HW) equilibrium p values, and the logistic regression p values for the allele frequency are also shown in table 2. Interestingly, the rs4657473 SNP is out of HW equilibrium, but the other SNPs are agreeing with it. These results indicate that the population itself is not at disequilibrium, but rather that this SNP may be a specific newer mutation, which is not inconsistent with it being associated with case status. The p value for the trend is also not significant, this is explained because the CT group is almost identical to the CC group, suggesting that only that the homozygous TT is necessary to be associated with the disease. No significant association was found for SNP rs4656461 or for any of the seven TMCO1 SNPs with the quantitative traits of POAG among the case subjects (online supplementary table 3).

EMSA binding studies

EMSA studies were performed to assess the binding of transcription factor TBX5 to its predicted binding motif encompassing 3 TMCO1 SNPs, rs4657475 (in LD with rs4656461), rs116461237 and rs61058666. EMSA showed that TBX5 formed complexes with the probes containing the three wildtype TMCO1 SNPs (figure 3). Introduction of variant rs116461237 (G>A) disrupted TBX5 binding. However, variants at rs4657475 and rs6105866 did not interfere with complex formation even when random alleles were introduced at these positions. Notably, TBX5 did not bind to the Alu consensus repeat motif, indicating that the affinity of this Alu region for TBX5 is an acquired trait.

Figure 3

(A) Three TMCO1 SNPs fall within a predicted TBX5 binding motif. (B) EMSA result for variants rs4657475, rs116461237 and rs61058666. Lane 1, Nppa promoter is known to bind to TBX5. Lane 2, wildtype SNP sequence forms complex with TBX5. Lanes 4–7, SNP rs116461237 variant disrupts TBX5 binding. Lanes 8–9, TBX5 does not form DNA-protein complex with Alu consensus sequence. Lanes 10–12, mutation of rs116461237 only disrupts TBX5-probe complex formation. Lanes 3, 7 and 9 contain cold competitor in addition to probe sequence. SNP, single nucleotide polymorphism.

Regulatory role of the TMCO1 SNP region

To test for a regulatory effect of SNP rs4656461 and its adjacent Alu region on gene transcription, a 3337 bp region surrounding SNP rs4656461 was cloned into a pGL4-derived vector containing a firefly luciferase reporter gene under the presence of a minimal promoter. Sequencing analysis revealed that Insert 2_703 had acquired the minor allele variant for SNP rs11430311, which was not in the expected haplotype. Sequence analysis of three inserts of all sizes (3_3337, 3_1274, 3_703 and 3_703Rev) also contained an undesired mutation 201 bp upstream of SNP rs11430311 (online supplementary figure 3).to

In HEK293T cells, inserts 1_3337, 2_3337, 3_3337 and 4_3337 repressed luciferase expression when compared to the empty pGL4.23 vector. Repression was significantly greater for Insert 2_3337 when compared to both Inserts 3_3337 (p=0.03) and 4_3337 (p=0.04) (figure 4A). Serial deletion of 3337 bp constructs demonstrated increasing levels of luciferase repression with 1274 bp and 703 bp constructs. ANOVA results comparing luciferase expression between the different sized inserts within each haplotype were all significant (Insert 1 p<0.0001, Insert 2 p=0.0006, Insert 3 p=0.003, Insert 4 p<0.0001) (figure 4B). In the reverse orientation, the 703 bp constructs showed enhancer activity, with Insert 4_703Rev showing significantly more reporter gene expression than wildtype Insert 1_703Rev (p=0.03) (figure 4C).

Figure 4

Luciferase assay in HEK293T Cells. (A) 3337 bp inserts repress luciferase expression. Insert with haplotype 2 (rs4657473 C>T) shows greater repression than Inserts 3 (rs116461237 G>A) (p=0.0261) and 4 (rs4656461 A>G) (p=0.0443). (B) Serial deletion of 3337 bp inserts identifies the core repressor 703 bp fragment surrounding SNP rs4656461. One-way ANOVA: Insert 1 (p<0.0001), Insert 2 (p=0.006), Insert 3 (p=0.0028), Insert 4 (p<0.001). (C) 703 bp insert acts as an enhancer when cloned in the reverse orientation. Greater enhancer activity is shown by Insert 4_703Rev (rs4656461 A>G) over wildtype sequence in Insert 1_703Rev (p=0.0252). SNP, single nucleotide polymorphism.

In primary TM cells, 3337 bp inserts did not have a regulatory effect on luciferase expression, while 1274 bp and 703 bp constructs demonstrated luciferase repression (online supplementary figure 4). Only one significant difference in repressor activity was found, which was between Insert 1_1274 and 3_1274. Serial deletion of the inserts did not demonstrate the same pattern of core repressor activity as in HEK293T cells, and changing the orientation of the 703 bp inserts did not result in enhancer activity .(online supplementary figure 5).

Discussion

The SNPs downstream from and within TMCO1 have been implicated in the genetic risk of POAG. The TMCO1 SNP rs4656461 was associated with POAG in Caucasians and then later in Asian populations with high IOP.7 15 Fewer studies have investigated the role of this variant in populations of African ancestry.16 The association of SNP rs4656461 with POAG was not replicated in our AA cohort. Our results are consistent with a prior case-control study, in which SNPs associated with POAG in populations of non-African ancestry (including SNP rs4656461) had no significance in those of African descent.10 However, we report a novel association between POAG and SNP rs4657473 (located within ~50 bp of rs4656461). Those with the homozygous TT genotype are at a lower risk of developing POAG. The allele frequencies of the SNP rs4657473 variant (C>T) are lower in Africans (7%) and African-Americans (15%) compared with European populations (25%) (1000Genomes). If the TT genotype is indeed protective, perhaps this lower prevalence contributes in some part to the greater prevalence of POAG in those of African ancestry.1 An additional association between a SNP in the TMCO1 region with POAG further strengthens the evidence that this region is involved in glaucoma pathogenesis. TMCO1 SNPs have previously been associated with high IOP and younger age at diagnosis.7 15 17 However, no associations were found between TMCO1 SNPs and quantitative glaucoma parameters such as IOP, CCT, CDR or RNFL thickness in our cohort and did not appear to define a glaucoma subtype or increase severity of disease.

It has largely been assumed that SNP rs4656461 is the functional POAG SNP and that it acts on its most proximal gene, TMCO1, that encodes a highly conserved transmembrane protein with a coiled-coil domain.18 19 However, the pathway through which rs4656461 is involved in POAG pathogenesis remains to be known. Recently, Wang et al described TMCO1 as an endoplasmic reticulum (ER) transmembrane protein that prevents ER calcium ion (Ca2+) stores from overfilling by forming a Ca2+-selective ion channel in response to Ca2+ overloading.18 TMCO1 protein expression has been found in several types of human ocular tissues, including the iris, ciliary body, retina and optic nerve.15 The mRNA expression levels of TMCO1 in human ocular tissue were shown to be highest in the TM and retina.17 TMCO1 protein may play a role in tumour suppression and/or cell cycle regulation, disrupting cellular functions and causing oxidative stress in ocular tissues.4 15

The intergenic region containing POAG-associated TMCO1 SNPs harbours an Alu element. Alu sequences belong to the Short Interspersed Element (SINE) family of human repetitive DNA elements.20 They are commonly located in introns, 3′ untranslated gene regions and intergenic genomic regions.20 Alu elements are transposable and after insertion into a new region of the genome acquire random mutations over time.20 The Alu element adjacent to SNP rs4656461 is from the AluJb subfamily, one of the oldest subfamilies, dating from 80 to 100 million years ago.21 Alus tend to be located near and within genes, which may be due to positive selection exerted as a result of functional roles, such as regulation of gene transcription.21 Alus have been shown to act as enhancers or silencers.22 They have an ability to introduce new TF binding sites, promote nucleosome assembly, induce methylation or form secondary structures.20 22 Alu consensus sequences were shown to be enriched for transcription factor binding motifs, including those involved in pathways of eye, brain and central nervous system development.21 The influence on POAG risk of TMCO1 SNPs may be due to the presence of an Alu sequence and the SNPs it contains, which could modulate differential transcriptional regulation and variable transcription factor binding.

To investigate a possible regulatory role of TMCO1 SNPs in this Alu region, we interrogated their transcription factor binding potential. A RegulomeDB query predicted that TMCO1 SNPs rs4657475, rs116461237 and rs61058666 are located in a TBX5 binding motif (figure 3A). TBX5 is a member of the T-box transcription factor gene family, whose genes regulate many developmental processes in vertebrates and invertebrates.23 TBX5 has been implicated in the dorsal-ventral patterning of the eye and the establishment of visual projections between the retina and tectum.24 Using EMSA, we confirmed that TBX5 forms complexes with an oligonucleotide sequence containing this predicted DNA binding motif, encompassing the three wildtype TMCO1 SNPs. The SNP variant rs116461237 disrupted complex formation, suggesting that variants at this locus may lead to differential transcription factor binding and thereby affect gene regulation. Although the SNP rs116461237 is not associated with POAG in our cohort, we infer that the previously associated SNP rs4656561 along with the LD-associated SNPs rs4657474 and rs4657475 (figure 2B) do not influence TBX5 binding (figure 3B).5 Using EMSA, we demonstrated that the AluJb sequence downstream from TMCO1 has affinity for TBX5, whereas the ancestral Alu sequence did not. This indicates that the TBX5-binding is a new trait, consistent with the proposal that this particular repetitive element has acquired a regulatory function.

Luciferase assays in HEK293T cells demonstrated that the TMCO1 SNP region is able to act as a transcriptional repressor in vitro. Serial deletion of this region revealed a core repressor element of 703 bp surrounding the AluJb sequence. As serial deletions removed the non-Alu parts of the sequence, it appears that the Alu region effected the gene repression. Alu sequences have been shown previously to silence gene transcription.25 26 They have also been shown to act as either negative or positive regulators of transcription based on their orientation.27 In our study, we confirmed this relationship between Alu orientation and differential effects on gene transcription. The 703 bp sequence when cloned into our reporter vector in the reverse orientation acted to enhance gene expression. The variying levels of repression and/or enhancer activity by the minimal regulatory region in our study need additional studies to establish if these AluJb haplotypes can modulate gene transcription.

When we tested our luciferase reporter constructs in an ocular cell line, primary TM cells, we did not observe the same core repressor pattern or orientation dependent regulatory activity. A different set of transcription factors may be expressed in TM than in HEK293T cells and perhaps TMCO1 SNPs are not relevant in the TM but rather a different ocular cell type. Further testing of other ocular tissues is required. In the luciferase reporter study, we compared the TMCO1 SNP region of interest to an empty vector rather than one with a similar sized random sequence. It has been demonstrated that luciferase reporter activity is inversely correlated with construct size.28 However, we do not believe that this confounded our results because we saw decreased luciferase expression with serial deletion (shortening) of our construct and saw both enhancement and repression of the 703 bp construct depending on its orientation without a change in its size.

If this Alu region is indeed regulating the transcription of TMCO1, it still remains to be shown that the TMCO1 protein is involved in POAG pathogenesis. In the TMCO1 defect syndrome, a distinct combination of craniofacial dysmorphism, skeletal anomalies and mental retardation are observed.29 There have been no eye-related abnormalities reported other than strabismus.29 However, the recently described function of TMCO1 as a regulator of cellular calcium homeostasis suggests a mechanism for its involvement in POAG.30 Dysregulation of cellular Ca2+ homeostasis, ER stress and Ca2+ overload have been implicated in the pathogenesis of neurodegenerative disorders such as Alzheimer’s disease, Parkinson’s disease, Huntington’s disease and amyotrophic lateral sclerosis.31 As suggested in these diseases, it is possible to speculate that perturbation of Ca2+ homeostasis through TMCO1 dysfunction could contribute to the degeneration of retinal ganglion cells in the optic nerve observed in POAG. Alterations in calcium flux may affect aqueous outflow pathways through TM cell volume regulation, possibly explaining why TMCO1 variants have been associated with high IOP.32 Alternatively, it is possible that these SNPs may influence the activity of the downstream gene, ALDH9A1, which has been implicated as a regulator of IOP by eQTL analaysis.33

In conclusion, this study supports the involvement of genetic variation at this TMCO1 SNP locus in the genetic risk of POAG in an AA population and suggests a possible regulatory role for this region. Future studies, however, are needed to determine if this gene region regulates the TMCO1 gene in vivo, and knockout TMCO1 in other ocular cell lines and in animal models are planned to determine the precise role of TMCO1 in the eye.