- Research
- Open access
- Published:
PARP1 expression drives the synergistic antitumor activity of trabectedin and PARP1 inhibitors in sarcoma preclinical models
Molecular Cancer volume 16, Article number: 86 (2017)
Abstract
Background
Enhancing the antitumor activity of the DNA-damaging drugs is an attractive strategy to improve current treatment options. Trabectedin is an isoquinoline alkylating agent with a peculiar mechanism of action. It binds to minor groove of DNA inducing single- and double-strand-breaks. These kinds of damage lead to the activation of PARP1, a first-line enzyme in DNA-damage response pathways. We hypothesized that PARP1 targeting could perpetuate trabectedin-induced DNA damage in tumor cells leading finally to cell death.
Methods
We investigated trabectedin and PARP1 inhibitor synergism in several tumor histotypes both in vitro and in vivo (subcutaneous and orthotopic tumor xenografts in mice). We searched for key determinants of drug synergism by comparative genomic hybridization (aCGH) and gene expression profiling (GEP) and validated their functional role.
Results
Trabectedin activated PARP1 enzyme and the combination with PARP1 inhibitors potentiated DNA damage, cell cycle arrest at G2/M checkpoint and apoptosis, if compared to single agents. Olaparib was the most active PARP1 inhibitor to combine with trabectedin and we confirmed the antitumor and antimetastatic activity of trabectedin/olaparib combination in mice models. However, we observed different degree of trabectedin/olaparib synergism among different cell lines. Namely, in DMR leiomyosarcoma models the combination was significantly more active than single agents, while in SJSA-1 osteosarcoma models no further advantage was obtained if compared to trabectedin alone. aCGH and GEP revealed that key components of DNA-repair pathways were involved in trabectedin/olaparib synergism. In particular, PARP1 expression dictated the degree of the synergism. Indeed, trabectedin/olaparib synergism was increased after PARP1 overexpression and reduced after PARP1 silencing.
Conclusions
PARP1 inhibition potentiated trabectedin activity in a PARP1-dependent manner and PARP1 expression in tumor cells might be a useful predictive biomarker that deserves clinical evaluation.
Background
An attractive strategy to improve current treatment options is to inflict cytotoxic DNA damage with chemotherapy, and then impede DNA repair by molecular targeting. Poly-ADP-ribosyl-transferase-1 (PARP1) is a key sensor of DNA damage and initiates recruitment of the DNA-repair machinery to the site of damage [1,2,3]. The development of PARP1 inhibitors has drawn closer the goal of combining these compounds with current therapies [2, 4,5,6,7,8,9]. Unfortunately, despite several preclinical data confirming an increased antitumor activity by combining PARP1 inhibitors with either chemotherapy or radiotherapy [10,11,12,13,14,15,16], dose escalation in phase 1 combination studies has been greatly hampered by the observed hematologic toxicities. These adverse events limited the possibility to exploit these combinations into clinical practice [4, 5, 9, 11, 17,18,19,20,21,22,23,24,25,26,27]. As a consequence, PARP1 inhibitors are today registered as monotherapy in cancers bearing DNA-repair deficiencies [28,29,30,31,32,33,34] and the strong rational to combine selected cytotoxics (especially alkylators) with PARP1 inhibitors has to face the risk of myelotoxicity.
Among chemotherapeutics, trabectedin has some peculiarities that point out this drug as an ideal candidate to be combined with PARP1 inhibitors: its favorable safety profile on the hematologic side and its unique mechanism of action [35, 36]. In particular, trabectedin traps enzymes belonging to the transcription-coupled nucleotide excision repair (TC-NER) system that in the attempt to remove trabectedin adducts, generates DNA single- and double-strand breaks. Understandably, trabectedin displays a greater clinical benefit in BRCA1/2-deficient tumors than in proficient ones [37, 38]. These characteristics prompted us to explore the combination of trabectedin and PARP1 inhibitors taking advantage of a large set of bone and soft tissue sarcoma (BSTS) cell lines. Despite a common mesenchymal origin, BSTS are characterized by a great degree of heterogeneity and, indeed, trabectedin displays a spectrum of activity in various sarcoma histotypes [39,40,41]. This heterogeneity allowed us to test cells with different intrinsic sensitivity to trabectedin, in order to explore if trabectedin could efficiently activate PARP1, and if PARP1-specific inhibition could be exploited sufficiently to cause irreversible DNA damage, and eventually tumor cell death. The observed results were subsequently validated in other tumor types and by combining PARP1 inhibitors with other cytotoxics characterized by different mechanisms of action.
Methods
Cell line characterization, cell viability and western blot
Cell lines characteristics are depicted in Additional file 1: Table S1 and S2. Short tandem repeat (STR) profile was checked and the genomic status of DNA-repair key components (BRCA1, BRCA2, ATM, CHEK2, PTEN) was analyzed by Multiplex ligation-dependent probe amplification (MLPA) and Denaturing high performance liquid chromatography (DHPLC, Wave 3500HT DNA Fragment Analysis System, Transgenomic Inc.) followed by direct sequencing (ABI PRISM3100 DNA Sequencer, Applied Biosystem, Forster City, CA). PARP1, RAD51, and BRCA1 copy number variations was confirmed by real-time polymerase chain reaction (PCR) on genomic DNA (TaqMan Assay, ABI PRISM 7900HT System, Applied Biosystem).
In order to evaluate the proliferation rate of each cell line, cells were seeded at a density of 2 × 104 cells/cm2, and the doubling-time (DT) of the harvested cells was calculated during the exponential growth phase (48 h of culture) by using the algorithm provided by http://www.doubling-time.com/compute.php: DT = t x lg2/(lgNt - lgN0) where N0 is the initial concentration of cells, Nt is final concentration of cells and t is the culture time in hours. Cell viability was determined with Cell Titer-Glo (Promega) after 72-h treatment with scalar doses (2–0.125 nM) of trabectedin (PharmaMar), as single agent or in constant combination with olaparib (20–1.25 μM) or veliparib (80–5 μM) (Sequoia research products). These ranges of concentrations were chosen on the results of previous studies testing the sensitivity to each single agent in sarcoma and non-sarcoma cell lines [42,43,44,45]. Protein extracts were obtained after 24-h treatment and resolved by western blot using primary antibodies from Cell Signaling Technology except for anti- poly(ADP-ribose) (PAR, Trevigen), phospho-histone H2AX (Ser139) (Millipore).
Quantitative real-time polymerase chain reaction (PCR), gene expression profiling and Gene Signature enrichment analysis (GSEA)
2.5 × 106 cells were plated in 150 mm diameter and grown for 24 h in complete medium. Cells were then treated for additional 24 h with trabectedin (0.125 nM) and olaparib (1.25 μM) as single agent and in combination. In each condition, 500.000 cells were fixed in 70% ethanol for cell cycle analysis and the remaining fraction was lysed with Qiazol reagent (Qiagen) for RNA extraction by means of RNeasy Kit (Qiagen) following manufacturer’s instructions. Gene expression profiling was performed with Human HT-12 v4.0 Expression BeadChip Kit (Illumina) Real-time PCR was performed by TaqMan Gene Expression Assay using ABI PRISM 7900HT System (Applied Biosystem). Fluorescence data were automatically converted into CT (cycle threshold) values. To export data, the threshold was 0.20. Raw data were analyzed by Microsoft Office Excel. Additional file 2: Figure S1 outlined the design of the expression profile project analyzed by GSEA (http://www.broadinstitute.org/gsea/msigdb/index.jsp) [46]. To highlight gene patterns associated and involved with trabectedin and olaparib synergism, we designed a specific expression profile project comparing cell lines displaying high synergism (HS-C: TC-106, 402.91, DMR) and cell lines displaying low/no synergism (LS-C: SJSA-1, HT1080, SW684). First, we compared gene expression profiles of HS-C and LS-C by GSEA to identify gene sets associated with high synergism. Then we selected the cores of the emerging gene sets matching the criterion of false discovery rate (FDR) < 0.05 and we challenged them with GSEA comparing each treated HS-C against untreated. GSEA analysis was performed with default parameters. Probes were collapsed into gene symbols. Gene set size for inclusion was set between 15 and 500, and the genes were permutated 1,000 times. Gene sets that met the FDR < 0.05 were considered significant. The catalog of gene sets was downloaded from Molecular Signature Database (C2, C5, C6, Hallmarks, MSigDB, version 5.0, http://www.broadinstitute.org/gsea/msigdb/index.jsp) for a total of 4965 curated gene sets.
Array comparative genomic hybridization (aCGH) analysis
Comparative genomic hybridization using aCGH microarrays was carried out using the enzymatic labeling method. Digestion, labeling, hybridization, washing and slide scanning were performed following the manufacturer protocols (Agilent Technologies). Briefly, 750 ng of samples and control DNA in a total volume of 10.1 μL were digested with restriction enzymes and labeled with Cy3 and Cy5. Subsequently, labeled DNA was cleaned up and hybridized using the Agilent-030587 CCMC CGH plus SNP 180 k Microarray platform. Each test DNA sample (Cy5) was hybridized with a reference DNA (Cy3) from Homo sapiens (Agilent). After hybridization for 24 h at 20 rpm and 65 °C, slides were washed following Agilent procedure and scanned with the dual-laser microarray scanner version C (G2505C, Agilent Technologies). Images were analyzed using Feature Extraction software version 10.7 (Agilent Technologies). Raw data were processed using the Agilent Genomic Workbench version 7. Aberrant regions were detected using ADM-2 algorithm with threshold = 6. To avoid false positive calls, the minimum number of consecutive probes for amplifications/deletions was set to 3, together with a minimum average absolute Log Ratio for aberrations > = 0.25.
Gene ontology (GO) enrichment analysis
Genes mapped in the regions with differential aberrations between HS-C and LS-C were uploaded into DAVID bioinformatics tool for functional enrichment analysis and GO biological processes (level 5). Enrichment p-values less than 0.01 were considered as statistically significant.
Mice xenograft models
Nonobese diabetic/severe combined immunodeficient (NOD/SCID) mice (Charles River) were injected with: a) 1–2.5 × 106 SJSA-1 subcutaneously (s.c.) into the right flank or 105 SJSA-1 intravenously (i.v.) into the tail vein to originate subcutaneous primary tumors and lung metastases, respectively, or b) 106 DMR s.c. or 4 × 105 DMR orthotopically into the uterine wall by laparotomy. Six mice per group (DMR s.c. model) were randomized to receive: a) 25 or b) 50 mg/kg/day intraperitoneal (i.p) injection of olaparib (5 days/week); c) 0.150 or d) 0.1 or e) 0.05 mg/kg weekly i.v. injection of trabectedin; f) combination of 25 mg/kg/day olaparib and 0.050 mg/kg trabectedin or g) left untreated for 21 days and then sacrificed for histological and molecular assays. SJSA-1 s.c. models were treated with a) 0.05 mg/kg or b) 0.025 mg/kg weekly i.v. injection of trabectedin; c) 25 mg/kg olaparib i.p; d) combination of 25 mg/kg/day olaparib and 0.050 mg/kg trabectedin; e) combination of 25 mg/kg/day olaparib and 0.025 mg/kg trabectedin or g) left untreated for 17 days. DMR orthotopic model was treated with 0.050 mg/kg trabectedin and 25 mg/kg/day olaparib as single agent, in combination, or left untreated. These in vivo experiments were conducted in accordance with the protocol approved by the Institutional Ethics Review Board (IRB) and by the Italian Ministry of Health (Aut. Min. 178/201S-PR). For in vivo imaging, each mouse received an intra-peritoneal dose of 150 mg luciferin/kg body weight, and ventral and dorsal images were obtained by IVIS Lumina II and quantified using Living Image software (PerkinElmer).
Flow cytometry assays: cell cycle, P-H2AX, and apoptosis analysis
The effects of trabectedin and olaparib on the cell cycle and the expression of the DNA-damage marker P-H2AX were determined by staining DNA content with propidium iodide (PI, Sigma Aldrich) and FITC-conjugated primary antibody anti-P-H2AX (Ser139) (Millipore), respectively, after treatment with trabectedin (0.125 nM) as single agent or in combination with olaparib (1.25 μM) for 24 h.
For apoptosis evaluation, after 96-h treatment APC-labelled Annexin V (eBioscence) and PI (0.5 μg/ml) staining was done. All samples were acquired by Cyan ADP Flow cytometer (Beckman Coulter) and analyzed by the Summit v4.3 (Dako) and FlowJo (Tree Star) softwares.
Comet assay (Single Cell Gel Electrophoresis assay)
For determination of DNA single- and double-strand breaks, a single-cell gel electrophoresis assay was used (Comet Assay TM, Trevigen) per the manufacturer’s instructions after 48-h treatment with trabectedin (0.125 nM) and olaparib (1.25 μM), as single agents or in combination. Quantization of the DNA in the tails of the comets was performed with Quantity One software (Bio-Rad Laboratories) and ImageJ software version 1.49 (http://rsbweb.nih.gov/ij/index.html).
Lentiviral vector production
pCCL.sin.cPPT.polyA.CTE.eGFP.minhCMV.hPGK.Luciferase.Wprep or lentiviral plasmid vector PLKO-puro (Sigma-Aldrich), the Precision LentiORF PARP1 (Id:PLOHS_100004266, Thermo Scientific) and packaging vectors pMDLg/pRRE pRSV-REV and pMD2.VSVG were used for lentiviral preparation. Efficiency of transduction was confirmed by western blot analysis and flow cytometry analysis by Cyan ADP Flow cytometer (Beckman Coulter). The luciferase activity was tested by in vitro bioluminescent assay (Caliper Life Sciences, Inc.) and measured by IVIS Lumina II (PerkinElmer).
Archival tumor samples, immunostaining, and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assay
Archival formalin-fixed, paraffin-embedded (FFPE) tumor samples were collected in accordance with an IRB approved protocol (202/2014). Four μm slices were cut and stained with primary antibodies against PARP1, BRCA1 and RAD51 (Abcam) following standard immunohistochemistry procedure. The same antibodies and one antibody against PAR (Calbiochem) were used for immunocytochemistry on cells grown as monolayer in chamber slides. TUNEL assays were used to evaluate the number of apoptotic cells in tumors explanted from mice, using the ApopTag kit (Millipore). Immunostaining on mice xenografts was performed according to standard protocols with primary antibodies from Sigma Aldrich (proliferating cell nuclear antigen, PCNA), Abcam (p-H2AX Serine 139). Nuclei were counterstained with hematoxylin. Immunofluorescence was done with FITC-conjugated antibody against P-H2AX (Millipore). Nuclei were counterstained with DAPI. Fluorescence in situ hybridization (FISH) analysis of PARP1 and centromeres was performed following standard procedures.
Visible images were acquired with a Leica DM1000 microscope equipped with a color 3.1 M Pixel CMOS digital camera. Fluorescent images were acquired with a Leica confocal laser-scanning microscope (TCS SP5 AOBS). For signal quantification, 5 images/sample were acquired by maintaining constant parameters. Image quantification was performed using ImageJ software and Leica Confocal Software.
Statistical and pharmacological combination analyses
All in vitro experiments were performed at least three times. Differences between treatment groups were analyzed by the two-tailed Student’s t test and the two-way ANOVA with post hoc Bonferroni’s correction for multiple tests using GraphPad Prism 5 (GraphPad Software Inc.). The concentration inhibiting 50% of the cell growth (IC50) with its 95% confidence intervals (95% CI), and drug synergism, expressed as combination index (CI) calculated at IC50 with its estimated standard deviation were obtained by using CalcuSyn software (Biosoft). The correlation analyses between ΔCT values (mRNA levels) or protein expression levels and combination indexes were performed by calculating Pearson correlation coefficient, t-distribution, and P values by Microsoft Excel.
Results
Trabectedin and olaparib synergism is related to PARP1 expression
Initially, we demonstrated that in bone and soft tissue sarcoma (BSTS) cell lines trabectedin treatment significantly increased phosphorylated histone H2AX (P-H2AX), the marker of DNA double-strand breaks (p < 0.001, Fig. 1a-c). Subsequently, we showed that after 24-h treatment with the combination of trabectedin (0.125 nM) and olaparib (1.25 μM) there was a variable but significant increase in P-H2AX positive cells if compared to single agents and untreated controls (Fig. 1a-c). Thereafter, we studied the role of PARP1 enzymatic activity (PARylation) after 4-h treatment with trabectedin (10 nM). We observed that in 3/6 (50%) bone and soft tissue sarcoma (BSTS) cell lines PARylation was highly increased, as compared to untreated controls (Fig. 1d). In these cells, the PARP1 inhibitor olaparib inhibited both basal and trabectedin-induced PARP1 activation remarkably well (Fig. 1d). To elucidate why trabectedin failed to activate PARP1 equally in all cell lines, we studied the correlation between PARP1 basal expression, activity and gene status. We found that after trabectedin exposure, PARylation was significantly increased in high- vs. low-PARP1-expressing cells (Fig. 1d). To further explain differences in PARylation activity, we performed mutational analysis of PARP1 gene and we found the Val762Ala single nucleotide polymorphism in SJSA-1 and SW684 cells (Additional file 2: Figure S2).
These results caused us to explore the effect of PARP1 inhibition on trabectedin-induced DNA damage and we found that 24-h treatment with the combination of trabectedin (0.125 nM) and olaparib (1.25 μM) significantly increased P-H2AX signaling in all cell lines tested (p < 0.001, Fig. 1a-c) if compared to both single agents and controls. Irreparable DNA fragmentation was significantly increased after 48-h treatment with the combination if compared to both single agents and controls in high-PARP1-expressing cells, as shown by COMET assay (p < 0.001; Fig. 1e-f). On the contrary, in low-PARP1-expressing cells, P-H2AX activation led to DNA repair. The intensity of the combination-induced DNA damage caused a cell line-dependent cell cycle arrest at G2/M checkpoint and the accumulation of cells into the late S phase (Fig. 2a), resulting in apoptosis after 96 h (Fig. 2b). This effect was consistently stronger in high-PARP1-expressing cells (TC-106, DMR, 402.91) and led to a significant inhibition of tumor cell colony growth (p < 0.001; Fig. 2c-d).
To expand our findings, we tested the antitumor activity of serial dilutions of two PARP1 inhibitors, olaparib (20–1.25 μM) and veliparib (80–5 μM) in a panel of 20 different BSTS cell lines, both as single agents and in constant combination with serial dilution of trabectedin (2–0.125 nM). We observed a strong synergism of trabectedin + olaparib in 18/20 cell lines while in the trabectedin + veliparib combination, we observed less synergism in 15/19 cell lines (Fig. 2e, Additional file 1: Table S1). In high-PARP1-expressing cells the trabectedin IC50s were significantly reduced if compared to low-PARP1-expressing cells both in combination with olaparib and veliparib (p = 0.01 and p = 0.02, respectively; Additional file 2: Figure S3).
Based on these in vitro experiments, we chose olaparib and trabectedin for further drug combination studies. As an example, Additional file 2: Figure S4 shows dose-effect curves obtained in TC-106, 402.91, DMR; SJSA-1, HT1080, and SW694 cells. Synergy was distributed along a spectrum of intensity; it was strongest in Ewing’s sarcoma cells and decreased progressively from liposarcoma (myxoid > dedifferentiated) to osteosarcoma to undifferentiated pleomorphic sarcoma to fibrosarcoma (Fig. 2e, Additional file 1: Table S1). Of note is that the doxorubicin-resistant leiomyosarcoma cell line MES-SA-DX5 also displayed cross-resistance to trabectedin, but addition of PARP1 inhibitors restored sensitivity to a clinically achievable trabectedin concentration (Additional file 1: Table S1 and Additional file 2: Figure S2e). Interestingly, MES-SA-DX5 cells displayed higher basal PARP1 expression and activity (Additional file 2: Figure S5 if compared with their parental counterparts, due to PARP1 gene amplification as shown by fluorescent in situ hybridization (FISH, Additional file 2: Figure S5 and Additional file 1: Table S4), real- time PCR (Additional file 1: Table S5), and by array comparative genomic hybridization (aCGH) analysis (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE77175).
Cell cycle control and DNA-repair pathways are involved in the activity of trabectedin and olaparib combination
To dissect the molecular mechanism behind trabectedin + olaparib synergism, we compared the gene expression profiles (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE76981) of three cell lines with high synergism (HS-C) and three cell lines with low/no synergism (LS-C). First, we compared gene expression profiles of HS-C and LS-C by gene set enrichment analysis (GSEA) to identify gene sets associated with the sensitive phenotype and we found 437 gene sets matching the criteria of false discovery rate (FDR) <0.05 and we challenged them with GSEA, comparing each treated cells against untreated. Of these ones, 57 gene sets were specifically enriched in combination treatment in HS-C, but not in single treatments. Of note, the expression of specific gene sets involved in DNA-damage response, DNA-repair pathways, and cell cycle control were higher in HS-C than in LS-C (Fig. 3a). Notably, these gene sets were specifically down-modulated by the drug combination-not by single agents- in HS-C, but were either unmodified or upregulated in LS-C (Fig. 3a). Taken together, these results indicated that specific gene sets involved in the DNA-damage response, G2/M cell cycle checkpoint, and DNA-repair pathways were implicated in the mechanisms behind trabectedin/olaparib synergism.
Consistent evidences were shown also by comparing genomic profiles of HS-C and LS-C as obtained by aCGH (Additional file 2: Figure S6, Additional file 3 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE77175). Differential aberration analysis demonstrated four groups of genomic alterations that clearly distinguished HS-C from LS-C (Additional file 3). Namely, regions of chromosomes enriched in genes involved in DNA-damage response/repair and G2/M cell cycle checkpoint were deleted and amplified in HS- and LS-cells, respectively (as highlighted in Additional file 3).
To validate these findings, we confirmed the genomic status of selected genes (Additional file 1: Table S4-S6) with different approaches (FISH, real-time PCR on genomic DNA, multiplex ligation-dependent probes amplification, denaturing high performance liquid chromatography, and sequencing) and analyzed expression of DNA-damage response/repair pathways key components in 19 cell lines using quantitative real-time PCR (Additional file 1: Table S3). Our results show that mRNA expression of PARP1, RAD51, and BRCA1 (but not their genomic status) is directly related to combination synergism (Fig. 3b, Additional file 1: Table S2-S5). Furthermore, in selected experiments, we confirmed these data at protein level by western blot and immunocytochemistry, and demonstrated that the expression of PARP1, BRCA1, and RAD51 was higher in HS-C than in LS-C (Fig. 3c-d).
Next, we evaluated the expression of these putative biomarkers in patients who might benefit from trabectedin + olaparib treatment. We employed immunohistochemistry (IHC) to measure PARP1, BRCA1, and RAD51 expression in a panel of 54 patient-derived sarcoma specimens. PARP1, BRCA1, and RAD51 were overexpressed in 34/53 (64%), 18/54 (33%), and 35/54 (65%) samples, respectively (Fig. 4, Additional file 1: Table S7). For all, IHC expression was characterized by high concordance rates (Additional file 1: Table S8) with an 85% (95% CI: 73–97%) and 86% (95% CI: 75–97%) probability of finding PARP1 overexpressed in BRCA1- and RAD51-positive samples, respectively.
Antitumor activity of trabectedin and olaparib combination in HS-C and LS-C mice models
Thereafter, we studied in vivo antitumor activity of the combination. We subcutaneously injected (s.c.) one tumorigenic HS-C (DMR) and one LS-C (SJSA-1) in nonobese diabetic/sever combined immunodeficient (NOD/SCID) mice. Mice with tumors ≥32 mm3 were randomized to receive drug treatment.
In DMR s.c. model, mean tumor volume was significantly reduced by the treatment with two doses of trabectedin (0.15 and 0.1 mg/kg, p < 0.0001) and not by the lowest dose (0.05 mg/kg; ns, non-significant, p > 0.05). Olaparib as single agent (two doses, 50 and 25 mg/kg/day) did not show antitumor activity. Interestingly, the combination of trabectedin and olaparib at the lowest concentrations (25 mg/kg/day olaparib + 0.05 mg/kg trabectedin) significantly reduced tumor growth if compared to untreated controls and each single agent (Fig. 5a). In SJSA-1 s.c. models the treatment with 0.05 mg/kg of trabectedin significantly reduced tumor growth and no further advantage was obtained by the combination with 25 mg/kg/day of olaparib (ns, p > 0.05) (Fig. 5b). At the lowest dosage (0.025 mg/kg) trabectedin did not impinge tumor growth both as single agent and in combination with olaparib 25 mg/kg (Fig. 5b). Tumor cell proliferation was evaluated by PCNA staining (Fig. 5c), DNA damage by P-H2Ax foci formation and apoptosis by TUNEL staining in combination-treated mice and compared to both controls and single agents (0.05 mg/kg trabectedin and 25 mg/kg olaparib). In DMR s.c. models, cell proliferation was significantly reduced (PCNA, Fig. 5c), and consistently both DNA damage (P-H2AX, Fig. 5c) and apoptosis (TUNEL, Fig. 5c) were significantly increased by the combination if compared to single agents and controls. Evaluation of the same parameters in LS-C SJSA-1 xenografts showed that trabectedin was effective as single agent, and only a marginal advantage was obtained by the combination (Fig. 5c-d). The combination of trabectedin and olaparib did not cause evident signs of toxicity after gross necropsy or weight loss during treatment.
To assess combination activity in the metastatic setting, we developed an orthotopic model of uterine leiomyosarcoma with HS-C DMR cells and an experimental metastatic model of osteosarcoma with LS-C SJSA-1 cells. DMR cells injected into uterine wall tissue developed local and distant metastases within 3 weeks in the untreated control group. Treatment with single agents delayed, but did not impede, tumor progression, whereas trabectedin + olaparib not only significantly reduced primary tumor growth, but also prevented metastatic spread (p < 0.001 vs. single agents and control; Fig. 5e-f). On the contrary, in the metastatic osteosarcoma model, the combination showed non-significant advantage over trabectedin alone (p > 0.05) (Fig. 5g-h). These in vivo experiments confirm data observed in vitro, wherein low-PARP1-expressing cells consistently demonstrated lower synergy of the combination.
PARP1 expression dictates the response to trabectedin and olaparib combination
To generalize to other tumor settings our observations on PARP1 inhibition, we tested trabectedin + olaparib activity in an independent panel of 11 tumor cell lines of different origins (prostate, lung, bile duct, breast), and verified a direct correlation between PARP1 protein expression and combination synergy (Fig. 6a, Additional file 1: Table S2). Finally, the functional role of PARP1 to determine trabectedin + olaparib synergy was further proved in silencing and overexpression experiments. Namely, PARP1 expression was downregulated with specific transient and stable silencing (Fig. 6b) in two HS-C (TC-106, 402.91), which resulted in significantly reduced trabectedin + olaparib activity compared to mock-treated cells (p < 0.001; Fig. 6c). On the other hand, the overexpression of PARP1, determined by stable lentiviral transduction (Fig. 6b) in two LS-C (SJSA-1, SW684), consistently and significantly increased combination activity vs. mock-treated controls (p < 0.001; Fig. 6d). Indeed, activity rose to levels observed in HS-C (Fig. 6c-d).
PARP1 expression determines the cellular response to different DNA-damaging agents
Finally, to generalize the key role played by PARP1 basal expression in synergism, we asked whether enzyme basal expression might drive tumor cell response to other cytotoxics used in combination with PARP1 inhibitors. We investigated PARylation in both low- vs. high-PARP1-expressing cells and PARP1-silenced vs. wild-type counterparts (Additional file 2: Figure S7A) after the treatment with different chemotherapeutic agents: cisplatin, dacarbazine, actinomycin-D, etoposide, doxorubicin, and gemcitabine characterized by different mechanism of action: crosslinking, methylation, transcription inhibition, topoisomerase II inhibition, DNA intercalation, and nucleotide mimicking, respectively. Consistently with trabectedin experiments, we found that PARylation correlated directly with PARP1 basal expression (higher in high- vs. low-PARP1-expressing cells), regardless of the studied cytotoxic or its mechanism of action (Additional file 2: Figure S7B). Indeed, in MSTO-H211 parental cells cisplatin, etoposide, gemcitabine and trabectedin were potent activator of PARP1 activity, while doxorubicin, actinomycin-D and dacarbazine (as pro-drug) less efficiently induced PAR (Additional file 2: Figure S7B). After PARP1 silencing in MSTO-H211 cells, PARylation was invariably decreased if compared to parental cells after each treatment.
Discussion
We identified and demonstrated that PARP1 basal expression is crucial in the mechanism behind the synergy between PARP1 inhibitors and trabectedin independent of BRCA1/2 status, in several robust models across different histotypes. The chemo-sensitization observed in high-PARP1-expressing cells suggests that PARP1 inhibitors may extend and improve the clinical application of DNA-repair targeting when combined with other cytotoxics also.
Our findings support a reappraisal of PARP1 inhibitors as chemo-sensitizers. In general, after chemotherapy exposure, tumor cell fate is highly dependent on the efficiency of the same DNA-damage response pathways that restore normal cell DNA integrity. This observation fueled the idea to target key points in DNA-repair machinery to increase cytotoxicity [2, 4] and brought about the development of several compounds, among which PARP1 inhibitors eventually demonstrated clinical efficacy [4]. However, so far, their clinical use is limited to tumors harboring hereditary or acquired genetic defects in homologous recombination (HR) exploiting an intrinsic cell weakness through so-called “synthetic lethality” [28,29,30,31,32,33,34]. Interestingly, patients affected by HR-defective tumors have the greatest benefit from treatment with DNA-damaging agents, such as platinum compounds and trabectedin [37, 38, 47].
We pursued to induce a sort of “chemical synthetic lethality” also in HR-proficient tumors, by combining trabectedin with PARP1inhibitors. Our hypothesis stemmed from the observation that trabectedin not only induced a peculiar DNA damage but also caused impairment in the DNA-repair machinery that might become lethal in presence of PARP1 inhibitors [35].
Demonstration that trabectedin strongly activates PARP1 enzyme strengthened our idea to combine this chemotherapy with PARP1 inhibitors. However, after trabectedin treatment, we observed significantly different PARylation across various histotypes. For example, single nucleotide polymorphism in the PARP1 gene (Val762Ala) is known to result in a less efficient PARP1 variant [48, 49], which may partly explain differences in the degree of PARylation. To understand the impact of PARylation differences, we turned to PARP1 basal expression and proved that PARP1 protein basal level was not only directly proportional to chemotherapy-induced PARP1 activation, but also that it dictated the synergism with trabectedin and several other chemotherapy compounds as well. Indeed, PARP1-depleted cells were tolerant to PARP1 inhibitors [50, 51], and therefore, PARP1 activity is the prerequisite to induce a significant amount of complexes formed by PARP1, damaged DNA and PARP1 inhibitor that are plausibly more cytotoxic than unrepaired single-strand breaks alone [52,53,54]. Clearly, trabectedin synergized with both tested PARP1 inhibitors, but we found olaparib was significantly more potent than veliparib. Our result reflects the fact that veliparib is a pure catalytic inhibitor of PARP1 activation, whereas olaparib is a poisoning drug that blocks activated PARP1 enzyme at the site of damage (PARP1-trapping activity) impeding the subsequent recruitment of repairing machinery [11, 52, 55]. We studied the synergism of PARP1 inhibitors and trabectedin in a large set of different tumor cell lines and in in vivo models, demonstrating that PARP1 inhibition improved the antitumor activity of trabectedin in a cell-line dependent intensity, as previously shown in breast cancer and in Ewing’s sarcoma cell lines [44, 55]. Among the tested BSTS histotypes, Ewing’s sarcoma cells were the most sensitive to the combination, a result that might be explained by the known exquisite sensitivity of the pathognomonic fusion protein EWS/FLI1-expressing cells to both trabectedin and PARP1 inhibition [10, 56,57,58,59].
Another important datum emerged from study of the MES-SA-DX5 model. This leiomyosarcoma cell line had been previously made resistant to doxorubicin [60], and interestingly, we found it also shared cross-resistance to trabectedin. Of note is that the addition of PARP1 inhibitors restored sensitivity to clinically achievable concentration of trabectedin. The increased activity in MES-SA-DX5 can be explained by comparing the genomic and protein status of this cell line with its parental counterpart. Indeed in MES-SA-DX5, the PARP1 gene was amplified and consequently caused higher PARP1 expression and activity that confirmed the mechanistic role played by PARP1 in determining synergy between trabectedin and PARP1 inhibitors. These consistent observations further prove the key role played by PARP1 expression.
Given the redundant complexity of DNA repair machinery, we took advantage of gene array analyses to delve deeper into the mechanisms behind trabectedin + olaparib synergism. We identified specific gene sets involved in the combination synergism. These genes belong to cell cycle control such as normally activated G2/M cell cycle checkpoints and DNA-damage response and repair pathways (i.e., Timeless, AURKA, MCM2, RAD54L, POLD1, MRE11A) [61,62,63,64]. Moreover, specific genomic aberrations of these genes differentially characterized HS-C and LS-C. We showed that enzymes involved in DNA-damage response and repair were actually expressed more in cell lines displaying higher PARP1 basal expression and high trabectedin + olaparib synergism. The cell cycle phase is crucial for trabectedin and olaparib response and these findings confirmed our data on cell cycle perturbation. Undeniably, trabectedin-induced DNA damage activates the DNA-damage response leading to G2/M cell cycle arrest [35]. However, in this phase, repairing enzymes may still fix the damage caused by either single agent, such that tumor cells survive after initial cell cycle arrest. On the contrary, when the combination was used, DNA repair and cell survival were impeded because of irreparable DNA fragmentation and apoptosis.
Again the intensity of this event directly correlated with PARP1 expression. We carefully validated the role of identified candidate genes in two independent panels of 20 and 11 different cell lines. We selected key enzymes of DNA-damage response and repair pathways and investigated their expression by real-time PCR in the first 20 BSTS cell-line panel. We confirmed a statistically significant direct correlation between PARP1 expression and trabectedin + olaparib synergism. Interestingly, we also observed a significant correlation between drug synergism and RAD51 and BRCA1 expression, which might be explained by the fact that in HR-proficient cells, PARP1 expression drives the transcription of BRCA1 and RAD51 genes co-activating E2F transcription factor [65,66,67,68,69]. If so, then the statistical correlation between combination synergism and RAD51 or BRCA1 expression might simply reflect the level of PARP1 expression. These data are consistent with previous observations of correlation between DNA repair protein expression and sensitivity to PARP1 inhibition [70]. Nonetheless, in HR-defective cells PARP1 activity was found increased [54]. However, we did not find any correlation between the genomic status of key DNA repair genes and the synergism of the combination. To strengthen our observation and generalize our hypothesis on the role of PARP1 expression also in non-sarcoma cells, we replicated our experiments in a second independent panel (11 cell lines of different tumor origin) and confirmed the correlation between trabectedin + olaparib synergy and PARP1 expression also at the protein level. The functional role of PARP1 was further validated in silencing and overexpression experiments that confirmed the data reported above. Noteworthy, the modulation of PARP1 expression or the inhibition of its activity might modulate the nuclear factor Kappa-b pathway that is a crucial determinant of drug resistance [71,72,73]. This intriguing aspect warrants further investigations.
Conclusions
Our data demonstrated that trabectedin is an ideal partner of PARP1 inhibitors because it potently activates PARP1. Furthermore, PARP1 basal expression drives the synergism intensity between PARP1 inhibitors and trabectedin. Indeed, the crucial role of PARP1 expression was confirmed in several BSTS histotypes and in other cancer types as well. Further studies of combinations of PARP1 inhibitors and other cytotoxics should consider basal PARP1 expression and activation after each drug exposure. Finally, olaparib and trabectedin combination is particularly attractive in tumors harboring high PARP1 expression and a specific DNA-damage response/repair and cell cycle control gene signatures that might become predictive biomarkers of response. Our findings pave the way for the future clinical study of trabectedin + olaparib combination independent of BRCA1/2 status.
Abbreviations
- aCGH:
-
Array comparative genomic hybridization
- Ala:
-
Alanine
- ATM:
-
Ataxia telangiectasia mutated
- AURKA:
-
Aurora kinase A
- BRCA1/2:
-
Breast Cancer Type 1/2 susceptibility protein
- BSTS:
-
Bone and soft tissue sarcoma
- CHEK2:
-
Checkpoint kinase 2
- DHPLC:
-
Denaturing high performance liquid chromatography
- DNA:
-
Deoxyribonucleic acid
- EWS:
-
Ewing sarcoma protein
- FDR:
-
False discovery rate
- FISH:
-
Fluorescent in situ hybridization
- FLI1:
-
Friend leukemia integration 1 transcription factor
- GO:
-
Gene ontology
- GSEA:
-
Gene set enrichment analysis
- HS-C:
-
High synergism- cells
- i.v.:
-
Intravenous
- LS-C:
-
Low synergism cells
- MCM2:
-
Minichromosome maintenance protein 2
- MLPA:
-
Multiplex ligation-dependent probes amplification
- MRE11A:
-
Homolog of meiotic recombination 11, S. cerevisiae
- ns:
-
Non-significant
- PARP1:
-
Poly-ADP-ribosyl-transferase-1
- PARylation:
-
Poly- ADP- ribosylation
- PCNA:
-
Proliferating cell nuclear antigen
- PCR:
-
Polymerase chain reaction
- P-H2AX:
-
Phosphorylated histone H2AX
- PI:
-
Propidium iodide
- POLD1:
-
Polymerase (DNA-directed) delta 1
- PTEN:
-
Phosphatase and tensin homolog
- RAD54L:
-
RAD54-Like protein
- s.c.:
-
Subcutaneous
- STR:
-
Short tandem repeat
- TC-NER:
-
transcription-coupled nucleotide excision repair
- TUNEL:
-
Terminal deoxynucleotidyl transferase dUTP nick end labeling
- Val:
-
Valine
References
Satoh MS, Lindahl T. Role of poly(ADP-ribose) formation in DNA repair. Nature. 1992;356:356–8.
Lord CJ, Ashworth A. The DNA damage response and cancer therapy. Nature. 2012;481:287–94.
Ying S, Chen Z, Medhurst AL, Neal JA, Bao Z, Mortusewicz O, McGouran J, Song X, Shen H, Hamdy FC, et al. DNA-PKcs and PARP1 bind to unresected stalled DNA replication forks where they recruit XRCC1 to mediate repair. Cancer Res. 2015;76(5):1078–88.
Scott CL, Swisher EM, Kaufmann SH. Poly (ADP-ribose) polymerase inhibitors: recent advances and future development. J Clin Oncol. 2015;33:1397–406.
Plummer R, Lorigan P, Steven N, Scott L, Middleton MR, Wilson RH, Mulligan E, Curtin N, Wang D, Dewji R, et al. A phase II study of the potent PARP inhibitor, Rucaparib (PF-01367338, AG014699), with temozolomide in patients with metastatic melanoma demonstrating evidence of chemopotentiation. Cancer Chemother Pharmacol. 2013;71:1191–9.
Helleday T. Putting poly (ADP-ribose) polymerase and other DNA repair inhibitors into clinical practice. Curr Opin Oncol. 2013;25:609–14.
LoRusso PM, Li J, Burger A, Heilbrun LK, Sausville E, Boerner S, Smith D, Pilat MJ, Zhang J, Tolaney S, et al. Phase I safety, pharmacokinetic and pharmacodynamic study of the poly (ADP-ribose) polymerase inhibitor veliparib with irinotecan in patients with advanced tumors. Clin Cancer Res. 2016;22(13):3227–37.
Rodler ET, Kurland BF, Griffin M, Gralow JR, Porter P, Yeh RF, Gadi VK, Guenthoer J, Beumer JH, Korde L, et al. Phase I Study of Veliparib (ABT-888) Combined with Cisplatin and Vinorelbine in Advanced Triple-Negative Breast Cancer and/or BRCA Mutation-Associated Breast Cancer. Clin Cancer Res. 2016;22(12):2855–64.
Lee JM, Hays JL, Annunziata CM, Noonan AM, Minasian L, Zujewski JA, Yu M, Gordon N, Ji J, Sissung TM, et al. Phase I/Ib study of olaparib and carboplatin in BRCA1 or BRCA2 mutation-associated breast or ovarian cancer with biomarker analyses. J Natl Cancer Inst. 2014;106:dju089.
Lee HJ, Yoon C, Schmidt B, Park DJ, Zhang AY, Erkizan HV, Toretsky JA, Kirsch DG, Yoon SS. Combining PARP-1 inhibition and radiation in Ewing sarcoma results in lethal DNA damage. Mol Cancer Ther. 2013;12:2591–600.
Murai J, Zhang Y, Morris J, Ji J, Takeda S, Doroshow JH, Pommier Y. Rationale for poly(ADP-ribose) polymerase (PARP) inhibitors in combination therapy with camptothecins or temozolomide based on PARP trapping versus catalytic inhibition. J Pharmacol Exp Ther. 2014;349:408–16.
Engert F, Schneider C, Weiβ LM, Probst M, Fulda S. PARP inhibitors sensitize Ewing sarcoma cells to temozolomide-induced apoptosis via the mitochondrial pathway. Mol Cancer Ther. 2015;14:2818–30.
Gill SJ, Travers J, Pshenichnaya I, Kogera FA, Barthorpe S, Mironenko T, Richardson L, Benes CH, Stratton MR, McDermott U, et al. Combinations of PARP inhibitors with temozolomide drive PARP1 trapping and apoptosis in Ewing’s sarcoma. PLoS One. 2015;10:e0140988.
Curtin NJ. DNA repair dysregulation from cancer driver to therapeutic target. Nat Rev Cancer. 2012;12:801–17.
Rouleau M, Patel A, Hendzel MJ, Kaufmann SH, Poirier GG. PARP inhibition: PARP1 and beyond. Nat Rev Cancer. 2010;10:293–301.
Verhagen CV, de Haan R, Hageman F, Oostendorp TP, Carli AL, O’Connor MJ, Jonkers J, Verheij M, van den Brekel MW, Vens C. Extent of radiosensitization by the PARP inhibitor olaparib depends on its dose, the radiation dose and the integrity of the homologous recombination pathway of tumor cells. Radiother Oncol. 2015;116:358–65.
Rajan A, Carter CA, Kelly RJ, Gutierrez M, Kummar S, Szabo E, Yancey MA, Ji J, Mannargudi B, Woo S, et al. A phase I combination study of olaparib with cisplatin and gemcitabine in adults with solid tumors. Clin Cancer Res. 2012;18:2344–51.
Bendell J, O’Reilly EM, Middleton MR, Chau I, Hochster H, Fielding A, Burke W, Burris H. Phase I study of olaparib plus gemcitabine in patients with advanced solid tumours and comparison with gemcitabine alone in patients with locally advanced/metastatic pancreatic cancer. Ann Oncol. 2015;26:804–11.
Liu JF, Barry WT, Birrer M, Lee JM, Buckanovich RJ, Fleming GF, Rimel B, Buss MK, Nattam S, Hurteau J, et al. Combination cediranib and olaparib versus olaparib alone for women with recurrent platinum-sensitive ovarian cancer: a randomised phase 2 study. Lancet Oncol. 2014;15:1207–14.
Khan OA, Gore M, Lorigan P, Stone J, Greystoke A, Burke W, Carmichael J, Watson AJ, McGown G, Thorncroft M, et al. A phase I study of the safety and tolerability of olaparib (AZD2281, KU0059436) and dacarbazine in patients with advanced solid tumours. Br J Cancer. 2011;104:750–5.
Menear KA, Adcock C, Boulter R, Cockcroft XL, Copsey L, Cranston A, Dillon KJ, Drzewiecki J, Garman S, Gomez S, et al. 4-[3-(4-cyclopropanecarbonylpiperazine-1-carbonyl)-4-fluorobenzyl]-2H-phthalazin-1-one: a novel bioavailable inhibitor of poly(ADP-ribose) polymerase-1. J Med Chem. 2008;51:6581–91.
Farmer H, McCabe N, Lord CJ, Tutt AN, Johnson DA, Richardson TB, Santarosa M, Dillon KJ, Hickson I, Knights C, et al. Targeting the DNA repair defect in BRCA mutant cells as a therapeutic strategy. Nature. 2005;434:917–21.
McCabe N, Turner NC, Lord CJ, Kluzek K, Bialkowska A, Swift S, Giavara S, O’Connor MJ, Tutt AN, Zdzienicka MZ, et al. Deficiency in the repair of DNA damage by homologous recombination and sensitivity to poly(ADP-ribose) polymerase inhibition. Cancer Res. 2006;66:8109–15.
Plummer R, Jones C, Middleton M, Wilson R, Evans J, Olsen A, Curtin N, Boddy A, McHugh P, Newell D, et al. Phase I study of the poly(ADP-ribose) polymerase inhibitor, AG014699, in combination with temozolomide in patients with advanced solid tumors. Clin Cancer Res. 2008;14:7917–23.
Kummar S, Chen A, Ji J, Zhang Y, Reid JM, Ames M, Jia L, Weil M, Speranza G, Murgo AJ, et al. Phase I study of PARP inhibitor ABT-888 in combination with topotecan in adults with refractory solid tumors and lymphomas. Cancer Res. 2011;71:5626–34.
Hussain M, Carducci MA, Slovin S, Cetnar J, Qian J, McKeegan EM, Refici-Buhr M, Chyla B, Shepherd SP, Giranda VL, Alumkal JJ. Targeting DNA repair with combination veliparib (ABT-888) and temozolomide in patients with metastatic castration-resistant prostate cancer. Invest New Drugs. 2014;32:904–12.
Kummar S, Ji J, Morgan R, Lenz HJ, Puhalla SL, Belani CP, Gandara DR, Allen D, Kiesel B, Beumer JH, et al. A phase I study of veliparib in combination with metronomic cyclophosphamide in adults with refractory solid tumors and lymphomas. Clin Cancer Res. 2012;18:1726–34.
Fong PC, Boss DS, Yap TA, Tutt A, Wu P, Mergui-Roelvink M, Mortimer P, Swaisland H, Lau A, O’Connor MJ, et al. Inhibition of poly(ADP-ribose) polymerase in tumors from BRCA mutation carriers. N Engl J Med. 2009;361:123–34.
Fong PC, Yap TA, Boss DS, Carden CP, Mergui-Roelvink M, Gourley C, De Greve J, Lubinski J, Shanley S, Messiou C, et al. Poly(ADP)-ribose polymerase inhibition: frequent durable responses in BRCA carrier ovarian cancer correlating with platinum-free interval. J Clin Oncol. 2010;28:2512–9.
Ledermann J, Harter P, Gourley C, Friedlander M, Vergote I, Rustin G, Scott CL, Meier W, Shapira-Frommer R, Safra T, et al. Olaparib maintenance therapy in patients with platinum-sensitive relapsed serous ovarian cancer: a preplanned retrospective analysis of outcomes by BRCA status in a randomised phase 2 trial. Lancet Oncol. 2014;15:852–61.
Ledermann J, Harter P, Gourley C, Friedlander M, Vergote I, Rustin G, Scott C, Meier W, Shapira-Frommer R, Safra T, et al. Olaparib maintenance therapy in platinum-sensitive relapsed ovarian cancer. N Engl J Med. 2012;366:1382–92.
Kaufman B, Shapira-Frommer R, Schmutzler RK, Audeh MW, Friedlander M, Balmaña J, Mitchell G, Fried G, Stemmer SM, Hubert A, et al. Olaparib monotherapy in patients with advanced cancer and a germline BRCA1/2 mutation. J Clin Oncol. 2015;33:244–50.
Mateo J, Carreira S, Sandhu S, Miranda S, Mossop H, Perez-Lopez R, Nava Rodrigues D, Robinson D, Omlin A, Tunariu N, et al. DNA-repair defects and olaparib in metastatic prostate cancer. N Engl J Med. 2015;373:1697–708.
Bryant HE, Schultz N, Thomas HD, Parker KM, Flower D, Lopez E, Kyle S, Meuth M, Curtin NJ, Helleday T. Specific killing of BRCA2-deficient tumours with inhibitors of poly(ADP-ribose) polymerase. Nature. 2005;434:913–7.
Larsen AK, Galmarini CM, D’Incalci M. Unique features of trabectedin mechanism of action. Cancer Chemother Pharmacol. 2015;77(4):663–71.
Le Cesne A, Blay JY, Domont J, Tresch-Bruneel E, Chevreau C, Bertucci F, Delcambre C, Saada-Bouzid E, Piperno-Neumann S, Bay JO, et al. Interruption versus continuation of trabectedin in patients with soft-tissue sarcoma (T-DIS): a randomised phase 2 trial. Lancet Oncol. 2015;16:312–9.
Schöffski P, Taron M, Jimeno J, Grosso F, Sanfilipio R, Casali PG, Le Cesne A, Jones RL, Blay JY, Poveda A, et al. Predictive impact of DNA repair functionality on clinical outcome of advanced sarcoma patients treated with trabectedin: a retrospective multicentric study. Eur J Cancer. 2011;47:1006–12.
Italiano A, Laurand A, Laroche A, Casali P, Sanfilippo R, Le Cesne A, Judson I, Blay JY, Ray-Coquard I, Bui B, et al. ERCC5/XPG, ERCC1, and BRCA1 gene status and clinical benefit of trabectedin in patients with soft tissue sarcoma. Cancer. 2011;117:3445–56.
Demetri GD, von Mehren M, Jones RL, Hensley ML, Schuetze SM, Staddon A, Milhem M, Elias A, Ganjoo K, Tawbi H, et al. Efficacy and safety of trabectedin or dacarbazine for metastatic liposarcoma or leiomyosarcoma after failure of conventional chemotherapy: results of a phase III randomized multicenter clinical trial. J Clin Oncol. 2016;34:786–93.
Moneo V, Serelde BG, Fominaya J, Leal JF, Blanco-Aparicio C, Romero L, Sánchez-Beato M, Cigudosa JC, Tercero JC, Piris MA, et al. Extreme sensitivity to Yondelis (Trabectedin, ET-743) in low passaged sarcoma cell lines correlates with mutated p53. J Cell Biochem. 2007;100:339–48.
Blay JY, Leahy MG, Nguyen BB, Patel SR, Hohenberger P, Santoro A, Staddon AP, Penel N, Piperno-Neumann S, Hendifar A, et al. Randomised phase III trial of trabectedin versus doxorubicin-based chemotherapy as first-line therapy in translocation-related sarcomas. Eur J Cancer. 2014;50:1137–47.
Germano G, Frapolli R, Simone M, Tavecchio M, Erba E, Pesce S, Pasqualini F, Grosso F, Sanfilippo R, Casali PG, et al. Antitumor and anti-inflammatory effects of trabectedin on human myxoid liposarcoma cells. Cancer Res. 2010;70:2235–44.
Mabuchi S, Hisamatsu T, Kawase C, Hayashi M, Sawada K, Mimura K, Takahashi K, Takahashi T, Kurachi H, Kimura T. The activity of trabectedin as a single agent or in combination with everolimus for clear cell carcinoma of the ovary. Clin Cancer Res. 2011;17:4462–73.
Ordóñez JL, Amaral AT, Carcaboso AM, Herrero-Martín D, Del Carmen G-MM, Sevillano V, Alonso D, Pascual-Pasto G, San-Segundo L, Vila-Ubach M, et al. The PARP inhibitor olaparib enhances the sensitivity of Ewing sarcoma to trabectedin. Oncotarget. 2015;6:18875–90.
Vilar E, Bartnik CM, Stenzel SL, Raskin L, Ahn J, Moreno V, Mukherjee B, Iniesta MD, Morgan MA, Rennert G, Gruber SB. MRE11 deficiency increases sensitivity to poly(ADP-ribose) polymerase inhibition in microsatellite unstable colorectal cancers. Cancer Res. 2011;71:2632–42.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Konstantinopoulos PA, Ceccaldi R, Shapiro GI, D’Andrea AD. Homologous recombination deficiency: exploiting the fundamental vulnerability of ovarian cancer. Cancer Discov. 2015;5:1137–54.
Lockett KL, Hall MC, Xu J, Zheng SL, Berwick M, Chuang SC, Clark PE, Cramer SD, Lohman K, Hu JJ. The ADPRT V762A genetic variant contributes to prostate cancer susceptibility and deficient enzyme function. Cancer Res. 2004;64:6344–8.
Smith TR, Levine EA, Freimanis RI, Akman SA, Allen GO, Hoang KN, Liu-Mares W, Hu JJ. Polygenic model of DNA repair genetic polymorphisms in human breast cancer risk. Carcinogenesis. 2008;29:2132–8.
Pettitt SJ, Rehman FL, Bajrami I, Brough R, Wallberg F, Kozarewa I, Fenwick K, Assiotis I, Chen L, Campbell J, et al. A genetic screen using the PiggyBac transposon in haploid cells identifies Parp1 as a mediator of olaparib toxicity. PLoS One. 2013;8:e61520.
Liu X, Han EK, Anderson M, Shi Y, Semizarov D, Wang G, McGonigal T, Roberts L, Lasko L, Palma J, et al. Acquired resistance to combination treatment with temozolomide and ABT-888 is mediated by both base excision repair and homologous recombination DNA repair pathways. Mol Cancer Res. 2009;7:1686–92.
Murai J, Huang SY, Das BB, Renaud A, Zhang Y, Doroshow JH, Ji J, Takeda S, Pommier Y. Trapping of PARP1 and PARP2 by Clinical PARP Inhibitors. Cancer Res. 2012;72:5588–99.
Murai J, Huang SY, Renaud A, Zhang Y, Ji J, Takeda S, Morris J, Teicher B, Doroshow JH, Pommier Y. Stereospecific PARP trapping by BMN 673 and comparison with olaparib and rucaparib. Mol Cancer Ther. 2014;13:433–43.
Gottipati P, Vischioni B, Schultz N, Solomons J, Bryant HE, Djureinovic T, Issaeva N, Sleeth K, Sharma RA, Helleday T. Poly(ADP-ribose) polymerase is hyperactivated in homologous recombination-defective cells. Cancer Res. 2010;70:5389–98.
Ávila-Arroyo S, Nuñez GS, García-Fernández LF, Galmarini CM. Synergistic effect of trabectedin and olaparib combination regimen in breast cancer cell lines. J Breast Cancer. 2015;18:329–38.
Garnett MJ, Edelman EJ, Heidorn SJ, Greenman CD, Dastur A, Lau KW, Greninger P, Thompson IR, Luo X, Soares J, et al. Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature. 2012;483:570–5.
Brenner JC, Feng FY, Han S, Patel S, Goyal SV, Bou-Maroun LM, Liu M, Lonigro R, Prensner JR, Tomlins SA, Chinnaiyan AM. PARP-1 inhibition as a targeted strategy to treat Ewing’s sarcoma. Cancer Res. 2012;72:1608–13.
Grohar PJ, Griffin LB, Yeung C, Chen QR, Pommier Y, Khanna C, Khan J, Helman LJ. Ecteinascidin 743 interferes with the activity of EWS-FLI1 in Ewing sarcoma cells. Neoplasia. 2011;13:145–53.
Grohar PJ, Segars LE, Yeung C, Pommier Y, D’Incalci M, Mendoza A, Helman LJ. Dual targeting of EWS-FLI1 activity and the associated DNA damage response with trabectedin and SN38 synergistically inhibits Ewing sarcoma cell growth. Clin Cancer Res. 2014;20:1190–203.
Harker WG, Sikic BI. Multidrug (pleiotropic) resistance in doxorubicin-selected variants of the human sarcoma cell line MES-SA. Cancer Res. 1985;45:4091–6.
Xie S, Mortusewicz O, Ma HT, Herr P, Poon RR, Helleday T, Qian C. Timeless interacts with PARP-1 to promote homologous recombination repair. Mol Cell. 2015;60:163–76.
Dedes KJ, Wilkerson PM, Wetterskog D, Weigelt B, Ashworth A, Reis-Filho JS. Synthetic lethality of PARP inhibition in cancers lacking BRCA1 and BRCA2 mutations. Cell Cycle. 2011;10:1192–9.
Bajrami I, Frankum JR, Konde A, Miller RE, Rehman FL, Brough R, Campbell J, Sims D, Rafiq R, Hooper S, et al. Genome-wide profiling of genetic synthetic lethality identifies CDK12 as a novel determinant of PARP1/2 inhibitor sensitivity. Cancer Res. 2014;74:287–97.
Young LM, Marzio A, Perez-Duran P, Reid DA, Meredith DN, Roberti D, Star A, Rothenberg E, Ueberheide B, Pagano M. TIMELESS forms a complex with PARP1 distinct from its complex with TIPIN and plays a role in the DNA damage response. Cell Rep. 2015;13:451–9.
Simbulan-Rosenthal CM, Rosenthal DS, Luo R, Samara R, Espinoza LA, Hassa PO, Hottiger MO, Smulson ME. PARP-1 binds E2F-1 independently of its DNA binding and catalytic domains, and acts as a novel coactivator of E2F-1-mediated transcription during re-entry of quiescent cells into S phase. Oncogene. 2003;22:8460–71.
Hegan DC, Lu Y, Stachelek GC, Crosby ME, Bindra RS, Glazer PM. Inhibition of poly(ADP-ribose) polymerase down-regulates BRCA1 and RAD51 in a pathway mediated by E2F4 and p130. Proc Natl Acad Sci U S A. 2010;107:2201–6.
Byers LA, Wang J, Nilsson MB, Fujimoto J, Saintigny P, Yordy J, Giri U, Peyton M, Fan YH, Diao L, et al. Proteomic profiling identifies dysregulated pathways in small cell lung cancer and novel therapeutic targets including PARP1. Cancer Discov. 2012;2:798–811.
Wang A, Schneider-Broussard R, Kumar AP, MacLeod MC, Johnson DG. Regulation of BRCA1 expression by the Rb-E2F pathway. J Biol Chem. 2000;275:4532–6.
Green AR, Caracappa D, Benhasouna AA, Alshareeda A, Nolan CC, Macmillan RD, Madhusudan S, Ellis IO, Rakha EA. Biological and clinical significance of PARP1 protein expression in breast cancer. Breast Cancer Res Treat. 2015;149:353–62.
Cardnell RJ, Feng Y, Diao L, Fan YH, Masrorpour F, Wang J, Shen Y, Mills GB, Minna JD, Heymach JV, Byers LA. Proteomic markers of DNA repair and PI3K pathway activation predict response to the PARP inhibitor BMN 673 in small cell lung cancer. Clin Cancer Res. 2013;19:6322–8.
Stilmann M, Hinz M, Arslan SC, Zimmer A, Schreiber V, Scheidereit C. A nuclear poly(ADP-ribose)-dependent signalosome confers DNA damage-induced IkappaB kinase activation. Mol Cell. 2009;36:365–78.
Godwin P, Baird AM, Heavey S, Barr MP, O’Byrne KJ, Gately K. Targeting nuclear factor-kappa B to overcome resistance to chemotherapy. Front Oncol. 2013;3:120.
Nakagawa Y, Sedukhina AS, Okamoto N, Nagasawa S, Suzuki N, Ohta T, Hattori H, Roche-Molina M, Narváez AJ, Jeyasekharan AD, Bernal JA, Sato K. NF-κB signaling mediates acquired resistance after PARP inhibition. Oncotarget. 2015;6:3825–39.
Acknowledgements
The authors would like to thank Joan C. Leonard for language editing; Roberta Porporato, Daniela Cantarella, Barbara Martinoglio and Stefania Giove for technical support; Claudio Isella for helpful discussion. Trabectedin was kindly provided by PharmaMar (Madrid, Spain).
Funding
This work was supported by the Italian Association for Cancer Research “Associazione Italiana per la Ricerca sul Cancro” AIRC: IG 2015 Id.17226 (GG); AIRC IG Id 11515 (MA); MFAG N.15731 (DS); “Fondazione per la Ricerca sui Tumori dell’Apparato muscoloscheletrico e Rari Onlus” CRT RF = 2016–0917 (GG); “Ricerca Finalizzata-Ministero della Salute” GR-2011-02349197 (DS); FPRC ONLUS-5 × 1000 "Ministero della Salute" 2012 (DS); University of Torino-“Progetti di Ateneo” 2011 grant Rethe- ORTO11RKTW (MA).
Availability of data and materials
The datasets generated and analyzed during the current study are available in the GEO repository, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE77175, and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE76981.
Authors’ contributions
YP, conceived, designed, performed experiments and wrote manuscript, FC, CDA, MB, FV, AL, GC and MC performed experiments; LDA wrote the manuscript; LDA, LG, DS edited the manuscript; BP supervised MLPA/DHPLC experiments and analysis, MSB provided cell lines; DG, SM, PB, LDA, SA, EP, acquired and managed patients; MA Supervised the study and edited the manuscript; and GG conceived and supervised the study and wrote the manuscript. All authors read and approved the final manuscript.
Competing interests
GG declares Board participation for Bayer, Novartis, PharmaMar, Eisai and reimbursement of travel expenses for meeting participation from PharmaMar. All other authors declare no competing financial interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
This work was done in accordance with IRB approval (202/2014). The archival sample were collected after signed informed consent obtainment. The in vivo experiments in mice were conducted in accordance with the protocol approved by the Institutional Ethics Review Board and by the Italian Ministry of Health (Aut. Min. 178/201S-PR).
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Authors and Affiliations
Corresponding authors
Additional files
Additional file 1: Table S1.
Concentrations inhibiting 50% of cell viability (IC50) after 72- h treatment with serial dilutions of trabectedin (2–0.125 nM), olaparib (20–1.25 μM), and veliparib (80–5 μM) as single agents (a) or in constant combination (b, c), and their 95% confidence intervals (95% CI). Population doubling time. Combination index ± estimated standard deviation (est st dev) calculated at IC50 by the Chou-Talalay method for trabectedin and olaparib in combination and trabectedin and veliparib in combination in cell lines of different histotypes: leiomyosarcoma (LMS), undifferentiated pleomorphic sarcoma (UPS), myxoid liposarcoma (LMS), dedifferentiated liposarcoma (DLS), fibrosarcoma (FSA), synovial sarcoma (SS), Ewing’s sarcoma (ES), and osteosarcoma (OS). Supplier and culture conditions are reported for each cell line. Table S2. Direct correlation between combination index (calculated at IC50) and PARP1 protein intensity expressed as Pearson score (r Pearson). IC50 and 95% confidence intervals (95% CI) were calculated after 72-h treatment with serial dilutions of trabectedin (2–0.125 nM), olaparib (20–1.25μM), and their constant combination. Cell line characteristics, population doubling time, purchasers and culture conditions were included. Table S3. Gene expression (ΔCT) of DNA-damage response and repair key components and drug synergism expressed by combination index (CI). The correlation between each gene expression and the CI was evaluated by Pearson score (r); t distribution and their relative P value were shown. Yellow cells highlight significant direct correlation. Table S4. PARP1 gene (chromosome 1 q42.12d) copy number obtained by FISH. Table S5. PARP1, BRCA1, RAD51 gene copy number obtained by real- time PCR on genomic DNA. The gene copy number of PARP1, RAD51, and BRCA1 did not correlate with the Combination index (CI) as shown by Pearson score. Table S6. Genomic status of selected genes analyzed by MLPA and DHPLC /Sequencing. Red cells indicate increased copy number, while blue cells indicate reduced copy number as obtained by DHPLC analysis. Table S7. Immunohistochemistry score of intensity for PARP1, BRCA1, and RAD51 protein expression in formalin-fixed paraffin-embedded sarcoma samples. Table S8. 2 × 2 contingency tables of immunohistochemistry (IHC) expression of PARP1, BRCA1, and RAD51 in patient-derived soft tissue and bone sarcoma specimens (a, b, c) and related concordance rates (d). (DOCX 2086 kb)
Additional file 2: Figure S1.
Overview of gene expression analysis. GSEA, gene signature enrichment analysis. Figure S2. DNA sequences of single nucleotide polymorphism at codon 762 of PARP1 gene in HT1080, SJSA-1, and SW684 cells. Figure S3. Distribution of trabectedin IC50 as single agent (TR alone) and in combination with veliparib (TR + VEL) or olaparib (TR + OL) among high-PARP1-expressing cells (red triangle) and low-PARP1-expressing cells (blue triangle). Figure S4. Dose- response curve obtained after 72-h treatment with trabectedin (2–0.125nM), olaparib (20–1.25 μM) as single agents and in constant combination. Figure S5. A, western blot analysis of PARylation and PARP1 expression in MES-SA and MES-SA-DX5 leiomyosarcoma cells; B, FISH analysis of PARP1 gene (red) and centromere of chromosome 1 (green) in MESSA and MESSA-DX5. Figure S6. Genomic status as obtained by aCGH analysis of TC-106, 402.91, DMR, SJSA-1, HT1080, SW684: gain (red) and loss (green) of chromosome regions. Figure S7. A, Western blot analysis of PARylation and PARP1 expression and B, quantitation of PAR in MSTO-H211, and PARP1-silenced MSTO-H211 untreated or treated with 10nM trabectedin, 20 μM cisplatin (Sandoz), 20 μM gemcitabine (Sandoz), 20 μM doxorubicin (Pfizer), 20 μM dacarbazine (Medac), 20 μM etoposide (Teva), 50 mM actinomycin-D (Thermo Fisher Scientific), β-actin was done as loading control. (DOCX 5982 kb)
Additional file 3:
Analysis of differential genomic aberrations in HS-C and LS-C cells. (PDF 255 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Pignochino, Y., Capozzi, F., D’Ambrosio, L. et al. PARP1 expression drives the synergistic antitumor activity of trabectedin and PARP1 inhibitors in sarcoma preclinical models. Mol Cancer 16, 86 (2017). https://doi.org/10.1186/s12943-017-0652-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12943-017-0652-5