- 1Stem Cell Institute, Department of Development and Regeneration, Katholieke Universiteit (KU) Leuven, Leuven, Belgium
- 2Leuven Cancer Institute (LKI), Department of Oncology, Gynecological Oncology Lab 3000, KU Leuven, Leuven, Belgium
- 3Stem Cell Institute, Department of Development and Regeneration, Laboratory for Stem Cell and Developmental Epigenetics, KU Leuven, Leuven, Belgium
- 4Translational Cell and Tissue Research – Department of Imaging & Pathology, KU Leuven, Leuven, Belgium
- 5Leuven Kanker Instituut (LKI) Leuven Cancer Institute - Trace, KU Leuven, Leuven, Belgium
- 6Centre for Gynecologic Oncology Amsterdam (CGOA), Antoni Van Leeuwenhoek-Netherlands Cancer Institute (AvL-NKI), University Medical Center (UMC), Amsterdam, Netherlands
- 7Division of Oncogenomics, Oncode Institute, The Netherlands Cancer Institute, Amsterdam, Netherlands
Triple-Negative Breast Cancer (TNBC) is the most aggressive breast cancer subtype, characterized by limited treatment options and higher relapse rates than hormone-receptor-positive breast cancers. Chemotherapy remains the mainstay treatment for TNBC, and platinum salts have been explored as a therapeutic alternative in neo-adjuvant and metastatic settings. However, primary and acquired resistance to chemotherapy in general and platinum-based regimens specifically strongly hampers TNBC management. In this study, we used carboplatin-resistant in vivo patient-derived xenograft and isogenic TNBC cell-line models and detected enhanced Wnt/β-catenin activity correlating with an induced expression of stem cell markers in both resistant models. In accordance, the activation of canonical Wnt signaling in parental TNBC cell lines increases stem cell markers’ expression, formation of tumorspheres and promotes carboplatin resistance. Finally, we prove that Wnt signaling inhibition resensitizes resistant models to carboplatin both in vitro and in vivo, suggesting the synergistic use of Wnt inhibitors and carboplatin as a therapeutic option in TNBC. Here we provide evidence for a prominent role of Wnt signaling in mediating resistance to carboplatin, and we establish that combinatorial targeting of Wnt signaling overcomes carboplatin resistance enhancing chemotherapeutic drug efficacy.
Introduction
Triple-negative breast cancer (TNBC) is a molecular subtype of breast cancer characterized by the lack of expression of estrogen-receptor, progesterone-receptor, and human epidermal growth factor receptor type 2 (1). TNBC accounts for 10-20% of all breast cancer cases, occurring with higher frequency in younger women, presenting with higher grade and mitotic counts than non-TNBCs, low differentiation, and frequent lymph node involvement, ultimately contributing to poor prognosis (1, 2).
The lack of hormone- and growth factor-receptors render chemotherapy the primary systemic treatment for TNBC. Interestingly, TNBC patients have high response rates to neoadjuvant chemotherapy, achieving pathological complete response (pCR) more frequently than those bearing non-TNBCs (3). Nonetheless, TNBC patients experience lower progression-free- and overall survival rates and higher distant metastatic relapse frequency than non-TNBC patients, highlighting the critical need for alternative therapeutic approaches (3).
The Food and Drugs Administration of the United States of America first approved platinum salts, namely cisplatin, to treat metastatic testicular cancer, ovarian cancer, and bladder cancer between 1978 and 1979 (4). Since then, the use of platinum-based chemotherapy has grown and is now applied in many other cancer types. Preclinical studies have highlighted that TNBC is particularly sensitive to DNA damaging agents (5, 6). For that reason, platinum salts – DNA-crosslinking agents – have gained traction as potential additions to the therapeutic toolbox for TNBC. Phase-II and Phase-III clinical trials have demonstrated the benefits of including carboplatin (CAR) in neoadjuvant regimens for TNBC (7–10). Importantly, pCR with neoadjuvant treatment is a robust predictor of survival in TNBC (3). However, chemotherapy-treated TNBC patients are likely to acquire resistance, and patients with residual disease (RD) have worsened prognosis and experience low survival rates, particularly within the first three years ensuing treatment (3).
During the last decades, our knowledge of platinum’s mechanism of action has increased significantly. Nonetheless, how cancer overcomes platinum-mediated cytotoxicity still holds unanswered questions. Several studies have shed light on how cancer cells adapt to platinum-based treatment by restoring DNA damage repair, increasing tolerance to DNA damage, decreasing its intracellular uptake and accumulation, and regulating apoptosis and autophagy (11). Chemotherapy resistance is also known to be induced and maintained by adaptations in pro-survival and anti-apoptotic signaling pathways. Like other chemotherapeutic agents, alterations of such cellular dynamics also affect platinum-based treatments. Several studies have demonstrated the involvement of NOTCH (12), MEK (13), Hedgehog (14), EGFR (15), among others, in mediating resistance to platinum in different cancer types. Also, cancer cells with stem cell-like properties have been described to significantly influence the response to different chemotherapeutic agents, including platinum compounds (16–19). Cancer stem cells constitute a subpopulation of cancer cells with tumorigenic and self-renewal capacities and are considered desirable therapeutic targets since their intrinsic cellular properties contribute extensively to treatment failure (20–22). Breast cancer stem cells were first isolated in 2003 based on cell surface markers CD44 and CD24 (23). Since then, many studies have demonstrated their tumorigenic and drug-resistance capacities, highlighting the need to develop therapeutic approaches that deplete this population (6, 24–26).
The Wnt/β-catenin signaling pathway is a developmental signaling cascade with a prominent role in cancer (27). It is initiated when Wnt ligands (secreted lipid-modified signaling molecules) bind the receptor complex at the cell membrane. A series of events ensues, culminating in the inhibition of glycogen synthase kinase 3 beta (GSK3β), and the subsequent cytoplasmic accumulation of β-catenin, the critical mediator protein of Wnt signaling. This accumulation leads to the nuclear translocation of β-catenin, eliciting Wnt target genes’ expression by interacting with different transcription factors. In the absence of Wnt ligands, β-catenin is constitutively phosphorylated by GSK3β and targeted for proteasomal degradation (28). Notably, Wnt is known to govern several cellular functions with the potential to contribute to chemotherapy resistance. Such functions include the control and regulation of proliferation (29), DNA damage repair (30), inhibition of apoptosis (31), and maintenance and regulation of embryonic, somatic, and cancer stem cell properties (32). Different studies have demonstrated Wnt pathway involvement in the mediation of platinum resistance in various cancer types, including squamous cell carcinoma (33, 34) and ovarian cancer (35). However, its involvement in platinum resistance in TNBC is not known.
To study how TNBCs acquire resistance to carboplatin treatment, we used a pair of isogenic carboplatin-sensitive and -resistant TNBC cell lines. Transcriptomic analysis was performed to gain insight into biological signaling pathways underlying acquired carboplatin resistance in vitro, leading to the identification of Wnt signaling as a candidate resistant-mediating pathway. Additionally, the resistant cell line displayed enhanced expression of pluripotency markers and stem cell features compared to the parental, carboplatin-sensitive cells.
In vitro pharmacological and genetic manipulation of Wnt signaling was employed to assess drug response alterations and stem cell potential functionally. Inducing Wnt signaling in parental non-resistant TNBC cell lines elicited the expression of pluripotency markers observed in isogenic resistant cells and enhanced stem cell features in vitro. Moreover, inhibition of Wnt activity in resistant cells resensitized them to carboplatin and hindered tumorsphere formation. Finally, carboplatin-resistant patient-derived xenograft (PDX) models were used to test the effect of in vivo Wnt inhibition on platinum-response. Similar to what we observed in vitro, inhibition of Wnt reduced expression of cancer stem cell markers and drastically reduced resistance to carboplatin treatment in vivo.
Altogether, our results suggest the potential for Wnt signaling inhibition in combination with carboplatin as a strategy to prevent or overcome platinum resistance in TNBC patients.
Materials and Methods
Cell Lines, Cell Culture, and Treatments
MDA-MB-468 (ATCC-HTB-132) and MDA-MB-231 (ATCC HTB-26) were maintained in DMEM high glucose (Gibco 41965039) supplemented with 10% fetal bovine serum, 1mM sodium pyruvate (Gibco, 11360070), 1X non-essential amino acids (Gibco, 11140035), 100 μg/mL penicillin-streptomycin (Gibco, 15140163) and 0.01mM 2-mercaptoethanol (Gibco, 31350010).
Unless otherwise specified in the text, all carboplatin treatments (CAR, Hospira UK, Ltd) were done at 2 μM for MDA-MB-468 cells and 35 μM for MDA-MB-231 cells. Small molecule Wnt activator CHIR99021 (CHIR, Sigma, SML1046) was administered to cells in DMSO (Sigma, D2650) at 4 μM. Small molecule Wnt inhibitor LGK-974 (LGK, Selleckchem, S7143) was, unless otherwise specified, administered at 200 nM in DMSO.
Carboplatin resistance was induced in MDA-MB-468 cells by continued maintenance in carboplatin containing medium, starting at a concentration of 0.4 μM. The concentration of carboplatin was increased in 9 increments until reaching 2 μM, once unhindered cell growth was obtained at each concentration level, allowing a 48h carboplatin-free recovery period with each splitting (35).
Lentiviral Particle Production and Transduction
Lentiviruses were produced according to the RNAi Consortium (TRC) protocol available from the Broad Institute (https://portals.broadinstitute.org/gpp/public/resources/protocols). In brief, 5x105 HEK293T cells seeded per well in 6-well plates and transfected the following day with 750 μg pCMV-dR8.91, 250 μg pCMV-VSV-G, and 1 μg of the specific lentiviral expression or silencing constructs using FugeneHD (Promega, E2311) in Opti-MEM (Gibco, 31985070). One day after, the culture medium was replaced. The same day, lentivirus-recipient cells were plated in 6-well plates at a density of 5x104 cells per well. Lentivirus-containing medium was collected from HEK293T cells 48h and 72h after transfection and added to recipient cancer cells after being filtered. Two days after infection, cells were washed thoroughly with PBS, medium refreshed, and appropriate selection antibiotics applied.
For overexpression of ΔN90-β-catenin, we used pLV-beta-catenin ΔN90 (Addgene, #36985) and pPRIME-CMV-NEO-recipient (CTRL, Addgene, #11659). For β-catenin shRNA mediated silencing, we used pXL002-ishRNA-beta-catenin-1 (Addgene, #36297) and pXL004-ishRNA-scramble (Addgene, #36311). For Wnt fluorescent reporter assay, we used 7TGP (Addgene, #24305).
In Vitro Carboplatin Response
For IC50 experiments, we seeded 2.5x104 cells per well in 12-well plates. Cells were treated with increasing concentrations of CAR (0.02 to 200 μM for MDA-MB-468 and 0.35 to 3500 μM for MDA-MB-231) for 72 hours. Viability was assessed by manual cell counting using a Neubauer hemocytometer using trypan blue for dead cell exclusion. Cell viability was determined as a percentage of untreated cells, and non-linear regressions of [CAR] vs. normalized-response were fitted using GraphPad Prism v.8.0.1. to mathematically determine the IC50.
Flow Cytometry
For annexin V apoptosis analysis, cells were detached and resuspended in annexin V binding buffer (BD Pharmigen, 51-66121E) and incubated for 15 minutes at room temperature with APC-conjugated AnnexinV (Thermo-eBioscience, BMS306APC-100). After incubation, cells were diluted in binding buffer containing 100 nM of 4’,6-diamidino-2-phenylindole (DAPI). Unstained and single-stained (annexin V and DAPI) were used as gating controls.
For ALDH activity assays, cells were detached, washed in PBS, and stained using the AldeRed ALDH detection assay kit (Merck SRC150) according to manufacturer specifications.
For immunolabeling of CD44 and CD24, cells were detached, washed twice in PBS with 4% FBS, and incubated with CD44-PE (BD Pharmigen, 555479) and CD24-APC (Invitrogen, 17-4714-81) antibodies according to manufacturer specifications at room temperature. After incubations, cells were washed twice in PBS with FBS and resuspended in PBS containing 4% FBS and 100 nM of DAPI. Cells incubated with PE- and APC- conjugated isotype-antibodies and single-stained cells were used as gating controls.
All data were collected on a BD FACS Canto II at the KU Leuven Flow Cytometry Core and analyzed using FlowJo v.10.6.2.
SDS-PAGE and Western Blot
For western blot, cells were collected and washed in PBS before being pelleted. Then, cells were lysed on ice with RIPA buffer (150 mM NaCl, 1% Nonidet P40, 0.5% sodium deoxycholate, 0.1% dodecyl sulfate, 50 mM Tris-HCL, pH 8.0) containing a cocktail of protease and phosphatase inhibitors (Sigma, #P5726, #P0044, #P8340). Lysates were centrifuged at 16,000x g for 10 minutes at 4°C to discard insoluble material, and protein concentration was determined using the Bradford method. For SDS-PAGE, 30 μg of protein were mixed with 4x Laemmli buffer (240 mM Tris/HCL pH 6.8, 8% SDS, 0.04% bromophenol blue, 5% 2-mercaptoethanol, 40% glycerol) and denatured for 5 minutes at 96°C prior to electrophoretic protein separation. Resolved protein extracts were transferred to PVDF membranes (BIORAD, 162-0177). Transfer success was assessed with Ponceau S solution, and membranes were blocked with 5% non-fat milk or 5% BSA in TBS-T (0.1% Tween-20®) for 60 minutes. After blocking, membranes were incubated with primary antibodies at 4°C overnight. The day after, membranes were washed 3 times with PBS-T for 10 minutes and incubated with secondary HRP-conjugated antibodies. Immunolabeled proteins were detected with Supersignal West Pico chemiluminescent kit (Fisher Scientific, 34077) on autoradiography film (Santa Cruz, SC-201697). The primary antibodies used were active rabbit anti-non-phosphorylated β-catenin (CellSignaling Technologies, #19807S), mouse anti-total β-catenin (BD, #610154), mouse anti-β-actin (Santa Cruz, #47778).
Next-Generation mRNA Sequencing
Total RNA was obtained from cells using the GenElute mammalian total RNA miniprep kit (Sigma, RTN350-1KT) for mRNA sequencing. According to the manufacturer’s specifications, libraries were prepared from 250 ng of total RNA using the KAPA stranded mRNA-seq kit (Roche, KK8421). KAPA-single index adapters (Roche, KK8700) were added to A-tailed cDNA, and libraries were amplified by 12 cycles. Finally, libraries were purified on Agencourt AMPure XP beads (Beckman Coulter, A63881). Libraries were controlled for fragment size using the High Sensitivity DNA analysis kit (Agilent, 5067-4626) on an Agilent Bioanalyzer 2100. Each library was diluted to 4 nM and pooled for single-end 50-bp sequencing on an Illumina Hiseq4000 20 – 27 million reads per sample (22 million reads on average).
Adapters, polyA tails, and bad quality reads (Phred score > 20) were trimmed using Trim Galore! (v0.6.4_dev) with default parameters. Reads were aligned to the transcriptome and quantified using Salmon (v0.14.1) (36) with default parameters using GENCODE release 36 of the human reference transcriptome sequences and the comprehensive gene annotation. Subsequently, the counts were imported into R (v4.0.2) using tximport (v1.18.0) and differentially expressed genes were defined using DEseq2 (v1.30.0) (37) and log fold changes corrected using “ashr” method (38) (FDR adjusted p.val < 0.05 & |log2(fold change)| > 1.5). TPM values were also calculated using tximport.
Functional Enrichment Analysis and Enrichment Maps
Datasets GSE103668 and E-MTAB-7083 were downloaded from the GeneExpression Omnibus and ArrayExpress public repositories, respectively. Differentially expressed genes with |log2(fold-change)| >1 and p-value < 0.05 were obtained using limma (v3.26.8) R package in R (v4.02) and by using the limma method on NetworkAnalyst (39). Differentially expressed genes were ranked by fold-change for Gene Set Enrichment Analysis (GSEA v4.1.0) using weighted enrichment statistic and KEGG, Hallmarks, and Wikipathways gene sets. Additionally, we used custom gene sets comprised of human embryonic stem cell-related genes (M1871: BENPORATH_ES1 and M4241: BENPORATH_ES2), pluripotency transcription factor target genes (M14573: BENPORATH_NOS_TARGETS), and cancer progenitor genes (ENGELMANN_CANCER_PROGENITORS_UP) obtained from www.gsea-msigdb.org. The statistical significance threshold was set at FDR<0.1 or (p<0.5 ∧ FDR<0.25). Additionally, gProfiler (https://biit.cs.ut.ee/gprofiler/gost) was used to assess the function of ranked DEGs using the ranked query mode and Benjamini-Hochberg FDR thresholding. GSEA and gProfiler analysis outputs were fed to the EnrichmentMap app on Cytoscape (v3.8.1) to generate visualizations of enriched biological features and pathways following published protocols (40). Differentially expressed genes from RNA-sequencing were processed for functional analysis and visualization in the same way, except for GSEA, differentially expressed genes were ranked by the absolute value of fold change (41).
Real-Time Quantitative Polymerase Chain Reaction and Gene Expression Analysis
For RT-qPCR, total RNA was extracted using the GenElute mammalian total RNA miniprep kit from Sigma (Sigma, RTN350-1KT) according to the manufacturer’s instructions, and DNA was digested during RNA extraction using on-column DNAse (Sigma, On-Column DNAse I digestion set, DNASE70). cDNA was synthesized from 500 ng of total RNA using the BIORAD iScript cDNA synthesis kit (BIORAD, CAT#1708891). Quantitative real-time PCR reactions were set up in technical triplicates with Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen, 11733-046) on a ViiA7 Real-Time PCR System (Thermo Scientific). Expression levels were normalized to housekeeping genes (HKG) GAPDH and RPL19. Statistical testing of differences in expression between samples was carried out on relative-expression values (2-ΔCT). In some figures, mRNA expression values are represented as fold-change for convenience of interpretation, although statistical testing was performed on relative expression values (2-ΔCT).
Tumorsphere Formation Assays
For tumorsphere formation assays, cells were collected as described above, washed, counted, and resuspended in serum-free tumorsphere assay medium containing DMEM/F12, 1x B27 (Thermo, 12587010), 10ng/mL bFGF, (Peprotech, 100-18b) 20 ng/mL EGF (Peprotech, AF-100-15), and 2% growth-factor reduced matrigel (Corning, 734-0268). Cells were seeded at a density of 1000 cells/mL in ultra-low attachment 6-well plates and allowed seven days to grow. On the 7th day, spheres were collected and centrifuged at 50g for 10 minutes, resuspended, and transferred to 96-well plates. Plates were briefly centrifuged at 50g for 1 minute to pull down larger spheroids (>60 μm) which were counted using a tally counter.
Immunohistochemistry
Tumor samples were dissected, washed in saline, and either snap-frozen in OCT compound (VWR 361603E) or fixed in 4% formalin for 24 hours. Frozen tissue was cut at 10 μm thickness using a cryostat and mounted on superfrost microscope slides (Thermo Scientific, °J1800AMNZ). Formalin-fixed tissue was embedded in paraffin and sectioned at 4 μm thickness using a microtome. Frozen and FFPE sections were stained with hematoxylin and eosin (HE), using a Leica Autostainer XL (Leica Microsystems), or stained in immunohistochemistry (IHC). In brief, frozen sections were fixed in acetone and preserved at -80°C until use. Slides were thawed at room temperature for 15 minutes and rehydrated in PBS. FFPE sections were deparaffinized in Leica Autostainer XL (Leica Microsystems) and pre-treated in citrate buffer (EnVision FLEX Target Retrieval Solution Low pH, Agilent-Dako, K8005) using a PT Link module (Dako), according to manufacturer’s instructions. For IHC, tumor sections were incubated with Envision Flex Peroxidase-Blocking Reagent (Dako, S202386-2) for at least 5 minutes, rinsed three times in wash buffer (Dako, K800721-2), and blocked with 5% bovine serum albumin (BSA) for 45 min at room temperature. After overnight incubation with anti-human Ki67 antibody (Abcam, EPR3610 clone, 1:1500) at 4°C, slides were incubated with HPR-conjugated secondary antibodies (Agilent, K400311-2) for 30 min. Tissue sections were stained with 3,3’-diaminobenzidine solution (DAB; Liquid DAB+ Substrate Chromogen System, K346889-2, Dako), counterstained with hematoxylin, and mounted using an automated coverslipper machine (Leica CV5030, Leica Biosystems). Pictures were acquired using a Zeiss Axiovision microscope. Quantification of Ki67+ cells was performed in at least 5 random 20x fields per sample using QuPath 0.2.3 (42).
Tunel Staining and Pan-Cytokeratin Immunofluorescent Staining
For TUNEL staining, cryopreserved tumor samples were cryosectioned at 10 μM thick and mounted on superfrost microscope slides. Slides were stored at -80°C and stained using the Click-iT Plus TUNEL (ThermoFisher C10617) according to the manufacturer’s instructions. After TUNEL staining, slides were blocked with 5% normal donkey serum in PBS (Gibco, 10010-050) for one hour and incubated overnight at 4°C with rabbit anti-pan-cytokeratin polyclonal antibody (Abcam ab217916 1:400). The following day, slides were washed three times in PBS containing 0.01% Triton X-100 and incubated for 2 hours with AlexaFluor conjugated donkey anti-rabbit secondary antibody (Abcam ab150073 1:1000). After secondary antibody incubation, slides were washed three times with PBS containing 0.01% Triton X-100 and mounted with ProLong™ Gold Antifade Mountant with DAPI (P36931, Thermo Fisher Scientific). Images were acquired using a Leica Sp8x confocal microscope. TUNEL positive/pan-cytokeratin positive cells were quantified in at least eight randomly sampled 10x fields per sample using QuPath 0.2.3.
PDX Models
BRC016 (primary, grade III, TNBC) was established at the University Hospital UZ Leuven and is available from Trace, the KU Leuven/LKI PDX Platform (https://www.uzleuven-kuleuven.be/lki/trace/trace-leuven-pdx-platform). C4O was previously obtained from a carboplatin treatment-refractory BRC016 tumor. Briefly, 10 BRC016 TNBC tumor-bearing NMRI-Fox1nu nude mice (Taconic) were treated intraperitoneally with 50 mg/kg carboplatin, once weekly for three weeks. The carboplatin dose chosen was based on previous studies of carboplatin-response in our lab and other groups (35, 43–45). Nine out of ten mice had a complete response to treatment. One mouse had delayed and incomplete response, bearing a residual tumor mass that became unresponsive to further carboplatin treatment and eventually resumed growth. In a previously published study, the regrown tumor was harvested and implanted on nude mice for propagation, re-testing, and confirmation of carboplatin resistance (35). Treatment experiments included 24 mice implanted with C4O tumor fragments. When tumors reached a volume of approximately 300 mm3, mice were randomly assigned to placebo (vehicle), CAR (50 mg/kg), LGK (2 mg/kg), or CAR+LGK (50 mg/kg + 2 mg/kg) treatment groups. Carboplatin was administered once a week intraperitoneally, and LGK974 was administered daily by oral gavage. Treatments were carried out for three weeks, and body weight was closely monitored throughout treatment Tumors were measured every 48 hours with digital calipers, and volume was estimated as V = L × W2 × π/6 (L: length, W: width). All animals were euthanized at the end of three weeks of treatment, and tumors were collected for downstream analysis.
Statistical Analysis
All data were analyzed using GraphPad Prism 8, except for transcriptomic datasets. Unless otherwise specified, comparisons between two groups were tested for statistical significance using unpaired t-tests with Welch correction. Comparisons between multiple groups were performed using a one-way analysis of variance (ANOVA). Comparisons between multiple groups across multiple time points were performed using two-way ANOVA. All statistical testing was corrected for multiple comparisons, using the Holm-Sidak method when comparing samples based on experimental design or the Tukey method when testing the comparison between all means in a dataset. For the reader’s convenience, all statistical tests and sample sizes are indicated in the figure legends.
Results
Carboplatin-Resistant TNBC Cells Are Characterized by Enhanced WNT/β-Catenin Pathway Activity, Stem Cell Marker Expression, and Tumorsphere Formation Capacity
To explore the mechanisms of in vitro carboplatin resistance in TNBC, we generated an isogenic carboplatin-resistant cell line (468’R) (Figure 1A) by exposing MDA-MB-468 cells to carboplatin treatment in incremental cycles.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/705384/fonc-11-705384-HTML-r1/image_m/fonc-11-705384-g001.jpg)
Figure 1 Carboplatin-resistant TNBC cells are characterized by enhanced WNT/β-catenin pathway activity, stem cell marker expression, and tumorsphere formation capacity. (A) Phase contrast microscope images of 468’P and 468’R cells. Scale bar: 100 μm. (B) Non-linear fit model of [CAR] μM vs normalized response for IC50 determination. (468’P: n=6, R2 = 0.92) (468’R: n=4, R2 = 0.95). (C) Growth curve and statistical analysis for 468’P and 468’R cells treated with VEH or CAR using two-way ANOVA. Statistical significance is reported for day 6 of treatment (n=4). (D) Representative flow cytometry scatterplots of annexin V staining of cells treated with 2 μM CAR for 72 hours (left) and respective statistical analysis (right) using multiple t-tests corrected for multiple comparisons with the Holms-Sidak method (n=3). (E) Enriched gene sets from Hallmarks and KEGG databases by one-tailed GSEA ranked by Normalized enrichment score (NES), illustrating pathways most significantly deregulated between 468’P and 468’R. (F) Enrichment map of one-tailed GSEA hits from Wikipathways database. Rectangles highlight clusters of gene sets with significant overlap and are labeled using AutoAnnotate on Cytoscape. (G) GSEA of hESCs (26) (left) and Cancer Progenitor (46) (right) gene sets in 468’R vs 468’P cells. (H) GSEA of NANOG, OCT4, and SOX2 target genes determined by ChIP-SEQ in hESC’s (26) in 468’R vs. 468’P cells. (I) Western blot of active non-phosphorylated β-catenin in 468’P and 468’R. β-actin was used as the loading control. (J) Representative scatterplots of flow cytometric analysis of aldehyde dehydrogenase activity (left). DEAB panels refer to internal controls in which ALDH activity is blocked with diethylaminobenzaldehyde to determine the background signal generated by unconverted ALDH substrate. TEST panels refer to the experimental samples where substrate for fluorimetric determination of ALDH activity is supplied. TEST samples are normalized to background fluorescence measured in DEAB internal controls and presented as the mean + standard error of the mean percentage of ALDH+ cells in 468’R (n=5) and 468’P (n=7) (right). Welch’s t-test. (K) Representative scatterplots of flow cytometric analysis of CD44-PE and CD24-APC immunolabeling (left) and corresponding statistical analysis of the mean percentage of CD44+/CD24- cells (right; n=3). Welch’s t-test. (L) qRT-PCR of Wnt target AXIN2 and stem cell markers in 468’R cells vs. 468’P (n=4). Multiple t-tests. (M) Representative brightfield images of tumorspheres generated from 468’P and 468’R cells (left, scale bar: 50 μm) and statistical analysis of mean tumorsphere forming units (number of spheres/number of seeded single cells) (right; n=3). Welch’s t-test. (Barplots represent mean + SEM. *p < 0.05; **p < 0.01; ****p < 0.0001; ns, non significant).
Half-maximal inhibitory concentration (IC50) profiles for carboplatin were determined for both parental (468’P) and 468’R cells (Figure 1B). 468’R cells displayed a 5-fold increase in IC50, thereby functionally confirming a significant increase in carboplatin resistance. When treated with vehicle, we did not observe differences in proliferation between resistant and parental cells. However, when treated with CAR, 468’R cells displayed unhindered proliferation capacity in striking contrast to parental cells (Figure 1C). Additionally, we did not find differences in Ki67 expression in standard culture conditions between cell lines (Supplementary Figure 1A). Altogether, these results suggest that differences in proliferation rate cannot justify resistance in 468’R cells compared to the parental, sensitive cell line. Flow-cytometric analysis of apoptosis further corroborated the establishment of the carboplatin-resistant phenotype. When exposed to 2 μM for 72 hours, 468’P stained positively and significantly for the apoptosis marker annexin V (47), whereas no significant increase in apoptotic cells was observed in 468’R (Figure 1D).
To study the underlying mechanisms of carboplatin resistance, we performed transcriptome analysis by next-generation mRNA sequencing in 468’P and 468’R cells. We used the ranked differentially expressed genes (Supplementary Table 1) followed by one-tailed Gene Set Enrichment Analysis (GSEA) (41) to identify changes in signaling pathways (Hallmarks, KEGG and Wikipathways). GSEA analysis identified alterations in several key cancer-related processes, such as epithelial-mesenchymal-transition (Hallmarks) and PPAR and P53 signaling (KEGG) (Supplementary Table 1). However, Wnt signaling was consistently enriched across the two databases (Figure 1E) with a clear differential expression pattern across the two cell lines (Supplementary Figure 1B). Moreover, enrichment maps of Wikipathway database terms highlighted a cluster of gene sets comprising Wnt signaling and pluripotency regulation (Figure 1F), suggesting a potential acquisition or enrichment of stem cell features in carboplatin resistant cells.
To further understand the differences in stem cell transcriptional features between 468’R and 468'P cells, we compared our transcriptomic data with curated gene sets comprised of genes found overexpressed in human embryonic stem cells (hESC) and cancer stem cells (26, 46). Interestingly, 468’R seems to be transcriptionally closer to embryonic and cancer stem cells than 468’P as determined by GSEA (Figure 1G). In addition, we compared our transcriptomic data with a gene set comprised of targets of pluripotency transcription factors NANOG, OCT4, and SOX2 in hESCs determined by chromatin immunoprecipitation followed by DNA sequencing (26). GSEA revealed enrichment of pluripotency transcription factor target genes in 468’R cells, supporting the putative acquisition of stem cell features in this cell line (Figure 1H).
To functionally validate the observed differences in Wnt/β-catenin signaling between 468’R and 468’P cells, we analyzed protein levels of non-phosphorylated (active) β-catenin. Importantly, western-blot analysis of total protein extracts in baseline untreated conditions revealed strong enrichment of active β-catenin in 468’R cells compared to the parental counterpart, thereby confirming the functional activation of Wnt signaling in drug-resistant cells (Figure 1I). Moreover, the accumulation of active β-catenin was accompanied by transcriptional activation of Wnt-reporter activity (Supplementary Figure 1C).
To quantify differences in frequency of putative cancer stem cell populations in both cell lines, we used flow cytometry to assess the enzymatic activity of aldehyde dehydrogenases (ALDH) and the expression level of the cell surface markers CD44 and CD24. Both methods have been used to identify, quantify, and isolate putative cancer stem cells from different cancer types. High ALDH activity and CD44/CD24 expression ratio in TNBC have been shown to correlate with enhanced tumorigenesis and metastatic potential as wells as radio- and chemotherapy resistance (23, 25, 48). Flow cytometric analysis showed significant differences in ALDH positive (Figure 1J) and CD44+/CD24- cells (Figure 1K), with 468’R cells expressing higher levels of both markers. Gene expression analysis also revealed an enrichment of the core pluripotency regulators OCT4, NANOG, and cMYC, as well as cancer stem cell marker LGR5 (49) (Figure 1L).
We performed a tumorsphere formation assay to functionally evaluate differences in cancer stem cell properties and in vitro tumor-initiating capacity. In vitro growth in non-adherent conditions has been described as an exclusive capability of cancer stem cells, thereby functioning as a surrogate measure of in vitro tumor-initiating capacity and as a method to enrich cancer stem cells (50). Importantly, in line with our stemness-related gene expression data and flow cytometry analysis of ALDH activity and CD44/CD24 expression, we observed a significantly higher tumorsphere formation frequency in 468’R cells compared to 468’P (Figure 1M).
Altogether, the transcriptomic evidence for alterations in Wnt signaling, and presumably stem cell features, between 468’R and 468’P cells suggests the possibility of its involvement in mediating carboplatin resistance in our isogenic cell line.
Pharmacological Activation of Wnt Signaling in 468'P Cells, Disrupts Carboplatin-Response and Enhances Stemness and Pluripotency Marker Expression
We hypothesized that modulation of Wnt signaling in the parental 468’P cell line could recapitulate the carboplatin-resistant phenotype and increase the expression level of pluripotency and stem cell-related genes. To test this hypothesis, we treated 468’P with a small molecule inhibitor of GSK3, CHIR99021 (CHIR), thereby preventing β-catenin degradation and activating the Wnt/β-catenin pathway. To confirm Wnt signaling activation, we used a lentiviral fluorescent reporter of canonical Wnt transcriptional activity (TOPGFP) (51). 468’P displayed low basal levels of Wnt-reporter activity but promptly induced the reporter upon GSK3β inhibition, with almost 100% of cells becoming GFP positive within 12 hours of treatment (Figure 2A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/705384/fonc-11-705384-HTML-r1/image_m/fonc-11-705384-g002.jpg)
Figure 2 Pharmacological in vitro Wnt induction prevents carboplatin-induced apoptosis and upregulates stem cell marker expression. (A) Representative flow cytometry scatterplots (left) of Wnt-reporter MDA-MB-468 TOPGFP cells treated with vehicle (DMSO) or GSK3 inhibitor (CHIR, 4 μM) for 12 hours and statistical analysis of the mean frequency of GFP+ cells using Welch’s t-test (right, n=3). (B) Phase-contrast microscopy images of 468’P cells treated with or without carboplatin in the presence of CHIR or DMSO (left, scale bar: 100 μm) and statistical analysis of absolute cell numbers after 72 hours of each treatment using One-way ANOVA corrected for multiple comparisons using the Holm-Sidak method (n=4). (C) Representative flow cytometry scatterplots of annexin V staining of 468’P cells (left) treated with or without carboplatin in the presence of DMSO or CHIR (4 μM) for 72 hours and statistical analysis of the mean frequency of annexin V positive cells (right) using one-way ANOVA corrected for multiple comparisons using the Holm-Sidak method (n=3). (D) Relative mRNA expression of Wnt target and stem cell markers upon 72-hour treatment with DMSO or CHIR (4 μM) in 468’P cells (n=3). Multiple t-tests with Holm-Sidak correction for multiple comparisons. (Barplots represent mean + SEM. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, non significant).
We observed a significant rescue of survival when treating 468’P with CHIR combined with carboplatin (Figure 2B). Carboplatin-induced apoptosis was also significantly reduced when cells were co-treated with CHIR and CAR, as determined by flow cytometric analysis of Annexin V positivity (Figure 2C). Also, GSK3β inhibition led to the upregulation of OCT4 and NANOG pluripotency markers (Figure 2D).
Similar results were observed by inhibiting GSK3β on a second TNBC cell line (MDA-MB-231) (Supplementary Figures 2A, B), suggesting this effect is not cell line-specific.
β-Catenin Overexpression Induces Carboplatin Resistance in 468’P and Enhances Stem Cell Features
GSK3 is a multi-substrate serine-threonine kinase that regulates a multitude of signaling pathways (52). As such, we proceeded to investigate whether the effects of pharmacological activation of Wnt signaling on stemness and carboplatin resistance observed upon GSK3 inhibition could also be induced by direct overexpression of β-catenin. To achieve this, we transduced 468’P cells with a lentiviral vector encoding a truncated, constitutively active mutant β-catenin (ΔN90-β-catenin) (53), generating β-catenin overexpressing cells (468’OE) (Figure 3A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/705384/fonc-11-705384-HTML-r1/image_m/fonc-11-705384-g003.jpg)
Figure 3 β-catenin overexpression in 468’P induces carboplatin-resistance, pluripotency-related gene expression, and cancer stem cell features. (A) Western blot (top) of total β-catenin in MDA-MB-468 cells transduced with an empty vector or truncated, constitutively active β-catenin isoform ΔN90 and phase-contrast microscopy (down). (B) Enriched gene sets from Wikipathways database by one-tailed GSEA of ranked DEGs between 468’OE and 468’P sorted by normalized enrichment score (left) and enrichment map illustrating pathways most significantly different between 468’OE and 468’P (right). (C) Non-linear fit model of [CAR] vs. normalized response for IC50 determination. (468’OE: n=6, R2 = 0.92; 468’CTRL: n=6, R2 = 0.95). (D) Representative flow cytometry scatterplots of annexin V staining (left) of 468’CTRL and 468’OE cells treated with carboplatin 2 μM for 72h and statistical analysis of the mean frequency of annexin V positive cells using one-way ANOVA corrected for multiple comparisons using the Holm-Sidak method (right, n=3). (E) mRNA level fold change (Log2) of CTNNB1 (β-catenin), Wnt target AXIN2, and stem cell markers in 468’OE cells vs. 468’CTRL (n=4). Multiple t-tests with Holms-Sidak correction for multiple comparisons. (F) Representative scatterplots of flow cytometric analysis of aldehyde dehydrogenase activity (left) and statistical analysis of the mean percentage of ALDH+ cells in 468’OE (n=5) and 468’CTRL (n=5) using Welch’s t-test (right). (G) Representative scatterplots of flow cytometric analysis of CD44-PE and CD24-APC immunolabeling (left) and corresponding statistical analysis of the mean percentage of CD44+/CD24- cells using Welch’s t-test (right; n=3). (H) Representative brightfield images of tumorspheres generated from 468’CTRL and 468’OE cells (left, scale bar: 50 μm) and statistical analysis of mean tumorsphere forming units (number of spheres/number of seeded single cells) using Welch’s t-test (right; n=3). (Barplots represent mean + SEM. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, non significant).
We then performed mRNA-sequencing to investigate biological processes altered by overexpression of β-catenin in TNBC cells. Differentially expressed genes (Supplementary Table 2) were ranked by fold-change and analyzed by GSEA using the Wikipathways database. We detected an interesting enrichment of gene sets related to both pluripotency and differentiation, namely Embryonic Stem Cell Pluripotency Pathways (M39530: WP_ESC_PLURIPOTENCY_PATHWAY) and Wnt signaling in pluripotency (M39387: WP_WNT_SIGNALING_PATHWAY_AND_PLURIPOTENCY) (Figure 3B). In addition, GSEA revealed a significant correlation between the transcriptional profile of 468’OE cells and hESC and pluripotency transcription factor gene sets (Supplementary Figure 3A).
468’OE displayed increased carboplatin resistance in comparison with 468’CTRL cells (empty lentiviral vector) (IC50 468’OE: ~7 μM vs. IC50 468’CTRL: ~1.5 μM) (Figure 3C). Accordingly, flow cytometric analysis revealed that upon overexpression of active β-catenin, 468’OE cells fail to induce apoptosis when challenged with carboplatin, displaying Annexin V positivity frequencies at the same level of untreated cells (Figure 3D). These data confirm the direct involvement of β-catenin in the acquisition of in vitro carboplatin resistance in TNBC.
Next, we probed gene expression of Wnt targets and pluripotency and cancer stem cell markers by RT-qPCR. As observed with pharmacological activation of Wnt, β-catenin overexpression also induced a significant upregulation of pluripotency markers OCT4 and NANOG but also cMYC and LGR5 (Figure 3E).
Further, in accordance with the increased expression of pluripotency and stemness-related genes, flow cytometric analysis of ALDH activity (Figure 3F) and CD44/CD24 expression (Figure 3G) corroborated the presumptive induction of an enhanced cancer stem cell phenotype upon overexpression of β-catenin in TNBC cells.
Finally, we investigated whether overexpression of β-catenin functionally endows TNBC cells with enhanced tumorsphere formation capacity. Indeed, 468’OE cells displayed a significantly higher sphere-forming efficiency when grown in 3D suspension culture conditions, confirming that β-catenin overexpression functionally enhances in vitro stem cell properties in TNBC cell lines (Figure 3H).
The overexpression of ΔN90-β-catenin on MDA-MB-231 cells yielded the same effect with a roughly 10-fold increase in the IC50 of carboplatin, as well as increased ALDH activity and tumorsphere formation capacity (Supplementary Figures 3B–F).
Wnt Inhibition Disrupts Carboplatin Resistance in 468’R Cells and Downregulates Cancer Stem Cell and Pluripotency Marker Expression
Wnt signaling is deregulated in 468’R cells, and β-catenin overexpression on parental cells confirmed its role in mediating the carboplatin-resistant stem-like phenotype. Given these observations, we hypothesized that inhibition of Wnt signaling could restore sensitivity in 468’R cells. To that end, we used LGK974 (LGK), a small molecule inhibitor of the endoplasmic reticulum palmitoyltransferase porcupine (PORCN). This enzyme is responsible for processing Wnt ligands for secretion, mediating a crucial step of Wnt-dependent signaling (54).
Combinatorial treatment of 468’R cells with 2 μM of CAR and increasing concentrations of LGK resulted in carboplatin-sensitivity in a dose-dependent manner, starting at 200 nM of the Wnt inhibitor (Figure 4A). Interestingly, LGK treatment alone (200 nM) was insufficient to induce apoptosis in both 468’R and 468’P. However, when added to carboplatin, LGK induced strong annexin V positivity in 468’R cells, indicating rescue of carboplatin sensitivity in the resistant cell line (Figure 4B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/705384/fonc-11-705384-HTML-r1/image_m/fonc-11-705384-g004.jpg)
Figure 4 Wnt inhibitor LGK974 disrupts carboplatin resistance and pluripotency gene expression. (A) Absolute cell number of 468’R cells treated for 72 hours with increasing concentrations of LGK974 in the presence or absence of 2 μM carboplatin. Multiple t-tests (n=3). (B) Phase-contrast microscopy of 468’R cells treated for 72 hours with 200 nM LGK974 or DMSO in the presence or absence of 2 μM carboplatin (left). Mean frequency of annexin V positive cells in 468’P and 468’R cells treated with or without 200 nM LGK974 in presence or absence of 2 μM carboplatin showing the resensitization of 468’R cells to carboplatin when co-treated with Wnt inhibitor (right, n=3). One-way Anova with correction for multiple comparisons using the Holm-Sidak method. (C) Relative mRNA expression of Wnt target and stem cell markers upon 72-hour treatment with DMSO or 200 nM LGK974 with or without 2 μM carboplatin in 468’R cells (n=3). Multiple t-tests: Unt vs. Unt+LGK & CAR vs. CAR+LGK. (Barplots represent mean + SEM. *p < 0.05; **p < 0.01; ****p < 0.0001; ns, non significant).
Gene expression analysis by qPCR of LGK-treated 468’R cells confirmed downregulation of Wnt signaling by reduced expression of Wnt target AXIN2. More importantly, LGK severely reduced transcript levels of pluripotency markers OCT4, NANOG, and cMYC, both in the presence or absence of carboplatin (Figure 4C). These results corroborate the hypothesis that Wnt primes TNBC cells for carboplatin resistance by maintaining a stem-cell-like phenotype.
Inducible β-Catenin Knockdown Resensitizes 468’R Cells to Carboplatin and Disrupts Expression of Stem Cell Markers and Tumorsphere Formation Capacity
Next, we investigated whether cancer stem cell markers and function could be manipulated by directly disrupting β-catenin making use of doxycycline (DOX) inducible short-hairpin RNA targeting CTNNB1 (β-catenin) transcripts in 468’R cells (iCTNNB1-KD) (55).
Gene expression analysis by RT-qPCR confirmed a reduction of roughly 90% in CTNNB1 transcripts upon DOX treatment of iCTNNB1-KD cells whereas, as expected, cells expressing inducible scrambled shRNA (iSCRMBL) maintained basal β-catenin transcript levels (Figure 5A). Notably, DOX-induced β-catenin knockdown led to a robust reduction of pluripotency and cancer stem cell markers, confirming the role of β-catenin in maintaining their high expression (Figure 5A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/705384/fonc-11-705384-HTML-r1/image_m/fonc-11-705384-g005.jpg)
Figure 5 Inducible β-catenin shRNA disrupts carboplatin-resistance and stem cell function in 468’R cells. (A) Relative mRNA expression level of CTNNB1 (β-catenin) and Wnt target and pluripotency markers in 468’R cells transduced with inducible CTNNB1-targeting or SCRMBL shRNAs, in the presence or absence of doxycycline (n=3). Welch’s t-test (Dox vs. no Dox). (B) Non-linear fit model of [CAR] vs normalized response for IC50 determination in iCTNNB1-KD cells in presence or absence of doxycyclin. (n=3, R2 + DOX: 0.93, R2 - DOX: 0.95). (C) Representative scatterplots of flow cytometric analysis of apoptosis by annexin V staining of 468’R iCTNNB1-KD and 468’R iSCRBML cells treated with or without 2 μM carboplatin for 72h in presence or absence of doxycycline (left) and corresponding statistical analysis of the mean frequency of annexin V positive cells (right, n=6). One-way Anova with correction for multiple comparisons using the Holm-Sidak method. (D) Representative brightfield images of tumorspheres generated from 468’R iSCRMBL and 468’R iCTNNB1-KD cells (left, scale bar: 50 μm) and statistical analysis of mean tumorsphere forming units (number of spheres/number of seeded single cells) (right; n=3). Welch’s t-test (Dox vs. no Dox). (Barplots represent mean + SEM. *p < 0.05; ***p < 0.001; ****p < 0.0001; ns, non significant).
Next, we investigated the effect of β-catenin knockdown on in vitro carboplatin response by determining the IC50 for 468’R iCTNNB1-KD cells in the presence or absence of DOX. Suppressing the expression of β-catenin had a remarkable effect on carboplatin resistance. DOX-treated iCTNNB1-KD cells displayed a substantial reduction in measured IC50 compared to non-induced cells (IC50 -DOX: ~11 μM vs. IC50 +DOX: ~1 μM) (Figure 5B). Importantly, 468’R iSCRMBL cells displayed no significant changes in IC50 (-DOX: ~13 μM vs. +DOX: ~11 μM) (Supplementary Figure 4). Additionally, when exposed to the IC50 of 468’P cells, 468’R iCTNNB1-KD cells strongly induced apoptosis in the presence of DOX, while no difference was observed in iSCRMBL cells (Figure 5C).
Finally, we assessed the effect of β-catenin suppression on tumorsphere formation as a functional readout for stem cell activity. In line with the downregulation of stem cell marker expression in the presence of DOX, iCTNNB1-KD cells displayed a significantly lower tumorsphere forming frequency when β-catenin shRNA was induced. On the other hand, differences upon the induction of SCRMBL shRNA were negligible (Figure 5D).
Wnt Inhibition Disrupts In Vivo Carboplatin Resistance in a Carboplatin-Resistant TNBC Patient-Derived Xenograft
Patient-derived xenografts are essential in vivo models of human neoplasms. Moreover, PDX models retain with excellent fidelity histological and molecular features of originating tumors, representing essential tools for assessing drug resistance and response (56–58).
To study whether, like in 468’R cells, in vivo PORCN inhibition with LGK re-establishes carboplatin sensitivity, we used an isogenic carboplatin-resistant TNBC PDX (C4O) obtained from a previously chemotherapy-sensitive model (BRC016) (35). In brief, ten mice bearing BRC016 tumors were treated with 50 mg/kg of carboplatin once weekly for three weeks. Nine out of ten BRC016-bearing mice achieved a complete response to treatment with no tumors detectable. However, one tumor displayed a late response and became resistant to further treatment despite undergoing a substantial volume reduction. The non-responder xenograft eventually regrew from the post-treatment residual tumor (Figure 6A). Material from the regrown tumor was collected to establish an isogenic model of carboplatin-resistant TNBC (C4O) and gene expression analysis revealed changes in Wnt signaling target AXIN2 and pluripotency and stem markers NANOG, OCT4, SOX2, and LGR5 (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/705384/fonc-11-705384-HTML-r1/image_m/fonc-11-705384-g006.jpg)
Figure 6 WNT inhibition disrupts in vivo carboplatin-resistance in a carboplatin-resistant TNBC Patient-Derived Xenograft. (A) Tumor growth curves of carboplatin-sensitive BRC016 TNBC PDX model. Nine out of 10 mice show a complete response to treatment. One animal (red line) had a very delayed response and still had residual tumor mass after three weeks of treatment. The residual xenograft resumed growth after carboplatin treatment was stopped. This tumor was collected to establish a carboplatin-resistant model (C4O). (B) Comparative gene expression analysis by qRT-PCR of Wnt target AXIN2 and stem cell markers in BRC016 carboplatin-sensitive PDX and the C4O carboplatin-resistant isogenic PDX (n=4). Welch’s t-test. (C) Tumor growth curves of C4O carboplatin-resistant PDX treated with VEH, LGK974, CAR, or CAR+LGK showing reduced tumor growth in the combinatorial treatment arm (VEH, CAR, CAR+LGK n=6 and LGK n=5) (left). Two-way ANOVA with Tukey correction. The green-shadowed area under the curve represents highlights the time points in which the difference between CAR and CAR+LGK is statistically significant. Representative photographs of tumors in each treatment arm at day 21 of treatment (right). (D) mRNA level fold change (Log2) vs. VEH treatment (no carboplatin and no Wnt inhibitor) in tumors dissected at treatment endpoint (21 days) (n=3). Multiple t-tests (VEH vs. LGK and CAR vs. CAR+LGK). (E) Representative confocal microscopy images (left) of TUNEL staining in green and human pan-cytokeratin immunolabeling in red and respective quantification and statistical analysis of TUNEL positive cells (VEH n=4, LGK n=4, CAR & CAR+LGK n=6). One-way ANOVA with correction for multiple comparisons using the Holm-Sidak method. (Barplots represent mean + SEM. *p < 0.05; **p < 0.01; ****p < 0.0001; ns, non significant).
In vivo carboplatin resistance was maintained in subsequent generations of transplanted C4O PDX models. We detected no significant differences in mean final tumor volumes in animals treated with vehicle or 50 mg/kg carboplatin, administered once weekly intraperitoneally, for three weeks (Supplementary Figure 5A). However, combinatorial treatment with daily dosing of PORCN inhibitor LGK drastically reduced C4O tumor growth (Figure 6C and Supplementary Figure 5A). LGK was well tolerated both as mono-therapy and in combination with carboplatin, and no significant changes in mouse body weight were observed during treatment (Supplementary Figure 5B). Interestingly, LGK or CAR alone could not reduce C4O growth (Figure 6C), and no significant differences in mean final tumor volumes between VEH, CAR, and LGK-treated animals were observed. Importantly, we did not find differences in Ki67 positivity between any treatment arms, indicating that reduced tumor growth in CAR+LGK treated animals was not due to differences in proliferation (Supplementary Figure 5C).
Gene expression analysis by RT-PCR revealed that similar to what we observed upon treating 468’R cells with LGK, inhibition of PORCN in C4O led to the depletion of pluripotency marker expression both in the presence or absence of carboplatin (Figure 6D).
Finally, given the absence of alterations in the expression of Ki67, we sought out to understand whether the drastic reduction in tumor growth in animals treated with the combination of CAR and LGK could be due to increased apoptotic cell death. For this, we performed fluorescent terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) to detect DNA fragmentation as a readout of apoptosis. We combined TUNEL staining with human pan-cytokeratin immunolabeling to enable quantification of apoptotic signal specifically in cancer cells. No differences in TUNEL positivity were measured in VEH, LGK, and CAR-treated animals. However, animals treated with the combination of CAR and LGK displayed a significant increase in apoptotic TUNEL signal (Figure 6E).
Wnt Signaling Is Deregulated in Patients With Platinum-Resistant TNBCs and High-Grade Serous Ovarian Cancer
To understand whether alterations in Wnt signaling are prevalent in platinum-resistant human TNBCs, we analyzed a public RNA microarray dataset (GSE103668) comprised of 21 pre-treatment samples from TNBC patients treated with cisplatin and bevacizumab (59). Due to the scarcity of additional platinum-treated TNBC transcriptomic data available in public repositories, we decided to include a supplementary dataset from high-grade serous ovarian cancer (HGSOC) (E-MTAB-7083) (Figure 7A). The reason for this choice lies in the extensive use of platinum-based chemotherapy in this type of cancer and the striking overlap in clinical and molecular features between HGSOC and TNBC (61).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/705384/fonc-11-705384-HTML-r1/image_m/fonc-11-705384-g007.jpg)
Figure 7 Wnt signaling is deregulated in patients with platinum-resistant TNBCs, high-grade serous ovarian cancer, and isogenic cisplatin-resistant ovarian cancer cell lines. (A) Summary of datasets and analysis methodology used. Functional enrichment and mapping as previously reported (40, 41, 60). (B) Enriched KEGG gene sets in patients with platinum-resistant TNBC (left) and HGSOC (right). (C) Enrichment maps for visualization of enriched KEGG gene sets in patients with platinum-resistant TNBC (left) and HGSOC (right). (D) Enrichment maps for visualization of gProfiler functional enrichment analysis of ranked, upregulated DEGs in patients with platinum-resistant TNBC (left) and HGSOC (right).
The available clinical metadata was analyzed in both datasets to classify patients based on response to treatment, and differential gene expression was calculated between responders and non-responders. Interestingly, GSEA analysis on ranked differentially expressed genes using the KEGG database retrieved Wnt signaling as one of the top enriched terms with False Discovery Rate (FDR) <10% in TNBC and HGSOC patients with no response to platinum therapy (Figure 7B). Other enriched KEGG terms comprised biological processes such as focal adhesion, extracellular matrix interaction, and other signaling pathways such as TGFβ and Hedgehog. Enrichment maps of GSEA hits for both datasets contained distinctive Wnt-related clusters involving gene sets with overlapping enriched genes such as “Melanogenesis” and “Basal Cell Carcinoma” for TNBC patients and HGSOC with the addition of “Hedgehog Signaling Pathway” in the latter (Figure 7C). To obtain a broader perspective of the function of differentially expressed genes, we performed functional enrichment analysis for both datasets using gProfiler to retrieve enriched gene ontology, Reactome, and Wikipathway gene sets to build enrichment maps. We obtained a distinctive cluster of Wnt-related terms in both datasets, including regulating both canonical and non-canonical Wnt signaling in TNBC and regulating pluripotency and Wnt ligand biogenesis and secretion in HGSOC (Figure 7D).
Altogether, these results highlight Wnt signaling’s importance in mediating platinum resistance in human TNBC and suggest transversal resistance mechanisms across TNBC and HGSOC.
Discussion
Primary and acquired resistance to chemotherapy poses a critical hurdle in the treatment of cancer. This is particularly important in TNBC due to the relatively limited therapeutic toolbox available and the daunting clinical characteristics of this disease. In the continued absence of targeted molecular therapies, we must strive to improve response to current therapeutic options given the high probability of shorter survival when pCR is not achieved. The use of platinum compounds in combination with other agents or as a standalone treatment in TNBC is still under intense investigation but already shows the potential to improve pCR rates in this breast cancer subtype. However, how TNBCs specifically develop resistance to platinum-based treatment is still a poorly understood process.
In this study, we used an isogenic carboplatin-resistant TNBC cell line and used next-generation mRNA sequencing to identify transcriptomic differences between sensitive and resistant cells. Functional enrichment analysis indicated, among others, the existence of profound differences in transcription of Wnt signaling and pluripotency-related genes. Active β-catenin protein levels, together with pluripotency transcription factor and cancer stem cell marker expression, confirmed enrichment of Wnt signaling and stem cell features in resistant cells. Importantly, resistant cells displayed significantly higher in vitro tumorsphere formation efficiency, a strong indicator of the functional acquisition of cancer stem cell features which could underly chemotherapy resistance. However, further in vivo characterization of this model of carboplatin resistance should be performed to validate enrichment tumor-initiating capacity thoroughly.
Nonetheless, we deemed the observed in vitro enrichment in Wnt signaling, and stem cell features significant since there is a link between this pathway and cancer stem cell biology and tumorigenesis. In addition, several studies have demonstrated the involvement of this pathway in mediating resistance to chemotherapy and radiation in different types of cancer, including breast (30, 34, 62–64). Wnt signaling is also known to specifically mediate platinum resistance in endometrial cancer (65), ovarian (66), and oral squamous cell carcinoma (67). However, little is known regarding the role of Wnt signaling in mediating resistance to platinum in TNBC. Interestingly, Wnt is often deregulated in breast cancers, particularly TNBCs, despite the negligible frequency of mutations in Wnt pathway components. Significantly, Wnt activation is associated with poor clinical outcomes in TNBC (68).
Based on the transcriptomic data herein generated, we hypothesized that stem-cell gene expression and carboplatin resistance could be induced on parental MDA-MB-468 cells by manipulating Wnt signaling. For this, we first used CHIR, a small molecule inhibitor of GSK3β, thereby activating Wnt. Our results showed a significant increase in pluripotency marker expression and reduced apoptosis upon concomitant treatment with Wnt agonist and carboplatin. Past studies regarding the role of this multi-substrate kinase in treatment resistance are rather intriguing. While our results confirm other studies that report the enrichment of stem cell features upon GSK3β inhibition (69, 70), they seem to contradict reports of GSK3 inhibition leading to reduction of tumor growth and apoptosis (71–73). However, it is essential to note that GSK3 regulates several pathways and a myriad of cellular functions. In addition, its’ role in cancer is not only complex but also controversial since it has both tumor-suppressing and tumor-promoting functions (74). One example of its’ tumor-suppressor capabilities is through its’ inhibitory role on the Wnt signaling pathway, with the expression of kinase-inactive mutants leading to oncogenic conversion of epidermal cells (75). On the other hand, GSK3β is overexpressed in certain cancers such as colon, liver, and pancreas. In the latter, depletion of GSK3β was shown to drastically reduce tumor growth, highlighting the pro-tumoral effects of GSK3β (76). Given this, it is highly likely that the different outcomes of GSK3β inhibition may be specific to different tissues and cancer types.
Nonetheless, by overexpressing β-catenin in parental cells, we could replicate the phenotype of the isogenic carboplatin-resistant cells and GSK3 inhibition extensively. Specifically, β-catenin overexpression induced a significant increase in expression of pluripotency markers and induced carboplatin resistance, highlighting the role of Wnt pathway activation on drug resistance in TNBC.
To further validate the role of Wnt signaling and specifically β-catenin in mediating carboplatin resistance in TNBC, we made use of another TNBC cell line: MDA-MB-231. Several studies report the enrichment of cancer stem cell features such as high expression of the CD44+/CD24- immunophenotype in these cells (25, 77). Both pharmacological activation of Wnt and β-catenin overexpression successfully induced pluripotency marker expression and enhanced resistance to carboplatin. Additionally, β-catenin overexpression significantly enhanced ALDH activity and tumorsphere formation. These are interesting observations given the pre-existing enrichment in CD44+/CD24- in this cell line.
Nonetheless, another study has established a clear correlation between enhanced Wnt activity and enriched cancer stem cell potential in MDA-MB-231 cells (78). Importantly, in standard culture conditions, these cells display very low Wnt transcriptional activity but respond robustly to exogenous Wnt activation. This suggests some level of phenotypic plasticity that can be further skewed towards a more stem-like state by high levels of Wnt activation.
Inhibition of Wnt ligand secretion using LGK in isogenic carboplatin-resistant cells disrupted stem cell markers’ expression and reversed resistance. Interestingly, inhibition of Wnt ligand secretion alone was insufficient to induce apoptosis in either carboplatin sensitive or resistant cells, despite severely downregulating stem cell marker expression. This indicates that Wnt secretory signals are not necessarily essential for the survival of either sensitive or resistant cells, but rather that these signals prime the latter for survival upon challenge by carboplatin. Secreted Wnt-ligands originating from cells of the tumor microenvironment have been shown to significantly contribute to resistance by mediating immune evasion, cancer progression, and cancer stem cell maintenance (79). We here provide in vitro evidence for cancer cells themselves being an additional source of resistant-mediating Wnt-ligands in a cell-autonomous manner. Understanding which cancer cell-derived ligands are being secreted and their influence over other cell types associated with the tumor microenvironment would prove beneficial for developing novel therapeutic approaches and discovering new biomarkers for drug response.
LGK prevents the secretion of all Wnt ligands by inhibiting the palmitoyl acyltransferase PORCN. In vertebrates, the Wnt family of lipid-modified secreted signaling proteins comprises 19 members, conferring a great deal of complexity to this pathway. Wnt signaling includes a canonical or Wnt/β-catenin dependent pathway and the non-canonical or β-catenin-independent pathway (80). For that reason, it was essential to determine whether LGK-induced carboplatin sensitivity was mediated directly by β-catenin. Silencing β-catenin downregulated stem cell marker expression and dramatically reversed resistance to carboplatin, phenocopying LGK effects, and directly implicating canonical Wnt signaling in the maintenance of carboplatin resistance.
Inhibition of Wnt ligand secretion, namely through the inhibition of PORCN, has been under scrutiny during the last decade as a potential therapeutic approach for different types of cancer. LGK specifically has shown excellent preclinical efficacy in Wnt-addicted models (81) and is under examination in phase I clinical trial for several Wnt-dependent solid malignancies, including TNBC (NCT01351103) (82). We evaluated whether this molecule could resensitize an in vivo model of carboplatin-resistant TNBC. Previous studies reported that LGK alone significantly reduces the growth of a murine mouse mammary tumor-Wnt3 model (81). To our surprise, daily dosing of LGK alone had no impact on tumor growth. This is analogous to what we observed when treating carboplatin-resistant cells with LGK alone. Importantly, we observed a significant reduction of tumor growth in animals treated with a combination of daily LGK and weekly carboplatin.
The resistance we observe in our models is likely multifactorial, and Wnt is probably one of several deregulated processes which ultimately contribute to it. Nonetheless, it is crucial to consider that Wnt signaling has an overarching set of functions that lead to chemotherapy resistance. Fischer and colleagues demonstrated that in colorectal cancer PDX models, monoclonal antibodies targeting RSPO3 (a potent Wnt-signaling enhancer that binds LGR5 to potentiate Wnt activation by reducing Wnt-receptor turnover) had a strong synergistic effect against tumor growth when administered in combination with paclitaxel (83). However, when anti-RSPO3 antibodies were administered alone, they failed to reduce tumor progression. Interestingly, anti-RSPO3 administration alone significantly decreased the frequency of cancer stem cells, despite the absence of deleterious effects on tumor viability. The combinatorial treatment further decreased cancer stem cell frequency, and a surge in differentiated paclitaxel-sensitive cells was observed. An earlier study reports a very similar effect in breast cancer xenografts treated with antibodies targeting Frizzled (Wnt-ligand receptors) and taxanes, where antibody treatment alone had a limited effect on tumor growth and combined with taxanes resulted in a potent synergism (84). We report a similar therapeutic dynamic with LGK in our resistant models of TNBC, wherein LGK or β-catenin knockdown disturb stem cell functions without affecting in vitro tumorsphere formation and proliferation and in vivo tumor growth.
In addition to regulating and maintaining cancer stem cells, Wnt is known to mediate chemotherapy and radiotherapy resistance in non-stem cancer cells. One example is the enhancement of DNA damage repair in ovarian cancers (85, 86). Moreover, Wnt directly regulates the transcription of several genes encoding proteins associated with resistance to chemotherapy in diverse cancer types, such as MDR1 (Multi-drug resistance protein 1, also known as ABCB1), Survivin, and MMP7 (87). However, lower expression levels of these proteins could likely have no deleterious effects in cancer cells in the absence of cytotoxic aggression by chemotherapeutic agents.
We thus hypothesize that the synergistic effect of LGK and carboplatin is likely mediated by the combination of a shift in population dynamics from a cancer stem cell-like state to a generally more differentiated and drug-sensitive phenotype. This effect compounds with the repression of other resistance-inducing features that are non-exclusive of stem cells, thereby enhancing the response to treatment upon exposure to exogenous stressors.
Finally, we were able to identify similarities in transcriptional profiles of patients with platinum-refractory ovarian and triple-negative breast neoplasms. Wnt signaling deregulation is a common denominator in both datasets we analyzed. Given the overlapping clinical and molecular features of both cancer types, it would be interesting to investigate whether transversal resistance mechanisms exist for other therapies and other cancers. There is growing evidence that different cancers harbor rare stem cell populations, significantly contributing to therapy resistance and relapse. In addition, in many cancer types, Wnt/β-catenin is enriched in cancer stem cells.
Interestingly, Wnt activation is also coupled with a poor prognosis in several cancers, possibly indicating the existence of a common role for Wnt signaling in cancer stem cell regulation across different neoplastic diseases and suggesting that a common therapeutic approach focusing on cancer stem cell eradication could be effective in a wide variety of cancers. This warrants the necessity for developing and characterizing isogenic in vitro and in vivo models of resistance and the preclinical evaluation of Wnt-targeting molecules and their effect on cancer stem cell-mediated treatment resistance in other cancer types.
Altogether, this study demonstrates that response to platinum can be improved, and stable in vitro and in vivo resistance can be reversed by a combinatorial approach to TNBC treatment which leverages the inhibition of Wnt signaling to disrupt resistance-inducing cancer stem cell functions.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ebi.ac.uk/arrayexpress/E-MTAB-10337.
Ethics Statement
Tumor tissue for PDX implantation was acquired from a patient who provided written informed consent, and the procedure was approved by the Commission of Medical Ethics of the University Hospitals Leuven (approval numbers S54185 and ML8713). The PDX models herein used are now banked at the Trace Patient Derived Tumor Xenograft Platform of UZ Leuven/KU Leuven. All animal experiments were performed at the Trace PDX platform and were approved by the Ethics Committee Research UZ/KU Leuven (approval number P038/2015). Experiments with cell lines were approved by the Ethische Comissie Onderzoek UZ/KU Leuven (S65166).
Author Contributions
WO and FL conceived and designed the study. SM generated in-vitro models of carboplatin resistance and the carboplatin-resistant C4O PDX under the supervision of DA and FA. YL generated the transgenic cell lines used in this study. WO and SM conducted all in vivo experimental work. EC conducted immunohistochemistry experiments. BV performed RNA-sequencing differential expression analysis, and WO performed downstream functional enrichment and visualizations. WO and partially YL carried out all flow-cytometry experiments. PA performed Western blots. WO carried out the remainder of the experimental work. Data analysis and figure preparation were performed by WO and reviewed by SM, DA, and FL. The manuscript was written by WO and reviewed, and approved by all authors. FL secured funding and supervised and guided experimental work and manuscript preparation. BV, KK, J-JV, MB and FA made comments on formatting, editing, and data analysis. All authors contributed to the article and approved the submitted version. WO and SM are co-first authors. DA and FL are co-last authors.
Funding
WO and BV are funded by the Research Foundation-Flanders with Ph.D. fellowships 1155619N and 11E7920N, respectively. EC holds a fellowship from Stichting Tegen Kanker – the Belgian Foundation against Cancer (FAF-F/2016/822, fellowship number ZKD2498-00-W02). The authors would like to extend their gratitude to the Research Foundation – Flanders (FWO) for the G091521N grant (FL and DA) and to Jo Van Biesbroeck, the breast cancer research fund Nadine de Beauffort and the KU Leuven Research Fund (C24/17/073) for funding the development of the models used in this research study (FA). FA is a senior researcher for the Research Foundation-Flanders (FWO). Trace staff is supported by the Stichting Tegen Kanker grant 2016-054.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We want to show our gratitude to Anchel de Jaime Soguero for his critical review of the manuscript. We also sincerely thank the Trace Leuven PDX platform for its help in setting up mouse experiments. Trace is part of the EuroPDX consortium.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.705384/full#supplementary-material
Supplementary Table 1 | Differentially expressed genes and GSEA of 468’R vs. 468’P.
Supplementary Table 2 | Differentially expressed genes and GSEA of 468’OE vs. 468’P.
Supplementary Table 3 | Primer List.
Supplementary Figure 1 | (A) Representative histogram of flow cytometry Ki67 frequency in 468’P and 468’R cells (left) and statistical analysis (n=3). (B) Expression heatmap of Wnt-related genes from the WNT Pathway Pluripotency cluster from Figure 1F. (C) Flow cytometry analysis of TOPGFP-Wnt reporter activity in 468’P and 468’R cells (n=3). (Barplots represent mean + SEM. *p <0.05, **p<0.01, ***p<0.001, ****p<0.0001, ns= non significant).
Supplementary Figure 2 | (A) Phase-contrast microscopy (left, scale bar: 100 μm) of MDA-MB-231 cells treated for 72 hours with 35 μM Carboplatin in the presence or absence of 4 μM CHIR and statistical analysis of mean absolute cell numbers (n=4). One-way ANOVA with correction for multiple comparisons using the Holm-Sidak method. (B) mRNA level fold change (Log2) Wnt target AXIN2, and stem cell markers in MDA-MB-231 cells treated with 4 μM CHIR vs. DMSO (n=1, average of technical replicates). (Barplots represent mean + SEM. *p <0.05, **p<0.01, ***p<0.001, ****p<0.0001, ns= non significant).
Supplementary Figure 3 | (A) GSEA of hESCs (26) (left), Cancer Progenitor (46) (center), and pluripotency transcription target gene sets (26) (right) in 468’OE vs 468’P cells. (B) Phase-contrast microscopy (scale bar: 100 μm) of control MDA-MB-23 and 231’OE (Δn90 β-catenin overexpression). (C) Western blot of total β-catenin in MDA-MB-231 cells transduced with an empty vector or truncated, constitutively active β-catenin isoform ΔN90. β-actin was used as a loading control. (D) Non-linear fit model of [CAR] vs. normalized response for IC50 determination (right). (n=2) (E) Representative scatterplots of flow cytometric analysis of aldehyde dehydrogenase activity (left) and statistical analysis of the mean percentage of ALDH+ cells in 231’OE and 231’Ctrl cells using Welch’s t-test (n=5) (right). (F) Representative brightfield images of tumorspheres generated from 231’Ctrl and 231’OE cells (left, scale bar: 100 μm) and statistical analysis of mean tumorsphere forming units (number of spheres/number of seeded single cells) (right; n=3). Welch’s t-test. (Barplots represent mean + SEM. *p <0.05, **p<0.01, ***p<0.001, ****p<0.0001, ns = non significant).
Supplementary Figure 4 | Non-linear fit model of [CAR] vs. normalized response for IC50 determination (right) in 468’R cells transduced with inducible scrambled non-targeting shRNAs in the presence or absence of DOX. (n=2).
Supplementary Figure 5 | (A) Final tumor volumes of VEH (n=6), LGK (n=5), CAR (n=6) and CAR+LGK (n=6) treated mice. One way-ANOVA with Holm-Sidak correction for multiple comparisons. (B) Bodyweight change in percentage of initial treatment bodyweight for VEH (n=6), LGK (n=5), CAR (n=6) and CAR+LGK (n=6) treated mice. No statistically significant changes detected by two way ANOVA. (C) Representative brightfield microscopy (20x magnification) of tumor sections labeled with anti-human Ki67 (left) and corresponding statistical analysis of the mean frequency of Ki67 positive cells. (VEH, CAR, CAR+LGK n=6 and LGK n=4). One-way ANOVA with Holm-Sidak correction for multiple comparisons. (Barplots represent mean + SEM. *p <0.05, **p<0.01, ***p<0.001, ****p<0.0001, ns = non significant).
References
1. Carey L, Winer E, Viale G, Cameron D, Gianni L. Triple-Negative Breast Cancer: Disease Entity or Title of Convenience? Nat Rev Clin Oncol (2010) 7:683–92. doi: 10.1038/nrclinonc.2010.154
2. Marra A, Trapani D, Viale G, Criscitiello C, Curigliano G. Practical Classification of Triple-Negative Breast Cancer: Intratumoral Heterogeneity, Mechanisms of Drug Resistance, and Novel Therapies. NPJ Breast Cancer (2020) 6:1–16. doi: 10.1038/s41523-020-00197-2
3. Liedtke C, Mazouni C, Hess KR, André F, Tordai A, Mejia JA, et al. Response to Neoadjuvant Therapy and Long-Term Survival in Patients With Triple-Negative Breast Cancer. J Clin Oncol (2008) 26:1275–81. doi: 10.1200/JCO.2007.14.4147
4. Rottenberg S, Disler C, Perego P. The Rediscovery of Platinum-Based Cancer Therapy. Nat Rev Cancer (2021) 21:37–50. doi: 10.1038/s41568-020-00308-y
5. Turner NC, Reis-Filho JS, Russell AM, Springall RJ, Ryder K, Steele D, et al. BRCA1 Dysfunction in Sporadic Basal-Like Breast Cancer. Oncogene (2007) 26:2126–32. doi: 10.1038/sj.onc.1210014
6. Shafee N, Smith CR, Wei S, Kim Y, Mills GB, Hortobagyi GN, et al. Cancer Stem Cells Contribute to Cisplatin Resistance in Brca1/p53-Mediated Mouse Mammary Tumors. Cancer Res (2008) 68:3243–50. doi: 10.1158/0008-5472.CAN-07-5480
7. Von Minckwitz G, Schneeweiss A, Loibl S, Salat C, Denkert C, Rezai M, et al. Neoadjuvant Carboplatin in Patients With Triple-Negative and HER2-Positive Early Breast Cancer (GeparSixto; GBG 66): A Randomised Phase 2 Trial. Lancet Oncol (2014) 15:747–56. doi: 10.1016/S1470-2045(14)70160-3
8. Sikov WM, Berry DA, Perou CM, Singh B, Cirrincione CT, Tolaney SM, et al. Impact of the Addition of Carboplatin and/or Bevacizumab to Neoadjuvant Once-Per-Week Paclitaxel Followed by Dose-Dense Doxorubicin and Cyclophosphamide on Pathologic Complete Response Rates in Stage II to III Triple-Negative Breast Cancer: CALGB 40603 (A. J Clin Oncol (2015) 33:13–21. doi: 10.1200/JCO.2014.57.0572
9. Loibl S, O’Shaughnessy J, Untch M, Sikov WM, Rugo HS, McKee MD, et al. Addition of the PARP Inhibitor Veliparib Plus Carboplatin or Carboplatin Alone to Standard Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer (BrighTNess): A Randomised, Phase 3 Trial. Lancet Oncol (2018) 19:497–509. doi: 10.1016/S1470-2045(18)30111-6
10. Iwase M, Ando M, Aogi K, Aruga T, Inoue K, Shimomura A, et al. Long-Term Survival Analysis of Addition of Carboplatin to Neoadjuvant Chemotherapy in HER2-Negative Breast Cancer. Breast Cancer Res Treat (2020) 180:687–94. doi: 10.1007/s10549-020-05580-y
11. Zhou J, Kang Y, Chen L, Wang H, Liu J, Zeng S, et al. The Drug-Resistance Mechanisms of Five Platinum-Based Antitumor Agents. Front Pharmacol (2020) 11:343. doi: 10.3389/fphar.2020.00343
12. Wang L, Jin F, Qin A, Hao Y, Dong Y, Ge S, et al. Targeting Notch1 Signaling Pathway Positively Affects the Sensitivity of Osteosarcoma to Cisplatin by Regulating the Expression and/or Activity of Caspase Family. Mol Cancer (2014) 13:139. doi: 10.1186/1476-4598-13-139
13. Jin L, Chun J, Pan C, Li D, Lin R, Alesi GN, et al. MAST1 Drives Cisplatin Resistance in Human Cancers by Rewiring Craf-Independent MEK Activation. Cancer Cell (2018) 34:315–30.e7. doi: 10.1016/j.ccell.2018.06.012
14. Seidl C, Panzitt K, Bertsch A, Brcic L, Schein S, Mack M, et al. MicroRNA-182-5p Regulates Hedgehog Signaling Pathway and Chemosensitivity of Cisplatin-Resistant Lung Adenocarcinoma Cells via Targeting GLI2. Cancer Lett (2020) 469:266–76. doi: 10.1016/j.canlet.2019.10.044
15. Atkinson CJ, Kawamata F, Liu C, Ham S, Győrffy B, Munn AL, et al. EGFR and Prion Protein Promote Signaling via FOXO3a-KLF5 Resulting in Clinical Resistance to Platinum Agents in Colorectal Cancer. Mol Oncol (2019) 13:725–37. doi: 10.1002/1878-0261.12411
16. Akrap N, Andersson D, Bom E, Gregersson P, Ståhlberg A, Landberg G. Identification of Distinct Breast Cancer Stem Cell Populations Based on Single-Cell Analyses of Functionally Enriched Stem and Progenitor Pools. Stem Cell Rep (2016) 6:121–36. doi: 10.1016/j.stemcr.2015.12.006
17. Pattabiraman DR, Weinberg RA. Tackling the Cancer Stem Cells — What Challenges Do They Pose? Nat Rev Drug Discov (2014) 13:497–512. doi: 10.1038/nrd4253
18. Moore N, Lyle S. Quiescent, Slow-Cycling Stem Cell Populations in Cancer: A Review of the Evidence and Discussion of Significance. J Oncol (2011) 2011:1–11. doi: 10.1155/2011/396076
19. Koury J, Zhong L, Hao J. Targeting Signaling Pathways in Cancer Stem Cells for Cancer Treatment. Stem Cells Int (2017) 2017:1–10. doi: 10.1155/2017/2925869
20. Badve S, Nakshatri H. Breast-Cancer Stem Cells-Beyond Semantics. Lancet Oncol (2012) 13:e43–8. doi: 10.1016/S1470-2045(11)70191-7
21. Clevers H. The Cancer Stem Cell: Premises, Promises and Challenges. Nat Med (2011) 17:313–9. doi: 10.1038/nm.2304
22. Yang F, Xu J, Tang L, Guan X. Breast Cancer Stem Cell: The Roles and Therapeutic Implications. Cell Mol Life Sci (2017) 74:951–66. doi: 10.1007/s00018-016-2334-7
23. Al-Hajj M, Wicha MS, Benito-Hernandez A, Morrison SJ, Clarke MF. Prospective Identification of Tumorigenic Breast Cancer Cells. Proc Natl Acad Sci (2003) 100:3983–8. doi: 10.1073/pnas.0530291100
24. Wei W, Lewis MT. Identifying and Targeting Tumor-Initiating Cells in the Treatment of Breast Cancer. Endocr Relat Cancer (2015) 22:R135–55. doi: 10.1530/ERC-14-0447
25. Li W, Ma H, Zhang J, Zhu L, Wang C, Yang Y. Unraveling the Roles of CD44/CD24 and ALDH1 as Cancer Stem Cell Markers in Tumorigenesis and Metastasis. Sci Rep (2017) 7:1–15. doi: 10.1038/s41598-017-14364-2
26. Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, et al. An Embryonic Stem Cell–Like Gene Expression Signature in Poorly Differentiated Aggressive Human Tumors. Nat Genet (2008) 40:499–507. doi: 10.1038/ng.127
27. Zhan T, Rindtorff N, Boutros M. Wnt Signaling in Cancer. Oncogene (2017) 36:1461–73. doi: 10.1038/onc.2016.304
28. Van Der Wal T, Van Amerongen R. Walking the Tight Wire Between Cell Adhesion and WNT Signalling: A Balancing Act for β-Catenin: A Balancing Act for CTNNB1. Open Biol (2020) 10:1–15. doi: 10.1098/rsob.200267rsob200267
29. Niehrs C, Acebron SP. Mitotic and Mitogenic Wnt Signalling. EMBO J (2012) 31:2705–13. doi: 10.1038/emboj.2012.124
30. Jun S, Jung YS, Suh HN, Wang W, Kim MJ, Oh YS, et al. LIG4 Mediates Wnt Signalling-Induced Radioresistance. Nat Commun (2016) 7:1–13. doi: 10.1038/ncomms10994
31. Chen S, Guttridge DC, You Z, Zhang Z, Fribley A, Mayo MW, et al. Wnt-1 Signaling Inhibits Apoptosis by Activating β-Catenin/T Cell Factor-Mediated Transcription. J Cell Biol (2001) 152:87–96. doi: 10.1083/jcb.152.1.87
32. Malanchi I, Santamaria-Martínez A, Susanto E, Peng H, Lehr HA, Delaloye JF, et al. Interactions Between Cancer Stem Cells and Their Niche Govern Metastatic Colonization. Nature (2012) 481:85–91. doi: 10.1038/nature10694
33. Roy S, Kar M, Roy S, Padhi S, Kumar A, Thakur S, et al. Inhibition of CD44 Sensitizes Cisplatin-Resistance and Affects Wnt/β-Catenin Signaling in HNSCC Cells. Int J Biol Macromol (2020) 149:501–12. doi: 10.1016/j.ijbiomac.2020.01.131
34. Xie S-L, Fan S, Zhang S-Y, Chen W-X, Li Q-X, Pan G-K, et al. SOX8 Regulates Cancer Stem-Like Properties and Cisplatin-Induced EMT in Tongue Squamous Cell Carcinoma by Acting on the Wnt/β-Catenin Pathway. Int J Cancer (2018) 142:1252–65. doi: 10.1002/ijc.31134
35. Moens S, Zhao P, Baietti MF, Marinelli O, Van Haver D, Impens F, et al. The Mitotic Checkpoint Is a Targetable Vulnerability of Carboplatin-Resistant Triple Negative Breast Cancers. Sci Rep (2021) 11:1–13. doi: 10.1038/s41598-021-82780-6
36. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon Provides Fast and Bias-Aware Quantification of Transcript Expression. Nat Methods (2017) 14:417–9. doi: 10.1038/nmeth.4197
37. Love MI, Huber W, Anders S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol (2014) 15:1–21. doi: 10.1186/s13059-014-0550-8
38. Stephens M. False Discovery Rates: A New Deal. Biostatistics (2017) 18:275–94. doi: 10.1093/biostatistics/kxw041
39. Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: A Visual Analytics Platform for Comprehensive Gene Expression Profiling and Meta-Analysis. Nucleic Acids Res (2019) 47:W234–41. doi: 10.1093/nar/gkz240
40. Müller AC, Giambruno R, Weißer J, Májek P, Hofer A, Bigenzahn JW, et al. Pathway Enrichment Analysis and Visualization of Omics Data Using G:Profiler, GSEA, Cytoscape and EnrichmentMap. Nat Protoc (2019) 22:924–34. doi: 10.1038/s41596-018-0103-9
41. Yoon S, Kim SY, Nam D. Improving Gene-Set Enrichment Analysis of RNA-Seq Data With Small Replicates. PloS One (2016) 11:1–16. doi: 10.1371/journal.pone.0165919
42. Bankhead P, Loughrey MB, Fernández JA, Dombrowski Y, McArt DG, Dunne PD, et al. QuPath: Open Source Software for Digital Pathology Image Analysis. Sci Rep (2017) 7:1–7. doi: 10.1038/s41598-017-17204-5
43. Schrauwen S, Coenegrachts L, Cattaneo A, Hermans E, Lambrechts D, Amant F. The Antitumor Effect of Metformin With and Without Carboplatin on Primary Endometrioid Endometrial Carcinoma In Vivo. Gynecol Oncol (2015) 138:378–82. doi: 10.1016/j.ygyno.2015.06.006
44. Depreeuw J, Hermans E, Schrauwen S, Annibali D, Coenegrachts L, Thomas D, et al. Characterization of Patient-Derived Tumor Xenograft Models of Endometrial Cancer for Preclinical Evaluation of Targeted Therapies. Gynecol Oncol (2015) 139:118–26. doi: 10.1016/j.ygyno.2015.07.104
45. Herath NI, Berthault N, Thierry S, Jdey W, Lienafa MC, Bono F, et al. Preclinical Studies Comparing Efficacy and Toxicity of DNA Repair Inhibitors, Olaparib, and AsiDNA, in the Treatment of Carboplatin-Resistant Tumors. Front Oncol (2019) 9:1097. doi: 10.3389/fonc.2019.01097
46. Engelmann K, Shen H, Finn OJ. MCF7 Side Population Cells With Characteristics of Cancer Stem/Progenitor Cells Express the Tumor Antigen MUC1. Cancer Res (2008) 68:2419–26. doi: 10.1158/0008-5472.CAN-07-2249
47. Koopman G, Reutelingsperger CP, Kuijten GA, Keehnen RM, Pals ST, van Oers MH. Annexin V for Flow Cytometric Detection of Phosphatidylserine Expression on B Cells Undergoing Apoptosis. Blood (1994) 84:1415–20. doi: 10.1182/blood.V84.5.1415.1415
48. Croker AK, Allan AL. Inhibition of Aldehyde Dehydrogenase (ALDH) Activity Reduces Chemotherapy and Radiation Resistance of Stem-Like ALDH Hicd44 + Human Breast Cancer Cells. Breast Cancer Res Treat (2012) 133:75–87. doi: 10.1007/s10549-011-1692-y
49. Yang L, Xie X, Tang H, Kong Y, Xie X, Chen J, et al. LGR5 Promotes Breast Cancer Progression and Maintains Stem-Like Cells Through Activation of Wnt/β-Catenin Signaling. Stem Cells (2015) 33:2913–24. doi: 10.1002/stem.2083
50. Ponti D, Costa A, Zaffaroni N, Pratesi G, Petrangolini G, Coradini D, et al. Isolation and In Vitro Propagation of Tumorigenic Breast Cancer Cells With Stem/Progenitor Cell Properties. Cancer Res (2005) 65:5506–11. doi: 10.1158/0008-5472.CAN-05-0626
51. Fuerer C, Nusse R. Lentiviral Vectors to Probe and Manipulate the Wnt Signaling Pathway. PloS One (2010) 5:e9370. doi: 10.1371/journal.pone.0009370
52. Doble BW, Woodgett JR. GSK-3: Tricks of the Trade for a Multi-Tasking Kinase. J Cell Sci (2003) 116:1175–86. doi: 10.1242/jcs.00384
53. Guo W, Keckesova Z, Donaher JL, Shibue T, Tischler V, Reinhardt F, et al. Slug and Sox9 Cooperatively Determine the Mammary Stem Cell State. Cell (2012) 148:1015–28. doi: 10.1016/j.cell.2012.02.008
54. Barrott JJ, Cash GM, Smith AP, Barrow JR, Murtaugh LC. Deletion of Mouse Porcn Blocks Wnt Ligand Secretion and Reveals an Ectodermal Etiology of Human Focal Dermal Hypoplasia/Goltz Syndrome. Proc Natl Acad Sci USA (2011) 108:12752–7. doi: 10.1073/pnas.1006437108
55. Lian X, Hsiao C, Wilson G, Zhu K, Hazeltine LB, Azarin SM, et al. Robust Cardiomyocyte Differentiation From Human Pluripotent Stem Cells via Temporal Modulation of Canonical Wnt Signaling. Proc Natl Acad Sci (2012) 109:E1848–57. doi: 10.1073/pnas.1200250109
56. Bruna A, Rueda OM, Greenwood W, Serra V, Garnett MJ. Caldas Correspondence C. A Biobank of Breast Cancer Explants With Preserved Intra-Tumor Heterogeneity to Screen Anticancer Compounds In Brief. Cell (2016) 167:260–74. doi: 10.1016/j.cell.2016.08.041
57. Gao H, Korn JM, Ferretti S, Monahan JE, Wang Y, Singh M, et al. High-Throughput Screening Using Patient-Derived Tumor Xenografts to Predict Clinical Trial Drug Response. Nat Med (2015) 21:1318–25. doi: 10.1038/nm.3954
58. Woo XY, Giordano J, Srivastava A, Zhao ZM, Lloyd MW, de Bruijn R, et al. Conservation of Copy Number Profiles During Engraftment and Passaging of Patient-Derived Cancer Xenografts. Nat Genet (2021) 53:86–99. doi: 10.1038/s41588-020-00750-6
59. Birkbak NJ, Li Y, Pathania S, Greene-Colozzi A, Dreze M, Bowman-Colin C, et al. Overexpression of BLM Promotes DNA Damage and Increased Sensitivity to Platinum Salts in Triple-Negative Breast and Serous Ovarian Cancers. Ann Oncol (2018) 29:903–9. doi: 10.1093/annonc/mdy049
60. Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, et al. G:Profiler: A Web Server for Functional Enrichment Analysis and Conversions of Gene Lists (2019 Update). Nucleic Acids Res (2019) 47:W191–8. doi: 10.1093/nar/gkz369
61. Jönsson JM, Johansson I, Dominguez-Valentin M, Kimbung S, Jönsson M, Bonde JH, et al. Molecular Subtyping of Serous Ovarian Tumors Reveals Multiple Connections to Intrinsic Breast Cancer Subtypes. PloS One (2014) 9:1–11. doi: 10.1371/journal.pone.0107643
62. Woodward WA, Chen MS, Behbod F, Alfaro MP, Buchholz TA, Rosen JM. WNT/beta-Catenin Mediates Radiation Resistance of Mouse Mammary Progenitor Cells. Proc Natl Acad Sci USA (2007) 104:618–23. doi: 10.1073/pnas.0606599104
63. Martins-Neves SR, Paiva-Oliveira DI, Wijers-Koster PM, Abrunhosa AJ, Fontes-Ribeiro C, Bovée JVMG, et al. Chemotherapy Induces Stemness in Osteosarcoma Cells Through Activation of Wnt/β-Catenin Signaling. Cancer Lett (2016) 370:286–95. doi: 10.1016/j.canlet.2015.11.013
64. Jang G-B, Kim J-Y, Cho S-D, Park K-S, Jung J-Y, Lee H-Y, et al. Blockade of Wnt/β-Catenin Signaling Suppresses Breast Cancer Metastasis by Inhibiting CSC-Like Phenotype. Sci Rep (2015) 5:12465. doi: 10.1038/srep12465
65. Xie W, Liu N, Wang X, Wei L, Xie W, Sheng X. Wilms’ Tumor 1-Associated Protein Contributes to Chemo-Resistance to Cisplatin Through the Wnt/b -Catenin Pathway in Endometrial Cancer. BioMed Res Inter (2021) 11:1–13. doi: 10.3389/fonc.2021.598344
66. He S, Wang W, Wan Z, Shen H, Zhao Y, You Z, et al. FAM83B Inhibits Ovarian Cancer Cisplatin Resistance Through Inhibiting Wnt Pathway. Oncogenesis (2021) 10. doi: 10.1038/s41389-020-00301-y
67. Mohapatra P, Shriwas O, Mohanty S, Ghosh A, Smita S, Kaushik SR, et al. CMTM6 Drives Cisplatin Resistance by Regulating Wnt Signaling Through ENO-1/AKT/Gsk3β Axis. JCI Insight (2021) 6:1–21. doi: 10.1172/jci.insight.143643
68. Geyer FC, Lacroix-Triki M, Savage K, Arnedos M, Lambros MB, MacKay A, et al. β-Catenin Pathway Activation in Breast Cancer Is Associated With Triple-Negative Phenotype But Not With CTNNB1 Mutation. Mod Pathol (2011) 24:209–31. doi: 10.1038/modpathol.2010.205
69. Liao S, Gan L, Qin W, Liu C, Mei Z. Inhibition of GSK3 and MEK Induced Cancer Stem Cell Generation via the Wnt and MEK Signaling Pathways. Oncol Rep (2018) 40:2005–13. doi: 10.3892/or.2018.6600
70. Yang Y, Wang QQ, Bozinov O, Xu RX, Sun YL, Wang SS. GSK–3 Inhibitor CHIR99021 Enriches Glioma Stem–Like Cells. Oncol Rep (2020) 43:1479–90. doi: 10.3892/or.2020.7525
71. Yoshino Y, Ishioka C. Inhibition of Glycogen Synthase Kinase-3 Beta Induces Apoptosis and Mitotic Catastrophe by Disrupting Centrosome Regulation in Cancer Cells. Sci Rep (2015) 5:1–14. doi: 10.1038/srep13249
72. Dewi FRP, Domoto T, Hazawa M, Kobayashi A, Douwaki T, Minamoto T, et al. Colorectal Cancer Cells Require Glycogen Synthase Kinase-3β for Sustaining Mitosis via Translocated Promoter Region (TPR)- Dynein Interaction. Oncotarget (2018) 9:13337–52. doi: 10.18632/oncotarget.24344
73. O’Flaherty L, Shnyder SD, Cooper PA, Cross SJ, Wakefield JG, Pardo OE, et al. Tumor Growth Suppression Using a Combination of Taxol-Based Therapy and GSK3 Inhibition in Non-Small Cell Lung Cancer. PloS One (2019) 14:1–17. doi: 10.1371/journal.pone.0214610
74. Glibo M, Serman A, Karin-Kujundzic V, Vlatkovic IB, Miskovic B, Vranic S, et al. The Role of Glycogen Synthase Kinase 3 (GSK3) in Cancer With Emphasis on Ovarian Cancer Development and Progression: A Comprehensive Review. Bosn J Basic Med Sci (2021) 21:5–18. doi: 10.17305/bjbms.2020.5036
75. Ma C, Wang J, Gao Y, Gao TW, Chen G, Bower KA, et al. The Role of Glycogen Synthase Kinase 3β in the Transformation of Epidermal Cells. Cancer Res (2007) 67:7756–64. doi: 10.1158/0008-5472.CAN-06-4665
76. Zhou W, Wang L, Gou Sm, Wang Tl, Zhang M, Liu T, et al. ShRNA Silencing Glycogen Synthase Kinase-3 Beta Inhibits Tumor Growth and Angiogenesis in Pancreatic Cancer. Cancer Lett (2012) 316:178–86. doi: 10.1016/j.canlet.2011.10.033
77. Fillmore CM, Kuperwasser C. Human Breast Cancer Cell Lines Contain Stem-Like Cells That Self-Renew, Give Rise to Phenotypically Diverse Progeny and Survive Chemotherapy. Breast Cancer Res (2008) 10:1–13. doi: 10.1186/bcr1982
78. Su Y-J, Chang Y-W, Lin W-H, Liang C-L, Lee J-L. An Aberrant Nuclear Localization of E-Cadherin is a Potent Inhibitor of Wnt/β-Catenin-Elicited Promotion of the Cancer Stem Cell Phenotype. Oncogenesis (2015) 4:e157–7. doi: 10.1038/oncsis.2015.17
79. Martin-Orozco E, Sanchez-Fernandez A, Ortiz-Parra I, Ayala-San Nicolas M. WNT Signaling in Tumors: The Way to Evade Drugs and Immunity. Front Immunol (2019) 10:2854. doi: 10.3389/fimmu.2019.02854
80. Bejsovec A. Wingless Signaling: A Genetic Journey From Morphogenesis to Metastasis. Genetics (2018) 208:1311–36. doi: 10.1534/genetics.117.300157
81. Liu J, Pan S, Hsieh MH, Ng N, Sun F, Wang T, et al. Targeting Wnt-Driven Cancer Through the Inhibition of Porcupine by LGK974. Proc Nat Acad Sci (2013) 110:20224–9. doi: 10.1073/pnas.1314239110
82. Novartis. A Study of LGK974 in Patients With Malignancies Dependent on Wnt Ligands. ClinicalTrials.gov identifier: NCT01351103 (2021), U.S. National Library of Medicine. Available at: https://clinicaltrials.gov/ct2/show/NCT01351103?id=NCT01351103&draw=2&rank=1&load=cart.
83. Fischer MM, Yeung VP, Cattaruzza F, Hussein R, Yen WC, Murriel C, et al. RSPO3 Antagonism Inhibits Growth and Tumorigenicity in Colorectal Tumors Harboring Common Wnt Pathway Mutations. Sci Rep (2017) 7:1–9. doi: 10.1038/s41598-017-15704-y
84. Gurney A, Axelrod F, Bond CJ, Cain J, Chartier C, Donigan L, et al. Wnt Pathway Inhibition via the Targeting of Frizzled Receptors Results in Decreased Growth and Tumorigenicity of Human Tumors. Proc Natl Acad Sci USA (2012) 109:11717–22. doi: 10.1073/pnas.1120068109
85. Fukumoto T, Zhu H, Nacarelli T, Karakashev S, Fatkhutdinov N, Wu S, et al. N6-Methylation of Adenosine of FZD10 mRNA Contributes to PARP Inhibitor Resistance. Cancer Res (2019) 79:2812–20. doi: 10.1158/0008-5472.CAN-18-3592
86. Yamamoto TM, McMellen A, Watson ZL, Aguilera J, Ferguson R, Nurmemmedov E, et al. Activation of Wnt Signaling Promotes Olaparib Resistant Ovarian Cancer. Mol Carcinog (2019) 58:1770–82. doi: 10.1002/mc.23064
Keywords: triple negative breast cancer, WNT pathway, platinum-resistance, cancer stem cells, patient-derived xenograft models
Citation: Abreu de Oliveira WA, Moens S, El Laithy Y, van der Veer BK, Athanasouli P, Cortesi EE, Baietti MF, Koh KP, Ventura J-J, Amant F, Annibali D and Lluis F (2021) Wnt/β-Catenin Inhibition Disrupts Carboplatin Resistance in Isogenic Models of Triple-Negative Breast Cancer. Front. Oncol. 11:705384. doi: 10.3389/fonc.2021.705384
Received: 05 May 2021; Accepted: 28 June 2021;
Published: 22 July 2021.
Edited by:
Daniela Spano, Institute of Biochemistry and Cell Biology, (CNR), ItalyReviewed by:
Laura Cerchia, Istituto per l’Endocrinologia e l’oncologia “Gaetano Salvatore (CNR), ItalyRadhika Nair, Rajiv Gandhi Centre for Biotechnology, India
Copyright © 2021 Abreu de Oliveira, Moens, El Laithy, van der Veer, Athanasouli, Cortesi, Baietti, Koh, Ventura, Amant, Annibali and Lluis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Frederic Lluis, frederic.lluisvinas@kuleuven.be
†These authors have contributed equally to this work
‡These authors share last authorship