Molecular Vision 2015; 21:548-554 <http://www.molvis.org/molvis/v21/548>
Received 11 February 2015 | Accepted 08 May 2015 | Published 11 May 2015

Expression-associated polymorphisms of CAV1-CAV2 affect intraocular pressure and high-tension glaucoma risk

Sewon Kim,1 Kyunglan Kim,1 Dong Won Heo,2 Jong-Sung Kim,3 Chan Kee Park,4 Chang-sik Kim,2 Changwon Kang1

1Department of Biological Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Korea; 2Department of Ophthalmology, Chungnam National University Hospital, Daejeon, Korea; 3Department of Family Medicine, Chungnam National University Hospital, Daejeon, Korea; 4Department of Ophthalmology and Visual Science, Seoul St. Mary’s Hospital, College of Medicine, the Catholic University of Korea, Seoul, Korea

Correspondence to: Changwon Kang, Department of Biological Sciences, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 305-701, Korea; Phone: +82-42-350-2628; FAX: +82-42-350-5628; email: ckang@kaist.ac.kr

Abstract

Purpose: The human CAV1-CAV2 locus has been associated with susceptibility to primary open-angle glaucoma in four studies of Caucasian, Chinese, and Pakistani populations, although not in several other studies of non-Korean populations. In this study with Korean participants, the CAV1-CAV2 locus was investigated for associations with susceptibility to primary open-angle glaucoma accompanied by elevated intraocular pressure (IOP), namely, high-tension glaucoma (HTG), as well as with IOP elevation, which is a strong risk factor for glaucoma.

Methods: Two single nucleotide polymorphisms (SNPs) were genotyped in 1,161 Korean participants including 229 patients with HTG and 932 healthy controls and statistically examined for association with HTG susceptibility and IOP. One SNP was rs4236601 G>A, which had been reported in the original study, and the other SNP was rs17588172 T>G, which was perfectly correlated (r2=1) with another reported SNP rs1052990. Expression quantitative trait loci (eQTL) analysis was performed using GENe Expression VARiation (Genevar) data.

Results: Both SNPs were associated with HTG susceptibility, but the rs4236601 association disappeared when adjusted for the rs17588172 genotype and not vice versa. The minor allele G of rs17588172 was associated significantly with 1.5-fold increased susceptibility to HTG (p=0.0069) and marginally with IOP elevation (p=0.043) versus the major allele T. This minor allele was also associated with decreased CAV1 and CAV2 mRNA in skin and adipose according to the Genevar eQTL analysis.

Conclusions: The minor allele G of rs17588172 in the CAV1-CAV2 locus is associated with decreased expression of CAV1 and CAV2 in some tissues, marginally with IOP elevation, and consequently with increased susceptibility to HTG.

Introduction

Glaucoma, a major eye disease that causes visual blindness, is characterized by optic nerve degeneration and visual defects. Primary open-angle glaucoma (POAG), the most common type of glaucoma, is accompanied by elevated intraocular pressure (IOP >21 mmHg) in high-tension glaucoma (HTG) or by normal IOP (≤21 mmHg) in normal-tension glaucoma (NTG). Although elevated IOP is a major risk factor for glaucoma, increased IOP is neither necessary nor sufficient to develop POAG [1-3]. Several years ago, an intergenic single nucleotide polymorphism (SNP), rs4236601, located between CAV2 (gene ID: 858; OMIM 601048; upstream) and CAV1 (gene ID: 857; OMIM 601047; downstream) on human chromosome 7q31.2 was identified as associated with susceptibility to POAG in a genome-wide association study (GWAS); the association was discovered in an Icelandic Caucasian population and replicated in additional European Caucasian, Australian Caucasian, and Chinese Asian populations [4].

Among eight subsequent replication studies using Iowa U.S. Caucasian [5], Afro-Caribbean [6], Saudi Arabian [7], U.S. Caucasian [8], Japanese Asian [9], African [10], U.S. Caucasian [11], and Pakistani [12] populations, only three studies using U.S. Caucasian and Pakistani populations successfully replicated the POAG association [8,11,12], as well as an additional meta-analysis [13]. In addition, the original study [4] and the first positive replication study [8] reported preferential association with susceptibility to NTG rather than to HTG.

Interestingly, the POAG risk-associated SNP alleles in the CAV1-CAV2 locus were also associated with IOP elevation [14-17]. Furthermore, the HTG risk association was observed in the Chinese [4] and Pakistani [12] populations. Because it seemed odd that the alleles associated with IOP elevation were preferentially associated with increased susceptibility to NTG with normal IOP rather than HTG with high IOP [4,8], these variants were examined to replicate the associations with HTG susceptibility, as well as with IOP, and CAV1 and CAV2 expression.

Methods

Study samples

Adhering to the ARVO statement on human subjects, this study followed the tenets of the Declaration of Helsinki. A total of 1,161 unrelated Korean participants including 229 patients with HTG (men 68.6%; mean age 56.1±15.9 years) and 932 healthy controls (men 50.1%; mean age: 57.6±13.2 years) were recruited at Chungnam National University Hospital and Seoul St. Mary’s Hospital under approval from the Institutional Review Board of each hospital, and each participant provided written informed consent.

In the case group, only patients diagnosed with POAG with high IOP or HTG by two glaucoma specialists (C.-S.K. and C.K.P.) at the glaucoma clinic of the two tertiary referral hospitals as previously described were included [18,19]. Among POAG patients with (1) characteristic glaucomatous optic nerve head defect, (2) visual field defect corresponding to the optic nerve head defect, (3) retinal nerve fiber layer defect and (4) normal open angle or deep anterior chamber angle, only those with IOP of >21 mmHg were included in this HTG case group.

All 932 controls were free of glaucoma and other eye diseases. Participants who were free of glaucoma but had a family history of glaucoma and a history of steroid treatment were excluded in the control group. Intraocular pressure was measured with a Goldman applanation tonometer for each visit, and the average of two measurements was recorded.

SNP genotyping and association tests

Under the approval of the KAIST Institutional Review Board, genomic DNA was extracted from peripheral venous blood collected in Vacutainer tubes using the Gentra Puregene blood kit (Qiagen, Valencia, CA) and stored at -80°C deep freezer until genotyped using the MassARRAY hME assay (Sequenom, San Diego, CA). The call rate was 99.9% for rs4236601 and 98.7% for rs17588172, and the genotype data were under Hardy–Weinberg equilibrium (p=1.00 for rs4236601 and p=0.54 for rs17588172).

SNP association with HTG susceptibility was assessed in a comparison of 229 patients with HTG and 932 controls using logistic regression analysis, with adjustment for age, gender, and recruitment site, unless described otherwise. Association was considered statistically significant when p≤α and marginal when α<p≤0.05, with a significance level of Bonferroni correction for multiple testing of α=0.05/2=0.025, as two SNPs were tested.

SNP association of IOP levels was additionally tested among 1,095 participants using nonparametric tests suitable for the non-normal distribution of IOP shown in Appendix 1. There were no pretreatment IOP measurement records for 66 patients with HTG, although their IOP levels were still high after treatment; therefore, these samples were excluded from the IOP association tests. The Jonckheere trend and Mann–Whitney U tests were used for rs17588172 of three genotypes and rs4236601 of two genotypes (with no minor-allele homozygotes), respectively.

All statistical tests were performed using PLINK v.1.07 or IBM SPSS Statistics v.18 (IBM Corporation, Chicago, IL). Statistical power was calculated using the CaTS interface.

Expression quantitative trait loci analysis

Two data sets from the GENe Expression VARiation (Genevar) database [20] were used for the cis-expression quantitative trait loci (eQTL) analysis. The CAV1 and CAV2 mRNA levels in adipose, skin, and the lymphoblastoid cell line (LCL) were retrieved from a data set of tissue-dependent mRNA level measurements from 856 British Caucasians by Grundberg et al. [21] and used to calculate standardized regression coefficient β values. To obtain Asian data, the CAV1 and CAV2 mRNA levels of 82 Japanese and 80 Chinese participants in the International HapMap Project were retrieved from another data set of population-dependent mRNA level measurements in LCLs from Stranger et al. [22] and used to calculate Spearman’s rank-order correlation coefficient ρ values.

Results

HTG susceptibility association of CAV1-CAV2 SNPs

A total of 1,161 case-control subjects of Korean ancestry, whose demographic characteristics are shown in Table 1, were recruited to this HTG association study. The rs4236601 (G>A) SNP had been reported in the original association study [4] and was the first SNP genotyped for replication in this study. Another SNP rs1052990 had been reported in the second positive-association study [8], but the SNP’s perfectly correlated (r2=1) rs17588172 (T>G) was the second SNP genotyped in this study. These two SNPs were tested for HTG susceptibility by comparing genotype distributions in 229 patients with HTG and 932 healthy controls (Table 2).

HTG susceptibility was significantly associated with the second SNP rs17588172 (odds ratio [OR]=1.5, p=0.0069) after Bonferroni correction for multiple testing of two SNPs (α=0.025). This association was maintained after adjusted for the genotype of the other, first SNP rs4236601 (pcon=0.010). However, the marginal HTG susceptibility association of the first SNP rs4236601 (p=0.026) disappeared when adjusted for the genotype of the second SNP rs17588172 (pcon=0.15). Thus, HTG susceptibility association of the CAV1-CAV2 locus was found with rs17588172 in this Korean population.

IOP association of CAV1-CAV2 SNPs

Next, the two SNPs were examined for association with IOP. All 932 healthy controls had IOP measurements, and only 163 patients with HTG had pretreatment IOP measurements. The other 66 patients with HTG had only post-treatment IOP measurements and therefore were excluded from this analysis to exclude the confounding treatment effect on IOP.

In a total of 1,095 samples, IOP was not normally distributed (Appendix 1), and nonparametric tests were used for this analysis. IOP was associated significantly with rs4236601 (p=0.010) but marginally with rs17588172 (p=0.043), while the HTG risk-associated alleles in both SNPs were associated with IOP elevation.

Because IOP elevation is a strong risk factor for glaucoma, we next examined whether the increased HTG susceptibility was derived from the increase in IOP (Table 2). After the adjustment for the IOP, the HTG association of rs17588172 disappeared (padj=0.40). Accordingly, the HTG risk-associated allele appeared to affect the IOP, which could mostly account for its effect on susceptibility to HTG.

Gene expression association of CAV1 and CAV2 SNPs

The HTG risk-associated SNP rs17588172 was highly correlated with 15 SNPs (r2≥0.8) in the Korean HapMap population, suggesting their possible association with HTG susceptibility. Because these SNPs from rs1018859 in chromosome position 116,067,347 to rs7801950 in position 116,154,783 spanned an 87 kb region encompassing CAV1 and CAV2, the possible effects on the two gene expression levels were examined using cis-eQTL analysis with two data sets from the Genevar database [20].

The CAV1 and CAV2 mRNA levels were associated with the two SNPs, but more strongly with the second SNP rs17588172 than the first SNP rs4236601 (Table 3). These eQTL results were unlikely to be derived from perturbed probe binding because no SNPs in the probe binding sites were in high linkage disequilibrium (LD) with either SNP. The risk allele G of rs1758872 was significantly correlated with decreased CAV1 mRNA in adipose (β=−0.32, p=3.9×10−23) and skin (β=−0.22, p=3.0×10−18) from the Caucasian subjects, although the effect was minimal in the LCLs derived from the Caucasian subjects (β=−0.027, p=0.00013), and no association was found in the LCLs from the Chinese (p=0.33) and Japanese (p=0.92) subjects.

The same risk allele was also correlated with the CAV2 mRNA level, but the direction of the effect was opposite depending on the tissue or cell line (Table 3). The risk allele was associated with increased CAV2 mRNA in the LCLs derived from the Japanese (ρ=0.82, p=9.5×10−21), Chinese (ρ=0.68, p=2.8×10−12), and Caucasian (β=0.045, p=4.9×10−8) subjects but decreased mRNA in adipose (β=−0.11, p=0.00027) and skin (β=−0.11, p=0.00014) from the Caucasian subjects. Accordingly, the risk-associated minor alleles, G for rs17518872 and A for rs4236601, decrease CAV1 and CAV2 mRNA in adipose and skin, although in the LCLs the alleles increase CAV2 mRNA but tend to decrease CAV1 mRNA.

Discussion

This report describes an SNP allele in the human CAV1-CAV2 locus associated with not only increased susceptibility to HTG but also increased IOP. The HTG risk-associated, minor allele G of rs17588172 was also associated with decreased CAV1 and CAV2 mRNA in skin and adipose. Thus, the rs17588172 minor allele appears to increase IOP and consequently increase HTG susceptibility possibly by decreasing the expression of CAV1, CAV2, or both.

The CAV1 and CAV2 genes are expressed in glaucoma- and IOP-related tissues [5], and the respective gene products, caveolin 1 and 2, involved in the formation of caveolae, are known to negatively regulate transforming growth factor beta (TGF-β) signaling [23,24], which has been implicated to cause IOP elevation and optic nerve degeneration [25]. If so, decreased expression of CAV1 or CAV2 would be associated with increased IOP and HTG risk. This hypothesis is consistent with our findings that the HTG risk-associated allele of rs17588172 confers decreased expression of CAV1 and CAV2 in adipose and skin, although not in LCLs, according to our eQTL analyses using the Genevar data for the mRNA levels.

Recently, Cav1-knockout mice were found to have higher IOPs than age- and sex-matched control mice [26]. In perfusion culture of human and porcine anterior segments, however, CAV2 knockdown rather than CAV1 knockdown decreased the outflow rate possibly leading to IOP elevation [27]. Thus, downregulation of caveolin 1 or 2 in the eye’s aqueous outflow pathway could enhance TGF-β signaling to increase IOP and consequently increase HTG susceptibility.

The HTG association was more evident with rs17588172 than rs4236601, as association of rs4236601 with HTG susceptibility disappeared after adjusting for rs17588172, possibly because the minor allele of rs4236601 was too rare in this Korean population to provide sufficient statistical power. This hierarchy of SNP association significance was not consistent with the original GWAS results in which rs1052990 (p=1.1×10−9), which was perfectly correlated with rs17588172, showed a less significant POAG association than rs4236601 (p=5.0×10−10) in a Caucasian population [4].

However, this hierarchy was consistent with a later finding that rs1052990 (p=3.0×10−4) showed a more significant POAG association than rs4236601 (p=1.4×10−3) in another U.S. Caucasian population [8]. Interestingly in a recent study, rs4236601 was the top SNP in association with susceptibility to all POAG, but rs17588172 was the top in association with susceptibility to POAG with early paracentral visual field loss [11].

The statistical power of this study to detect the originally reported association of rs1052990 (OR=1.32) or the perfectly correlated, second SNP rs17588172 was 61% with a minor allele frequency of 21% in our samples. In contrast, the power to detect the original Caucasian association of the first SNP rs4236601 (OR=1.36) was too low (6%), possibly because the minor allele was too rare (0.8%) in Koreans compared with 28% in Caucasians. Therefore, the marginal association of rs4236601 could have been caused by insufficient statistical power.

In summary, in a Korean population, HTG susceptibility and IOP were associated with SNP rs17588172 (and presumably its perfectly correlated SNP rs1052990 as well) in the human CAV1-CAV2 locus, which appeared to affect CAV1 and CAV2 expression. These findings suggest that changes in the expression of CAV1, CAV2, or both could affect IOP and consequently susceptibility to HTG or POAG with elevated IOP.

Appendix 1. Distribution of intraocular pressure (IOP).

Acknowledgments

The authors thank Yong-jun Yun for assistance with the blood sample collection and Yuri Chae for technical and administrative assistance. This work was supported in part by a grant (2014006849) from the National Research Foundation of Korea. Three authors, Changwon Kang (ckang@kaist.ac.kr), Chang-sik Kim (kcs61@cnu.ac.kr) and Chan Kee Park (ckpark@catholic.ac.kr) are co-corresponding authors. The authors declare no conflict of interest.

References

  1. Anderson DR. Collaborative normal tension glaucoma study. Curr Opin Ophthalmol. 2003; 14:86-90. [PMID: 12698048]
  2. Armaly MF, Krueger DE, Maunder L, Becker B, Hetherington J, , Jr Kolker AE, Levene RZ, Maumenee AE, Pollack IP, Shaffer RN. Biostatistical analysis of the collaborative glaucoma study. I. Summary report of the risk factors for glaucomatous visual-field defects. Arch Ophthalmol. 1980; 98:2163-71. [PMID: 7447768]
  3. Heijl A, Leske MC, Bengtsson B, Hyman L, Hussein M. Reduction of intraocular pressure and glaucoma progression: results from the Early Manifest Glaucoma Trial. Arch Ophthalmol. 2002; 120:1268-79. [PMID: 12365904]
  4. Thorleifsson G, Walters GB, Hewitt AW, Masson G, Helgason A, DeWan A, Sigurdsson A, Jonasdottir A, Gudjonsson SA, Magnusson KP, Stefansson H, Lam DS, Tam PO, Gudmundsdottir GJ, Southgate L, Burdon KP, Gottfredsdottir MS, Aldred MA, Mitchell P, St Clair D, Collier DA, Tang N, Sveinsson O, Macgregor S, Martin NG, Cree AJ, Gibson J, Macleod A, Jacob A, Ennis S, Young TL, Chan JC, Karwatowski WS, Hammond CJ, Thordarson K, Zhang M, Wadelius C, Lotery AJ, Trembath RC, Pang CP, Hoh J, Craig JE, Kong A, Mackey DA, Jonasson F, Thorsteinsdottir U, Stefansson K. Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma. Nat Genet. 2010; 42:906-9. [PMID: 20835238]
  5. Kuehn MH, Wang K, Roos B, Stone EM, Kwon YH, Alward WL, Mullins RF, Fingert JH. Chromosome 7q31 POAG locus: ocular expression of caveolins and lack of association with POAG in a US cohort. Mol Vis. 2011; 17:430-5. [PMID: 21321670]
  6. 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:e39278 [PMID: 22761751]
  7. 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. [PMID: 22876122]
  8. Wiggs JL, Kang JH, Yaspan BL, Mirel DB, Laurie C, Crenshaw A, Brodeur W, Gogarten S, Olson LM, Abdrabou W, DelBono E, Loomis S, Haines JL, Pasquale LR. Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma in Caucasians from the USA. Hum Mol Genet. 2011; 20:4707-13. [PMID: 21873608]
  9. Kato T, Meguro A, Nomura E, Uemoto R, Nomura N, Ota M, Kashiwagi K, Mabuchi F, Iijima H, Kawase K, Yamamoto T, Nakamura M, Negi A, Sagara T, Nishida T, Inatani M, Tanihara H, Aihara M, Araie M, Fukuchi T, Abe H, Higashide T, Sugiyama K, Kanamoto T, Kiuchi Y, Iwase A, Chin S, Ohno S, Inoko H, Mizuki N. Association study of genetic variants on chromosome 7q31 with susceptibility to normal tension glaucoma in a Japanese population. Eye (Lond). 2013; 27:979-83. [PMID: 23743525]
  10. Liu Y, Hauser MA, Akafo SK, Qin X, Miura S, Gibson JR, Wheeler J, Gaasterland DE, Challa P, Herndon LW, Ritch R, Moroi SE, Pasquale LR, Girkin CA, Budenz DL, Wiggs JL, Richards JE, Ashley-Koch AE, Allingham RR. 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. [PMID: 23963167]
  11. Loomis SJ, Kang JH, Weinreb RN, Yaspan BL, Cooke Bailey JN, Gaasterland D, Gaasterland T, Lee RK, Lichter PR, Budenz DL, Liu Y, Realini T, Friedman DS, McCarty CA, Moroi SE, Olson L, Schuman JS, Singh K, Vollrath D, Wollstein G, Zack DJ, Brilliant M, Sit AJ, Christen WG, Fingert J, Kraft P, Zhang K, Allingham RR, Pericak-Vance MA, Richards JE, Hauser MA, Haines JL, Pasquale LR, Wiggs JL. Association of CAV1/CAV2 genomic variants with primary open-angle glaucoma overall and by gender and pattern of visual field loss. Ophthalmology. 2014; 121:508-16. [PMID: 24572674]
  12. Micheal S, Ayub H, Khan MI, Bakker B, Schoenmaker-Koller FE, Ali M, Akhtar F, Khan WA, Qamar R, den Hollander AI. Association of known common genetic variants with primary open angle, primary angle closure, and pseudoexfoliation glaucoma in Pakistani cohorts. Mol Vis. 2014; 20:1471-9. [PMID: 25489222]
  13. Huang W, Wang W, Zhou M, Zhang X. Association of single-nucleotide polymorphism rs4236601 near caveolin 1 and 2 with primary open-angle glaucoma: a meta-analysis. Clin Experiment Ophthalmol. 2014; 42:515-21. [PMID: 24034151]
  14. van Koolwijk LM, Ramdas WD, Ikram MK, Jansonius NM, Pasutto F, Hysi PG, Macgregor S, Janssen SF, Hewitt AW, Viswanathan AC, ten Brink JB, Hosseini SM, Amin N, Despriet DD, Willemse-Assink JJ, Kramer R, Rivadeneira F, Struchalin M, Aulchenko YS, Weisschuh N, Zenkel M, Mardin CY, Gramer E, Welge-Lussen U, Montgomery GW, Carbonaro F, Young TL, Bellenguez C, McGuffin P, Foster PJ, Topouzis F, Mitchell P, Wang JJ, Wong TY, Czudowska MA, Hofman A, Uitterlinden AG, Wolfs RC, de Jong PT, Oostra BA, Paterson AD, Mackey DA, Bergen AA, Reis A, Hammond CJ, Vingerling JR, Lemij HG, Klaver CC, van Duijn CM. Common genetic determinants of intraocular pressure and primary open-angle glaucoma. PLoS Genet. 2012; 8:e1002611 [PMID: 22570627]
  15. Chen F, Klein AP, Klein BE, Lee KE, Truitt B, Klein R, Iyengar SK, Duggal P. Exome Array Analysis Identifies CAV1/CAV2 as a Susceptibility Locus for Intraocular Pressure. Invest Ophthalmol Vis Sci. 2015; 56:544-51. [PMID: 25525164]
  16. Hysi PG, Cheng CY, Springelkamp H, Macgregor S, Bailey JN, Wojciechowski R, Vitart V, Nag A, Hewitt AW, Hohn R, Venturini C, Mirshahi A, Ramdas WD, Thorleifsson G, Vithana E, Khor CC, Stefansson AB, Liao J, Haines JL, Amin N, Wang YX, Wild PS, Ozel AB, Li JZ, Fleck BW, Zeller T, Staffieri SE, Teo YY, Cuellar-Partida G, Luo X, Allingham RR, Richards JE, Senft A, Karssen LC, Zheng Y, Bellenguez C, Xu L, Iglesias AI, Wilson JF, Kang JH, van Leeuwen EM, Jonsson V, Thorsteinsdottir U, Despriet DD, Ennis S, Moroi SE, Martin NG, Jansonius NM, Yazar S, Tai ES, Amouyel P, Kirwan J, van Koolwijk LM, Hauser MA, Jonasson F, Leo P, Loomis SJ, Fogarty R, Rivadeneira F, Kearns L, Lackner KJ, de Jong PT, Simpson CL, Pennell CE, Oostra BA, Uitterlinden AG, Saw SM, Lotery AJ, Bailey-Wilson JE, Hofman A, Vingerling JR, Maubaret C, Pfeiffer N, Wolfs RC, Lemij HG, Young TL, Pasquale LR, Delcourt C, Spector TD, Klaver CC, Small KS, Burdon KP, Stefansson K, Wong TY, Viswanathan A, Mackey DA, Craig JE, Wiggs JL, van Duijn CM, Hammond CJ, Aung T. Genome-wide analysis of multi-ancestry cohorts identifies new loci influencing intraocular pressure and susceptibility to glaucoma. Nat Genet. 2014; 46:1126-30. [PMID: 25173106]
  17. Ozel AB, Moroi SE, Reed DM, Nika M, Schmidt CM, Akbari S, Scott K, Rozsa F, Pawar H, Musch DC, Lichter PR, Gaasterland D, Branham K, Gilbert J, Garnai SJ, Chen W, Othman M, Heckenlively J, Swaroop A, Abecasis G, Friedman DS, Zack D, Ashley-Koch A, Ulmer M, Kang JH, Liu Y, Yaspan BL, Haines J, Allingham RR, Hauser MA, Pasquale L, Wiggs J, Richards JE, Li JZ. Genome-wide association study and meta-analysis of intraocular pressure. Hum Genet. 2014; 133:41-57. [PMID: 24002674]
  18. Kim K, Yun YJ, Kim S, Kim JS, Kim CS, Kang C. Analysis of an extended chromosome locus 2p14–21 for replication of the 2p16.3 association with glaucoma susceptibility. Mol Vis. 2011; 17:1136-43. [PMID: 21552472]
  19. 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:409-13. [PMID: 23838595]
  20. Yang TP, Beazley C, Montgomery SB, Dimas AS, Gutierrez-Arcelus M, Stranger BE, Deloukas P, Dermitzakis ET. Genevar: a database and Java application for the analysis and visualization of SNP-gene associations in eQTL studies. Bioinformatics. 2010; 26:2474-6. [PMID: 20702402]
  21. Grundberg E, Small KS, Hedman AK, Nica AC, Buil A, Keildson S, Bell JT, Yang TP, Meduri E, Barrett A, Nisbett J, Sekowska M, Wilk A, Shin SY, Glass D, Travers M, Min JL, Ring S, Ho K, Thorleifsson G, Kong A, Thorsteindottir U, Ainali C, Dimas AS, Hassanali N, Ingle C, Knowles D, Krestyaninova M, Lowe CE, Di Meglio P, Montgomery SB, Parts L, Potter S, Surdulescu G, Tsaprouni L, Tsoka S, Bataille V, Durbin R, Nestle FO, O'Rahilly S, Soranzo N, Lindgren CM, Zondervan KT, Ahmadi KR, Schadt EE, Stefansson K, Smith GD, McCarthy MI, Deloukas P, Dermitzakis ET, Spector TD. Mapping cis- and trans-regulatory effects across multiple tissues in twins. Nat Genet. 2012; 44:1084-9. [PMID: 22941192]
  22. Stranger BE, Montgomery SB, Dimas AS, Parts L, Stegle O, Ingle CE, Sekowska M, Smith GD, Evans D, Gutierrez-Arcelus M, Price A, Raj T, Nisbett J, Nica AC, Beazley C, Durbin R, Deloukas P, Dermitzakis ET. Patterns of cis regulatory variation in diverse human populations. PLoS Genet. 2012; 8:e1002639 [PMID: 22532805]
  23. Razani B, Zhang XL, Bitzer M, von Gersdorff G, Bottinger EP, Lisanti MP. Caveolin-1 regulates transforming growth factor (TGF)-beta/SMAD signaling through an interaction with the TGF-beta type I receptor. J Biol Chem. 2001; 276:6727-38. [PMID: 11102446]
  24. Xie L, Vo-Ransdell C, Abel B, Willoughby C, Jang S, Sowa G. Caveolin-2 is a negative regulator of anti-proliferative function and signaling of transforming growth factor-beta in endothelial cells. Am J Physiol Cell Physiol. 2011; 301:C1161-74. [PMID: 21832243]
  25. Fuchshofer R, Tamm ER. The role of TGF-beta in the pathogenesis of primary open-angle glaucoma. Cell Tissue Res. 2012; 347:279-90. [PMID: 22101332]
  26. Elliott MH, Gu X, Ashpole NE, Griffith GL, Boyce TM, Tanito M, Tamm ER, Stamer WD. Role of caveolin-1 in intraocular pressure and conventional outflow regulation. ARVO 2014, May 4-8.
  27. Aga M, Bradley JM, Wanchu R, Yang YF, Acott TS, Keller KE. Differential effects of caveolin-1 and −2 knockdown on aqueous outflow and altered extracellular matrix turnover in caveolin-silenced trabecular meshwork cells. Invest Ophthalmol Vis Sci. 2014; 55:5497-509. [PMID: 25103269]