Next Article in Journal
Circulating Tumor DNA: The Dawn of a New Era in the Optimization of Chemotherapeutic Strategies for Metastatic Colo-Rectal Cancer Focusing on RAS Mutation
Next Article in Special Issue
Knockdown of UBQLN1 Functions as a Strategy to Inhibit CRC Progression through the ERK-c-Myc Pathway
Previous Article in Journal
Patient Preferences for Unresectable Hepatocellular Carcinoma Treatments: A Discrete-Choice Experiment
Previous Article in Special Issue
Microbiome Profiling from Fecal Immunochemical Test Reveals Microbial Signatures with Potential for Colorectal Cancer Screening
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hedgehog-GLI and Notch Pathways Sustain Chemoresistance and Invasiveness in Colorectal Cancer and Their Inhibition Restores Chemotherapy Efficacy

1
Department of Experimental Medicine, Sapienza University of Rome, Viale Regina Elena 324, 00161 Rome, Italy
2
Department of Molecular Medicine, Sapienza University of Rome, Viale Regina Elena 291, 00161 Rome, Italy
3
Faculty of Medicine and Surgery, Kore University of Enna, Cittadella Universitaria, 94100 Enna, Italy
4
Soft and Living Matter Laboratory, Institute of Nanotechnology, Consiglio Nazionale Delle Ricerche, Piazzale Aldo Moro 5, 00185 Rome, Italy
5
Center for Life Nano- and Neuro-Science, Istituto Italiano di Tecnologia, Viale Regina Elena 291, 00161 Rome, Italy
6
D-TAILS srl, Viale Regina Elena 291, 00161 Rome, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2023, 15(5), 1471; https://doi.org/10.3390/cancers15051471
Submission received: 22 January 2023 / Revised: 16 February 2023 / Accepted: 22 February 2023 / Published: 25 February 2023
(This article belongs to the Special Issue Molecular Biology of Colorectal Cancers)

Abstract

:

Simple Summary

Colorectal cancer is a leading cause of cancer-related deaths, mainly caused by resistance to therapy and metastatic spread, in turn sustained by the activation of mechanisms such as the epithelial-to-mesenchymal transition (EMT). We investigate here the role of the Hedgehog-GLI and NOTCH signaling pathways, already associated with poor prognosis in CRC, in the mechanism of chemoresistance and EMT, using monolayer and organoids from two models of common mutations in CRC: KRAS and BRAF. Our results show that treatment with the chemotherapeutic drug 5-fluorouracil activated both pathways in the investigated contexts. However, we observed a different behavior in the investigated models: in KRAS-mutated CRC, the inhibition of both the HH-GLI and NOTCH pathways is necessary to enhance chemosensitivity, while in BRAF-mutated CRC the inhibition of HH-GLI is sufficient to impair both signaling pathways and promote chemosensitivity.

Abstract

Colorectal cancer (CRC) is a leading cause of cancer-related mortality and chemoresistance is a major medical issue. The epithelial-to-mesenchymal transition (EMT) is the primary step in the emergence of the invasive phenotype and the Hedgehog-GLI (HH-GLI) and NOTCH signaling pathways are associated with poor prognosis and EMT in CRC. CRC cell lines harboring KRAS or BRAF mutations, grown as monolayers and organoids, were treated with the chemotherapeutic agent 5-Fluorouracil (5-FU) alone or combined with HH-GLI and NOTCH pathway inhibitors GANT61 and DAPT, or arsenic trioxide (ATO) to inhibit both pathways. Treatment with 5-FU led to the activation of HH-GLI and NOTCH pathways in both models. In KRAS mutant CRC, HH-GLI and NOTCH signaling activation co-operate to enhance chemoresistance and cell motility, while in BRAF mutant CRC, the HH-GLI pathway drives the chemoresistant and motile phenotype. We then showed that 5-FU promotes the mesenchymal and thus invasive phenotype in KRAS and BRAF mutant organoids and that chemosensitivity could be restored by targeting the HH-GLI pathway in BRAF mutant CRC or both HH-GLI and NOTCH pathways in KRAS mutant CRC. We suggest that in KRAS-driven CRC, the FDA-approved ATO acts as a chemotherapeutic sensitizer, whereas GANT61 is a promising chemotherapeutic sensitizer in BRAF-driven CRC.

Graphical Abstract

1. Introduction

Colorectal cancer (CRC) is the third most frequent cancer and the second cause of cancer-related death worldwide [1]. Mutations in KRAS and BRAF oncogenes represent the most common genetic drivers in CRC. Indeed, KRAS and BRAF mutations occur in 40% and 10% of CRC, respectively [2], and they are both associated with a poor outcome [3]. Even though they both belong to the MAPK pathway, KRAS and BRAF mutations are mutually exclusive in CRC and these two types of cancer are characterized by distinct clinical and molecular features. BRAF-mutant CRC often displays genome-wide hypermethylation, high microsatellite instability and mutation rates, while KRAS mutant CRC is associated with lower levels of microsatellite instability and gene methylation [2].
First-line and palliative treatments for metastatic CRC, bearing KRAS or BRAF mutations, include the cytotoxic chemotherapeutic agent 5-fluorouracil (5-FU) [1,4]; however, patients often present disease recurrence after 5-FU therapy [5]. Chemoresistance is conferred by a plethora of mechanisms, including the modulation of signaling pathways involved in the emergence of the cancer stem features and epithelial-to-mesenchymal transition (EMT) [6]. Other mechanisms for resistance to therapy include the inhibition of apoptosis driven by upregulation of autophagy [7], metabolic reprogramming [8], upregulation of molecules involved in drug efflux and drug metabolism and activation of alternative pathways [9].
Hedgehog-GLI (HH-GLI) and NOTCH signaling are pivotal developmental pathways involved in the regulation of multiple biological and pathological processes. The canonical HH-GLI pathway is activated upon the interaction between the extracellular ligands Shh, Ihh and Dhh and the receptor Patched (PTCH), which in turn derepresses Smoothened (Smo), thus activating the transcription factors GLI1, GLI2 and GLI3. Activated GLI translocate into the nucleus where they bind to DNA and activate the transcription of target genes [10]. In cancers, GLI1 can also be activated in a non-canonical way by the “oncogenic load” of the cancer cell [11]. NOTCH cascade is activated upon binding of ligands Jag1, Jag2, Dll1, Dll3 and Dll4 to NOTCH receptors (from 1 to 4). The binding leads to proteolytic cleavages of the NOTCH receptor, releasing the NOTCH intracellular domain (ID) into the cytoplasm. Then, NOTCH ID migrates into the nucleus where, in complex with CBF1 (also known as RPBJ), it activates its transcriptional program [12]. Downstream target genes include HES1, which is involved in EMT and transcriptionally regulates ATP-binding cassettes transporters (ABC transporters), involved in multidrug resistance [13]. Interestingly, deregulation of the NOTCH pathway was described in numerous cancerous and non-cancerous diseases, with its role being highly context-dependent [14].
Deregulated HH-GLI is involved in the development and maintenance of numerous cancers [10] and, together with NOTCH signaling, plays a crucial role in the maintenance of stem cells of the intestinal epithelia [15]. The crosstalk of HH-GLI and NOTCH signaling is fundamental for spinal cord patterning [16], and several previous reports highlighted how several molecules belonging to the NOTCH pathway regulate the key components of the HH-GLI pathway and vice versa, as reviewed Kumar et al. [17].
Both HH-GLI and NOTCH pathways were described as deregulated and associated with poor prognosis in CRC [18,19]. In this context, our previous work has described a chemoresistance mechanism operated by the HH-GLI signaling in CRC, where chemotherapy treatment resulted in aberrant activation of the HH-GLI pathway which in turn led to the transcription of ATP-binding cassette transporters (ABC transporters), involved in multidrug resistance [20].
Therefore, our current work aimed to evaluate the role of HH-GLI and NOTCH signaling pathways as regulatory molecular mechanisms responsible for chemotherapy resistance in models of KRAS- or BRAF-driven CRC.

2. Materials and Methods

2.1. Cell Cultures and Treatments

HCT116 (KRAS G13D mutant) and HT29 (BRAF V600E mutant) were obtained from American Type Culture Collection (ATCC) and grown in DMEM high glucose (supplemented with 10% (v/v) fetal bovine serum, 1% (v/v) penicillin (50 U mL−1)—streptomycin (50 U mL−1)—and 2 mM L-glutamine. Cells were routinely checked for mycoplasma contamination by testing with PCR Mycoplasma Detection Kit (Cat. G238, ABM, Richmond, BC, Canada).
Cells were treated with 10 μM GANT61 (ENZO Lifesciences, New York, NY, USA), 10 μM DAPT (Merk Life Science S.r.l., Milan, Italy), 10 μM Arsenic Trioxide (ATO) (Merck, Merk Life Science S.r.l., Milan, Italy) and 10 μM 5-Fluorouracil (5-FU).
For combined treatments, GANT61 and DAPT or Arsenic Trioxide (ATO) were administered to the cells 24 h before 5-FU.

2.2. Cell Viability by Trypan Blue Exclusion Assay

Cell proliferation was assessed by trypan blue dye exclusion test using 0.4% (w/v) Trypan Blue solution (Merk Life Science S.r.l., Milan, Italy). Blue-stained cells were scored as non-viable and unstained cells were scored as viable cells. The percentage of viable cells was obtained as the ratio between the percentage of viable cells in treated cells versus control.

2.3. Transwell Invasion Assay

Transwell invasion assay was performed using Corning® Transwell® chambers (8 μm pore size, Corning®). HCT116 and HT29 cells (2.5 × 104 in each well) were seeded in the upper chambers of the 48-well plates (Corning, Somerville, MA, USA) while lower chambers were filled with 1 mL of medium with indicated treatments. Cells in the lower chambers were fixed with 95% ethanol for 10 min, stained with crystal violet and counted.

2.4. Western Blot

Cells were lysed as previously described [20]. Lysates were separated on 8% acrylamide gel and immunoblotted using standard procedures [21]. Primary antibodies were Anti-GLI1 (L42B10, Cell Signalling Technology Inc., Boston, MA, USA), anti-PARP p85 Fragment (G7341, Promega, Madison, WI, USA) and anti-Cleaved NOTCH1 (D3B8, Cell Signalling Technology Inc., Boston, MA, USA). HRP-conjugated secondary antisera (Santa Cruz Biotechnology, Shanghai, China) were used, followed by enhanced chemiluminescence (ECL Amersham, Merk Life Science S.r.l., Milan, Italy).

2.5. RNA Isolation and Real Time qPCR

cDNA was obtained as described earlier [20]. RNA expression was analyzed on cDNAs using the ViiA™ 7 Real-Time PCR System, SensiFAST™ Probe Lo-ROX (Bioline, Memphis, TN, USA), TaqMan gene expression assay according to the manufacturer’s instructions (Life Technologies, Waltham, MA, USA). mRNA quantification was expressed in arbitrary units, as the ratio of the sample quantity to the calibrator or to the mean values of control samples. All values were normalized to three endogenous controls: HPRT, GAPDH and β-ACTIN.
Primers for gene expression are listed in Supplementary Table S1. Gene expression of GLI1, HES1, c-MET, ABCG2, CD133, KRAS, BRAF, HPRT, GAPDH and β-ACTIN was assessed using Life technologies “best coverage” assays (Life Technologies, Waltham, MA, USA).

2.6. Organoids

Organoids were produced by seeding 1500 cells per well. Cells were mixed with 33% growth-factor-reduced phenol red-free Matrigel (Corning, Somerville, MA, USA). Cultures were grown using a flat-bottom 24-well microplate in advanced DMEM-F12 (Cat. 12634010, Gibco, Waltham, MA, USA) supplemented with Epidermal growth factor and Fibroblast growth factor both at final concentrations of 20 ng/µL.
For in vivo live imaging experiments, GFP-labelled organoids were obtained by transducing HCT116 with PLKO lentiviral particles carrying pTWEEN-GFP vector.
Transduced green fluorescent cells were selected by cell sorting and used for organoids production.

2.7. Whole Mount Immunofluorescence

Organoids were fixed with 4% paraformaldehyde and permeabilized with Triton X-100 in PBS (Sigma-Aldrich, St. Louis, MO, USA). Organoids were stained with anti-vimentin (ab11256, ABCAM, Cambridge, UK) antibody. Nuclei were DAPI-counterstained. Phalloidin was used for f-actin staining. Images were acquired using an LSM 900 (Zeiss, Milan, Italy) laser scanning confocal microscope with 40×/0.75 NA objective. Images were analyzed by using the program Zeiss ZEN 2.3 blue edition (https://www.zeiss.com/microscopy/int/products/microscope-software/zen-lite.html (accessed on 10 September 2022)).

2.8. Datasets and In Silico Analyses

Datasets available on R2 platform (https://hgserver1.amc.nl/cgi-bin/r2/main.cgi (accessed on 15 December 2021)) were interrogated to evaluate GLI1 and NOTCH1 correlation in patients carrying BRAF or KRAS mutations. In detail, Tumor Colon Mutation status (Core Exon)—Sieber—211—rma_sketch—huex10p investigated gene correlation between GENE/REPORTER1: GLI1 and GENE/REPORTER2: NOTCH1, in 29 samples of CRC-carrying braf_v600e mutation; Tumor Colon (after surgery)—Beissbarth—363—custom—4hm44k investigated gene correlation between GENE/REPORTER1: GLI1 and GENE/REPORTER2: NOTCH1, in 32 samples of CRC carrying kras_g13d mutation.

2.9. Single-Particle Tracking Analysis

Single-particle tracking (SPT) diffusibility analysis was performed in five steps. In the first step, single cells were detected from the time-lapse movies of oHCT116-treated and control group organoids using Imaris spot model. Spots were taken in each frame and were linked to the spots corresponding to the same cell in the successive frame. Frame-to-frame tracking was implemented using the linear assignment problem (LAP) method [22,23]. In the third step, MSD, the mean square distance travelled by a cell given a certain time interval (see Supplementary Figure S3), was computed from single trajectories, as described by Michalet X [24]. In the fourth step, the diffusion parameter D was calculated for each tracked cell. To this end, the MSD was plotted for different time intervals (Δt) for each cell trajectory and the slope was computed using the Least Squares Method. In the fifth step, the Kolmogorov–Smirnov test was applied on the diffusion parameters D, obtained from 5-FU-treated oHCT116 and control group organoids.

2.10. Statistical Analysis

Results are representative of at least three independent experiments and are expressed as means +/− SD. Differences were analyzed using One-way ANOVA and Two-way ANOVA tests where appropriate, using the GraphPad Prism software Version 8.0. Adjusted p-values of less than 0.05 were considered as statistically significant.

3. Results

3.1. HH-GLI and NOTCH Signaling Pathways Sustain Resistance to 5-FU in KRAS Mutant CRC Cells

5-fluorouracil (5-FU) is a chemotherapeutic agent used for adjuvant and palliative treatment of CRC; however, patients often present disease recurrence [5]. Therefore, we evaluated the role of HH-GLI and NOTCH signaling pathways as molecular mechanisms responsible for chemotherapy resistance.
KRAS mutant HCT116 CRC cells were treated with 5-FU, alone or in combination with the HH-GLI inhibitor GANT61 and/or the NOTCH inhibitor DAPT. Our results showed that GLI1 and NOTCH1 ID were significantly upregulated after 5-FU treatment (Figure 1A). HH-GLI inhibition by GANT61 resulted in the downregulation of GLI1 and, interestingly, in the upregulation of NOTCH1 ID; vice versa, NOTCH inhibition by DAPT resulted in the downregulation of NOTCH1 ID and the upregulation of GLI1 (Figure 1A).
We then determined the effects of treatments with the combination of HH-GLI inhibitor, NOTCH inhibitor and 5-FU. NOTCH1 ID expression was impaired in all combined treatments that included DAPT (DAPT+GANT61, DAPT+5-FU and DAPT+GANT61+5FU), while it was unaffected by the combination of 5-FU+GANT61.
On the other hand, GLI1 was downregulated by GANT61 alone, as well as when combined with 5-FU and 5-FU+DAPT, while the combination of NOTCH inhibition and 5-FU failed to inhibit GLI1.
Overall, our results show that GLI1 and NOTCH1 ID levels were concomitantly significantly downregulated only after the combined treatment of the chemotherapeutic agent 5-FU together with the inhibition of HH-GLI and NOTCH.
In addition, we observed that the combination of HH-GLI and NOTCH pathway inhibition prevents the GLI1 upregulation and NOTCH1 activation induced by 5-FU.
To determine the effects of treatments on apoptosis, levels of cleaved PARP (c-PARP) were evaluated; our results show that c-PARP was significantly induced by the combination of 5-FU and GANT61 or DAPT, and the three drugs combined (Figure 1A).
We further analyzed the effects of treatments on cell viability, and we found it significantly impaired in cells treated with the combination of GANT61 and DAPT and with the combination of 5-FU with either GANT61 or DAPT or the combination of the three drugs (Figure 1B).
To discern potential interdependence between the HH-GLI and NOTCH signaling pathways in mutant KRAS CRC cells, we analyzed GLI1 and NOTCH1 levels in an available cohort of CRC patients carrying this mutation (Tumor Colon (after surgery)—Beissbarth—363—custom—4hm44k; https://hgserver1.amc.nl/cgi-bin/r2/main.cgi, accessed on 15 December 2021) and no correlation was found (Figure 1C)
We therefore envisioned a model where the oncogenic force of the driver gene KRASG13D sustains both HH-GLI and NOTCH pathways and both pathways need to be targeted to achieve a successful impairment of cells after chemotherapy.
Hence, to clarify if the KRASG13D driver mutation sustained expression of GLI1 and NOTCH, we performed silencing of KRAS in HCT116 (Supplementary Figure S1A), which resulted in the significant downregulation of GLI1 and NOTCH1 ID protein levels (Supplementary Figure S1B). KRAS silencing was also accompanied by a significant downregulation of ABCG2 and HES1, target genes of HH-GLI and NOTCH1 ID, respectively (Supplementary Figure S1C).
Arsenic Trioxide (ATO) is an organic compound approved for the therapy of adult patients with acute promyelocytic leukemia [25] and was shown to successfully inhibit both HH-GLI and NOTCH pathways [26]. ATO’sability to inhibit both GLI1 and NOTCH ID levels was confirmed in the KRASG13D-driven CRC model (Figure 1D). As previously shown, 5-FU alone was able to upregulate both GLI1 and NOTCH1 ID, while the combination with ATO impaired both signaling pathways (Figure 1D). Cleaved-PARP levels showed that apoptosis was significantly increased by the combination of ATO and 5-FU, while we observed a non-significant trend in ATO-treated cells (Figure 1D).
A pivotal feature of CRC aggressiveness relies on the epithelial-to-mesenchymal transition (EMT), a process that includes the acquisition by cancer cells of properties including motility and migration, early steps in cancer invasion and metastasis.
Therefore, we investigated whether the targeting of HH-GLI and NOTCH could impair KRAS mutant CRC’s migratory ability.
We investigated the effects of the combined treatment of 5-FU and ATO on the migration ability of HCT116 cells. We observed that the migration was unaffected by 5-FU treatment, while it was impaired with ATO treatment and was completely abrogated after ATO plus 5-FU combined treatment (Figure 1E). Then, we evaluated epithelial differentiation through E-cadherin levels, which increased after the combined treatment of ATO and 5-FU (Supplementary Figure S1D).

3.2. HH-GLI Signaling Pathway Sustains Resistance to 5-FU in BRAF Mutant CRC Cells

BRAF V600E is the activating driving mutation in 10% of CRC and correlates with poor prognosis, however targeted therapy against the mutation was proven ineffective and first-line treatment includes cytotoxic chemotherapy [1]; thus, we investigated the role of the HH-GLI and NOTCH signaling pathways in 5-FU chemotherapy resistance in BRAF mutant HT29 cells.
HT29 cells were treated with 5-FU alone or in combination with the HH-GLI inhibitor GANT61 and the NOTCH1 inhibitor DAPT (Figure 2A). We observed that 5-FU induced upregulation of GLI1 and NOTCH1. GANT61 treatment resulted in the downregulation of both GLI1 and NOTCH1 ID, while DAPT treatment caused the downregulation only of NOTCH1 ID, without exerting any effect on GLI1 levels compared to control cells. The combination of GANT61 and DAPT successfully targeted both GLI1 and NOTCH1 ID. The combined treatment of GANT61 plus 5-FU was able to revert the 5-FU-induced upregulation of GLI1 and NOTCH1 ID, and the combined treatment of DAPT and 5-FU was able to revert the 5-FU-induced upregulation of NOTCH1 ID and partially of GLI1. Only when both HH-GLI and NOTCH pathways were inhibited together with 5-FU treatment were both GLI1 and NOTCH1 ID significantly downregulated (Figure 2A).
Apoptosis was evaluated through c-PARP levels; treatment with 5-FU and single inhibition of HH-GLI and NOTCH1 failed to induce apoptosis; c-PARP levels indeed increased only when cells were treated with GANT61 in combination with 5-FU, or with the combination of the three drugs (Figure 2A).
We then investigated cell viability and our results showed a significant impairment after GANT61 treatment, alone or in combination with 5-FU (Figure 2B).
Based on these results, chemotherapy resistance to apoptosis in BRAF V600E mutated cells seems to be driven by the HH-GLI signaling, which in turn sustains the activation of the NOTCH pathway. To gain more insight into the interdependence between the HH-GLI and NOTCH pathways, we interrogated GLI1 and NOTCH1 levels in a cohort of CRC patients carrying BRAFV600E mutation (Mutation status (Core Exon)—Sieber—211—rma_sketch—huex10p; https://hgserver1.amc.nl/cgi-bin/r2/main.cgi, accessed on 15 December 2021) and found a positive and significant correlation between GLI1 and NOTCH1 (Figure 2C). The above presented data suggest an upstream role of HH-GLI in the regulation of NOTCH signaling in the BRAF-driven CRC model.
To investigate whether BRAFV600E acted as a driver on the regulation of HH-GLI and NOTCH, we performed BRAF silencing (Supplementary Figure S2A). BRAF silencing resulted in decreased GLI1 and NOTCH1 ID protein levels (Supplementary Figure S2B).
We also evaluated mRNA levels of HH-GLI and NOTCH1 ID readout, ABCG2 and HES1, respectively, and both were significantly decreased after BRAF silencing (Supplementary Figure S2C).
The above-reported data demonstrate that chemotherapy stress induced increased levels of both HH-GLI and NOTCH1 pathways in the BRAF-driven CRC model. Interestingly, we observed that the GLI1 inhibitor GANT61 was also able to decrease NOTCH1 ID levels; conversely, the NOTCH1 inhibitor DAPT did not affect GLI1 levels.
Since we observed that the targeting of HH-GLI was able to indirectly also target the NOTCH pathway, we wondered if the combination of 5-FU and GANT61 could affect cell motility, a key feature of EMT and therefore of CRC aggressiveness.
Our experiments showed that 5-FU did not affect cell motility, while GANT61 resulted in decreased cell motility, which was further impaired by the combination of GANT61 with 5-FU (Figure 2D).
Then, we investigated the expression of two HT29 cell-specific epithelial differentiation markers, Axin and Muc2. We observed upregulation of Axin only after combined treatment, while Muc2 was affected by both 5-FU and GANT61 alone and by their combination (Supplementary Figure S2D), suggesting that treatments enhance the differentiated phenotype.

3.3. 5-FU Increases Motility of CRC Organoids

The previous set of experiments allowed us to point out the role of HH-GLI and NOTCH pathways as regulators of EMT in KRAS mutant and BRAF mutant CRC, a key feature of chemoresistance [27]. Organoid models in pre-clinical studies have become widespread due to their high reproducibility and high similarity to in vivo models [28,29]. Indeed, cell features and behavior depend on the architecture of the cell population, e.g., the cell–cell contact, the stiffness of the extracellular matrix and the interaction with the microenvironment. All these conditions concur with specific characteristics related to cell polarity, stemness and differentiation status.
Thus, to obtain CRC organoids, we seeded HCT116 and HT29 cells in Matrigel and after 7 days we observed organoid growth, as shown in Figure 3A,B. We compared basal levels of GLI1 and NOTCH1 in organoids and in 2D monolayer and our results reported higher GLI1 and HES1 expression levels in organoids, indicating that both pathways were more active in organoids compared to monolayer cellular models (Figure 3C). We then evaluated levels of the EMT marker c-MET in both organoids and monolayers and observed that c-MET was expressed at higher levels in organoids (Figure 3C). Since our results showed that 5-FU was not able to impair the migratory ability of CRC (Figure 1 and Figure 2), and that CRC patients often present disease progression despite chemotherapy, we wondered if 5-FU itself favored aggressiveness in organoids, unleashing the migratory potential.
Increased motility and migration capacity are features of EMT, thus we performed in vivo live cell imaging in the KRASG13D-driven CRC organoid model, the HCT116-derived organoids (oHCT116) at basal state and after 5-FU treatment (Supplementary Material Supplementary Video S1).
To investigate the behavior of CRC cells within organoids, we investigated the diffusion parameters that allow the motility of individual cells to be quantified. The diffusion parameters from the oHCT116 control or 5-FU-treated organoids are reported (Figure 3E) along with the single cell trajectories that were used for the calculation of the diffusion parameters (Supplementary Figure S3). Interestingly, 5-FU-treated oHCT116 cells mostly present lower diffusion parameters compared with CTRL (Figure 3D), with a long tail corresponding to a sub-group of cells presenting very high diffusion (Figure 3E). Based on these results, we believe that cells with augmented motility after chemotherapy represent a subset of aggressive cells able to initiate the metastatic process.

3.4. HH-GLI and NOTCH Inhibition Impairs 5-FU-Driven Mesenchymal Phenotype in KRASG13D-Driven CRC Organoids

We then proceeded to investigate the inhibition of HH-GLI and NOTCH by using ATO in combination with 5-FU in KRAS-driven CRC organoids, oHCT116.
Treatment with 5-FU alone did not affect organoid growth, while organoids treated with ATO were significantly smaller; the association of 5-FU and ATO further impaired organoid growth (Figure 4A). Expression levels of the EMT marker c-MET, cancer stemness markers ABCG2 and CD133a, which is both HH-GLI target and cancer stemness marker, were significantly decreased in the combined treatment of 5-FU and ATO (Figure 4B). Interestingly, ATO was able to counteract the 5-FU-driven upregulation of ABCG2.
Mesenchymal features were also investigated by the immunofluorescence of the EMT marker vimentin, whose levels increased after 5-FU treatment and were reduced when organoids were treated with ATO alone or in combination with 5-FU. Of note F-actin, revealed by phalloidin staining, underwent a marked rearrangement in 5-FU-treated oHCT116, where cells lost their pseudopodia, probably due to a modification in the cell polarity (Figure 4C).

3.5. HH-GLI Inhibition Impairs 5-FU-Driven Mesenchymal Phenotype in BRAFV600E-Driven CRC Organoids

We then investigated the effects of 5-FU alone or in combination with the HH-GLI blockade in the BRAFV600E-driven CRC organoids (oHT29).
Our results showed that the size of oHT29 treated with 5-FU did not differ from the control group, while organoids treated with GANT61 were smaller in size and the combination of 5-FU and GANT61 strongly impaired organoid growth (Figure 5A).
Gene expression analysis showed that the levels of cancer stem cell and EMT markers ABCG2, CD133 and c-MET significantly increased after chemotherapy treatment and were impaired by HH-GLI inhibition and the combination of 5-FU and GANT61 (Figure 5B).
To better investigate EMT, we performed whole-mount immunofluorescence staining for the mesenchymal marker vimentin and observed that vimentin levels were upregulated in 5-FU-treated organoids, they decreased with GANT61 and were strongly impaired in the combined treatment (Figure 5C).
Altogether, our experiments show that in KRAS-driven and BRAF-driven CRC, the HH-GLI and NOTCH pathways sustain the resistance to 5-FU through the activation of the EMT. Of note, ATO, the drug targeting both HH-GLI and NOTCH pathways, reverted the mesenchymal phenotype, therefore supporting the action of the chemotherapeutic drug.

4. Discussion

Despite recent advances in cancer therapy, CRC is still among the prevalent causes of cancer-related death [30]. Even though medical research has focused on identifying genetic mutations linked to CRC progression and tumor prognosis to improve patient treatment, drug resistance often occurs. One of the mechanisms conferring drug resistance is the misactivation of evolutionarily conserved pathways, such as Wingless (WNT) [31,32], phosphoinositide-3-kinase [33,34], extracellular signal-regulated kinase (ERK) [35,36], nuclear factor-κB (NF-κB) [37,38] and the Hedgehog-GLI (HH-GLI) signaling pathway [20]. The HH-GLI pathway has a crucial role in correct embryonic development and plays a role in the physiological maintenance of many tissues, including the colonic mucosa [39,40]. While canonical activation of the HH-GLI pathway transduces the signal through the Hedgehog/PTCH/SMO/GLI axis, non-canonical regulation of GLI is external to Hedgehog signaling. Of note, it was demonstrated that transforming growth factor-beta (TGF-β) [41], epidermal growth factor receptor (EGFR) [42], mitogen-activated protein kinases (MAPK) [11], β-arrestin [43] and WNT/β-catenin [44,45] were able to induce the expression of GLI, regardless of SMO activation. Since both canonical and non-canonical routes culminate with the activation of the GLI1 transcriptional program, GLI1 inhibition could be useful to prevent chemoresistance in cancer cells. Our group has previously demonstrated that HH-GLI signaling regulates the expression of ATP-binding cassette transporters (ABC transporters), which are correlated to multidrug resistance in cancer cells, providing a rationale for the consideration of the HH-GLI pathway as a therapeutic target in CRC [20]. NOTCH signaling has been reported to play a crucial role in the development of the normal mucosa [15] and its aberrant activation is related to carcinogenesis in CRC. HH-GLI and NOTCH signaling pathways together with the WNT and BMP pathways are responsible for the development of intestinal mucosa, which is the innermost layer of the colon. Stem cells, transit amplifying cells and terminally differentiated secretory cells or enterocytes, concur in the formation of the structural unit of the colon, known as the crypt of Lieberkuhn [46]. A recent paper showed that the HH-GLI blockade with GANT61 was able to inhibit NOTCH and WNT/β-catenin in cellular models of CRC [47]. Since the HH-GLI and NOTCH pathways play a fundamental role in the correct patterning of the colonic mucosa and HH-GLI is upregulated by chemotherapeutic stress, we wondered whether HH-GLI and NOTCH crosstalk could be involved in the resistance mechanism of CRC cells related to 5-FU chemotherapeutic stress.
The results of this study show how the HH-GLI and NOTCH pathways sustain CRC chemoresistance in different ways depending on the driver oncogene mutation. In detail, in KRASG13D-driven HCT116 cells we observed an upregulation of HH-GLI and NOTCH pathways after 5-FU and the inhibition of HH-GLI resulted in increased levels of NOTCH1 ID and vice versa (Figure 1A). These results, coupled with the interrogation of public datasets (Figure 1C) suggested that the HH-GLI and NOTCH signaling pathways are connected in a positive feedback loop aiming to escape apoptosis induced by 5-FU (Figure 6).
Importantly, the combined inhibition of HH-GLI and NOTCH was able to impair EMT, shown both as an impairment of transwell migration ability and with EMT markers in organoids (Figure 1E and Figure 4). ATO, which was used to target both HH-GLI and NOTCH pathways, has been approved by the FDA for the therapy of adult patients with acute promyelocytic leukemia (APL). A phase I trial investigating the co-administration of ATO and 5-FU/Leucovorin in patients with advanced/relapsed CRC showed that ATO was well tolerated and that in some patients it was associated with therapeutic response and increased survival; a later study investigated GLI1 levels in biopsies from the above-mentioned clinical trial and found that it resulted to be down-modulated after ATO administration. Of note, data on the mutational status of enrolled patients are not available [48,49]. In BRAFV600E-driven CRC, both pathways were upregulated after 5-FU treatment, but importantly GANT61 downregulated not only its specific target GLI1 but also NOTCH (Figure 2A), suggesting an upstream role of HH-GLI over the NOTCH pathway (Figure 6), thus explaining the positive correlation between these two signaling pathways (Figure 2C). Importantly, HH-GLI inhibition was able to impair EMT features, both in monolayer and organoids (Figure 2D and Figure 5).

5. Conclusions

In conclusion, our study describes for the first time two distinct models for KRAS- and BRAF-driven CRC where the HH-GLI and NOTCH signaling pathways play different roles in the chemoresistance and mesenchymal phenotype of CRC (Figure 6). Indeed, we described that in KRASG13D-driven CRC, chemotherapy resistance is directed by the concurrent activation of the HH-GLI and NOTCH pathways and the inhibition of both is crucial to revert the resistant phenotype. Conversely, in BRAFV600E-mutated CRC, the resistance to apoptosis induced by chemotherapy is mainly sustained by the HH-GLI signaling pathway. The implications of this novel information can be far-reaching if taken into consideration for the management of CRC patients, providing clinicians with further tools for the development of more effective treatment plans.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/cancers15051471/s1. Supplementary Table S1. List of gene expression primers for quantitative real-time PCR. Gene expression of GLI1, HES1, c-MET, ABCG2, CD133, KRAS, BRAF, HPRT, GAPDH and β-ACTIN was assessed using Life Technologies “best coverage” assays (Life Technologies). Supplementary Figure S1. (A) Quantitative real-time PCR of KRAS in HCT116 cells after KRAS silencing (siKRAS) and control group (siCTRL). (B) Western blot analysis of KRAS, GLI1 and NOTCH1 ID in HCT116 cells after KRAS silencing (siKRAS) and control group (siCTRL). (C) Quantitative real-time PCR of ABCG2 and HES1 in HCT116 cells after KRAS silencing (siKRAS) and control group (siCTRL). (D) mRNA levels of E-cadherin (ECAD) expressed in arbitrary units in HCT116 treated with 5-FU, ATO, combined treatment and control group (CTRL). Data are representative of three independent experiments * p  <  0.05 versus control; ** p < 0.01 versus 5-FU (Two-way ANOVA test). Supplementary Figure S2. (A) Quantitative real-time PCR of BRAF in HCT116 cells after BRAF silencing (siBRAF) and control group (siCTRL). (B) Western blot analysis of GLI1 and NOTCH1 ID in HT29 cells after BRAF silencing (siBRAF) and control group (siCTRL). (C) Quantitative real-time PCR of ABCG2 and HES1 in HT29 cells after BRAF silencing (siBRAF) and control group (siCTRL). (D) mRNA levels of AXIN and MUC2 expressed in arbitrary units, in HT29 treated with 5-FU, GANT61, the combined treatment and control group (CTRL). Data are representative of three independent experiments * p < 0.05; ** p < 0.01 (Two-way ANOVA test). Supplementary Figure S3. The plots show how the mean square displacement (MSD) changes for different time intervals (Δt) for each tracked single cell trajectory from either the HCT116 control group 3D organoids (CTRL) or the HCT116 5-FU-treated group (5-FU). Supplementary Figure S4. Uncropped full scan for Figure 1 (panels refer to main figure panel). Supplementary Figure S5. Uncropped full scan for Figure 2 (panels refer to main figure panel). Supplementary Video S1. time lapse of gfp transduced oHCT116 at basal state (oHCT116-gfp CTRL); Supplementary Video S2. treated with 5-FU (oHCT116-gfp 5-FU); length: 6 h.

Author Contributions

Conceptualization, A.C. and A.P.; methodology, G.G.; validation, G.C., S.T. and T.M.A.; formal analysis, Z.M.B., G.G., M.L. and L.M.; investigation, A.C., G.C., S.T., F.B., A.D.F., L.C., T.M.A., Z.S., M.A.V.; resources, M.L., E.D.S., M.A.V.; writing—original draft preparation, A.C., Z.M.B., A.P.; writing—review and editing, G.C., A.V., E.D.S., E.F. and L.M.; visualization, A.C., A.P. and L.M.; supervision, E.F., A.P.; project administration, A.P.; funding acquisition, M.L., F.B., A.V., E.F. and L.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Sapienza University of Rome (Ateneo Project 2021) (A.V.); Sapienza Ateneo 2021 RM12117A7F986146 (L.M.); Istituto Pasteur Italia—Fondazione Cenci Bolognetti (E.F.); PRIN PRIN2020 20203AMKTW (F.B.); Project LOCALSCENT, Grant PROT. A0375-2020-36549, Call POR-FESR “Gruppi di Ricerca 2020” (G.G., M.L.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data are included in the manuscript.

Acknowledgments

T.M.A. is fellow of Network Oncology and Precision Medicine, Dpt. Experimental Medicine, Sapienza University of Rome. Graphical abstract and Figure 6 were created with BioRender.com.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cervantes, A.; Adam, R.; Roselló, S.; Arnold, D.; Normanno, N.; Taïeb, J.; Seligmann, J.; De Baere, T.; Osterlund, P.; Yoshino, T.; et al. Metastatic colorectal cancer: ESMO Clinical Practice Guideline for diagnosis, treatment and follow-up†. Ann. Oncol. 2023, 34, 10–32. [Google Scholar] [CrossRef]
  2. Morkel, M.; Riemer, P.; Bläker, H.; Sers, C.; Morkel, M.; Riemer, P.; Bläker, H.; Sers, C. Similar but different: Distinct roles for KRAS and BRAF oncogenes in colorectal cancer development and therapy resistance. Oncotarget 2015, 6, 20785–20800. [Google Scholar] [CrossRef] [Green Version]
  3. Li, Z.N.; Zhao, L.; Yu, L.F.; Wei, M.J. BRAF and KRAS mutations in metastatic colorectal cancer: Future perspectives for personalized therapy. Gastroenterol. Rep. 2020, 8, 192–205. [Google Scholar] [CrossRef] [PubMed]
  4. Vodenkova, S.; Buchler, T.; Cervena, K.; Veskrnova, V.; Vodicka, P.; Vymetalkova, V. 5-fluorouracil and other fluoropyrimidines in colorectal cancer: Past, present and future. Pharmacol. Ther. 2020, 206, 107447. [Google Scholar] [CrossRef] [PubMed]
  5. Cho, Y.H.; Ro, E.J.; Yoon, J.S.; Mizutani, T.; Kang, D.W.; Park, J.C.; Il Kim, T.; Clevers, H.; Choi, K.Y. 5-FU promotes stemness of colorectal cancer via p53-mediated WNT/β-catenin pathway activation. Nat. Commun. 2020, 11, 5321. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, L.; Yang, F.; Chen, S.; Tai, J. Mechanisms on chemotherapy resistance of colorectal cancer stem cells and research progress of reverse transformation: A mini-review. Front. Med. 2022, 9, 2592. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, J.; Na, R.; Xiao, C.; Wang, X.; Wang, Y.; Yan, D.; Song, G.; Liu, X.; Chen, J.; Lu, H.; et al. The loss of SHMT2 mediates 5-fluorouracil chemoresistance in colorectal cancer by upregulating autophagy. Oncogene 2021, 40, 3974–3988. [Google Scholar] [CrossRef]
  8. Zhang, K.; Zhang, T.; Yang, Y.; Tu, W.; Huang, H.; Wang, Y.; Chen, Y.; Pan, K.; Chen, Z. N6-methyladenosine-mediated LDHA induction potentiates chemoresistance of colorectal cancer cells through metabolic reprogramming. Theranostics 2022, 12, 4802–4817. [Google Scholar] [CrossRef]
  9. Marin, J.J.G.; Macias, R.I.R.; Monte, M.J.; Herraez, E.; Peleteiro-vigil, A.; de Blas, B.S.; Sanchon-sanchez, P.; Temprano, A.G.; Espinosa-escudero, R.A.; Lozano, E.; et al. Cellular Mechanisms Accounting for the Refractoriness of Colorectal Carcinoma to Pharmacological Treatment. Cancers 2020, 12, 2605. [Google Scholar] [CrossRef]
  10. Pietrobono, S.; Gagliardi, S.; Stecca, B. Non-canonical hedgehog signaling pathway in cancer: Activation of GLI transcription factors beyond smoothened. Front. Genet. 2019, 10, 556. [Google Scholar] [CrossRef] [Green Version]
  11. Po, A.; Silvano, M.; Miele, E.; Capalbo, C.; Eramo, A.; Salvati, V.; Todaro, M.; Besharat, Z.M.; Catanzaro, G.; Cucchi, D.; et al. Noncanonical GLI1 signaling promotes stemness features and in vivo growth in lung adenocarcinoma. Oncogene 2017, 36, 4641–4652. [Google Scholar] [CrossRef] [Green Version]
  12. Artavanis-Tsakonas, S.; Rand, M.D.; Lake, R.J. Notch Signaling: Cell Fate Control and Signal Integration in Development. Science 1999, 284, 770–776. [Google Scholar] [CrossRef] [Green Version]
  13. Sun, L.; Ke, J.; He, Z.; Chen, Z.; Huang, Q.; Ai, W.; Wang, G.; Wei, Y.; Zou, X.; Zhang, S.; et al. HES1 Promotes Colorectal Cancer Cell Resistance To 5-Fu by Inducing Of EMT and ABC Transporter Proteins. J. Cancer 2017, 8, 2802–2808. [Google Scholar] [CrossRef] [PubMed]
  14. Zhou, B.; Lin, W.; Long, Y.; Yang, Y.; Zhang, H.; Wu, K.; Chu, Q. Notch signaling pathway: Architecture, disease, and therapeutics. Signal Transduct. Target. Ther. 2022, 7, 95. [Google Scholar] [CrossRef] [PubMed]
  15. Radtke, F.; Clevers, H.; Riccio, O. From Gut Homeostasis to Cancer. Curr. Mol. Med. 2006, 6, 275–289. [Google Scholar] [CrossRef]
  16. Jacobs, C.T.; Huang, P. Notch signalling maintains hedgehog responsiveness via a gli-dependent mechanism during spinal cord patterning in zebrafish. elife 2019, 8, e49252. [Google Scholar] [CrossRef] [PubMed]
  17. Kumar, V.; Vashishta, M.; Kong, L.; Wu, X.; Lu, J.J.; Guha, C.; Dwarakanath, B.S. The Role of Notch, Hedgehog, and Wnt Signaling Pathways in the Resistance of Tumors to Anticancer Therapies. Front. Cell Dev. Biol. 2021, 9, 857. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, J.; Fan, J.; Zeng, X.; Nie, M.; Luan, J.; Wang, Y.; Ju, D.; Yin, K. Hedgehog signaling in gastrointestinal carcinogenesis and the gastrointestinal tumor microenvironment. Acta Pharm. Sin. B 2021, 11, 609–620. [Google Scholar] [CrossRef]
  19. Tyagi, A.; Sharma, A.K.; Damodaran, C. A Review on Notch Signaling and Colorectal Cancer. Cells 2020, 9, 1549. [Google Scholar] [CrossRef] [PubMed]
  20. Po, A.; Citarella, A.; Catanzaro, G.; Besharat, Z.M.; Trocchianesi, S.; Gianno, F.; Sabato, C.; Moretti, M.; De Smaele, E.; Vacca, A.; et al. Hedgehog-GLI signalling promotes chemoresistance through the regulation of ABC transporters in colorectal cancer cells. Sci. Rep. 2020, 10, 13988. [Google Scholar] [CrossRef]
  21. Celano, M.; Schenone, S.; Cosco, D.; Navarra, M.; Puxeddu, E.; Racanicchi, L.; Brullo, C.; Varano, E.; Alcaro, S.; Ferretti, E.; et al. Cytotoxic effects of a novel pyrazolopyrimidine derivative entrapped in liposomes in anaplastic thyroid cancer cells in vitro and in xenograft tumors in vivo. Endocr. Relat. Cancer 2008, 15, 499–510. [Google Scholar] [CrossRef] [PubMed]
  22. Meijering, E.; Dzyubachyk, O.; Smal, I.; van Cappellen, W.A. Tracking in cell and developmental biology. Semin. Cell Dev. Biol. 2009, 20, 894–902. [Google Scholar] [CrossRef] [PubMed]
  23. Jaqaman, K.; Loerke, D.; Mettlen, M.; Kuwata, H.; Grinstein, S.; Schmid, S.L.; Danuser, G. Robust single-particle tracking in live-cell time-lapse sequences. Nat. Methods 2008, 5, 695–702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Michalet, X. Mean square displacement analysis of single-particle trajectories with localization error: Brownian motion in an isotropic medium. Phys. Rev. E-Stat. Nonlinear Soft Matter Phys. 2010, 82, 041914. [Google Scholar] [CrossRef] [Green Version]
  25. Soignet, S.L.; Frankel, S.R.; Douer, D.; Tallman, M.S.; Kantarjian, H.; Calleja, E.; Stone, R.M.; Kalaycio, M.; Scheinberg, D.A.; Steinherz, P.; et al. United States multicenter study of arsenic trioxide in relapsed acute promyelocytic leukemia. J. Clin. Oncol. 2001, 19, 3852–3860. [Google Scholar] [CrossRef]
  26. Ding, D.; Lim, K.S.; Eberhart, C.G. Arsenic trioxide inhibits Hedgehog, Notch and stem cell properties in glioblastoma neurospheres. Acta Neuropathol. Commun. 2014, 2, 31. [Google Scholar] [CrossRef] [Green Version]
  27. Yang, L.; Shi, P.; Zhao, G.; Xu, J.; Peng, W.; Zhang, J.; Zhang, G.; Wang, X.; Dong, Z.; Chen, F.; et al. Targeting cancer stem cell pathways for cancer therapy. Signal Transduct. Target. Ther. 2020, 5, 8. [Google Scholar] [CrossRef] [Green Version]
  28. Clevers, H. Modeling Development and Disease with Organoids. Cell 2016, 165, 1586–1597. [Google Scholar] [CrossRef] [Green Version]
  29. Lancaster, M.A.; Knoblich, J.A. Organogenesis in a dish: Modeling development and disease using organoid technologies. Science 2014, 345, 1247125. [Google Scholar] [CrossRef]
  30. Xi, Y.; Xu, P. Global colorectal cancer burden in 2020 and projections to 2040. Transl. Oncol. 2021, 14, 101174. [Google Scholar] [CrossRef]
  31. Kukcinaviciute, E.; Jonusiene, V.; Sasnauskiene, A.; Dabkeviciene, D.; Eidenaite, E.; Laurinavicius, A. Significance of Notch and Wnt signaling for chemoresistance of colorectal cancer cells HCT116. J. Cell. Biochem. 2018, 119, 5913–5920. [Google Scholar] [CrossRef] [PubMed]
  32. Yuan, S.; Tao, F.; Zhang, X.; Zhang, Y.; Sun, X.; Wu, D. Role of Wnt/ β-Catenin Signaling in the Chemoresistance Modulation of Colorectal Cancer. Biomed Res. Int. 2020, 2020, 9390878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. West, K.A.; Castillo, S.S.; Dennis, P.A. Activation of the PI3K/Akt pathway and chemotherapeutic resistance. Drug Resist. Updat. 2002, 5, 234–248. [Google Scholar] [CrossRef]
  34. Deng, J.; Bai, X.; Feng, X.; Ni, J.; Beretov, J.; Graham, P.; Li, Y. Inhibition of PI3K/Akt/mTOR signaling pathway alleviates ovarian cancer chemoresistance through reversing epithelial-mesenchymal transition and decreasing cancer stem cell marker expression. BMC Cancer 2019, 19, 618. [Google Scholar] [CrossRef] [PubMed]
  35. Zhao, Y.; Shen, S.; Guo, J.; Chen, H.; Yu Greenblatt, D.; Kleeff, J.; Liao, Q.; Chen, G.; Friess, H.; Sing Leung, P. Mitogen-Activated Protein Kinases and Chemoresistance in Pancreatic Cancer Cells. J. Surg. Res. 2006, 136, 325–335. [Google Scholar] [CrossRef] [Green Version]
  36. Salaroglio, I.C.; Mungo, E.; Gazzano, E.; Kopecka, J.; Riganti, C. ERK is a Pivotal Player of Chemo-Immune-Resistance in Cancer. Int. J. Mol. Sci. 2019, 20, 2505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Godwin, P.; Baird, A.M.; Heavey, S.; Barr, M.P.; O’Byrne, K.J.; Gately, K. Targeting nuclear factor-kappa B to overcome resistance to chemotherapy. Front. Oncol. 2013, 3, 120. [Google Scholar] [CrossRef] [Green Version]
  38. Li, Q.; Yang, G.; Feng, M.; Zheng, S.; Cao, Z.; Qiu, J.; You, L.; Zheng, L.; Hu, Y.; Zhang, T.; et al. NF-κB in pancreatic cancer: Its key role in chemoresistance. Cancer Lett. 2018, 421, 127–134. [Google Scholar] [CrossRef] [PubMed]
  39. Echelard, Y.; Epstein, D.J.; St-Jacques, B.; Shen, L.; Mohler, J.; McMahon, J.A.; McMahon, A.P. Sonic hedgehog, a member of a family of putative signaling molecules, is implicated in the regulation of CNS polarity. Cell 1993, 75, 1417–1430. [Google Scholar] [CrossRef] [PubMed]
  40. Bitgood, M.J.; McMahon, A.P. HedgehogandBmpGenes Are Coexpressed at Many Diverse Sites of Cell–Cell Interaction in the Mouse Embryo. Dev. Biol. 1995, 172, 126–138. [Google Scholar] [CrossRef] [Green Version]
  41. de Reyniès, A.; Javelaud, D.; Elarouci, N.; Marsaud, V.; Gilbert, C.; Mauviel, A. Large-scale pan-cancer analysis reveals broad prognostic association between TGF-β ligands, not Hedgehog, and GLI1/2 expression in tumors. Sci. Rep. 2020, 10, 14491. [Google Scholar] [CrossRef]
  42. Eberl, M.; Klingler, S.; Mangelberger, D.; Loipetzberger, A.; Damhofer, H.; Zoidl, K.; Schnidar, H.; Hache, H.; Bauer, H.C.; Solca, F.; et al. Hedgehog-EGFR cooperation response genes determine the oncogenic phenotype of basal cell carcinoma and tumour-initiating pancreatic cancer cells. EMBO Mol. Med. 2012, 4, 218–233. [Google Scholar] [CrossRef] [PubMed]
  43. Miele, E.; Po, A.; Begalli, F.; Antonucci, L.; Mastronuzzi, A.; Marras, C.E.; Carai, A.; Cucchi, D.; Abballe, L.; Besharat, Z.M.; et al. β-arrestin1-mediated acetylation of Gli1 regulates Hedgehog/Gli signaling and modulates self-renewal of SHH medulloblastoma cancer stem cells. BMC Cancer 2017, 17, 488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Varnat, F.; Siegl-Cachedenier, I.; Malerba, M.; Gervaz, P.; Ruiz, I. Altaba, A. Loss of WNT-TCF addiction and enhancement of HH-GLI1 signalling define the metastatic transition of human colon carcinomas. EMBO Mol. Med. 2010, 2, 440–457. [Google Scholar] [CrossRef] [PubMed]
  45. Noubissi, F.K.; Goswami, S.; Sanek, N.A.; Kawakami, K.; Minamoto, T.; Moser, A.; Grinblat, Y.; Spiegelman, V.S. Wnt signaling stimulates transcriptional outcome of the hedgehog pathway by stabilizing GLI1 mRNA. Cancer Res. 2009, 69, 8572–8578. [Google Scholar] [CrossRef] [Green Version]
  46. Medema, J.P.; Vermeulen, L. Microenvironmental regulation of stem cells in intestinal homeostasis and cancer. Nature 2011, 474, 318–326. [Google Scholar] [CrossRef]
  47. Si, Y.; Li, L.; Zhang, W.; Liu, Q.; Liu, B. GANT61 exerts anticancer cell and anticancer stem cell capacity in colorectal cancer by blocking the Wnt/β-catenin and Notch signalling pathways. Oncol. Rep. 2022, 48, 182. [Google Scholar] [CrossRef]
  48. Kerl, K.; Moreno, N.; Holsten, T.; Ahlfeld, J.; Mertins, J.; Hotfilder, M.; Kool, M.; Bartelheim, K.; Schleicher, S.; Handgretinger, R.; et al. Arsenic trioxide inhibits tumor cell growth in malignant rhabdoid tumors in vitro and in vivo by targeting overexpressed Gli1. Int. J. Cancer 2014, 135, 989–995. [Google Scholar] [CrossRef] [PubMed]
  49. Ardalan, B.; Subbarayan, P.R.; Ramos, Y.; Gonzalez, M.; Fernandez, A.; Mezentsev, D.; Reis, I.; Duncan, R.; Podolsky, L.; Lee, K.; et al. A phase I study of 5-fluorouracil/leucovorin and arsenic trioxide for patients with refractory/relapsed colorectal carcinoma. Clin. Cancer Res. 2010, 16, 3019–3027. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. HH-GLI and NOTCH signaling pathways sustain resistance to 5-FU in KRAS mutant CRC cells. (A) Western blot analysis of NOTCH1 Intracellular domain (NOTCH1 ID), GLI1 and cleaved-PARP (c-PARP) in HCT116 cells after 5-FU treatment in combination with GANT61 and DAPT. Numbers indicate intensity ratio of bands. Bar graphs show densitometric quantification of the band intensity values normalized to the loading control. * p  <  0.05; ** p  <  0.01; *** p < 0.001 versus control; § p < 0.05 versus 5-FU; §§ p  <  0.01; §§§ p < 0.001 (Two-way ANOVA test). Uncropped full scan in Supplementary Figure S4. (B) Evaluation of cell viability by Trypan Blue exclusion assay in HCT116 after 5-FU treatment with or without GANT61 and/or DAPT; * p  <  0.05, ** p  <  0.01 versus CTRL (One-way ANOVA test). (C) Correlation analysis between GLI1 and NOTCH1 expression from dataset interrogated on R2 platform, as indicated in main text. (D) Western blot analysis of GLI1, NOTCH1 ID and c-PARP in HCT116 cells treated with 5-FU and ATO. Numbers indicate intensity ratio of bands. Bar graphs show densitometric quantification of the band intensity values normalized to the loading control; * p  <  0.05; ** p  <  0.01; (Two-way ANOVA test). (E) Transwell invasion assay in HCT116 cells treated with 5-FU, ATO and the combined treatment and control group (CTRL). Scale bar 150 µm. ** p  < 0.01; *** p < 0.001; versus control; §§§ p < 0.001 (One-way ANOVA test).
Figure 1. HH-GLI and NOTCH signaling pathways sustain resistance to 5-FU in KRAS mutant CRC cells. (A) Western blot analysis of NOTCH1 Intracellular domain (NOTCH1 ID), GLI1 and cleaved-PARP (c-PARP) in HCT116 cells after 5-FU treatment in combination with GANT61 and DAPT. Numbers indicate intensity ratio of bands. Bar graphs show densitometric quantification of the band intensity values normalized to the loading control. * p  <  0.05; ** p  <  0.01; *** p < 0.001 versus control; § p < 0.05 versus 5-FU; §§ p  <  0.01; §§§ p < 0.001 (Two-way ANOVA test). Uncropped full scan in Supplementary Figure S4. (B) Evaluation of cell viability by Trypan Blue exclusion assay in HCT116 after 5-FU treatment with or without GANT61 and/or DAPT; * p  <  0.05, ** p  <  0.01 versus CTRL (One-way ANOVA test). (C) Correlation analysis between GLI1 and NOTCH1 expression from dataset interrogated on R2 platform, as indicated in main text. (D) Western blot analysis of GLI1, NOTCH1 ID and c-PARP in HCT116 cells treated with 5-FU and ATO. Numbers indicate intensity ratio of bands. Bar graphs show densitometric quantification of the band intensity values normalized to the loading control; * p  <  0.05; ** p  <  0.01; (Two-way ANOVA test). (E) Transwell invasion assay in HCT116 cells treated with 5-FU, ATO and the combined treatment and control group (CTRL). Scale bar 150 µm. ** p  < 0.01; *** p < 0.001; versus control; §§§ p < 0.001 (One-way ANOVA test).
Cancers 15 01471 g001
Figure 2. HH-GLI signaling pathway sustains resistance to 5-FU in BRAF mutant CRC cells. (A) Western blot analysis of NOTCH1 Intracellular domain (NOTCH1 ID), GLI1 and cleaved-PARP (c-PARP) in HT29 cells after 5-FU treatment in combination with GANT61 and DAPT. Numbers indicate intensity ratio of bands. Bar graphs show densitometrically quantified band intensity values normalized to the loading control; * p  <  0.05 ** p  <  0.01; *** p < 0.001; § p < 0.05 versus 5-FU; §§ p < 0.01 (Two-way ANOVA test). Uncropped full scan in Supplementary Figure S5. (B) Evaluation of cell viability by Tripan Blue exclusion assay in HT29 after 5-FU treatment in combination with GANT61 and DAPT; * p  <  0.05, ** p  <  0.01 versus CTRL (One-way ANOVA test). (C) Correlation analysis between GLI1 and NOTCH1 expression from dataset interrogated on R2 platform as indicated in main text. (D) Transwell invasion assay in HT29 cells treated with 5-FU, GANT61, the combined treatment and control group (CTRL); Scale bar 150 µm; * p  <  0.05, ** p  <  0.01 versus CTRL (One-way ANOVA test).
Figure 2. HH-GLI signaling pathway sustains resistance to 5-FU in BRAF mutant CRC cells. (A) Western blot analysis of NOTCH1 Intracellular domain (NOTCH1 ID), GLI1 and cleaved-PARP (c-PARP) in HT29 cells after 5-FU treatment in combination with GANT61 and DAPT. Numbers indicate intensity ratio of bands. Bar graphs show densitometrically quantified band intensity values normalized to the loading control; * p  <  0.05 ** p  <  0.01; *** p < 0.001; § p < 0.05 versus 5-FU; §§ p < 0.01 (Two-way ANOVA test). Uncropped full scan in Supplementary Figure S5. (B) Evaluation of cell viability by Tripan Blue exclusion assay in HT29 after 5-FU treatment in combination with GANT61 and DAPT; * p  <  0.05, ** p  <  0.01 versus CTRL (One-way ANOVA test). (C) Correlation analysis between GLI1 and NOTCH1 expression from dataset interrogated on R2 platform as indicated in main text. (D) Transwell invasion assay in HT29 cells treated with 5-FU, GANT61, the combined treatment and control group (CTRL); Scale bar 150 µm; * p  <  0.05, ** p  <  0.01 versus CTRL (One-way ANOVA test).
Cancers 15 01471 g002
Figure 3. 5-FU increases motility in CRC organoids. (A) Workflow for CRC organoid growth. (B) Brightfield image of HCT116 and HT29 organoids (respectively oHCT116 and oHT29) after 7 days of culture. (C) Quantitative real-time PCR of HES1, GLI1, c-MET expressed in HCT116 and HT29 cultured in monolayer and as organoids; * p  <  0.05 versus CTRL; ** p  <  0.01 (Two-way ANOVA test). (D) Fluorescent images of GFP-labeled HCT116 on sequential hours, scale bar 200 µm. Supplementary Videos S1 and S2 of time lapse experiments are available in Supplementary Figure S3A,B. (E) Violin plot of the diffusion parameters obtained from the single cell trajectories for CTRL and 5-FU-treated organoids. Kolmogorov–Smirnov test p-value 0.0032.
Figure 3. 5-FU increases motility in CRC organoids. (A) Workflow for CRC organoid growth. (B) Brightfield image of HCT116 and HT29 organoids (respectively oHCT116 and oHT29) after 7 days of culture. (C) Quantitative real-time PCR of HES1, GLI1, c-MET expressed in HCT116 and HT29 cultured in monolayer and as organoids; * p  <  0.05 versus CTRL; ** p  <  0.01 (Two-way ANOVA test). (D) Fluorescent images of GFP-labeled HCT116 on sequential hours, scale bar 200 µm. Supplementary Videos S1 and S2 of time lapse experiments are available in Supplementary Figure S3A,B. (E) Violin plot of the diffusion parameters obtained from the single cell trajectories for CTRL and 5-FU-treated organoids. Kolmogorov–Smirnov test p-value 0.0032.
Cancers 15 01471 g003
Figure 4. HH-GLI and NOTCH inhibition impairs 5-FU-driven mesenchymal phenotype in KRASG13D-driven CRC organoids. (A) Brightfield image of oHCT116 treated with 5-FU, ATO, their combination and the control group (CTRL); Scale bar 150 µm; * p  <  0.05 versus CTRL; §§ p < 0.01 versus 5-FU (One-way ANOVA test). (B) Quantitative real-time PCR of stem markers expressed in oHCT116 treated with 5-FU, ATO, their combination and the control group (CTRL). mRNA levels of ABCG2, CD133a and c-MET expressed in oHCT116 were expressed in arbitrary units; * p  <  0.05 versus CTRL; §§ p < 0.01 versus 5-FU (Two-way ANOVA test). (C) Whole-mount immunofluorescence staining of oHCT116 stained with phalloidin-594 (F-actin, red), vimentin (mesenchymal marker, green) and DNA (DAPI). Images were analyzed by using the program Zeiss ZEN 2.3 blue edition. Scale bar 200 µm.
Figure 4. HH-GLI and NOTCH inhibition impairs 5-FU-driven mesenchymal phenotype in KRASG13D-driven CRC organoids. (A) Brightfield image of oHCT116 treated with 5-FU, ATO, their combination and the control group (CTRL); Scale bar 150 µm; * p  <  0.05 versus CTRL; §§ p < 0.01 versus 5-FU (One-way ANOVA test). (B) Quantitative real-time PCR of stem markers expressed in oHCT116 treated with 5-FU, ATO, their combination and the control group (CTRL). mRNA levels of ABCG2, CD133a and c-MET expressed in oHCT116 were expressed in arbitrary units; * p  <  0.05 versus CTRL; §§ p < 0.01 versus 5-FU (Two-way ANOVA test). (C) Whole-mount immunofluorescence staining of oHCT116 stained with phalloidin-594 (F-actin, red), vimentin (mesenchymal marker, green) and DNA (DAPI). Images were analyzed by using the program Zeiss ZEN 2.3 blue edition. Scale bar 200 µm.
Cancers 15 01471 g004
Figure 5. HH-GLI inhibition impairs 5-FU-driven mesenchymal phenotype in BRAFV600E-driven CRC organoids. (A) Brightfield image of oHT29 treated with 5-FU, GANT61, their combination and the control group (CTRL); Scale bar 150µm; * p  <  0.05 versus Ctrl; §§ p < 0.01 versus 5-FU (One-way ANOVA test). (B) Quantitative real-time PCR of ABCG2, CD133a and c-MET expressed in oHT29 were expressed in arbitrary units. Data are representative of three independent experiments, * p  <  0.05 versus Ctrl; §§ p < 0.01 versus 5-FU (Two-way ANOVA test). (C) Whole-mount immunofluorescence of HT29 organoids stained using phalloidin-594 (f-actin, red), vimentin (mesenchymal marker, green) and DNA (DAPI). Images were analyzed by using the program Zeiss ZEN 2.3 blue edition. Scale bar 200 µm.
Figure 5. HH-GLI inhibition impairs 5-FU-driven mesenchymal phenotype in BRAFV600E-driven CRC organoids. (A) Brightfield image of oHT29 treated with 5-FU, GANT61, their combination and the control group (CTRL); Scale bar 150µm; * p  <  0.05 versus Ctrl; §§ p < 0.01 versus 5-FU (One-way ANOVA test). (B) Quantitative real-time PCR of ABCG2, CD133a and c-MET expressed in oHT29 were expressed in arbitrary units. Data are representative of three independent experiments, * p  <  0.05 versus Ctrl; §§ p < 0.01 versus 5-FU (Two-way ANOVA test). (C) Whole-mount immunofluorescence of HT29 organoids stained using phalloidin-594 (f-actin, red), vimentin (mesenchymal marker, green) and DNA (DAPI). Images were analyzed by using the program Zeiss ZEN 2.3 blue edition. Scale bar 200 µm.
Cancers 15 01471 g005
Figure 6. Hedgehog-Gli and NOTCH pathways sustain chemoresistance and the mesenchymal phenotype in CRC. Model of the activity of the Hedgehog-GLI and NOTCH pathways after chemotherapy stress in BRAFV600E and KRASG13D models. In KRAS-driven CRC, the chemotherapy stress activates both HH-GLI and NOTCH, which independently sustain the EMT program; in BRAF-driven CRC, chemotherapy stress induces the activation of the HH-GLI pathway, which in turn sustains the activation of NOTCH1 signaling, determining the acquisition of the EMT phenotype.
Figure 6. Hedgehog-Gli and NOTCH pathways sustain chemoresistance and the mesenchymal phenotype in CRC. Model of the activity of the Hedgehog-GLI and NOTCH pathways after chemotherapy stress in BRAFV600E and KRASG13D models. In KRAS-driven CRC, the chemotherapy stress activates both HH-GLI and NOTCH, which independently sustain the EMT program; in BRAF-driven CRC, chemotherapy stress induces the activation of the HH-GLI pathway, which in turn sustains the activation of NOTCH1 signaling, determining the acquisition of the EMT phenotype.
Cancers 15 01471 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Citarella, A.; Catanzaro, G.; Besharat, Z.M.; Trocchianesi, S.; Barbagallo, F.; Gosti, G.; Leonetti, M.; Di Fiore, A.; Coppola, L.; Autilio, T.M.; et al. Hedgehog-GLI and Notch Pathways Sustain Chemoresistance and Invasiveness in Colorectal Cancer and Their Inhibition Restores Chemotherapy Efficacy. Cancers 2023, 15, 1471. https://doi.org/10.3390/cancers15051471

AMA Style

Citarella A, Catanzaro G, Besharat ZM, Trocchianesi S, Barbagallo F, Gosti G, Leonetti M, Di Fiore A, Coppola L, Autilio TM, et al. Hedgehog-GLI and Notch Pathways Sustain Chemoresistance and Invasiveness in Colorectal Cancer and Their Inhibition Restores Chemotherapy Efficacy. Cancers. 2023; 15(5):1471. https://doi.org/10.3390/cancers15051471

Chicago/Turabian Style

Citarella, Anna, Giuseppina Catanzaro, Zein Mersini Besharat, Sofia Trocchianesi, Federica Barbagallo, Giorgio Gosti, Marco Leonetti, Annamaria Di Fiore, Lucia Coppola, Tanja Milena Autilio, and et al. 2023. "Hedgehog-GLI and Notch Pathways Sustain Chemoresistance and Invasiveness in Colorectal Cancer and Their Inhibition Restores Chemotherapy Efficacy" Cancers 15, no. 5: 1471. https://doi.org/10.3390/cancers15051471

APA Style

Citarella, A., Catanzaro, G., Besharat, Z. M., Trocchianesi, S., Barbagallo, F., Gosti, G., Leonetti, M., Di Fiore, A., Coppola, L., Autilio, T. M., Spinello, Z., Vacca, A., De Smaele, E., Venneri, M. A., Ferretti, E., Masuelli, L., & Po, A. (2023). Hedgehog-GLI and Notch Pathways Sustain Chemoresistance and Invasiveness in Colorectal Cancer and Their Inhibition Restores Chemotherapy Efficacy. Cancers, 15(5), 1471. https://doi.org/10.3390/cancers15051471

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop