Association between variant Y402H in age-related macular degeneration (AMD) susceptibility gene CFH and treatment response of AMD: a meta-analysis.

PURPOSE: To investigate the association between polymorphism rs1061170 (T1277C, Y402H) in age-related macular degeneration (AMD) susceptibility gene Complement Factor H (CFH) and treatment response of neovascular AMD.
METHODS: We performed a literature-based meta-analysis including 10 published association studies involving 1,510 patients. Treatments included anti-VEGF (bevacizumab and ranibizumab) or photodynamic therapy. Summary odds ratios (ORs) and 95% confidence intervals (CIs) were estimated using fixed- and random-effects models. Q-statistic test was used to assess heterogeneity.
RESULTS: Polymorphism rs1061170 showed a significant summary OR of 1.68 (95% CI, 1.09 to 2.60; P = 0.020; CC versus TT; random-effects) for treatment response of neovascular AMD with heterogeneity of 0.09. In subgroup analysis, rs1061170 was more likely to be a predictor of response to anti-VEGF therapy (P = 0.011). However, heterozygous TC genotype was not associated with altered treatment response (OR = 1.18, 95% CI, 0.95 to 1.47; P = 0.145; fixed-effects). Influence analysis indicated the robustness of our findings.
CONCLUSIONS: rs1061170 might be associated with treatment response of neovascular AMD, especially for the anti-VEGF agents. It might be the first meta-analytically confirmed genetic marker predictive for AMD treatment response though a further validation in larger studies is needed.

Introduction

Age-related macular degeneration (AMD) is the leading cause of irreversible blindness among elderly population. There are two types of AMD, one is non-exudative (dry or atropic) AMD, and the other is exudative (wet or neovascular). The later form, neovascular AMD, is characterized by the presence of choroidal neovascularization (CNV) beneath the fovea. [1], [2] CNV is a pathologic ocular occurrence in which aberrant blood vessels expand from the choroid to the retinal pigment epithelium, reaching the retina in a high proportion of patients. In addition, while traditionally neovascular AMD was characterized by CNV, advances in imaging technology have also identified retinal angiomatous proliferative (RAP) lesions, which in part originate in the retina in some patients. Neovascular AMD is caused by genetic and various environmental factors such as age, smoking, and serum level of cholesterol. [1] Over the past several years, some genes have been identified to be AMD susceptibility gene. [1], [2] Of these, the Complement Factor H (CFH) gene is one of the most important genes. [3], [4], [5] Variant rs1061170 (T1277C, Y402H) in the CFH gene located in the heparin and CRP-binding domain may cause complement dysregulation and lead to the pathogenesis of AMD. [2], [6] Many independent studies as well as a meta-analysis have indicated that harboring the risk C-allele results in an approximately 2-fold increased risk of AMD. [3], [4], [5], [7]

The vision loss due to CNV development in neovascular AMD can be mainly managed by two different therapeutic strategies. One is photodynamic therapy (PDT) with verteporfin, and the other is intravitreal administration of drugs acting against vascular endothelial growth factor (VEGF), known as anti-VEGF treatment. Anti-VEGF agents (i.e., ranibizumab and bevacizumab) have shown efficacy in counteracting the damage owing to choroidal neovessels. Many efforts are underway to identify clinical, genetic, and pharmacologic biomarkers that could predict response to therapy, thereby providing important information to clinical decision making and treatment option. New pharmacogenetic knowledge provides data regarding the role of several single nucleotide polymorphisms (SNPs) as genetic predictors of treatment responsiveness to PDT as well as to anti-VEGF therapy. Since immune factors and inflammation are relatively important concepts linked to AMD and the complement system is a key component in the pathogenesis of AMD, [8] they may play a major role in therapeutic interventions. Previous studies proposed that patients with the variant rs1061170 in the CFH gene have higher background levels of inflammation, may continue to affect the disease progression, and probably lead to more rapid recurrence of neovascularization. [9] Patients with the different rs1061170 genotypes may also respond differently to treatment, and even require additional injections of agents.

Thus far, several association studies regarding the predictive role of rs1061170 in treatment response of neovascular AMD have been reported, [9], [10], [11], [12], [13], [14], [15], [16], [17], [18] though the results were inconclusive yet. Therefore, we carried out a meta-analysis focusing on the relationship between polymorphism rs1061170 and treatment response of AMD in order to get a more convincing and precise conclusion.

Methods

Study identification and data extraction

Publication search of this meta-analysis was performed as described previously. [19], [20], [21] Briefly, relevant studies were searched in the PubMed, Medline, and Web of Science database (updated to February-10, 2012) using the following search terms: (“CFH” or “complement factor H”) and (age-related macular degeneration). There were 374 results, 155 of which regarding rs1061170 (T1277C, Y402H) of the CFH gene. Only those published studies in English language with full text articles were included in this meta-analysis; we did not define the minimum number of patients to be included for meta-analysis. The abstracts of those crudely identified 155 articles were reviewed. The inclusion criteria were: (i) evaluating the relationship between rs1061170 and treatment response of neovascular AMD, (ii) independent retrospective or prospective association study, and (iii) with sufficient available data to estimate an odds ratio (OR) with 95% confidence interval (CI). As a result, we identified 10 candidate studies for systematic review ( ). The following variables were extracted from each study if available: first author's surname, publication year, ethnicity, numbers of cases, OR with 95% CI of response to treatment, treatment modality, and information of comparison. The information was collected independently by the two authors (C.H. and X.G.Z.), and any discrepancy were resolved by discussion.

Figure 1

The literature search process.

Statistical methods

For each study, OR with 95% CI were recorded or calculated to assess the relation strength between rs1061170 genotype and treatment response. The pooled OR was calculated by a fixed-effects model (using the Mantel-Haenszel method) or a random-effects model (using the DerSimonian and Laird method) according to the heterogeneity. Heterogeneity assumption was checked by the Q test. If P-value was not more than 0.10, the inter-study heterogeneity was considered to be significant and we would choose the random-effects model to pool the ORs. Otherwise, the fixed-effects model was employed. Because of the limitation of original data, two types of OR were calculated in present meta-analysis: TC (YH) genotype versus TT (YY) genotype, and CC (HH) genotype versus TT (YY) genotype. The potential publication bias was examined visually in a funnel plot of ln[OR] against its standard error (SE), and the degree of asymmetry was tested using Egger's test. We also performed sub-population analysis by treatment and ethnicity. Influence analysis (sensitivity analysis) was conducted by omitting each study to find potential outliers. All of the statistical analyses were performed using Stata/SE version 10.0 (Stata Corporation, College Station, TX).

Results

Systematic review of eligible studies

We identified 10 eligible studies for the present systematic review ( ) [9], [10], [11], [12], [13], [14], [15], [16], [17], [18]. Other four candidate studies [22], [23], [24], [25], although also investigating the relationship between neovascular AMD treatment response and rs1061170 genotype, failed to provide sufficient original data for OR estimation and were thus excluded from meta-analysis. Most reported studies were of a relative small sample size, with mean subject number of 151 (range, 86–265). We recorded the genotypic risks. If some data in the original papers were not presented in an appropriate form for meta-analysis, we estimated the ORs with 95% CIs from the raw data. Regarding ethnicity, seven studies were performed in Caucasians, one in Asians, and the remainders were unknown. The frequencies of variant C-allele of rs1061170 (T is the ancestral allele) were similar among Caucasian studies, from 50% to 58%. In contrast, the Asian study reported a much lower C-allele frequency of 7%, [16] which however was consistent with that in Asian population according to the public dbSNP and HapMap databases. All the subjects included in the meta-analysis are patients, and we need to not check the departure from Hardy-Weinberg equilibrium (HWE). [26], [27] It is because that departure from HWE in patients might indicate a genetic association between studied locus and disease rather than a biased population selection. Testing for HWE in patients is not meaningful for quality control. Of note, due to heterogeneity between studies, some studies did not directly evaluate the treatment response by comparing responders with non-responders, but used the disease progression or recurrence as a surrogate of poor response. We clarified the OR of comparison and definition of good/poor response in . In terms of predictive role of rs1061170 in treatment response, when TC genotype versus TT genotype, six of the ten studies showed C-allele tended to be a predictor of poor response; when CC genotype versus TT genotype, seven of nine studies indicated a resistant role of C-allele in treatment.

Table 1

Characteristics of the studies included in the meta-analysis.

Year	Author	Treatment	Frq*	Total	Ethnicity	OR (95% CI)		Study endpoints	OR in MTA	Definition of good response
						TC vs.TT	CC vs.TT
2007	Brantley [10]	Bev.	0.55	86	CA	0.84 (0.17–4.09)	8.50 (0.93–105.50)	VA (Snellen)	Rspr vs. non-Rspr	Not clearly defined
2007	Seitsonen [11]	PDT	0.58	88	CA	0.66 (0.15–3.08)	0.62 (0.12–3.34)	VA (Snellen/ETDRS)	Rspr vs. non-Rspr	Not clearly defined
2008	Klein [12]	AO&Zn	0.53	239	CA	1.96 (0.71–6.24)	5.09 (1.79–16.45)	Progression	Progr. vs. no Progr.	No progression to advanced AMD after treatment
2009	Feng [13]	PDT	0.54	265	CA	1.62 (0.78–3.32)	1.77 (0.79–3.98)	VA (Snellen)	Rspr vs. non-Rspr	An improved/unchanged VA, or lost <3 lines of vision (final VA ≥6/60)
2009	Lee [9]	Ran.	0.54	156	CA	1.25 (0.94–1.67)	1.37 (1.01–1.87)	VA, recurrence	Recur. vs. no Recur.	No recurrence in follow-up time (the first 70 days excluded)
2011	Kloeckener -Gruissem [14]	Ran.	0.50	122	CA	0.68 (0.26–1.77)	2.67 (0.85–8.58)	VA (ETDRS)	Rspr vs. non-Rspr	Improve VA at ≥75th percentile
2011	Nischler [15]	Bev.	0.48	197	UD	1.26 (0.53–2.95)	2.39 (0.76–8.29)	VA	Rspr vs. non-Rspr	Gained ≥3 lines in reading VA
2011	Tsuchihashi [16]	PDT	0.07	110	AA	1.41 (0.38–5.28)	NA	VA, recurrence	Recur. vs. no Recur.	No recurrence in 1-year after PDT
2012	McKibbin [17]	Ran.	0.50	104	CA	0.41 (0.13–1.19)	0.42 (0.12–1.48)	VA (ETDRS)	Rspr vs. non-Rspr	More than 5 BCVA letter score gain after 6 months
2012	Orlin [18]	Ran/Bev	0.47	143	UD	1.07 (0.45–2.60)	1.42 (0.51–3.97)	VA (Snellen)	Rspr vs. non-Rspr	An improved/unchanged VA in ≥1 eye (most recent VA ≥20/200)

Abbreviations: AA, Asian; AMD, age-related macular degeneration; AO&Zn, antioxidants and Zinc; Bev., bevacizumab; CA, Caucasian; CI, confidence interval; ETDRS, Early Treatment Diabetic Retinopathy Study; Frq, frequency; MAF, minor allele frequency; MTA., meta-analysis; NA, not available; OR, odds ratio; PDT, photodynamic therapy; Progr., progression; Ran., ranibizumab; Recur., recurrence; Rspr, responders; UD, Undeclared; VA, visual acuity.

Frequency of the variant C-allele of rs1061170 (T1277C, Y402H). The ancestral allele of rs1061170 is T according to dbSNP (http://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=1061170).

Meta-analysis

We subsequently meta-analyzed the ten studies for the pooled association between treatment response of neovascular AMD and rs1061170 genotype. The combined sample size was 1,510. presents the results of meta-analysis (ten studies for TC versus TT and nine studies for CC versus TT). The patients harboring homozygous for variant C-allele (CC genotype) seemed to be associated with a reduced response to treatment of neovascular AMD (pooled OR = 1.68, 95% CI, 1.09 to 2.60, random-effects, ). However, heterozygous (TC genotype) was not associated with altered treatment response (pooled OR = 1.18, 95% CI, 0.95 to 1.47, fixed-effects). Due to limitation of original data, we could not do meta-analysis under the dominant or recessive genetic model. In subgroup analysis by treatment, we found that rs1061170 was more likely to be a predictor of anti-VEGF therapy rather than a powerful predictor for PDT effectiveness (P = 0.011; OR = 1.42, 95% CI, 1.08 to 1.86) ( ). When we divided the patients according to ethnicity (Caucasians or not), the association of treatment response with rs1061170 was also observed. CC genotype was associated with a reduced response to treatment of neovascular AMD with an OR of 1.73 (95% CI, 1.05 to 2.86) in Caucasians. Since only one study performed in Asians was reported, we could not conduct a sub-meta-analysis in Asians.

Table 2

Results of meta-analysis for the CFH rs1061170 (T1277C, Y402H) and treatment response of neovascular AMD.

	TC versus TT (n * = 10)				CC versus TT (n * = 9)
	OR	95% CI	P for test	P for heterogeneity	P	95% CI	P for test	P for heterogeneity
Overall
Fixed-effects model	1.18	0.95 to 1.47	0.145	0.597	1.52	1.19 to 1.94	0.001	0.087
Random-effects model	/	/	/		1.68	1.09 to 2.60	0.020	0.087
Treatment (Fixed-effects)
Anti-VEGF	1.12	0.88 to 1.42	0.370	0.419	1.42	1.09 to 1.86	0.011	0.171
PDT	1.38	0.77 to 2.47	0.283	0.576	1.45	0.70 to 3.00	0.318	0.266
Ethnicity (Random-effects)
Caucasians &	1.17	0.91 to 1.50	0.229	0.402	1.73	1.05 to 2.86	0.032	0.055

Abbreviations: CI, confidence interval; PDT, photodynamic therapy; OR, odds ratio; VEGF, vascular endothelial growth factor.

n represents the number of studies included.

including seven studies in Caucasians and one study reported by institution in Austria [15]. The study in Asians [16] is not included; one study reported by institutions in South Korea and the USA [18] is also excluded because of a concern of the mixture of Korean patients in that study.

Figure 2

Forest plot of meta-analysis of association between CFH rs1061170 (T1277C, Y402H) and treatment response of neovascular AMD.

A, Each study is shown by the point estimate of the odds ratio (OR) (the size of the square is proportional to the weight of each study) and 95% confidence interval (CI) for the OR (extending lines). Genotype comparison: CC versus TT. B, Forest plot of meta-analysis according to treatment modality. Anti-VEGF agents indicate bevacizumab and ranibizumab; PDT, photodynamic therapy. Genotype comparison: CC versus TT.

In addition, we evaluated the influence of any individual study on the overall OR for rs1061170 ( ). The results showed that no study fundamentally changed the positive relationship between rs1061170 and treatment response. Moreover, Begg's funnel plot and Egger's test were performed to evaluate the publication bias of literature, and no significant publication bias was observed either in the CC versus TT comparison ( , P = 0.27 for Egger's test) or the TC versus TT comparison (data not shown).

Figure 3

Influence analysis and publication bias plot for meta-analysis.

shows the influence of individual studies on the summary OR (CC genotype versus TT genotype). The horizontal axis represents the odds ratio. Every circle indicates the pooled OR when the left study is omitted in this meta-analysis. The two ends of every broken line represent the respective 95% CI. shows the Begg's funnel plot of studies included in the meta-analysis (genotype comparison: CC versus TT). The vertical axis represents ln[OR] and the horizontal axis means the standard error of ln[OR]. Horizontal line and sloping lines in funnel plot represent summary OR and expected 95% CI for a given standard error, respectively. Area of each circle represents the contribution of each study to the pooled OR.

Discussion

Because relatively small sample sizes in studies have hindered reliable assessment of the association between polymorphism rs1061170 (T1277C, Y402H) in the CFH gene and treatment response of neovascular AMD, we performed the current systematic meta-analysis of 10 relevant studies involving 1,510 patients. The results showed that polymorphism rs1061170 was a predictor of treatment response of neovascular AMD, especially for anti-VEGF agents. Individuals harboring homozygous for the variant risk C-allele corresponded to a decreased response to treatment by approximately 1.6-fold when compared with patients carrying homozygous for the ancestral T-allele. Sensitivity analysis indicated the robustness of our findings, and evidence of publication bias was not observed either. Our study for the first time provides evidence of an association between rs1061170 and treatment response of neovascular AMD.

Complement dysregulation has emerged as an important pathogenetic factor in AMD. CFH is a critical negative regulator of complement activation. [28] Previous biological studies have demonstrated that SNP rs1061170, which is located in the coding sequence, exerts allelic differences on the binding affinity to C-reactive protein (CRP), with the risk allele showing reduced affinity. Previous observations support the functional relevance of rs1061170 to AMD pathogenesis as well as treatment response. A recent meta-analysis and many original papers have provided substantial evidence that rs1061170 is significantly associated with AMD risk. [3], [4], [5], [7] In the present study, our results solidly expand the predictive role of rs1061170 in neovascular AMD treatment response, especially for anti-VEGF therapy.

However, considering the relative low OR of 1.7 with 95% CI of 1.1 to 2.6, using rs1061170 only to predict AMD treatment response is not accurate enough, probably coupling with high false-positivity and high false-negativity. Therefore, combining a series of predictive loci as well as clinical and pathological markers would be a better way. Thus far, some polymorphisms in the CRP (rs2808636, rs876538), VEGF (rs699947, rs2146323), pigment epithelium derived factor (PEDF, rs12603825), HTRA1 (rs11200638), and Factor XIII-A (rs5985) genes were also involved in prediction of neovascular AMD treatment response. Additionally, other polymorphisms (such as rs1410996, rs2274700) and haplotypes in the CFH gene may be also helpful to determine the treatment outcome. [16], [29] However, because of limited study number, those SNPs are not available for meta-analysis yet and further molecular epidemiological studies are needed. In addition, methodologic issues such as publication bias could have particular impact when combining results of observational studies. [30] Our results are purely based on literature that has been published and are not validated by our own data. Although we did not observe apparent publication bias by Eggers test or funnel plots, these analyses however cannot rule out this problem. Further validation in larger studies is needed.

In conclusion, our analysis might provide a piece of evidence that polymorphism rs1061170 is associated with treatment response of neovascular AMD, especially for the anti-VEGF agents (i.e., ranibizumab and bevacizumab). Our results for the first time meta-analytically confirmed a genetic marker predictive for AMD treatment response, and these results, if validated, might promote the potential for clinical treatment response prediction by combining it with other potential clinicopathological and genetic markers in the future.