HTRA1 promoter variant differentiates polypoidal choroidal vasculopathy from exudative age-related macular degeneration

Abstract

Exudative age-related macular degeneration (AMD) and polypoidal choroidal vasculopathy (PCV) share similar abnormal choroidal vasculature, but responses to treatments are different. In this study, we sequenced the whole HTRA1 gene and its promoter by direct sequencing in a Hong Kong Chinese PCV cohort. We identified rs11200638, c.34delCinsTCCT, c.59C>T, rs1049331 and rs2293870 significantly associated with PCV. Notably, rs2672598 was significantly associated with exudative AMD (p = 1.31 × 10^-4 ) than PCV (p = 0.11). Logistic regression indicated that rs2672598 (p = 2.27 × 10^-3 ) remained significant after adjusting for rs11200638 in exudative AMD. Moreover, the rs11200638-rs2672598 joint genotype AA-CC conferred higher risk to exudative AMD (43.11 folds) than PCV (3.68 folds). Promoter analysis showed that rs2672598 C-allele showed higher luciferase expression than wildtype T-allele (p = 0.026), independent of rs11200638 genotype (p = 0.621). Coherently, vitreous humor HTRA1 expression with rs2672598 CC genotype was significantly higher than that with TT genotype by 2.56 folds (p = 0.02). Furthermore, rs2672598 C-allele was predicted to alter the transcription factor binding sites, but not rs11200638 A-allele. Our results revealed that HTRA1 rs2672598 is more significantly associated with exudative AMD than PCV in ARMS2/HTRA1 region, and it is responsible for elevated HTRA1 transcriptional activity and HTRA1 protein expression.

promoter variant

R P

Ng,*, XiaoYing Liang,* Lai, Li Ma Tam, Jian XiongWang, Li Jia Pang

>

p= p=

p= p= p=

p=

region, and

Age-related macular degeneration (AMD) is a progressive macular disease and a major cause of irreversible blindness in the elderly population1. At the advanced stage, it can be classied into exudative and non-exudative AMD2. Exudative AMD, characterized by choroidal neovascularization (CNV), shares similar phenotype with polypoidal choroidal vasculopathy (PCV). Both of them involve abnormal vasculopathies arising from the choroidal vasculature leading to sub-retinal serous exudation and hemorrhage3,4. PCV is characterized by the polypoidal structures beneath the relatively intact Bruchs membrane, branching from abnormal inner choroidal vasculature5. The polypoidal structures can be visualized by indocyanine green angiography (ICGA), a diagnostic tool for dierentiating PCV lesion from CNV. PCV and CNV are dierent in their risk factors6, clinical manifestations7, natural courses5, responses to treatment8, and the overall visual prognosis9. PCV occurs predominately in Asian populations and tends to be younger6. Its disease course is relatively benign with less drusen but a tendency of relapse5. Without treatment, the visual outcome of PCV is usually guarded10. Patients with AMD, however, tend to present with a more aggressive course and uncontrollable visual deterioration6. In addition, exudative AMD responds better to anti-vascular endothelial growth factor (VEGF) treatment, whereas PCV better to photodynamic therapy (PDT)9,11,12.

Some reports propose that PCV could be a variant of typical CNV and regarded as one of the subtypes of exudative AMD6,13. This could be reected by their similar genetic and environmental factors. Multiple genetic variants have exhibited signicant association with both exudative AMD and PCV6. These include the complement factor H (CFH) gene, age-related maculopathy susceptibility 2 (ARMS2)/high temperature requirement factor A1 (HTRA1) locus and complement component 2 (C2)/factor B (CFB) genes14. Previously, we screened

Department of Ophthalmology and Visual Sciences, The Chinese University of Hong Kong, Hong Kong. The Joint Shantou International Eye Center of Shantou University and The Chinese University of Hong Kong, Shantou, China.

* These authors contributed equally to this work. Correspondence and requests for materials should be addressed to C.P.P. (email: [email protected])

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

Genotypic Frequencies

PCV (n=188) AMD (n=195)

Control

(n=183) pG pA

625G>A rs11200638 70 84 33 109 63 23 38 90 55 1.00104 7.49104 502C>T Reported 0 15 173 0 14 181 0 20 163 0.33 0.34 497C>T Reported 0 7 181 0 6 189 0 8 175 0.75 0.75 487T>C rs2672598 109 67 10 162 31 2 97 68 18 0.23 0.13c.34delCinsTCCT Novel L12insS 0 5 183 0 2 193 0 16 167 2.90102 3.10102c.59C>T Novel A20V 7 30 148 3 42 150 2 56 125 3.70103 0.01c.102C>T rs1049331 A34A 70 84 33 108 64 23 38 90 55 1.00104 9.43104c.108G>T rs2293870 G36G 70 84 33 108 64 23 38 90 55 1.00104 9.43104c.IVS1-176C>G rs12267142 0 17 170 0 19 176 0 25 158 0.17 0.18c.569G>A Novel R190H 0 1 186 0 0 195 0 0 183 1.00 1.00c.IVS2+100C>T rs80158665 0 4 183 0 4 191 0 8 175 0.23 0.23c.IVS2+216A>G Reported 0 1 186 0 0 195 0 1 182 1.00 1.00c.IVS2+317C>T Reported 0 2 185 0 2 193 0 1 182 1.00 1.00c.672C>T Novel N224N 0 1 186 0 0 195 0 0 183 1.00 1.00 IVS3+93C>T rs2239586 30 73 84 19 89 87 22 92 69 0.09 0.65 IVS3+167G>A rs2239587 30 73 84 19 89 87 22 92 69 0.09 0.65 IVS4+99C>T rs2672582 38 76 73 35 105 55 30 86 67 0.41 0.85 IVS5+76_79delGTTT Reported 7 56 124 0 69 126 0 48 135 IVS5+169G>A rs2672583 39 77 71 38 105 52 28 87 68 0.30 0.51 IVS6+58C>T Novel 0 1 186 0 0 195 0 0 183 1.00 .001 IVS6+111G>A rs79778361 2 42 144 5 39 151 1 57 125 IVS6+115C>G rs2672585 39 77 72 33 107 55 28 88 67 0.27 0.60 IVS7+149C>G rs76357476 2 43 143 5 40 150 1 57 125 c.1221C>T rs11538140 D407D 0 3 184 0 3 192 0 5 178 0.46 0.50 IVS8+14G>A rs2272599 71 76 40 53 107 35 67 85 31 0.42 0.67 IVS836C>T rs2293871 34 81 73 28 90 77 24 102 57

Table 1. The HTRA1 polymorphisms identied in PCV subjects. pG= genotypic p-value for PCV;pA= allelic p-value for PCV. The genotypic distribution was presented as homozygous/heterozygous/wild type. The rare variants in PCV (R190H) and exudative AMD (R59P, V417I) did not showed signicant association either in single marker or grouped analysis.

the entire ARMS2 gene by direct sequencing in our Hong Kong Chinese exudative AMD and PCV cohort15.

Although signicant dierences in the genotypic distributions of 11 polymorphisms in ARMS2 gene were found between exudative AMD and PCV, all of the associated polymorphisms showed similar trends. Since the ARMS2 and HTRA1 genes are in complete linkage disequilibrium (LD), it is also important to delineate the association pattern of HTRA1 in PCV patients. Three HTRA1 haplotypes are associated individually with typical AMD and PCV in a Japanese study although no dierence between them was reported16. Therefore, the dierential association between exudative AMD and PCV remains equivocal. The association of other HTRA1 variants in PCV should be investigated. In this study, 9 exons, intron-exon junctions and the promoter region of HTRA1 gene were sequenced in our Hong Kong Chinese PCV cohort with a view to clarify the association of HTRA1 with PCV.

Totally 26 polymorphisms in HTRA1 were identied in PCV patients (Table1). Four of them, IVS5+76_79delGTTT, rs79778361 (IVS6+111G>A), rs76357476 (IVS7+ 149C>G) and rs2293871 (IVS8-36C>T), were excluded from further analysis since they did not follow the Hardy-Weinberg equilibrium. Among the remaining 22 polymorphisms, 10 uncommon variants (minor allele frequency <3%) and 6 common variants did not showed signicant associations. They were excluded in subsequent analysis. The remaining 6 variants were further analyzed for their corresponding ORs under dominant model (Table2). Two SNPs, c.34delCinsTCCT and rs2672598 (487T>C), were also included in the analysis since they showed signicant association with exudative AMD. Four variants were associated with exudative AMD: rs11200638 (625G>A; p=1.16105, OR= 3.21, 95% condent intervals (CI): 1.885.50), c.34delCinsTCCT (L12insS; p = 0.001, OR = 0.13, 95% CI: 0.130.57), rs1049331 (c.102C>T, A34A; p = 1.16 105,

OR= 3.21, 95% CI: 1.885.50), and rs2293870 (c.108G>T, G36G; p=1.16105, OR= 3.21, 95% CI: 1.885.50) as well as PCV (rs11200638: p= 0.005, OR = 2.01, 95% CI: 1.233.28; c.34delCinsTCCT: p= 0.029, OR = 0.34, 95% CI: 0.130.90; rs1049331: p=0.009, OR= 1.91, 95% CI: 1.173.11; and rs2293870: p=0.009, OR=1.91, 95% CI: 1.173.11). None of these variants showed signicant dierences between AMD and PCV under dominant model (p > 0.05). Instead, a signicant dierential association (p= 0.015, OR = 5.48, 95% CI: 1.1925.37) was detected in rs2672598. It was strongly associated with exudative AMD (p = 1.31 104, OR = 10.53, 95% CI: 2.4146.04) but not associated with PCV (p= 0.106). Our result is coherent to a recent report in northern China

Polymorphism db SNP ID

Codon

change

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

Control

(n=183) p OR (95% CI) p OR (95% CI) p OR (95% CI)

625G>A rs11200638 109/63/23 70/84/33 38/90/55 1.16105 3.21 (1.885.50) 0.005 2.01 (1.233.28) 0.106 1.60 (0.902.85) 487T>C rs2672598 162/31/2 109/67/10 97/68/18 1.31104 10.53 (2.4146.04) 0.106 1.92 (0.864.28) 0.015 5.48 (1.1925.37)c.34delCinsTCCT 0/2/193 0/6/182 0/16/167 0.001 0.13 (0.030.57) 0.029 0.34 (0.130.90) 0.220 0.38 (0.081.89)c.59C>T 3/42/150 3/31/153 2/56/125 0.130 0.69 (0.431.11) 0.004 0.49 (0.300.79) 0.190 1.42 (0.842.31)c.102C>T rs1049331 108/64/23 70/81/34 38/90/55 1.16105 3.21 (1.885.50) 0.009 1.91 (1.173.11) 0.072 1.68 (0.952.99)c.108G>T rs2293870 108/64/23 70/81/34 38/90/55 1.16105 3.21 (1.885.50) 0.009 1.91 (1.173.11) 0.072 1.68 (0.952.99)

Table 2. The odds ratio analysis of associated HTRA1 polymorphisms. The genotypic distribution was presented as homozygous/heterozygous/wild type. The odds ratios were assessed under dominant model.

Polymorphism db SNP ID

Genotypic frequency AMD PCV AMD-PCV AMD

(n=195)

Figure 1. Meta-analysis of rs2672598 in 3 Chinese PCV cohorts. Three Chinese cohorts (Hong Kong, Shantou and Beijing) were included in the meta-analysis under dominant model. Random eect was chosen. Squares indicate study-specic OR. The size of the box is proportional to the weight of the study. Horizontal lines indicate 95% CI. Diamond indicates summary OR with its corresponding 95% CI.

that rs2672598 is associated with exudative AMD (p=4.60103), but not PCV (p=0.24)17. In order to further illustrate the association of rs2672598 with PCV, we genotyped another PCV cohort in southern China. Unlike the Hong Kong and Beijing cohort, rs2672598 was associated with PCV in Shantou cohort (p=0.031, OR=2.49, 95% CI: 1.065.83). The meta-analysis of the 3 Chinese cohorts at rs2672598 also showed a signicant association with PCV (p=0.02, OR= 1.94, 95% CI: 1.133.32; Fig.1). In addition, c.59C>T (A20V) was signicantly associated with PCV (p=0.004, OR= 0.49, 95% CI: 0.300.79), but not with exudative AMD (p=0.13).

LD and haplotype-based analyses of the 6 associated variants showed that rs11200638, rs1049331, and rs2293870 were in perfect LD (r2 > 0.95) in both exudative AMD and PCV. LD blocks across all variants were observed in exudative AMD (Fig.2A), PCV (Fig.2B) and in control subjects (Fig.2C), except c.34delCinsTCCT probably because of its low frequency. It was excluded for haplotype-based association analysis. Three associated haplotypes dened by these 5 variants were predicted (Table3). Haplotype ACCTT, dened by the risk alleles, was associated with increasing risk in both exudative AMD (p=8.701014, OR= 3.10, 95% CI: 2.294.19) and PCV (p=1.60104, OR= 1.66, 95% CI: 1.242.22). Haplotype GTCCG, dened by the wildtype alleles, was associated with decreasing risk in exudative AMD (p=5.351012, OR= 0.25, 95% CI: 0.160.38), but not in PCV (p= 0.17). These 2 haplotypes were dierent between AMD and PCV (p < 0.001). Another haplotype, GCTCG, was associated with PCV (p=0.04, OR= 0.57, 95% CI: 0.370.89), but not with exudative AMD (p=0.363).

The 4 associated variants between exudative AMD and PCV were examined by joint genotype analysis SNP rs11200638 was used as a proxy because it was in perfect LD with other two SNPs (rs1049331 and rs2293870). Therefore, only rs11200638 and rs2672598 were included in the analysis. The results of the joint genotypes included homozygous carriers of rs2672598-C risk allele. All genotypes of rs11200638 were associated with exudative AMD (pAA-CC = 4.31 108, pGA-CC = 3.00 103 and pGG-CC = 2.00 103,

respectively). Significant differences were found between exudative AMD and PCV (pAA-CC = 4.00 103,

OR = 11.70, 95% CI: 1.4594.53; pGA-CC = 0.01, OR = 10.16, 95% CI: 1.2284.81 and pGG-CC = 0.02, OR = 11.25, 95% CI: 1.17108.41, respectively). However, only the homozygous risk allele carriers were associated with PCV (p=3.00103). The highest OR in homozygous risk allele joint genotypes for exudative AMD (OR=43.11, 95% CI: 5.56334.49) has 11.7-fold higher risk (p=4.00103) than for PCV (OR= 3.68, 95% CI: 1.518.99) (Fig.3 and Supplementary Table 1). Logistic regression analysis showed that rs11200638 (p=1.48104) and rs2672598 (p=2.27103) were independently associated with exudative AMD aer adjusting for age and gender. On the contrary, in PCV, rs11200638 remained signicant (p=9.05103) but not rs2672598 (p=0.20).

Cloning the HTRA1 promoter identied three rs11200638-rs2672598 haplotypes: G-T (wildtype), G-C and A-C. No A-T haplotype was found. Luciferase expression analysis showed that the rs11200638-rs2672598 haplotypes aect dierential transcriptional activities of the HTRA1 promoter (Fig.4). Comparing to G-T haplotype, elevated transcriptional activity was observed in G-C and A-C haplotypes

PCV

(n=188)

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

Figure 2. Haplotype block structure for the associated HTRA1 polymorphisms. The haplotype analysis revealed linkage disequilibrium (LD) blocks lying across all of the associated HTRA1 polymorphisms except for the c.34delCinsTCCT in (A) exudative AMD, in (B) PCV and in (C) control subjects. The displayed value of the LD was D (the value of D100, threshold= 0.8), and the r2 threshold was 0.8.

Figure 3. The odds ratio plot of rs11200638rs2672598 joint genotypes. The rs11200638-A and rs2672598-C were the risk alleles. The odds ratio plot revealed the highest joint odds ratio of 43.11 from the AA-CC joint genotype in exudative AMD group. In addition, >10-fold higher risk was observed in exudative AMD, comparing to PCV.

AMD vs Control PCV vs Control AMD vs PCV

Ratio p OR (95% C.I.) Ratio p OR (95% C.I.) Ratio p OR (95% C.I.) ACCTT 0.718:0.451 <0.0001 3.10 (2.294.19) 0.577:0.450 0.0007 1.66 (1.242.22) 0.718:0.578 <0.0001 1.87 (1.382.52) GTCCG 0.090:0.284 <0.0001 0.25 (0.160.38) 0.219:0.283 0.0128 0.71 (0.511.00) 0.089:0.217 <0.0001 0.35 (0.230.54) GCTCG 0.121:0.164 0.1677 0.70 (0.461.06) 0.100:0.163 0.0500 0.57 (0.370.89) 0.120:0.100 1.0000 1.22 (0.781.92) GCCCG 0.069:0.096 0.3826 0.70 (0.421.19) 0.081:0.098 0.1232 0.65 (0.381.11) 0.070:0.066 1.0000 1.04 (0.591.84)

Table 3. Haplotype-based association analysis of HTRA1 polymorphisms. The haplotypes were denedby rs11200638, rs2672598, c.59C>T, rs1049331 and rs2293870. The p-value was corrected by permutation test (n=10,000).

(p = 0.009, one-way ANOVA; Post-hoc tests adjusted by Tukey HSD: G-C vs G-T, p = 0.026, and A-C vs G-T, p= 0.009) with 1.78 0.35 and 1.99 0.30 fold increase, respectively. However, there was no dierence in lucif-erase expression level between G-C and A-C haplotypes (p= 0.621), suggesting that A-allele of rs11200638 might not alter the transcriptional activity of the HTRA1 promoter. Therefore, rs2672598 should be the variant that contributed to the altered transcriptional activity of the HTRA1 promoter.

Our luciferase-reporter analysis showed that transcription activity increased with rs2672598 C-allele but not rs11200638 A-allele. This phenomenon may be due to inuences of transcription factor binding sites within these regions. With G-T haplotype as a reference, one transcription factor E2F-1 binding site was predicted at rs11200638, and one each of STAT4, NFB, c-Ets-1, RelA, Elk-1 and WT1 binding sites were predicted at rs2672598 site (Fig.5). When rs11200638 G-allele was replaced by A-allele, the transcription factor E2F-1 binding site was unchanged. In contrast, when rs2672598 T-allele was replaced by C-allele, all the transcription factor binding sites were abolished and changed to GR-alpha, AP-2A and Sp1 binding sites. Bioinformatic analysis conrmed that rs2672598 contributed to the altered transcriptional activity of the HTRA1 promoter.

Haplotypes

SCIENTIFIC REPORTS

4

www.nature.com/scientificreports/

Figure 4. Luciferase expression driven by the HTRA1 promoter with dierent HTRA1 rs11200638-rs2672598 haplotypes. A 726-bp genomic fragment (from 658 to +68) containing the HTRA1 promoter was cloned into empty pGL3-Basic vector. Human retinal pigment epithelial cell line ARPE-19were transfected with dierent luciferase constructs for 72 hours. (A) Detection of luciferase expression was performed by immunoblotting. (B) Normalized luciferase intensities were calculated by dividing the quantied luciferase intensities by the housekeeping -actin intensities. Two dierent clones of each haplotype were used in the transfection experiment. The wildtype rs11200638-rs2672598 haplotype (G-T) was set as a reference. Recombinant rey luciferase was used as a positive control, whereas empty pGL3 was used as negative control.

In order to validate the elevated transcriptional activity of rs2672598 C-allele, we performed the genotype-expression correlation analysis in vitreous humor collected from our previous study18, which HTRA1 protein expression has been determined. In this study, the rs2672598 genotype of the vitreous humor samples was determined, and 5 samples were in TT genotype, 23 samples in TC genotype and 27 samples in CC genotypes. The HTRA1 expression was higher in the vitreous humor with rs2672598 CC genotype than that with TT genotype by 2.56 folds (p=0.009, one-way ANOVA; Post-hoc tests adjusted by Tukey HSD: CC vs TT, p= 0.02; Fig.6). HTRA1 level of TC geno-type is also higher than TT genotype by 2.03 folds but not statistically signicant (p= 0.32). As we have previously shown than rs11200638 is not correlated with HTRA1 expression in vitreous humor18, this result conrmed that rs2672598 should be the variant leading to the dierential HTRA1 expression in human ocular samples.

PCV has been regarded as a dierent disease entity from AMD, with a distinct set of demographic, pathological and clinical characteristics3. However, emerging evidences suggest that PCV could be a variant of typical CNV5, sharing common genetic determinants and environmental risk factors6. Moreover, elevated serum C reactive protein and VEGF levels in CNV were also found in PCV19,20. In clinical practice, PCV usually responds better to PDT8,11, whereas exudative AMD responds better to anti-VEGF treatment12,21. Meanwhile, they share some common genetic determinants. ARMS2 rs10490924 and HTRA1 rs11200638 were associated with both exudative AMD and PCV in a Japanese study, in which neither allelic nor genotypic frequencies showed signicant dierence22. Similar association of rs10490924 with advanced AMD and PCV has also been reported in Caucasians23. In contrast, our previous study identied signicantly dierent association of ARMS2 polymorphisms between exudative AMD and PCV in the Chinese population15.

Results in this study revealed a HTRA1 promoter variant (rs2672598) that is signicantly associated with exudative AMD, but not PCV (Tables1 and 2). Our result is coherent to a recent report in a northern Chinese population17. Although meta-analysis of the 3 Chinese cohorts at rs2672598 also showed a signicant association with PCV (p = 0.02, OR = 1.94, 95% CI: 1.133.32; Fig.1), rs2672598 is more signicantly associated with exudative AMD than PCV as shown in Hong Kong (exudative AMD: p = 1.31 104; PCV: p = 0.11) and Beijing cohorts (exudative AMD: p=4.60103; PCV: p= 0.24). Furthermore, the joint genotype constructed by homozygous carriers of risk allele (rs2672598-C and rs11200638-A) were prone to exudative AMD (OR=43.11, 95% CI: 5.56334.49), conferring 11.7-fold higher risk (p = 4.00 103) when compared to PCV (OR = 3.68, 95% CI: 1.518.99; Fig.3 and Supplementary Table 1). Moreover, joint genotypes constructed by rs2672598-CC and rs11200638-GA/GG conferred increasing AMD risk (OR= 12.86, 95% CI: 1.64100.86; OR= 18.00, 95% CI: 2.00161.83, respectively), but not associated with PCV. Logistic regression results further indicated dierential

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

Figure 5. Transcription factor binding site prediction at rs11200638 and rs2672598 in the HTRA1 promoter. DNA sequence of the cloned HTRA1 promoter was used to predict the transcription factor binding sites. PROMO was used to predict the transcription factors that would bind to the region of rs11200638 and rs2672598 in the HTRA1 promoter. Predictions with dierent alleles of rs11200638 and rs2672598 were performed. Transcription factors: 10: GR-alpha; 12: E2F-1; 13: AP-2A; 22: STAT4; 23: NFB; 24: c-Ets-1; 25: RelA; 26: Elk-1; 38: Sp1; 39: WT1.

Figure 6. Correlation analysis of rs2672598 genotypes and HTRA1 protein expression in human vitreous humor samples. HTRA1 expressions in 55 vitreous humor samples were determined by immunoblotting, and expressed in intensity (count/mm2). Genotyping of rs2672598 polymorphism revealed TT genotype in 5 samples, TC genotype in 23 samples and CC genotypes in 27 samples. The HTRA1 expression of CC genotype (homozygous) was signicantly higher than that of TT genotype (wildtype).

roles of rs2672598 in exudative AMD and PCV. This SNP might contribute independently of rs11200638 to exu-dative AMD but not PCV.

HTRA1 has its own role in both AMD and PCV24,25, governed by dierent molecular determinants. In this study, rs2672598, located in the promoter region of HTRA1, not only showed dierential associations between AMD and PCV, but also altered the transcriptional activity of the HTRA1 promoter. Higher transcriptional activity of the HTRA1 promoter was observed in the constructs with rs2672598 C-allele rather than that with

SCIENTIFIC REPORTS

6

www.nature.com/scientificreports/

rs11200638 A-allele (Fig.4), suggesting that the rs2672598 C-allele, not rs11200638 A-allele, is responsible for the alteration in the transcriptional activity of the HTRA1 promoter. Our luciferase-reporter analysis agreed with the transcription factor binding site prediction analysis (Fig.5). No alteration of transcription factor binding site was found in rs11200638 A-allele when compared to the wildtype G-allele. In contrast, the transcription factor binding sites in rs2672598 C-allele were totally dierent from that in the wildtype T-allele. This further conrmed that rs2672598 should be the variant contributed to the altered transcriptional activity of the HTRA1 promoter. In order to validate the results of in vitro transcriptional activity, we performed a correlation analysis on the HTRA1 protein expression in vitreous humor with the rs2672598 genotypes. Our results showed that HTRA1 expression was higher in the vitreous humor with rs2672598 CC genotype than that with TT genotype by 2.56 folds (p=0.02; Fig.6). As we previously showed than rs11200638 is not correlated with HTRA1 expression in vitreous humor18,

the vitreous humor expression results in this study not only validate our ndings in the promoter analysis, but also conrm that rs2672598 is responsible for the higher HTRA1 expression in human ocular samples. Our lucif-erase activity and vitreous humor expression results of rs11200638 are similar to previous reports that there is no signicant change of HTRA1 gene expression level in human retina among genotypes of HTRA1 rs11200638, rs1049331 and rs2293870, and ARMS2 rs10490924 and the del443ins54 variant26,27.

In summary, our results revealed that a HTRA1 promoter variant, rs2672598, is more signicantly associated with exudative AMD than PCV. It should be responsible for elevated transcriptional activity in the HTRA1 promoter.

We recruited 188 PCV patients from the Prince of Wales Hospital Eye Centre and the Hong Kong Eye Hospital in Hong Kong. Together with 195 exudative AMD patients and 183 normal controls reported in our previous study28, a total of 534 unrelated Han Chinese were included (Supplementary Table 2). For validation, we recruited 187 PCV patients and 190 control subjects from the Joint Shantou International Eye Center of Shantou University and The Chinese University of Hong Kong, Shantou, China. All study subjects, including patients and controls, were given comprehensive ophthalmic evaluations and examined by clinical ophthalmologists. Both uorescein angiography (FA) and ICGA were performed in the exudative AMD and PCV patients. AMD was graded according to an international classication and grading system29. Patients with exudative AMD had non-drusenoid RPE detachment, choroidal neovascularization, serous or hemorrhagic retinal detachments, sub-retinal or sub-RPE hemorrhage or brosis. PCV patients had sub-retinal red or orange nodules and hemorrhagic pigment epithelial detachment and characteristic sacculated vascular abnormalities in the inner choroid as visualized on ICGA. The diagnosis of PCV was distinguished from AMD by FA and ICGA staining30. PCV patients with geographic atrophy or early signs of AMD were excluded, together with patients who were difficult to be clearly dierentiated (at late stage of disease that possess brosis and disciform scar or at advanced stage that present with intensive hemorrhage). The control subjects were recruited from elderly people greater than 60 years. They did not have any identiable signs of AMD, PCV or other major eye diseases, except for mild senile cataract and slight refractive errors ranging from 1.5 to +1.5 diopters. The study protocol, approved by the Ethics Committee for Human Research at the Chinese University of Hong Kong, is in accordance with the tenets of the Declaration of Helsinki. Informed consent was obtained from each study subject.

Genomic DNA from whole blood was extracted (Qiagen QIAamp DNA Blood Mini kit, Qiagen, Hiden, Germany) according to the suppliers instructions. The 9 exons, intron-exon junctions of HTRA1 (ENSG00000166033) and its promoter region were screened by polymerase chain reaction (PCR) with specic primers followed by direct sequencing (BigDye Terminator Cycle Sequencing Reaction Kit, v3.1; Applied Biosystems, Foster City, CA) on a DNA sequencer (ABI 3130XL, Applied Biosystems)28.

All the identied polymorphisms were assessed for Hardy-Weinberg equilibrium using 2 analysis. Allelic and genotypic distributions for association between PCV and controls, between exudative AMD and controls, and between exudative AMD and PCV) were compared using the 2 test or Fishers exact test (SPSS, version 16.0; SPSS Science, Chicago, IL). Logistic regression was used to examine the association proles of the associated polymorphisms in HTRA1 between exudative AMD and PCV (SPSS). LD and haplotype-based association analyses were performed based on Haploview, version 4.2; http://www.broadinstitute.org/)31.

Multiple testing correction was performed by permutation test (n=10,000).

For the meta-analysis, the Review Manager soware (RevMan, version 5.2, The Cochrane Collaboration, Copenhagen, Denmark) was used to generate the Z scores. The exact P values with the Z scores was then calculated using R (v3.0.0, http://cran.r-project.org/). The I2 statistic was used to assess the heterogeneity among studies, which corresponds to no (<25%), low (2550%), moderate (5075%), and high heterogeneity (75%). We adopted the random eects model in the meta-analysis. Summary p< 0.05 was considered statistically signicant.

The genotype combinations of rs11200638 and rs2672598 in exudative AMD, PCV and control subjects were counted. Their corresponding odds ratios (OR) were calculated by the 2 test (SPSS). The ORs were compared to the baseline genotype of the two genes that showed the lowest frequency of the disease risk alleles (homozygous G carrier of rs11200638 with homozygous T carrier of rs2672598).

A 726-bp genomic fragment (from 658 to +68) containing the HTRA1 promoter was cloned into empty pGL3-Basic vector, pGL3 (Promega, Madison, WI) between the NheI and BglII sites (OriGene Technologies, Rockville, MD). Genomic DNA of the study subjects were amplified by PCR (Ex Taq DNA Polymerase, TakaRa, Japan): Forward primer (5-TAATGCTA GCTCTCTGCGAATACGGACACG), reverse

SCIENTIFIC REPORTS

7

www.nature.com/scientificreports/

primer (5-TAATAGATCTGGGAGAGTGCAG GAGGG). Constructs with dierent rs11200638-rs2672598 haplotypes were generated. The cloned sequence of all constructs was veried by direct sequencing.

Human retinal pigment epithelial cell line ARPE-19 (American Type Culture Collection, Manassas, VA)32

was cultured in Dulbeccos modied Eagles medium and F-12 nutrient mixture (DMEM/F-12) supplemented with 10% fetal bovine serum (Gibco BRL, Rockville, MD). Cells were plated in 60mm tissue culture dishes at a density of 23 105 cells per dish one day before transfection. Aer 24 hours, cells were transfected with 2.5g of luciferase constructs in 7.5 l FuGene HD transfection reagent (Roche, Indianapolis, IN) per dish. Empty pGL3 was used as negative control. At 72 hours aer transfection, cell lysates were extracted by RIPA reagent (Sigma-Aldrich, St. Louis, MO) for detection of luciferase expression.

Luciferase expression was detected by immunoblotting according to our published method33. Briey, the denatured cell lysates of the transfected cells were resolved on 10% SDS-polyacrylamide gel and electro-transferred to nitrocellulose membranes for probing with a mouse monoclonal primary antibody against rey luciferase (LifeSpan BioSciences, Seattle, WA) and a secondary antibody against mouse IgG conjugated with horseradish peroxidase (Santa Cruz Biotechnology, Dallas, TX). Signals were detected by the enhanced chemiluminescence (ECL) system (Amersham Pharmacia, Cleveland, OH) and quantied by ChemiDoc (BioRad, Hercules, CA). Normalized luciferase intensities were calculated by dividing the quantied luciferase intensities by the housekeeping -actin intensities. Two dierent clones of each haplotype were used in the transfection experiment. Triplicated experiments were performed. The wildtype rs11200638-rs2672598 haplotype (G-T) was set as a reference. One-way Analysis of Variance (ANOVA) with post-hoc tests adjusted by Tukey HSD (SPSS) was used to compare the means among dierent groups. Recombinant rey luciferase (Promega) was used as a positive control.

DNA sequence of the cloned HTRA1 promoter was used to predict the transcription factor binding sites. We predicted the transcription factors that would bind to the region of rs11200638 and rs2672598 in the HTRA1 promoter by PROMO (version 3.0.2; http://alggen.lsi.upc.es/cgi-bin/promo_v3/promo/promoinit.cgi?dirDB= TF_8.3)34. Predictions with dierent alleles of rs11200638 and rs2672598 were performed.

Fiy-ve unrelated Chinese patients underwent ocular surgeries at the Prince of Wales Hospital in Hong Kong and the Joint Shantou International Eye Center (JSIEC) of Shantou University and the Chinese University of Hong Kong were recruited and given complete ophthalmoscopic examinations18. In all patients, a standard three-port pars plana vitrectomy was performed as a part of the regular surgical procedures. Undiluted vitreous humor samples (0.51ml) were collected into sterile tubes at the time of surgery and aliquots were rapidly frozen at 80C until assay.

HTRA1 protein expression has been previously determined18. Briey, total protein concentrations in the vitreous samples were measured by Protein assay (BioRad, Hercules, CA). Equal amount of total protein (10ug) for each denatured vitreous humor samples were resolved on 12.5% SDS-polyacrylamide gel and electro-transferred to nitrocellulose membranes for probing with the mouse monoclonal antibody against HTRA1 (R&D Systems Inc., Minneapolis, MN), and secondary antibody against mouse IgG conjugated with horseradish peroxidase (Jackson Immuno. Res., West Grove, PA). The signals were detected by enhanced chemiluminescence (ECL) (Amersham Pharmacia, Cleveland, OH) and the band intensities quantied by Quantity One Image Analysis

soware (BioRad).

DNA was extracted from all of the vitreous humor samples using a commercially available DNA extraction kit (Qiagen, Germany). The HTRA1 rs2672598 genotypes were determined in all 55 samples by polymerase chain reaction (PCR) and direct DNA sequencing as described above28.

1. Pascolini, D. et al. 2002 global update of available data on visual impairment: a compilation of population-based prevalence studies. Ophthalmic Epidemiol 11, 67115 (2004).

2. Jager, R. D., Mieler, W. F. & Miller, J. W. Age-related macular degeneration. N Engl J Med 358, 26062617 (2008).3. Ciardella, A. P., Donso, I. M., Huang, S. J., Costa, D. L. & Yannuzzi, L. A. Polypoidal choroidal vasculopathy. Surv Ophthalmol 49, 2537 (2004).

4. Ambati, J., Ambati, B. K., Yoo, S. H., Ianchulev, S. & Adamis, A. P. Age-related macular degeneration: etiology, pathogenesis, and therapeutic strategies. Surv Ophthalmol 48, 257293 (2003).

5. Imamura, Y., Engelbert, M., Iida, T., Freund, K. B. & Yannuzzi, L. A. Polypoidal choroidal vasculopathy: a review. Surv Ophthalmol 55, 501515 (2010).

6. Laude, A. et al. Polypoidal choroidal vasculopathy and neovascular age-related macular degeneration: same or dierent disease? Prog Retin Eye Res 29, 1929 (2010).

7. Yannuzzi, L. A. et al. Polypoidal choroidal vasculopathy and neovascularized age-related macular degeneration. Arch Ophthalmol 117, 15031510 (1999).

8. Lai, T. Y., Chan, W. M., Liu, D. T., Luk, F. O. & Lam, D. S. Intravitreal bevacizumab (Avastin) with or without photodynamic therapy for the treatment of polypoidal choroidal vasculopathy. Br J Ophthalmol 92, 661666 (2008).

9. Gomi, F. & Tano, Y. Polypoidal choroidal vasculopathy and treatments. Curr Opin Ophthalmol 19, 208212 (2008).10. Kwok, A. K., Lai, T. Y., Chan, C. W., Neoh, E. L. & Lam, D. S. Polypoidal choroidal vasculopathy in Chinese patients. Br J Ophthalmol 86, 892897 (2002).

11. Honda, S. et al. Photodynamic therapy for typical age-related macular degeneration and polypoidal choroidal vasculopathy: a 30-month multicenter study in Hyogo, Japan. Jpn J Ophthalmol 53, 593597 (2009).

12. Lim, J. Y. et al. Intravitreal bevacizumab alone versus in combination with photodynamic therapy for the treatment of neovascular maculopathy in patients aged 50 years or older: 1-year results of a prospective clinical study. Acta Ophthalmol 90, 6167 (2012).13. Maruko, I., Iida, T., Saito, M., Nagayama, D. & Saito, K. Clinical characteristics of exudative age-related macular degeneration in Japanese patients. Am J Ophthalmol 144, 1522 (2007).

14. Nakata, I. et al. Signicance of C2/CFB Variants in Age-related Macular Degeneration and Polypoidal Choroidal Vasculopathy in a Japanese Population. Invest Ophthalmol Vis Sci 53, 794798 (2012).

SCIENTIFIC REPORTS

8

www.nature.com/scientificreports/

15. Liang, X. Y. et al. Dierentiation of exudative age-related macular degeneration and polypoidal choroidal vasculopathy in the ARMS2/HTRA1 locus. Invest Ophthalmol Vis Sci 53, 31753182 (2012).

16. Gotoh, N. et al. Haplotype analysis of the ARMS2/HTRA1 region in Japanese patients with typical neovascular age-related macular degeneration or polypoidal choroidal vasculopathy. Jpn J Ophthalmol 54, 609614 (2010).

17. Cheng, Y. et al. Genetic and functional dissection of ARMS2 in age-related macular degeneration and polypoidal choroidal vasculopathy. PLoS One 8, e53665 (2013).

18. Ng, T. K. et al. Interactive expressions of HtrA1 and VEGF in human vitreous humors and fetal RPE cells. Invest Ophthalmol Vis Sci 52, 37063712 (2011).

19. Tong, J. P. et al. Aqueous humor levels of vascular endothelial growth factor and pigment epithelium-derived factor in polypoidal choroidal vasculopathy and choroidal neovascularization. Am J Ophthalmol 141, 456462 (2006).

20. Kikuchi, M. et al. Elevated C-reactive protein levels in patients with polypoidal choroidal vasculopathy and patients with neovascular age-related macular degeneration. Ophthalmology 114, 17221727 (2007).

21. Matsumiya, W. et al. Early responses to intravitreal ranibizumab in typical neovascular age-related macular degeneration and polypoidal choroidal vasculopathy. J Ophthalmol 2011, 742020 (2011).

22. Kondo, N., Honda, S., Ishibashi, K., Tsukahara, Y. & Negi, A. LOC387715/HTRA1 variants in polypoidal choroidal vasculopathy and age-related macular degeneration in a Japanese population. Am J Ophthalmol 144, 608612 (2007).

23. Lima, L. H. et al. Three major loci involved in age-related macular degeneration are also associated with polypoidal choroidal vasculopathy. Ophthalmology 117, 15671570 (2010).

24. Jones, A. et al. Increased expression of multifunctional serine protease, HTRA1, in retinal pigment epithelium induces polypoidal choroidal vasculopathy in mice. Proc Natl Acad Sci USA 108, 1457814583 (2011).

25. Ng, T. K., Liang, X. Y. & Pang, C. P. HTRA1 in age-related macular degeneration. Asia Pac J Ophthalmol 1, 5163 (2012).26. Kanda, A. et al. Age-related macular degeneration-associated variants at chromosome 10q26 do not signicantly alter ARMS2 and HTRA1 transcript levels in the human retina. Mol Vis 16, 13171323 (2010).

27. Wang, G. et al. Variants at chromosome 10q26 locus and the expression of HTRA1 in the retina. Exp Eye Res 112, 102105 (2013).28. Tam, P. O. et al. HTRA1 variants in exudative age-related macular degeneration and interactions with smoking and CFH. Invest Ophthalmol Vis Sci 49, 23572365 (2008).

29. Bird, A. C. et al. An international classification and grading system for age-related maculopathy and age-related macular degeneration, The International ARM Epidemiological Study Group. Surv Ophthalmol 39, 367374 (1995).

30. Cackett, P., Wong, D. & Yeo, I. A classication system for polypoidal choroidal vasculopathy. Retina 29, 187191 (2009).31. Barrett, J. C., Fry, B., Maller, J. & Daly, M. J. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263265 (2005).

32. Dunn, K. C., Aotaki-Keen, A. E., Putkey, F. R. & Hjelmeland, L. M. ARPE-19, a human retinal pigment epithelial cell line with dierentiated properties. Exp Eye Res 62, 155169 (1996).

33. Ng, T. K. et al. AC and AG dinucleotide repeats in the PAX6 P1 promoter are associated with high myopia. Mol Vis 15, 22392248 (2009).

34. Messeguer, X. et al. PROMO: detection of known transcription regulatory elements using species-tailored searches. Bioinformatics 18, 333334 (2002).

We express our greatest appreciation to all the participants in the study. This study was supported by the Health and Medical Research Fund (project number: 12130791 to T.K.N.), the General Research Fund from the Research Grants Council (grant number: 473410 to C.P.P.), Hong Kong, and Direct Grants from the Medical Panel (grant number: 2041771 and 4054029 to C.P.P.), The Chinese University of Hong Kong.

T.K.N., X.Y.L. and L.M. mainly performed the experiments. P.O.S.T. and J.X.W. helped with the experiments. T.Y.Y.L. and H.C. recruited study subjects. C.P.P. and L.J.C. supervised the project. T.K.N., X.Y.L. and C.P.P. interpreted the results and wrote the manuscript.

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing nancial interests: The authors declare no competing nancial interests.

How to cite this article: Ng, T. K. et al. HTRA1 promoter variant dierentiates polypoidal choroidal vasculopathy from exudative age-related macular degeneration. Sci. Rep. 6, 28639; doi: 10.1038/srep28639 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

9