Age-related macular degeneration (AMD) is a complex disease with both genetic and environmental factors implicated in its development and progression. Twin studies have shown that 46% to 71% of AMD severity can be explained by genetic factors, while environmental variations are responsible for 19% to 37%.1 Furthermore, environmental factors such as smoking, body mass index, serum C-reactive protein, and dietary intake were found to modify AMD genetic susceptibility.2–5
Epigenetics is the study of covalent modifications of the genome that can change gene structure and function without altering its sequence. These mechanisms include DNA methylation, histone acetylation/deacetylation, chromatin remodeling and noncoding RNA-mediated gene silencing. Epigenetic mechanisms can explain the impact of gene-environment interactions on different disease phenotypes.6
The introduction of antiangiogenic agents has revolutionized the management of wet AMD.7 Despite good results obtained by these drugs, 10% to 15% of patients continue to lose vision.8–12 This variability in patient response to therapy has been attributed to several clinical, behavioral, and genetic factors, as well as to possible gene-gene or gene-environmental interactions.13–16
This study examines the effects of antiangiogenic drugs on the transcription of acetylation genes in immortalized human retinal pigment epithelium cells (ARPE-19) in vitro, to determine if epigenetic pathways are modulated by antiangiogenic treatments.
Materials and Methods
Immortalized ARPE-19 cells were cultured until confluent in 175 cm2-flasks containing DMEM/F-12 culture medium (Dulbecco's Modification of Eagle's Medium; Mediatech, Manassas, VA), 10% dialyzed fetal bovine serum, supplementary antibiotics, and glucose. The cells were then plated into six-well plates for 24 hours. Culture media were removed and replaced with the same media containing aflibercept (Eylea; Regeneron, Tarrytown, NY), ranibizumab (Lucentis; Genentech, South San Francisco, CA), or bevacizumab (Avastin; Genentech, South San Francisco, CA) at one (1×) or two times (2×) the concentrations of the clinical intravitreal dose (12.5 μl/ml and 25 μl/ml, respectively). Cells were treated with antivascular endothelial growth factor (VEGF) drugs for 24 hours, and then pelleted for RNA isolation using a PureLink RNA Mini Kit (Ambion; Thermo Fisher Scientific, Waltham, MA). For real-time quantitative reverse transcription polymerase chain reaction (qRTPCR) analyses, 100 ng of individual RNA samples were reverse transcribed into complementary DNA using SuperScript VILO Master Mix (Thermo Fisher Scientific, Waltham, MA).
To identify genes associated with epigenetic pathways, we used the TaqMan Array 96-well plate (Thermo Fisher Scientific, Waltham, MA). This plate contains 28 genes that were specifically related to DNA methylation and transcription repression associated genes along with four candidate endogenous control genes. TaqMan Assay and complementary DNA samples from ARPE-19 cell cultures treated with aflibercept or bevacizumab at 1× concentration were added to each well, and qRT-PCR was performed using a StepOnePlus Real-Time PCR system (Thermo Fisher Scientific, Waltham, MA). Untreated samples were used as control.
After identification of possible target genes with the TaqMan Array, the transcription levels were verified using qRT-PCR in triplicate. We used primers for genes associated with epigenetic acetylation pathways: histone acetyltransferase 1 (HAT1) and histone deacetylases 1, 6 and 11 (HDAC1, HDAC6, and HDAC11). The qRT-PCR was performed on individual samples using a QuantiFast SYBR Green PCR Kit (Qiagen, Germantown, MD) on a StepOnePlus Real-Time PCR system. For the various target genes, housekeeping genes that had comparable amplification efficiencies to the genes of interest were chosen to maximize the accuracy of our differences in cycle threshold (ΔΔCt) values. The housekeeper genes were either hypoxanthine phosphoribosyltransferase 1 (HPRT1) or hydroxymethylbilane synthase (HMBS). Untreated samples were used as control. ΔΔCts were obtained, and folds were calculated using the formula 2^ΔΔCt.
Statistical analyses of the data were performed by unpaired t test using GraphPad Prism, Version 5 (GraphPad Software, La Jolla, CA). A P value less than .05 was considered statistically significant. Untreated samples (controls) were normalized to a value of 100% for comparison to treated samples.
Aflibercept-treated (1×) ARPE-19 cells exhibited changes in expression of all markers of histone deacetylation: HAT1 (1.18-fold, P = .02), HDAC1 (0.74-fold, P < .001), HDAC6 (1.11-fold, P = .04), and HDAC11 (1.27-fold, P = .02). At 2× concentration of aflibercept, only HAT1 (1.21-fold, P = .01) and HDAC11 (1.48-fold, P = .02) were significantly altered.
Bevacizumab-treated (1×) ARPE-19 cells also had significant change in the expression of all four genes compared with untreated cells: HAT1 (1.26-fold, P = .003), HDAC1 (0.73-fold, P < .001), HDAC6 (1.3-fold, P = .004), and HDAC11 (1.52-fold, P = .009). In cultures treated with 2× bevacizumab, expression of only HAT1 (1.97-fold, P = .01) and HDAC11 (1.39-fold, P = .01) were significantly affected compared with controls.
Compared with untreated APRE-19 cells, 1× ranibizumab-treated cells expressed higher levels of HAT1 (1.55-fold, P = .007) and HDAC6 (1.26-fold, P = .02), but none of the other epigenetic acetylation pathway markers were changed. In ARPE-19 cells treated with 2× ranibizumab, HAT1 (4-fold, P = .02) and HDAC11 (1.32-fold, P = .006) were the only epigenetic markers expressed differently compared with untreated cells (Table, Figure).
Fold Expression of Histone Acetylation Genes in Untreated and Anti-VEGF treated ARPE-19 cultures
Expression profiles for histone acetylation genes in untreated and anti-vascular endothelial growth factor treated human retinal pigment epithelium cell (ARPE-19) cultures.
*P < .05; **P < .01; ***P < .001
CT = difference in cycle thresholds; A1X = aflibercept at one time the concentration of clinical intravitreal dose; A2X = aflibercept at two times the concentration of clinical intravitreal dose; B1X = bevacizumab at one time the concentration of clinical intravitreal dose; B2X = bevacizumab at two times the concentration of intravitreal clinical dose; HAT = histone acetyltransferase; HDAC = histone deacetylase; R1X = ranibizumab at one time the concentration of clinical intravitreal dose; R2X = ranibizumab at two times the concentration of clinical intravitreal dose; UNT = untreated.
VEGF plays a key role in the pathogenesis of neovascularization and has been the main target for wet AMD therapy.17 Many studies have attempted to explain the resistance of a subset of patients to anti-VEGF therapy, but to date, pharmacogenetic studies have yielded inconclusive results. Several studies reported a significant association between specific single nucleotide polymorphisms, or SNPs, in genes known to be associated with AMD risk and anti-VEGF treatment outcomes.13,14,18–25 Conversely, other studies, including analyses of the major CATT and IVAN trials, were not able to prove this association.26–30 This inconsistency may be the result of different study designs, or due to the modification of genetic susceptibility by environmental factors.
Epigenetic factors have been implicated in the pathogenesis of AMD by selective transcription of genes involved in angiogenesis, inflammation, and oxidative stress pathways.31,32 These findings led us to the hypothesis that anti-VEGF drugs might interact with the epigenome of retinal cells, leading to differential responses to therapy in different patients. We used immortalized ARPE-19 cells to investigate this potential effect in vitro. Our results show that anti-VEGF drugs can alter the expression profiles of genes regulating the histone acetylation status in ARPE-19 cell line. We demonstrated this effect with both the clinically used concentration (1×) and a twofold concentration (2×).
Histone acetylation and deacetylation affect DNA transcription by two mechanisms: controlling the accessibility of DNA for transcription factors by altering chromatin folding properties, and modifying lysine residues to enhance specific binding sites for recruitment of repressors and activators of gene activity.6 HDAC inhibitors are potent angiogenesis suppressors in vivo and in vitro. Treatment of ARPE-19 cells with HDAC inhibitors increased the expression of clusterin protein; a major regulator of the complement pathway that acts to prevent cytolysis. Furthermore, HDAC inhibitors can significantly reduce ischemic retinal injury.33–35 Based on these findings, regulation of acetylation in the retina might have an important role in neuroprotection.
Our findings suggest the so-called drug resistance or tachyphylaxis may be in part due to the treatment interaction with different epigenome profiles of different patients. Greater understanding of this event could possibly help us in the future to tailor therapies for patients based on an epigenetic analysis, hopefully by means of a simple blood test.
- Seddon JM, Cote J, Page WF, Aggen SH, Neale MC. The US twin study of age-related macular degeneration: relative roles of genetic and environmental influences. Arch Ophthalmol. 2005;123(3):321–327. doi:10.1001/archopht.123.3.321 [CrossRef]
- Seddon JM, George S, Rosner B, Klein ML. CFH gene variant, Y402H, and smoking, body mass index, environmental associations with advanced age-related macular degeneration. Hum Hered. 2006;61(3):157–165. doi:10.1159/000094141 [CrossRef]
- Seddon JM, Gensler G, Rosner B. C-reactive protein and CFH, ARMS2/HTRA1 gene variants are independently associated with risk of macular degeneration. Ophthalmology. 2010;117(8):1560–1566. doi:10.1016/j.ophtha.2009.11.020 [CrossRef]
- Ho L, van Leeuwen R, Witteman JM, et al. Reducing the genetic risk of age-related macular degeneration with dietary antioxidants, zinc, and -3 fatty acids: the Rotterdam Study. Arch Ophthalmol. 2011;129(6):758–766. doi:10.1001/archophthalmol.2011.141 [CrossRef]
- Reynolds R, Rosner B, Seddon JM. Dietary omega-3 fatty acids, other fat intake, genetic susceptibility, and progression to incident geographic atrophy. Ophthalmology. 2013;120(5):1020–1028. doi:10.1016/j.ophtha.2012.10.020 [CrossRef]
- Hjelmeland LM. Dark matters in AMD genetics: epigenetics and stochasticity. Invest Ophthalmol Vis Sci. 2011;52(3):1622–1631. doi:10.1167/iovs.10-6765 [CrossRef]
- Augustin AJ, Offermann I. Emerging drugs for age-related macular degeneration. Expert Opin Emerg Drugs. 2006;11(4):725–740. doi:10.1517/1472818.104.22.1685 [CrossRef]
- Brown DM, Kaiser PK, Michels M, et al. Ranibizumab versus verteporfin for neovascular age-related macular degeneration. N Engl J Med. 2006;355(14):1432–1444. doi:10.1056/NEJMoa062655 [CrossRef]
- Rosenfeld PJ, Brown DM, Heier JS, et al. Ranibizumab for neovascular age-related macular degeneration. N Engl J Med. 2006;355(14):1419–1431. doi:10.1056/NEJMoa054481 [CrossRef]
- Martin DF, Maguire MG, Ying GS, Grunwald JE, Fine SL, Jaffe GJCATT Research Group. Ranibizumab and bevacizumab for neovascular age-related macular degeneration. N Engl J Med. 2011;364(20):1897–1908. doi:10.1056/NEJMoa1102673 [CrossRef]
- Chakravarthy U, Harding SP, Rogers CA, et al. Ranibizumab versus bevacizumab to treat neovascular age-related macular degeneration: one-year findings from the IVAN randomized trial. Ophthalmology. 2012;119(7):1399–1411. doi:10.1016/j.ophtha.2012.04.015 [CrossRef]
- Heier JS, Brown DM, Chong V, et al. Intravitreal aflibercept (VEGF trap-eye) in wet age-related macular degeneration. Ophthalmology. 2012;119(12):2537–2548. doi:10.1016/j.ophtha.2012.09.006 [CrossRef]
- Shah AR, Williams S, Baumal CR, Rosner B, Duker JS, Seddon JM. Predictors of response to intravitreal anti-vascular endothelial growth factor treatment of age-related macular degeneration. Am J Ophthalmol. 2016;163:154–166.e8. doi:10.1016/j.ajo.2015.11.033 [CrossRef]
- Chen G, Tzekov R, Li W, Jiang F, Mao S, Tong Y. Pharmacogenetics of complement factor H Y402H polymorphism and treatment of neovascular AMD with anti-VEGF agents: a meta-analysis. Sci Rep. 2015;5:14517. doi:10.1038/srep14517 [CrossRef]
- Park UC, Shin JY, McCarthy LC, et al. Pharmacogenetic associations with long-term response to anti-vascular endothelial growth factor treatment in neovascular AMD patients. Mol Vis. 2014;20:1680–1694.
- Dedania VS, Grob S, Zhang K, Bakri SJ. Pharmacogenomics of response to anti-VEGF therapy in exudative age-related macular degeneration. Retina. 2015;35(3):381–391. doi:10.1097/IAE.0000000000000466 [CrossRef]
- Grisanti S, Tatar O. The role of vascular endothelial growth factor and other endogenous interplayers in age-related macular degeneration. Prog Retin Eye Res. 2008;27(4):372–390. doi:10.1016/j.preteyeres.2008.05.002 [CrossRef]
- Medina FM, Alves Lopes da Motta A, Takahashi WY, et al. Pharmacogenetic effect of complement factor H gene polymorphism in response to the initial intravitreal injection of bevacizumab for wet age-related macular degeneration. Ophthalmic Res. 2015;54(4):169–174. doi:10.1159/000439172 [CrossRef]
- Hermann MM, van Asten F, Muether PS, et al. Polymorphisms in vascular endothelial growth factor receptor 2 are associated with better response rates to ranibizumab treatment in age-related macular degeneration. Ophthalmoloy. 2014;121(4):905–910. doi:10.1016/j.ophtha.2013.10.047 [CrossRef]
- Zhao L, Grob S, Avery R, et al. Common variant in VEGFA and response to anti-VEGF therapy for neovascular age-related macular degeneration. Curr Mol Med. 2013;13(6):929–934. doi:10.2174/15665240113139990048 [CrossRef]
- Hautamaki A, Kivioja J, Vavuli S, et al. Interleukin 8 promoter polymorphism predicts the initial response to bevacizumab treatment for exudative age-related macular degeneration. Retina. 2013;33(9):1815–1827. doi:10.1097/IAE.0b013e318285cf92 [CrossRef]
- Abedi F, Wickremasinghe S, Richardson AJ, Islam AF, Guymer RH, Baird PN. Genetic influences on the outcome of anti-vascular endothelial growth factor treatment in neovascular age-related macular degeneration. Ophthalmology. 2013;120(8):1641–1648. doi:10.1016/j.ophtha.2013.01.014 [CrossRef]
- Francis PJ. The influence of genetics on response to treatment with ranibizumab (Lucentis) for age-related macular degeneration: the Lucentis genotype study (An American Ophthalmological Society Thesis). Trans Am Ophthalmol Soc. 2011;109:115–156.
- Brantley MA Jr, Fang AM, King JM, Tewari A, Kymes SM, Shiels A. Association of complement factor H and LOC387715 genotypes with response of exudative age-related macular degeneration to intravitreal bevacizumab. Ophthalmology. 2007;114(12):2168–2173. doi:10.1016/j.ophtha.2007.09.008 [CrossRef]
- Smailhodzic D, Muether PS, Chen J, et al. Cumulative effect of risk alleles in CFH, ARMS2, and VEGFA on the response to ranibizumab treatment in age-related macular degeneration. Ophthalmology. 2012;119(11):2304–2311. doi:10.1016/j.ophtha.2012.05.040 [CrossRef]
- Hagstrom SA, Ying G-s, Pauer GJT, et al. Pharmacogenetics for genes associated with age-related macular degeneration in the comparison of AMD treatments trials (CATT). Ophthalmology. 2013;120(3):593–599. doi:10.1016/j.ophtha.2012.11.037 [CrossRef]
- Lotery AJ, Gibson J, Cree AJ, et al. Pharmacogenetic associations with vascular endothelial growth factor inhibition in participants with neovascular age-related macular degeneration in the IVAN study. Ophthalmology. 2013;120(12):2637–2643. doi:10.1016/j.ophtha.2013.07.046 [CrossRef]
- Hagstrom SA, Ying G, Pauer GT, et al. VEGFA and VEGFR2 gene polymorphisms and response to anti-vascular endothelial growth factor therapy: comparison of age-related macular degeneration treatments trials (CATT). JAMA Ophthalmol. 2014;132(5):521–527. doi:10.1001/jamaophthalmol.2014.109 [CrossRef]
- Hagstrom SA, Ying GS, Pauer GJ, et al. Endothelial PAS domaincontaining protein 1 (EPAS1) gene polymorphisms and response to anti-VEGF therapy in the comparison of AMD treatments trials (CATT). Ophthalmology. 2014;121(8):1663–1664.e1. doi:10.1016/j.ophtha.2014.02.025 [CrossRef]
- Hagstrom SA, Ying GS, Maguire MG, et al. VEGFR2 gene polymorphisms and response to anti-vascular endothelial growth factor therapy in age-related macular degeneration. Ophthalmology. 2015;122(8):1563–1568. doi:10.1016/j.ophtha.2015.04.024 [CrossRef]
- Gemenetzi M, Lotery AJ. The role of epigenetics in age-related macular degeneration. Eye (Lond). 2014;28(12):1407–1417. doi:10.1038/eye.2014.225 [CrossRef]
- Oliver VF, Jaffe AE, Song J, et al. Differential DNA methylation identified in the blood and retina of AMD patients. Epigenetics. 2015;10(8):698–707. doi:10.1080/15592294.2015.1060388 [CrossRef]
- Yoshida M, Horinouchi S, Beppu T. Trichostatin A and trapoxin: novel chemical probes for the role of histone acetylation in chromatin structure and function. Bioessays. 1995;17(5):423–430. doi:10.1002/bies.950170510 [CrossRef]
- Suuronen T, Nuutinen T, Ryhänen T, Kaarniranta K, Salminen A. Epigenetic regulation of clusterin/apolipoprotein J expression in retinal pigment epithelial cells. Biochem Biophys Res Commun. 2007;357(2):397–401. doi:10.1016/j.bbrc.2007.03.135 [CrossRef]
- Crosson CE, Mani SK, Husain S, Alsarraf O, Menick DR. Inhibition of histone deacetylase protects the retina from ischemic injury. Invest Ophthalmol Vis Sci. 2010;51(7):3639–3645. doi:10.1167/iovs.09-4538 [CrossRef]
Fold Expression of Histone Acetylation Genes in Untreated and Anti-VEGF treated ARPE-19 culturesa
|Gene||Fold||P Value||Fold||P Value||Fold||P Value||Fold||P value||Fold||P Value||Fold||P Value|
|HDAC1||0.74||< .001||0.9||.25||0.73||< .001||0.97||.73||0.91||.54||1.34||.07|