Abstract
With 5-year survival rates below 5%, small cell lung carcinoma (SCLC) has very poor prognosis and requires improved therapies. Despite an excellent overall response to first-line therapy, relapses are frequent and further treatments are disappointing. The goal of the study was to improve second-line therapy of SCLC.
The effect of chemotherapeutic agents was evaluated in cell lines (apoptosis, reactive oxygen species, and RNA and protein expression) and in mouse models (tumour development).
We demonstrate here that valproic acid, a histone deacetylase inhibitor, improves the efficacy of a second-line regimen (vindesine, doxorubicin and cyclophosphamide) in SCLC cells and in mouse models.
Transcriptomic profiling integrating microRNA and mRNA data identifies key signalling pathways in the response of SCLC cells to valproic acid, opening new prospects for improved therapies.
Abstract
Valproic acid improves second-line regimen of SCLC response in preclinical models http://ow.ly/Rsyd8
Introduction
Lung cancer is the leading cause of cancer-related death worldwide. The outcome of small cell lung carcinoma (SCLC) patients is the poorest of any histological subtype, with 5-year survival rates of <25% and <5% for limited- and extensive-stage disease, respectively [1]. Despite overall first-line response rates ranging between 60% and 80% (extensive), and 80% and 90% (limited), most tumours relapse. The prognosis remains very poor, with median survival rates of only 8–13 months (extensive) and 14–20 months (limited) [2]. Although significant efforts to develop new therapeutic strategies have been made during the last decade, results are still disappointing [2–5]. Future improvements in outcomes will require clarification of the molecular basis of this disease [1].
Epigenetic errors contribute to the initiation, progression and response to therapy of cancer (reviewed by Barnes et al. [6] and Petta et al. [7]). We and others previously proposed a working hypothesis postulating that histone deacetylase (HDAC) inhibitors induce antitumor activity by reversing epigenetic errors [8–11]. In particular, valproic acid (VPA) is an inhibitor of HDACs displaying appropriate pharmacokinetic properties, and yielding only moderate toxicity that is acceptable in the context of an anticancer treatment [12–14]. By modulating a broad range of activities, including proliferation, apoptosis and differentiation, VPA has antitumoural properties in several cancers, including SCLC [15–21].
Although there is no standard second-line therapy for SCLC, possible treatments most often comprise a combination of three chemotherapeutic agents: a DNA crosslinking agent (e.g. cyclophosphamide), an inhibitor of topoisomerase II (e.g. doxorubicin) and a mitotic spindle poison (e.g. vindesine) (here referred to collectively as “VAC”). With the aim of improving the treatment of extensive SCLC, we evaluated the capacity of VPA to increase the anticancer effect of the VAC regimen in cell cultures and in xenograft mouse models. The mechanisms involved in chemotherapeutic response to VPA were then studied by transcriptomic analyses.
Materials and methods
Cell culture conditions
Human SCLC cell lines (H146, H526 and H69) were purchased from the ATCC (Manassas, VA, USA) and cultivated as detailed previously [19]. Cells were incubated with VPA (Sigma-Aldrich, Diegem, Belgium), mafosfamide (Baxter, Braine-l'Alleud, Belgium), cyclophosphamide (Baxter), doxorubicin (Pfizer, Elsene, Belgium) and vindesine (Lilly, Brussels, Belgium), alone or in combination. Since cyclophosphamide needs to be activated in vivo by the hepatic metabolism, its active form, mafosfamide, was used for in vitro experiments. Optimal drug concentrations were determined by MTS viability assays.
Detection of apoptosis
Apoptosis was quantified by flow cytometry after ethanol fixation and propidium iodide incorporation, as outlined previously [22]. A synergy index was calculated using the formula: The percentage of specific apoptosis was determined using the formula: When the synergy index was >1, 1 or <1, the effects were defined as synergistic, additive or antagonistic, respectively.
To assess the role of caspases in apoptotic pathways, 5×105 cells were incubated with or without: 20 µM Z-Val-Ala-Asp(OMe)-CH2F (Becton Dickinson, Erembodegem, Belgium), a total pan-caspase inhibitor; 20 µM negative control (Z-FA-fmk) (Becton Dickinson); 40 µM Z-Ile-Glu(OMe)-Thr-Asp(OMe)-CH2F (Calbiochem, Overijse, Belgium), a caspase-8 specific inhibitor; or 40 µM Z-Leu-Glu(OMe)-His-Asp(OMe)-CH2F (Calbiochem), caspase 9 specific inhibitor; all compounds being diluted in dimethylsulfoxide.
Quantification of reactive oxygen species
Reactive oxygen species (ROS) were detected using 5,6-chloromethyl-2′,7′-dichlorodihydrofluorescein diacetate acetyl ester (CM-H2DCFDA; InVitrogen, Ghent, Belgium). After 30 min of pre-incubation with 5 µM CM-H2DCFDA, the different drugs were added alone or in combination. After 24 h of culture, SCLC cell lines (5×105 cells per mL in 24-well plates) were harvested, washed with PBS and analysed by flow cytometry (FACS Aria; Becton Dickinson). ROS production was quantified using the fluorescence intensity of chloromethyldichlorofluorescein. 10 000 events were collected and analysed with the FACS Diva software (Becton Dickinson). Cells were also treated with 100 µM hydrogen peroxide or 10 mM N-acetylcysteine (Calbiochem), a free-radical scavenger, as positive and negative controls, respectively.
Immunoblot analysis
Protein expression levels and intracellular translocations were assessed using cytosolic and nucleic buffers, and standard protocols for western blotting as detailed previously [19]. The antibodies used were: anti-acetylated histone H3 (Upstate, Overijse, Belgium), anti-actin, anti-Erk (both Sigma-Aldrich), anti-Bax, anti-Bcl-2 (both Dako Cytomation, Heverlee, Belgium), anti-Bid, anti-cytochrome c (Becton Dickinson), anti-BclxL, anti-phospho-Erk, anti-caspase 8, anti-caspase 9 (Cell Signaling, Leiden, the Netherlands), anti-γH2AX and anti-VDAC1 (Abcam, Cambridge, UK).
Evaluation of regimen efficacy in severe combine immunodeficiency mice
The Institutional Animal Care and Usage Committee of the University of Pennsylvania (Philadelphia, PA, USA) and the University of Liege (Liege, Belgium) approved all animal protocols in compliance with the Guide for the Care and Use of Laboratory Animals, according to the Declaration of Helsinki. The severe combined immunodeficiency (SCID) mice (BALB/c HanHsd-Prkdc; Jackson Laboratories, Sacramento, CA, USA) or NOD/SCID mice received a standard research diet throughout the experiment. H146 and H69 cells (2×106), embedded in 50% Matrigel Basement Membrane Matrix High Concentration (BD Biosciences, Erembodegem, Belgium), were implanted subcutaneously into the flanks of 7-week-old female SCID mice. When tumours reached a volume of 300–400 mm³, mice were administered with daily intraperitoneal injections of VPA (400 mg·kg−1·day−1) or PBS as a control. 3 days after the first VPA administration, intraperitoneal injections of cyclophosphamide (40 mg·kg−1), vindesine (0.5 mg·kg−1) and doxorubicin (3 mg·kg−1) were performed (at days 14 and 17, and days 8, 13, 25 and 34 for H69 and H146, respectively). The schedule of drugs injections was first determined in preliminary dose response experiments. Tumour volumes were calculated twice a week using the formula: where D is the largest diameter and d the smallest. Groups of at least six mice were tested under each experimental condition.
RNA extraction
Cell lines were cultivated for 4 h in presence of 1 mM VPA. RNA was then extracted after two washes in PBS and lysis Trizol RNA Isolation Solution (Ambion, Ghent, Belgium), according to the manufacturer's protocol. For microarray analysis, a second separation between the aqueous and organic phases was performed on Phase Lock Gel (5 Prime, Leuven, Belgium) to optimise RNA recovery and purity. All the extracted RNAs were assessed for quantity and purity using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Temse, Belgium) and stored at −80°C until the amplification step.
Microarray analysis
The details of the microarray analysis have been described previously [19]. Briefly, reference RNAs were obtained by pooling equal amounts of total RNAs from three different SCLC cell lines (H146, H526 and H69). RNAs were reverse-transcribed with Moloney murine leukaemia virus reverse transcriptase using a polydeoxythymidine primer and amplified by PCR. The cDNAs were transcribed into copy RNAs (cRNAs) transcribed and labelled with Cy3 (reference) or Cy5 (sample) dyes using the Low Input RNA Linear Amplification Plus kit (Agilent Technologies, Diegem, Belgium). After quality checking, the labelled cRNAs and RNAs spiked in were hybridised on an Agilent oligonucleotide microarray (two-colour Whole Human Genome 4×44K arrays) and analysed according to manufacturer's protocol (GE2_v5_95_Feb07; Feature Extraction software, version 9.5.3.1). The statistical analyses were performed with Genespring GX, version 7.3.1 (Agilent Technologies). Additional normalisation steps were performed (per spot: division by the control signal; per array: normalisation to the 50th centile; per gene: normalisation to the median). Transcripts that were not present in at least one sample were excluded from additional analyses. Welch's t-test was used to assess the statistical significance, excluding genes whose expression varied by a factor inferior to 2 across the sample set of interest. The threshold for statistical significance after Benjamini's correction for multiple testing was 0.05.
Real-time reverse transcription PCR
mRNA levels were analysed by real-time quantitative reverse transcription PCR (qRT-PCR) as described previously [19], using the primers in table 1.
MicroRNA analysis
MicroRNA (miRNA) expression was determined in triplicates using the TaqMan Low Density Array method (DNAVision, Gosselies, Belgium). Fluorescence data (Bioanalyzer 2100 Nano chip; Agilent Technologies) were normalised by means of two endogenous controls (RNU44 and RNU48). Results were expressed as fold change (2−(mean(ΔCT(VAC+VPA))–mean(ΔCT(VAC)), where CT is the cycle threshold). The Benjamini and Hochberg method was used for multiple testing corrections.
Statistical analysis
All cell culture experiments were performed at least three times and data are presented as mean±sd. Statistical significance was calculated using Student's t-test and data were considered statistically significant, very statistically significant or highly statistically significant when p<0.05, p<0.01 and p<0.001, respectively.
Results
VPA increases apoptosis in SCLC cells treated with VAC in vitro and in vivo
To evaluate the synergism between VPA and the second-line regimen VAC, VPA, at a concentration achievable in patients (1 mM) [23], was combined with therapeutically relevant doses of mafosfamide (10 µM), doxorubicin (0.3 µM) and vindesine (20 nM). After 24 h of culture, the combination of VPA and VAC significantly increased apoptotic rates in three SCLC cell lines (H146, H69 and H526; p<0.01 by Student's t-test) (figure 1). In fact, the effect of VPA on the combined VAC treatment was synergistic (synergy index of 2.9, 2.4 and 1.6 in H146, H69 and H526 cells, respectively).
We next determined the involvement of specific caspases using pharmacological inhibitors of caspases 3, 8 and 9 (figure 2). To facilitate comparison between treatments, control results (Z-FA-fmk) were arbitrarily normalised to 1. All three inhibitors significantly reduced apoptosis of SCLC cells, indicating that the mechanism is caspase-dependent, and involves both extrinsic and intrinsic pathways.
To evaluate the role of reactive oxygen species in the apoptotic process, intracellular ROS levels were monitored by flow cytometry in H146, H69 and H526 cells cultivated for 24 h in presence of VPA and/or VAC. VPA, but not VAC, significantly increased ROS production in all three cell lines (p<0.05) (figure 3a, c and e). Despite significantly decreasing ROS levels, the free-radical scavenger N-acetylcysteine was not sufficient to inhibit apoptosis (figure 3b, d and f).
We next evaluated expression of a series of proteins involved in the cell cycle and apoptosis (figure 4). Western blot analysis demonstrated that VPA-induced apoptosis involved hyperacetylation of histone H3 and phosphorylation of Erk1/2. Double-stranded DNA breaks in the presence of VPA were revealed by γH2AX blotting. VPA created an imbalance between pro- and anti-apoptotic modulators: a decrease in cytoplasmic Bcl-2 and BclxL; cleavage of Bid (into t-Bid); processing of caspase 8 and caspase 9; the appearance of a cleaved form of Bax in the mitochondria; and release of cytochrome c into the cytoplasm.
Based on the proapoptotic synergy between VPA and VAC in SCLC cells, we investigated the ability of VPA to improve antitumor efficacy in mouse models. Therefore, SCLC cells (H69 and H146) were injected subcutaneously into the flanks of SCID mice, and the mice were treated with VPA, cyclophophamide (40 mg·kg−1), doxorubicin (3 mg·kg−1) and vindesine (0.5 mg·kg−1), as described in the materials and methods section. Although partial responses were observed with VPA or VAC, tumours eventually relapsed (figure 5). In contrast, tumour growth was restricted in both SCID mouse xenograft models when VPA was combined with VAC (p<0.001 and p<0.05 by Student's t-test, for H69 and H146, respectively).
We conclude that, in two preclinical models of SCLC, VPA improves the efficacy of VAC.
VPA modulates key cellular pathways in SCLC, including cell death and tumour invasion
To better characterise the molecular mechanisms involved, we analysed the transcriptome of H526 cells treated for 4 h with VPA and/or VAC using Agilent microarrays. Bioinformatic analyses revealed a list of genes that were significantly up- or down-regulated by a factor >2 in presence of VPA. Among these, 138 were also specifically identified in cells treated with VPA+VAC compared to VAC (online supplementary material). Expression of the most significantly modulated genes was confirmed by qRT-PCR in three SCLC cell lines (table 2).
With 54 genes involved, cell death was the top scored biological function affected, as revealed by Ingenuity (Qiagen, Venlo, the Netherlands) (figure 6). The network related to this function included several genes involved in apoptosis: BBC3 (Puma, a pro-apoptotic member of the Bcl2 family), TNFRSF19 and TNFRSF12A (upstream effectors of apoptosis). To extend this transcriptomic analysis, we next analysed expression of miRNAs [24, 25]. After correction for multiple testing, only two miRNAs were significantly modulated when adding VPA, both of them being involved in tumour invasion: miR-589 (fold increase of 441 with VPA, corrected p-value <0.01) and miR-575 (fold decrease of 0.018, corrected p-value <0.05) (table 3).
Discussion
The preclinical evidence provided in this report indicates that VPA synergistically increases the apoptotic rates of three different cell lines in response to a second-line therapy for SCLC (figure 1). The mechanism is caspase-dependent, and involves both extrinsic and intrinsic pathways (figure 2). VPA increases the level of ROS, and modifies the balance between pro- and antiapoptotic modulators (figures 3 and 4). The new regimen of VPA+VAC efficiently reduces tumour growth in SCID mouse models engrafted with human SCLC cells (figure 5). This evidence spurred the European Lung Cancer Work Party to launch a phase I/II clinical trial to assess the combination of VPA with VAC in patients presenting with relapsing or refractory SCLC (registered at www.elcwp.org with identifier number 01081). Results presented in the accompanying article by Berghmans et al. [26] indicate that the new regimen induces a significant response in relapsing/refractory SCLC patients. However, since the treatment does not translate into a significant progression-free survival, it is not recommended for second-line therapy.
Our transcriptomic data provide some evidence to better understand the mechanisms involved and to give insights for novel regimen. VPA modulates key genes modulating pathways linked to cell cycle and apoptosis. Among these, Puma (BBC3, Bcl-2 binding component 3), TNFRSF12A (TWEAK/Fn14), sodium/potassium ATPase and HEY1 (a mediator of Notch signalling) are significantly upregulated by VPA. In particular, the sodium/potassium ATPase is involved in sensitivity to platinium [27] and could be a major factor involved in VPA efficacy. Puma is able to interact with Bax and BclxL, promoting translocation of Bax to the mitochondria and competitive binding to BclxL, impeding the latter's ability to inhibit Bax proapoptotic activity [28]. Parallel analyses of miRNA expression identified two particular miRNAs whose expression is tightly correlated with VPA treatment. When the VAC+VPA regimen is compared to VAC alone, miR-589 is increased 441-fold while miR-575 expression is reduced to 0.018. miRNA-589 modulates epithelial–mesenchymal transition in human mesothelial cells [29] while miR-575 promotes growth and invasion of lung cancer cells [30]. Causal involvement of these microRNAs in treatment efficacy could be further investigated by antagomiRs and miRNA mimics.
VPA also regulates expression of modulators of DNA damage response such as Bcl-6, which inhibits TP53 and modulates DNA damage-induced apoptotic responses [31]. Most important is the concomitant modulation of several genes involved in the Wnt pathway, including FZD2, FZD5, FZD7, FZD10 and WNT11. In this context, our preliminary data show that VPA activates TCF/LEF transcription factors in SCLC cells (data not shown). Furthermore, two general Wnt inhibitors (quercetin and ethacrynic acid) impair VPA-induced apoptosis. At first glance, involvement of Wnt in the therapeutic response may appear surprising, since this pathway has been directly correlated with tumorigenicity. For example, overexpression of FZD7 was associated with aberrant activation of the Wnt canonical pathway in oesophagus and colorectal cancers [32, 33]. However, the complexity of the Wnt pathways is clearly less well-understood or at least heterogenous in lung cancer [34]. In fact, expression of WNT7A, which targets FZD7, has been shown to be downregulated in lung cancer [35]. However, SCLC is characterised by overexpression of Wnt inhibitors such as NLK, Sox11 and TCF-4 [36]. In contrast, a tumour suppressor role of the Wnt pathway in SCLC has also been reported [37]. In this perspective, our results indicate that Wnt signalling is associated with a chemotherapeutic response. Similarly, there is a relationship between hyperinduced canonical Wnt activity and enhanced apoptosis in HDAC inhibitor-treated colorectal tumour cells [38].
In conclusion, this report provides preclinical evidence for the use of VPA in second-line therapy of SCLC and opens new prospects for improvement using Wnt activators [39, 40].
Acknowledgements
We are grateful to Jean-Paul Sculier and Thierry Berghmans (Department of Intensive Care Unit and Thoracic Oncology, Institut Jules Bordet, Centre des Tumeurs de l'Université Libre de Bruxelles, Brussels, Belgium) for their collaboration in setting up the clinical trial reported in their companion article to this one [26]. We thank the GIGA technology platforms.
Footnotes
This article has supplementary material available from openres.ersjournals.com
Support statement: The work was supported by the Fonds National de la Recherche Scientifique (FNRS), the Télévie, the Interuniversity Attraction Poles programme BELVIR initiated by the Belgian Science Policy Office, the Action de Recherche Concertée Glyvir of the Communauté française de Belgique, the Belgian Foundation against Cancer, the Sixth Research Framework Programme of the European Union (project INCA LSHC-CT-2005-018704), the Neoangio excellence programme and the Partenariat Public Privé INCA of the Direction générale des Technologies, de la Recherche et de l'Energie/DG06 of the Walloon government, the Centre anticancéreux près ULg, the Synbiofor and Agricultureislife projects of GxABT, the ULg Fonds Spéciaux pour la Recherche, the Plan Cancer of the Service Public Fédéral. L. Willems (research director) is member of the FNRS. Funding information for this article has been deposited with FundRef.
Conflict of interest: Disclosures can be found alongside this article at openres.ersjournals.com
- Received June 2, 2015.
- Accepted August 24, 2015.
- The content of this work is ©the authors or their employers. Design and branding are ©ERS 2015
This article is open access and distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0.