2Adrem Data Lab, Department of Computer Science,
3Molecular Parasitology Unit, Department of Biomedical Sciences,
4Department of Pathology,
© The Author(s) 2020. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as 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.
Aim: The purpose of this manuscript is to study the potential characteristics of acquired nutlin-3 resistant non-small cell lung cancer cells (NSCLC). Nutlin-3 is an inhibitor of the murine-double minute 2 protein, the main negative regulator of wild type p53, of which several derivatives are currently in clinical development.
Methods: A549 NSCLC cells were exposed to increasing concentrations of nutlin-3 for a period of 18 weeks. Monoclonal derivates were cultured, and the most resistance subclone was selected for whole transcriptome analysis. Gene set enrichment analysis was performed on differentially expressed genes between A549 nutlin-3 resistant cancer cells and the parental A549 p53 wild type cancer cells. Relevant findings were validated at the gene, protein and/or functional level.
Results: All nutlin-3 resistant subclones acquired mutations in the TP53 gene, resulting in overexpression of the mutant p53 protein. The most resistant subclone was enriched for genes related to epithelial to mesenchymal transition (EMT), resulting in increased migratory and invasive potential. Furthermore, these cells were enriched in genes related to inflammation, tissue remodelling, and angiogenesis. Importantly, expression of several immune checkpoints, including PD-L1 and PD-L2, was significantly upregulated, and cisplatin-induced cell death was reduced.
Conclusion: Transcriptome analysis of a highly nutlin-3 resistant A549 subclone shows the relevance of studying (1) resistance to standard of care chemotherapy; (2) secretion of immunomodulating chemo- and cytokines; (3) immune checkpoint expression; and (4) EMT and invasion in nutlin-3 resistant cancer cells in addition to acquired mutations in the TP53 gene.
p53, murine-double minute 2, non-small cell lung cancer, acquired resistance, nutlin-3
The p53 protein is considered to be “the guardian of the genome” and protects healthy cells against malignant transformation. Consequently, its function is often disturbed in cancer cells by either inactivating or gain-of-function mutations in the TP53 gene. Alternatively, wild type p53 levels can be suppressed by overexpression of p53’s main negative regulator, i.e., murine double minute 2 (MDM2) protein[2,3]. For the latter, selective inhibitors of the interaction between MDM2 and p53 have been developed with nutlin-3 being the first-in-class small molecule. Increased levels of wild type p53 following nutlin-3 treatment lead to activation of downstream effector pathways related to cell cycle arrest, senescence, or apoptosis. This activation can be further enforced in combination with stress-inducing agents like cisplatin, inducing DNA damage and leading to synergistic anti-cancer effects in wild type p53 cancer cells.
Resistance to targeted therapies remains a major problem in the clinic and predictive biomarkers for therapy response are crucial. For nutlin-3, or its derivates currently in clinical development (RG7112 and RG7388), it is clear that wild type p53 is an important prerequisite for response. However, little is known about acquired resistance mechanisms and the phenotype of these resistant cancer cells. Therefore, we generated nutlin-3 resistant A549 non-small cell lung cancer cells (NSCLC). Following monoclonal expansion of these resistant cells, we selected the most resistant subclone for in-depth analysis of its characteristics by means of whole transcriptome analysis and provide an overview of characteristics worth studying in nutlin-3 resistant cancer cells.
The NSCLC cell lines A549 (p53 wild type) and NCI-H1975 (p53 mutant R273H) were obtained from ATCC. Cells were grown as monolayers and maintained in exponential growth in 5% CO2/95% air in a humidified incubator at 37 °C. A549 and NCI-H1975 cells were, respectively, cultured in DMEM or RPMI (Life Technologies) supplemented with 10% fetal bovine serum (FBS, Life Technologies), 1% penicillin/streptomycin (Life technologies) and 2 mmol/L L-glutamine (Life technologies). A549 cells were exposed to increasing concentrations of nutlin-3 (Tocris) over a period of 18 weeks [Figure 1A, Supplementary Table 1] to generate nutlin-3 resistant cancer (A549.R) cells. After 18 weeks, cells were diluted to 1000 cells/mL and 1 µL was added to a 96-well plate to generate monoclonal subclones. Four subclones (A549.R1 - R4) were successfully generated. A549.R2 was identified as the most resistant subclone and used for further characterization of nutlin-3 resistance mechanisms.
Figure 1. Acquired nutlin-3 resistance. A549 cells were exposed to increasing concentrations of nutlin-3 for a period of 18 weeks. From these nutlin-3 exposed A549 cells, four monoclonal subclones were cultured (A); dose response curve for nutlin-3 of A549 cells, the parental resistant A549.R cells and four monoclonal subclones using the SRB assay (A549.R1-R4) (B); corresponding nutlin-3 IC50 values (C); percentage survival (D); and percentage cell death of A549 and A549.R2 cells treated with idasanutlin (0-50 µmol/L) for 96 h determined by the IncuCyte ZOOM cytotoxicity assay (E). Data presented as mean +/- SD of at least 3 biological replicates
A separate batch of A549 and A549.R2 cells were stably transduced with the Incucyte® Nuclight red lentivirus reagent with a puromycin selection gene to allow cell tracking on the IncuCyte ZOOM system (Sartorius). All cell lines were confirmed as mycoplasma-free using the MycoAlert Mycoplasma Detection Kit.
A colometric sulforhodamine B (SRB) assay was used to measure treatment-induced cytotoxicity to nutlin-3. 1500 cells/well were plated and allowed to settle overnight. Afterwards they were exposed to 0-30 µmol/L nutlin-3 for 72 h. Cell monolayers were then fixed with 10% trichloroacetic acid for 1 h at 4 °C and stained with 100 µL 0.1% SRB, as previously described. Dose response curves were plotted and corresponding IC50 values (i.e., drug concentrations causing 50% growth inhibition) were calculated using GraphPad Prism software.
Nuclight red transduced A549 and A549.R2 cells were seeded at a concentration of 500 cells/well in a 384-well plate and allowed to settle overnight. A 7-point titration of idasanutlin (0-50 µmol/L, Tocris, DMSO), APR-246 (0-50 µmol/L, Tocris, DMSO), cisplatin (0-20 µmol/L, Tocris, 0.3%Tween20 in 0.9%NaCl) and 75 nM cytotox green reagent (Sartorius, DMSO) was added using the Tecan D300e digital drug dispenser. After 96 hours the plate was scanned with the IncuCyte ZOOM and % survival was determined with the formula:
and % cell death with the formula:
Dose response curves were plotted with GraphPad Prism software.
Nuclight red transduced A549 and A549.R2 cells were seeded at a concentration of 20.000 cells/well in an ImageLock 96-well plate (Sartorius) and allowed to settle overnight. Scratch wounds were made with the WoundMaker tool according to the manufacturer’s instruction. For invasion, scratched cells were overlayed with 6 mg/mL Matrigel matrix. Cell migration and invasion was monitored every 2 h using the IncuCyte Zoom system and analysed using the IncuCyte Scratch Wound Analysis Software module.
Total RNA was isolated from 3 biological replicates of the parental A549 and nutlin-3 resistant A549.R2 cell lines using the Qiagen RNA Plus Mini Kit, according to the manufacturer’s instructions. RNA concentration and purity were checked using the NanoDrop ND-1000 (ThermoFisher). Samples were frozen at -80 °C and sent to BGI Genomics for quality control and transcriptome sequencing (oligo dT, non-stranded DNBseq platform, PE150 30M reads/9Gb). Raw sequencing reads were aligned to the human reference genome v38, and per gene read counts generated with STAR v2.7.2b. Differential expression analysis was carried out in R with the DEseq2 package. Volcano plots were generated with EnhancedVolcano in R. Preranked Gene Set Enrichment Analysis (GSEA) was performed using the www.webgestalt.org analysis tool kit v2019 for WikiPathwaysCancer and KEGG databases and the Hallmark and Oncogenic Signatures gene sets downloaded from the Molecular Signatures Database v7.1 database. Gene sets with a false discovery rate (FDR) ≤ 0.05 and P-value ≤ 0.05 were considered.
DNA from A549 and its isogenic nutlin-3 resistant subclones was isolated using the GenElute Blood Genomic DNA kit (Sigma Aldrich) and DNA concentration and purity were checked using the NanoDrop ND-1000.
TP53 NGS was performed using the TP53 MASTRTM with MID for Illumina Miseq kit (Multiplicom) according to the manufacturer’s instructions. Sequencing was performed on a MiSeq system using the MiSeq Reagent Kit v2 (500 cycles, Illumina). Using an in-house annotation and filtering tool, VariantDB, single nucleotide variants were annotated to the TP53 NM_000546 transcript file. SNVs present in the dbSNP137 and 1000 Genomes Project databases were identified as single nucleotide polymorphisms and excluded. The reads were visualised using the Integrative Genomics Viewer software (IGV, version 2.3.67) and aligned to the human reference genome (hg19, NCBI build 37). Using the MUT-TP53 2.0 tool, inactivating TP53 mutations were identified.
Cells were lysed in lysis buffer (10 mmol/L TrisHCl, 400 mmol/L NaCl, 1 mM EDTA, 0.1% NP40 and protease inhibitor). After centrifugation (10 min, 13 000 rpm, 4 °C), cleared lysates containing the isolated proteins were harvested and kept at -80 °C. Protein concentrations were determined using the Pierce BCA protein kit, according to the manufacturer’s instructions. Western blot was performed as described previously. Blocking, primary and secondary antibody incubation was performed using the SNAP id 2.0 protein detection system, according to the manufacturer’s instructions. Membranes were incubated with rabbit anti-p53 (1:1000, Cell Signalling #9282) and mouse anti-b-actin (1:2500, Sigma-Aldrich). Goat anti-rabbit (1:10 000, Licor IRDye 800CW) and goat anti-mouse (1:10 000, Licor IRDye 680RD) fluorescently labelled secondary antibodies were used. Fluorescent detection was performed using the Odyssey imaging system (Li-Cor). Image Studio Lite software (Li-Cor) was used to perform pixel quantification of the images. Normalization against the internal actin control was performed for each sample.
Long term exposure to nutlin-3 resulted in a slight decrease in sensitivity of the A549 vs. A549.R cells (IC50: 7.7 µmol/L vs. 11.0 µmol/L; Figure 1B and C). Within the monoclonal subclones, a strong variability in sensitivity was observed, with the A549.R2 subclone being the most resistant clone (IC50: 22.2 µmol/L; Figure 1B and C). Reduced sensitivity is also retained for idasanutlin (RG7388, IC50: 8.45 µmol/L vs. 16.04 µmol/L), a more potent second-generation nutlin-3 currently in clinical development, which induced a predominant cytostatic response [Figure 1D and E]. We selected the A549.R2 clone for further downstream analysis of the molecular characteristics of acquired nutlin-3 resistance by means of transcriptome analysis. Supplementary Figure 1 represents the most differentially expressed genes between A549.R2 and the parental A549 cells. This differentially expressed gene sets was pre-ranked and subjected to GSEA enrichment analysis and the top gene signatures are presented in Figure 2 and Supplementary Table 2 (FDR ≤ 0.05; P-value ≤ 0.05).
Figure 2. Gene set enrichment analysis of A549.R2 nutlin-3 resistant cells compared to A549 cells. Normalized enrichment scores (NES) of gene expression signatures with a FDR ≤ 0.05 and P-value ≤ 0.05 are presented from the WikiPathwaysCancer, KEGG, Hallmark and Oncogenic Signatures gene sets. A detailed overview and gene set names are presented in Supplementary Table 2.
GSEA analysis indicates that A549.R2 cells have a gene signature related to cell lines with mutated p53 (gene set: M2697, Figure 3A). Targeted NGS further confirmed the presence of two heterogenous missense mutations, Y236N and R248W, in the A549.R2 cells [Figure 3B, Supplementary Table 3]. No mutations were detected in the parental A549 cell line, and similar mutations were detected in the other monoclonal subclones A549.R1 (Y236N), A549.R3 (Y236N, R248W) and A549.R4 (Y236N) [Supplementary Table 3]. At the protein level, p53 was overexpressed in the A549.R2 cell line, which is often observed in mutant p53 cancer cells (like NCI-H1975) due to the disturbance of the MDM2 regulated negative feedback loop [Figure 3C and D]. Finally, we determined if the presence of these TP53 mutations could sensitize A549.R2 cancer cells to APR-246, a reactivator of mutant p53 and/or inducer of oxidative stress. However, we did not observe any difference in cell survival [Figure 3E] or cell death [Figure 3F] between the wild type A549 and mutant A549.R2 cell lines.
Figure 3. Acquired TP53 mutations. Genes downregulated in NCI-60 cell line panel with mutated p53 (gene set: M2697) (A); mutations in the TP53 gene detected in the A549.R2 cell line (B); Western Blot analysis of p53 protein expression (C) and fluorescent intensity normalized to b-actin (D); percentage survival (E) and percentage cell death (F) of A549 and A549.R2 cells treated with APR-246 (0-50 µmol/L) for 96 h determined by the IncuCyte ZOOM cytotoxicity assay. Data presented as mean +/- SD of at least 3 biological replicates
A549.R2 cells have an increased enrichment of an EMT-like expression signature (gene set: M5930, Figure 4A). EMT is potentially mediated by activation of LEF1, a transcription factor which regulates transcription of hallmark EMT effectors, since A549.R2 cells are enriched for genes downregulated in LEF1 overexpressing cells (Figure 4B; gene set: M2903). Consistent with a mesenchymal-like phenotype, A549.R2 cells have an increased migratory and invasive potential compared to the parental A549 cells [Figure 4C-E].
Figure 4. Epithelial to mesenchymal transition (EMT). Genes defining EMT (gene set M5930) (A); genes down-regulated in DLD1 cells (colon carcinoma) over-expressing LEF1 (gene set: M2903) (B); representative images (24 h) of a scratch wound migration (0 mg/mL Matrigel) and invasion (6 mg/mL Matrigel) assay of the nutlin-3 resistant A549.R2 cell line and parental A549 cell line (C); real-time migration (D) and invasion (E). RWD: relative wound density presented as mean +/- SD of 3 biological replicates
A549.R2 cells are enriched in genes related to NF-kB signalling (gene set WP3617, Supplementary Figure 2), IL-17 signalling (gene set hsa04657, Supplementary Figure 3), TNF signalling (gene set hsa04668, Supplementary Figure 4) and cytokines and inflammatory response (gene set WP530, Supplementary Figure 5). The affected genes are mostly related to inflammation, tissue remodelling, angiogenesis, and consist mainly of chemo- and cytokines [Supplementary Figures 2-5]. When focusing on an independent chemo- and cytokine gene set, we also observed downregulated genes as presented in Figure 5A. In addition, we observed an increased expression of several immune checkpoints including PD-L1 (CD274), PD-L2 (PDCD1LG2), CD73 (NT5E) and galectin-3 (LGALS3) as presented in Figure 5B.
Figure 5. Volcano plot of differentially expressed genes in nutlin-3 resistant A549.R2 cells versus parental A549 cells. Chemokine and cytokine gene set derived from the Nanostring nCounter® PanCancer Immune Profiling Panel (A); immune checkpoints gene set based on literature (B). Log2 fold change cut-off ≤ -0.5 / ≥ 0.5, P-value cut-off ≤ 0.05
A549.R2 cells showed a reduced expression of genes defining KRAS dependency (gene set M2851; Supplementary Figure 6B) or related to overexpression of oncogenic KRAS (gene set M2892; Supplementary Figure 6B) compared to the parental A549 cancer cells (KRAS mutation p.G12S c.34G>A). Next, A549.R2 cells are enriched in gene sets related to knockdown of BMI-1 (gene sets: M2730, M2775, M2782, M2779; Supplementary Figure 6A), an oncogene related to chemoresistance, including cisplatin[17,18]. Therefore, we determined the sensitivity of A549.R2 and A549 cancer cells for cisplatin. Only a small difference in survival was observed (IC50: 1.07 µmol/L vs. 1.69 µmol/L), but the induction of cell death was strongly reduced in the A549.R2 cells (EC50: 10.88 µmol/L vs. 21.93 µmol/L) [Figure 6A and B]. Finally, A549.R2 cells showed gene expression characteristics related to NFE2L2 activation (gene set: M2903; Supplementary Figure 6C), PCGF2 knockdown (gene sets: M2775, M2783, M2784; Supplementary Figure 6A), activation of TGF-b receptor signalling (gene set WP560, Supplementary Figure 7), SUZ12 knockdown (gene set: M2740; Supplementary Figure 6D) and MEK signalling (gene set: M2725; Supplementary Figure 6E).
Figure 6. Cisplatin-induced cell death is reduced in A549.R2 cells. Percentage survival (A) and percentage cell death (B) of A549 and A549.R2 cells treated with cisplatin (0-50 µmol/L) for 96 hours determined by the IncuCyte ZOOM cytotoxicity assay. Data presented as mean +/- SD of at least 3 biological replicates
In this study, we successfully generated nutlin-3 resistant A549 NSCLC cancer cells to study the characteristics of acquired resistance. Consistent with previous studies[19,20], we observed that nutlin-3 selects for mutant p53 cancer cells since all subclones harboured at least one TP53 mutation. Michaelis et al. suggest that nutlin-3 induces de novo p53 mutations due to the high variability in mutations they detected in their sublines. Similarly, Aziz et al. report that four different p53 mutations were obtained in four separate experiments in which nutlin-3 resistant SJSA-1 cell were generated. Although we only observed two different mutations (R248W and Y236N), the fact that we obtained two subclones with a double heterogenous mutation also suggest that these mutations are de novo mutations that can accumulate rather than a selection of pre-existing mutant p53 A549 cells since all clones have the same Y236N mutation and two clones acquired an additional R248W mutation. This is further supported by the fact that no mutations were detected in the parental A549 cells following deep-sequencing analysis. Finally, Skalniak et al. report similar findings for Idasanutlin which generates p53-mutated drug-resistant subpopulations.
Despite the presence of p53 mutations in all the resistant subclones, a strong variability in nutlin-3 sensitivity was observed. Prolonged exposure to increasing concentrations of nutlin-3 would further select for the most resistant clones. Therefore, we selected the most resistant subclone A549.R2 for in-depth characterization by whole-transcriptome analysis to provide novel insights into the characteristics of nutlin-3 resistant cancer cells besides acquisition of p53 mutations and multidrug resistance as previously described[20-22]. Although our results cannot be generalized to all nutlin-3 resistant cells, our findings highlight interesting adaptations of cancer cells following continuous exposure to nutlin-3 that could be considered in other studies.
Acquiring a TP53 mutation did not sensitize A549.R2 cells for treatment with APR-246 (PRIMA-1Met), a first in class reactivator of mutant p53. It is important to note that increasing evidence suggests that the main mechanism of action of this compound is the induction of oxidative stress through its influence on thioredoxin and glutathione antioxidant systems and that restoring wild type p53 conformation is limited to specific mutations[24,25]. Nevertheless, this also suggests that A549.R2 cells are not sensitized to antioxidant inhibitors. In line with this observation, genes upregulated in cells with knockout of NFE2L2 were downregulated in the A549.R2 cells showing evidence of increased activity of the Nrf2 transcription factor, a master regulator of the cellular antioxidant response.
Several gene sets related to BMI-1 knockdown were enriched in the A549.R2 cells. BMI-1 plays a vital role in cancer cell proliferation, invasion/metastasis, chemo-sensitivity and patient survival. Wang et al. found that BMI-1 gene silencing can reduce intracellular glutathione levels in ovarian cancer cells, thus decreasing their resistance to the chemotherapeutic drug cisplatin[17,18]. However, we observed that A549.R2 cells were more resistant to cisplatin treatment, manifested as a reduced induction of cell death, consistent with the findings of Michaelis et al. in acquired nutlin-3 resistant neuroblastoma cells. Cisplatin acts by inducing DNA damage, thereby triggering the p53-regulated DNA-damage response. Our group recently showed that nutlin-3 synergistically enhances the cytotoxic effects of cisplatin through stimulation of wild type p53 leading to increased cell cycle arrest and apoptosis. Consequently, loss of wild type p53 functions due to acquired mutations can disturb this response leading to increased resistance.
Further related to NSCLC standard of care therapy, A549.R2 cells acquired several characteristics that can modulate the tumour microenvironment, which could impact the response to standard of care immunotherapies like anti-PD-1/PD-L1 blockers. Firstly, gene expression of several matrix metalloproteinases (MMPs) was upregulated as indicated by pathway analysis of enriched gene sets WP3617 and hsa04657. These MMPs play a role in tissue remodelling and neo-angiogenesis which could affect treatment response to immunotherapy since combination strategies of immunotherapy and antiangiogenic compounds show promise for synergistic treatment of advanced NSCLC[27-29]. Secondly, the expression of a variety of chemo- and cytokines was affected in the A549.R2 cells which includes genes related to leukocyte recruitment (CXCL1-3) and activation (CSF1-2) and Th2 response (IL-5, IL-13) as indicated by the pathway analysis of gene sets hsa04668 and hsa04675, respectively. IL-32 was the most significantly upregulated cytokine, which has been shown to play a crucial role in tumour initiation, proliferation and maintenance. Thirdly, we observed an increased expression of several immune checkpoints, including PD-L1 (CD274), PD-L2 (PDCD1LG2), CD73 (NT5E), and galectin-3 (LGALS3) in the nutlin-3 resistant cancer cells.
Finally, whole transcriptome analysis showed that nutlin-3 resistant cells acquire EMT-like properties. EMT-related transcription factors (Snail, Slug and ZEB1/2) suppress the expression of epithelial markers such as E-cadherin, while increasing the expression of mesenchymal markers N-cadherin and fibronectin. An important characteristic of these mesenchymal-like cells is their increased motility which was confirmed by the increased migration and invasion of A549.R2 cells. LEF1 is another transcription factor with an important role in EMT by activating transcription of hallmark EMT effectors such as N-cadherin, Vimentin and Snail. A549.R2 were enriched in genes related to overexpression of LEF1, indicating that EMT is induced through induction of LEF1, although this needs to be confirmed at the protein and functional level.
In conclusion, transcriptome analysis of a highly nutlin-3 resistant A549 subclone shows the relevance of studying (1) resistance to standard of care chemotherapy; (2) secretion of immunomodulating chemo- and cytokines; (3) immune checkpoint expression; and (4) EMT and invasion in nutlin-3 resistant cancer cells in addition to acquired mutations in the TP53 gene.
Willy Floren (Wijnegem, Antwerp) is greatly acknowledged for his private funding of some of the equipment used in this study.Authors’ contributions
Made substantial contributions to conception and design of the study and performed data analysis and interpretation: Deben C, Wouters A, Lardon F, Pauwels P
Performed data acquisition, as well as provided administrative, technical, and material support: Deben C, Boullosa LF, Domen A
Bioinformatic analysis: Deben C, Cuypers B, Laukens KAvailability of data and materials
RNA-sequencing data will be deposited into data repositories upon acceptance of the manuscript.Financial support and sponsorship
Boullosa LF was supported by the Research Foundation Flanders (11E7719N). Domen A was funded by the University Research Fund of the University of Antwerp.Conflicts of interest
All authors declared that there are no conflicts of interest.Ethical approval and consent to participate
Not applicable.Consent for publication
© The Author(s) 2020.
1. Deben C, Van den Bossche J, Van Der Steen N, et al. Deep sequencing of the TP53 gene reveals a potential risk allele for non-small cell lung cancer and supports the negative prognostic value of TP53 variants. Tumour Biol 2017;39:1010428317694327.DOIPubMed
2. Deben C, Op de Beeck K, Van den Bossche J, et al. MDM2 SNP309 and SNP285 act as negative prognostic markers for non-small cell lung cancer adenocarcinoma patients. J Cancer 2017;8:2154-62.DOIPubMedPMC
3. Deben C, Deschoolmeester V, Lardon F, Rolfo C, Pauwels P. TP53 and MDM2 genetic alterations in non-small cell lung cancer: evaluating their prognostic and predictive value. Crit Rev Oncol Hematol 2016;99:63-73.DOIPubMed
4. Vassilev LT. Small-molecule antagonists of p53-MDM2 binding: research tools and potential therapeutics. Cell Cycle 2004;3:419-21.PubMed
5. Deben C, Wouters A, Op de Beeck K, et al. The MDM2-inhibitor Nutlin-3 synergizes with cisplatin to induce p53 dependent tumor cell apoptosis in non-small cell lung cancer. Oncotarget 2015;6:22666-79.DOIPubMedPMC
6. Tisato V, Voltan R, Gonelli A, Secchiero P, Zauli G. MDM2/X inhibitors under clinical evaluation: perspectives for the management of hematological malignancies and pediatric cancer. J Hematol Oncol 2017;10:133.DOIPubMedPMC
7. Pauwels B, Korst AE, de Pooter CM, et al. Comparison of the sulforhodamine B assay and the clonogenic assay for in vitro chemoradiation studies. Cancer Chemother Pharmacol 2003;51:221-6.DOIPubMed
8. Dobin A, Davis CA, Schlesinger F, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013;29:15-21.DOIPubMedPMC
9. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 2014;15:550.DOIPubMedPMC
10. Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res 2019;47:W199-205.DOIPubMedPMC
11. Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The molecular signatures database (MSigDB) hallmark gene set collection. Cell Syst 2015;1:417-25.DOIPubMedPMC
12. Subramanian A, Tamayo P, Mootha VK, et al. 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.DOIPubMedPMC
13. Vandeweyer G, Van Laer L, Loeys B, Van den Bulcke T, Kooy RF. VariantDB: a flexible annotation and filtering portal for next generation sequencing data. Genome Med 2014;6:74.DOIPubMedPMC
14. Soussi T, Hamroun D, Hjortsberg L, Rubio-Nevado JM, Fournier JL, Beroud C. MUT-TP53 2.0: a novel versatile matrix for statistical analysis of TP53 mutations in human cancer. Hum Mutat 2010;31:1020-5.DOIPubMed
15. Pauwels B, Korst AE, Andriessen V, et al. Unraveling the mechanism of radiosensitization by gemcitabine: the role of TP53. Radiat Res 2005;164:642-50.DOIPubMed
16. Wang MC, Li CL, Cui J, et al. BMI-1, a promising therapeutic target for human cancer. Oncol Lett 2015;10:583-8.DOIPubMedPMC
17. Wu Z, Min L, Chen D, et al. Overexpression of BMI-1 promotes cell growth and resistance to cisplatin treatment in osteosarcoma. PLoS One 2011;6:e14648.DOIPubMedPMC
18. Wang E, Bhattacharyya S, Szabolcs A, et al. Enhancing chemotherapy response with Bmi-1 silencing in ovarian cancer. PLoS One 2011;6:e17918.DOIPubMedPMC
19. Kucab JE, Hollstein M, Arlt VM, Phillips DH. Nutlin-3a selects for cells harbouring TP53 mutations. Int J Cancer 2017;140:877-87.DOIPubMedPMC
20. Michaelis M, Rothweiler F, Barth S, et al. Adaptation of cancer cells from different entities to the MDM2 inhibitor nutlin-3 results in the emergence of p53-mutated multi-drug-resistant cancer cells. Cell Death Dis 2011;2:e243.DOIPubMedPMC
21. Aziz MH, Shen H, Maki CG. Acquisition of p53 mutations in response to the non-genotoxic p53 activator Nutlin-3. Oncogene 2011;30:4678-86.DOIPubMedPMC
22. Skalniak L, Kocik J, Polak J, et al. Prolonged idasanutlin (RG7388) treatment leads to the generation of p53-mutated cells. Cancers (Basel) 2018;10.DOIPubMedPMC
23. Lambert JM, Gorzov P, Veprintsev DB, et al. PRIMA-1 reactivates mutant p53 by covalent binding to the core domain. Cancer Cell 2009;15:376-88.DOIPubMed
24. Zhang Q, Bykov VJN, Wiman KG, Zawacka-Pankau J. APR-246 reactivates mutant p53 by targeting cysteines 124 and 277. Cell Death Dis 2018;9:439.DOIPubMedPMC
25. Liu DS, Duong CP, Haupt S, et al. Inhibiting the system xC-/glutathione axis selectively targets cancers with mutant-p53 accumulation. Nat Commun 2017;8:14844.DOIPubMedPMC
26. Bai X, Chen Y, Hou X, Huang M, Jin J. Emerging role of NRF2 in chemoresistance by regulating drug-metabolizing enzymes and efflux transporters. Drug Metab Rev 2016;48:541-67.DOIPubMed
27. Stetler-Stevenson WG. Matrix metalloproteinases in angiogenesis: a moving target for therapeutic intervention. J Clin Invest 1999;103:1237-41.DOIPubMedPMC
28. Page-McCaw A, Ewald AJ, Werb Z. Matrix metalloproteinases and the regulation of tissue remodelling. Nat Rev Mol Cell Biol 2007;8:221-33.DOIPubMedPMC
29. Manegold C, Dingemans AC, Gray JE, et al. The potential of combined immunotherapy and antiangiogenesis for the synergistic treatment of advanced NSCLC. J Thorac Oncol 2017;12:194-207.DOIPubMed
30. Yan H, He D, Huang X, et al. Role of interleukin-32 in cancer biology. Oncol Lett 2018;16:41-7.DOIPubMedPMC
31. Kobayashi W, Ozawa M. The epithelial-mesenchymal transition induced by transcription factor LEF-1 is independent of β-catenin. Biochem Biophys Rep 2018;15:13-8.DOIPubMedPMC
32. Son H, Moon A. Epithelial-mesenchymal transition and cell invasion. Toxicol Res 2010;26:245-52.DOIPubMedPMC
33. Santiago L, Daniels G, Wang D, Deng FM, Lee P. Wnt signaling pathway protein LEF1 in cancer, as a biomarker for prognosis and a target for treatment. Am J Cancer Res 2017;7:1389-406.PubMedPMC
Deben C, Boullosa LF, Domen A, Wouters A, Cuypers B, Laukens K, Lardon F, Pauwels P. Characterization of acquired nutlin-3 resistant non-small cell lung cancer cells. Cancer Drug Resist 2020;3:[Online First]. http://dx.doi.org/10.20517/cdr.2020.91
et al., Mini-invasive Surgery, 2021
et al., Mini-invasive Surgery, 2020
et al., Journal of Cancer Metastasis and Treatment, 2020
et al., Mini-invasive Surgery, 2020
et al., Mini-invasive Surgery, 2020
et al., Journal of Cancer Metastasis and Treatment, 2020
et al., Cancer Drug Resistance, 2020
et al., Mini-invasive Surgery, 2020
et al., Cancer Drug Resistance, 2020
et al., Mini-invasive Surgery, 2020