Next Article in Journal
The Release of Inflammatory Mediators from Acid-Stimulated Mesenchymal Stromal Cells Favours Tumour Invasiveness and Metastasis in Osteosarcoma
Next Article in Special Issue
CDK5RAP3, a New BRCA2 Partner That Regulates DNA Repair, Is Associated with Breast Cancer Survival
Previous Article in Journal
IGF1R/IR Mediates Resistance to BRAF and MEK Inhibitors in BRAF-Mutant Melanoma
Previous Article in Special Issue
Individualized Proteogenomics Reveals the Mutational Landscape of Melanoma Patients in Response to Immunotherapy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 22
10.3390/cancers13225858
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Hedgehog Pathway Inhibition on Nevoid Basal Cell Carcinoma Syndrome Fibroblasts and Basal Cell Carcinoma-Associated Fibroblasts: Are Vismodegib and Sonidegib Useful to Target Cancer-Prone Fibroblasts?

by
Laura Eibenschutz
1,
Silvia Caputo
2,
Emanuela Camera
2,
Anna Carbone
1,
Vitaliano Silipo
1,
Emilia Migliano
3,
Caterina Aurizi
4,
Carlo Cota
5,
Pasquale Frascione
1 and
Barbara Bellei
2,*
1
Oncologic and Preventative Dermatology, San Gallicano Dermatological Institute, IRCCS, 00144 Rome, Italy
2
Laboratory of Cutaneous Physiopathology and Integrated Center of Metabolomics Research, San Gallicano Dermatological Institute, IRCCS, 00144 Rome, Italy
3
Department of Plastic and Reconstructive Surgery, San Gallicano Dermatological Institute, IRCCS, 00144 Rome, Italy
4
Department of Dermatological Clinic, San Gallicano Dermatological Institute, IRCCS, 00144 Rome, Italy
5
Department of Dermophatology, San Gallicano Dermatological Institute, IRCCS, 00144 Rome, Italy
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(22), 5858; https://doi.org/10.3390/cancers13225858
Submission received: 13 October 2021 / Accepted: 21 October 2021 / Published: 22 November 2021
(This article belongs to the Special Issue Functional Genomics of Cancer)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Nevoid basal cell carcinoma syndrome (NBCCS) is a genetic disorder of autosomal dominant inheritance that dramatically predisposes a patient to the formation of basal cell carcinoma (BCC) due to causative mutations in several genes associated with the Hedgehog (Hh) pathway. Somatic mutations in this pathway are also associated with sporadic BCC, the most common form of skin cancer. Hh signalling extends its effect on tumorigenesis by modulating the tumour microenvironment in a paracrine fashion. Consistently, NBCCS fibroblasts could facilitate BCC occurrence. Here, we investigated vismodegib and sonidegib, two molecules currently used to target this pathway in cancer cells, as a therapeutic option against syndromic and BCC-associated fibroblasts.

Abstract

Activating mutations in the Hh pathway underlies the development of sporadic and familial skin BCC. For these oncogenic proliferations displaying ligand-independent activation of the intracellular pathway, two molecules have been approved for therapeutic purposes: vismodegib and sonidegib. Improper Hh signalling occurs in many human tumours also via a paracrine mechanism (ligand-dependent) in which the secretion of Hh ligands by stromal cells support tumour growth. On the other hand, the mobilization of neoplastic stroma by cancer cells is sustained by the activation of Hh signalling in surrounding fibroblasts suggesting a central role of this bidirectional crosstalk in carcinogenesis. Additionally, loss-of-function mutations in the PTCH1 gene in the context of NBCCS, an autosomal dominant disorder predisposing to multiple BCCs, determine tumour permissive phenotypes in dermal fibroblasts. Here, profiling syndromic and BCC-associated fibroblasts unveiled an extraordinary similarity characterized by overexpression of several Hh target genes and a marked pro-inflammatory outline. Both cell types exposed to Hh inhibitors displayed reversion of the tumour-prone phenotype. Under vismodegib and sonidegib treatment, the Wnt/β-catenin pathway, frequently over-active in tumour stroma, resulted down-regulated by pAKT-GSK3β axis and consequent increase of β-catenin turnover. Overall, this study demonstrated that vismodegib and sonidegib impacting on fibroblast tumour supportive functions might be considered in therapy for BCC independently to the mutation status of Hh components in neoplastic cells.

1. Introduction

Basal cell carcinoma (BCC) accounts for almost 80% of skin cancers [1,2] and its healthcare workload is principally within dermatology departments. Although most BCCs are small, well-defined tumours amenable to surgery or conservative procedures, in a small proportion of patients, BCCs can progress to an advanced stage including locally advanced BCC and metastatic BCC. Advanced BCC can cause disfigurement, morbidity, and lower a patient’s quality of life [3,4]. Thus, in these selected cases pharmacological therapies are preferable to surgical treatment. Non-surgical techniques include photodynamic therapy and topical imiquimod or 5-fluorouracil treatment whereas cytotoxic systemic chemotherapy has not been approved for the treatment of non-resectable BCCs. However, most of these approaches aim to obtain local control of slow-growing lesions, whereas advanced BCC can incompletely benefit from this non-resolutive care. Recently, vismodegib (GDC0449; trade name Erivedge) and sonidegib (LDE225; trade name Odomzo) two hedgehog pathway (Hh) inhibitors have been approved for oral treatment of patients with metastatic BCC or patients with recurrent, locally advanced BCC who are inappropriate for surgery and radiotherapy. The Hh signalling is a major regulator of cell differentiation, cell proliferation, and tissue polarity. The Hh pathway is largely inactive in adult tissue except for its function in tissue repair and regulation of adult stem cells but it is aberrantly active in several types of tumours [1,2,5]. In the canonical pathway, binding of Sonic Hedgehog (SHH), the most widely expressed hedgehog protein to its receptor, patched homolog 1 (PTCH1), releases inhibition on Smoothened (SMO), a 7-span transmembrane protein, decreasing the interaction between suppressor-of-fused homolog (SUFU) and the glioma-associated transcription factors (GLI1, GLI2, and GLI3), the terminal effectors of the pathway [6,7,8]. In mammals, Desert Hedgehog (DHH) and Indian Hedgehog (IHH) also function as ligands for PTCH1 resulting in intracellular signalling activation. The Hh target genes include proliferation and differentiation regulating genes (CyclinD1 and D2, E2F1, PDGFRα, IGFBP3, and IGFBP6), cell survival regulating gene (BCL-2), pro-angiogenic factors (VEGF, Cyr61) as well as feedback genes of Hh pathway (PTCH1 and 2, GLI1 and huntingtin interacting protein-1, Hip1) that further amplify signalling activation [8,9]. The Hh pathway is in conjunction with other important pathways, including epidermal growth factors (EGF/EGFR), transforming growth factor-β (TGF-β), Wnt/β-catenin signalling, PI3K/mTOR, NF-kappaB and Vitamin D [10,11,12]. Improper activation of Hh signalling is the hallmark of BCC pathogenesis [13,14]. The causative role of aberrant Hh signalling in BCC is underlined by the fact that approximatively 70–80%, 20%, and 8% have PTCH1, SMO, and SUFU driver mutations respectively [15,16,17]. Bi-allelic inactivating mutations are necessary for the loss of Hh pathway control that confers cellular growth deregulation, and, potentially, tumour development. Thus, pharmacological inhibition of the Hh pathway exerts an anti-tumor effect in several human cell lines, including renal carcinoma [18], oral squamous cell carcinoma [19], breast cancer cells [20] and skin basal cell carcinoma [21]. Another class of patients that benefices from Hh inhibitors are individuals with nevoid basal cell carcinoma syndrome (NBCCS) (also referred as nevoid BCC syndrome, Gorlin syndrome, or Gorlin Goltz syndrome, OMIM#109400), an autosomal dominant rare hereditary condition frequently due to germline mutation of PTCH1 gene lacking significant genotype-phenotype correlations [22,23,24]. Prevalence of NBCCS is estimated from 1 per 57,000 to 1 per 164,000 with an even sex distribution and approximately 20–30% de novo mutation [25,26,27]. Affected individuals have two major phenotypes: developmental defects associated with generalized overgrowth and an increased risk of developing cancers including multiple BCCs (with a median age of onset of 20 years) frequently occurring in unexposed body areas [26,28]. In NBCCS patients presenting systemic loss of one functional copy (haploinsufficiency) of PTCH1, inactivation of the normal homolog by environmental mutagenesis or random genetic rearrangement confers ligand-independent activation of the Hh pathway and tumour permissive phenotype [29]. Unlike sporadic BCC patients, target therapy for these individuals necessitates a dual approach, focusing on the long-term preservation of healthy skin and reversing/limiting the growth of invasive tumours. A clinical trial including 41 NBCCS patients demonstrated that vismodegib significantly reduced the number of new BCCs per year and lowered the rate of surgically eligible tumour lesions, compared to the placebo arm [30]. However, the use of vismodegib in syndromic patients is limited by medication side effects and the consequent discontinuation rate. Thus, based on the idea that long-term therapy of syndromic patient needs a safer profile treatment, an alternative of topically applied sonidegib has been tested in a proof-of-concept clinical study confirming the therapeutic utility in BCC [31]. The response rate of syndromic patients to orally administered vismodegib is higher compared to sporadic probably due to the reduced mutational load of NBCCS-BCCs and the lack of intrinsic pre-existing resistance to Hh inhibitors [32,33,34,35]. Both vismodegib and sonidegib block hedgehog signalling by selective inhibition of SMO, even if chemical structures are different and they bind different residues of the same functional pocket eventually selecting distinctive mutation-based resistance [36,37,38]. Although the two principal patient-based studies, ERIVANCE (vismodegib) and BOLT (sonidegib), validated the therapeutic use demonstrating similar efficacy and patterns of adverse events, some studies evidenced relevant pharmacological differences [39,40]. One of the most interesting discrepancies resides in molecular distribution. Sonidegib may achieve a significantly higher concentration in the skin and other tissue due to the major grade of lipophilicity [41] and effectively cross the blood-brain barrier [13]. Considerations regarding the bioavailability and tissue concentration are not restricted to cancer cells but involve exposure of tumour surrounding tissue especially in the case of topically administered therapy. The effect of SMO antagonists in tumour stroma is largely unknown. Frequently, tumour-associated stromal cells not only provide paracrine factors but also present a coordinated regulation of major oncogenic intracellular signalling [42,43,44]. A large number of tumours have been demonstrated associated with Hh peptide secretion by tumour cells themselves and by the stromal microenvironment [42,45,46]. Moreover, two different studies demonstrated that CAFs can promote epithelial-to-mesenchymal transition through the activation of the Hh pathway [47,48]. On the other hand, NBCCS fibroblasts carrying chronic active intracellular Hh signal transduction, functionally exhibit CAF features [49,50] underling the concept that Hh signalling deregulation may represent the joint point of similarity between CAFs and NBCCS fibroblasts.
This project, based on the idea that dermal-epidermal interaction plays a crucial role in BCC onset and progression, explored the possibility that inhibition of Hh signalling could be effective in targeting tumour-prone fibroblasts in the context of subjects predisposed to skin cancer and BCC-CAFs.

2. Materials and Methods

2.1. Isolation of Primary Fibroblastic Population from Human Tissue Skin

NBCCS patients were diagnosed by clinical examination and family history and then genetically characterized by genomic DNA sequencing. Tissue samples were collected from NBCCS (unaffected skin at the time of surgical operation, n = 10; age range 32–78, mean age 57.7 ± 16 years), BCCs (tumoral tissue, n = 7; age range 66–87, mean age 74.4 ± 6.9 years) and control healthy donors (n = 12, age range 29–96, mean age 61.2 ± 22.6 years). Skin fragments obtained from surgery, were catted into approximately 4 mm2 sized pieces and digested overnight at 4 °C with dispase (2.5 mg/mL) to separate epidermis from dermis. Dermis was digested with collagenase for 2 h at 37 °C and extracted normal human fibroblasts (NHFs), CAFs and NBCCS-human fibroblasts (NBCCS-HFs) were maintained in culture with DMEM (EuroClone S.p.A., Milan, Italy) supplemented with 10% FBS and antibiotics (Hyclone Laboratories, South Logan, UT, USA). NBCCS patient characteristics and corresponding PTCH1 mutation are reported in Table 1. BCC patient and healthy donor characteristics are reported in Supplementary Tables S1 and S2. Vismodegib (GDC0449) and sonidegib (LDE225) were both purchased from MedChemExpress (MedCheExpress LLC, NJ, USA) and Adooq Biosciences (Irvine, CA, USA) for comparison.

2.2. Nevoid Basal Cell Carcinoma Syndrome (NBCCS) Patients’ Genomic Characterization

PTCH1 (NM_000264.3) mutations were identified in the genome from non-tumoral tissue (fibroblasts or blood cells). Genomic DNA was extracted using Tissue Kit (Qiagen, Milan, Italy) following the manufacturer’s instructions. About 100–200 ng of genomic DNA was subject to PCR in a total volume of 50 μL containing 25 μL of 2× AmpliTaq GoldTM 360 Master Mix (Thermo Fisher Scientific, Monza, Italy) and 22 sets of primers (25 pmol each). Primer sets were designed to cover the entire coding sequences plus a few nucleotides into the intron sequences on both ends. PCR primer sequences are available upon request. DNA fragments were checked by electrophoresis in 2% agarose gel and purified before bidirectional direct Sanger sequencing and chromatograms inspection using ChromasPro software. Variants that changed their nucleotide sequence, to assess the pathogenicity data, were matched to https://varsome.com, https://clinvarminer.genetics.utah.edu and https://www.lovd.nl databases (accessed on 20 June 2021).

2.3. MMT Assays

Briefly, 0.8 × 104 fibroblasts were seeded into the 24-well plates for 24 h to adhere. Then, growth medium was changed with fresh medium containing treatments (or not for control cells) at the appropriate concentrations. At the experimental end point (72 h), cells were incubated with 3-(4,5 dimethylthiazol)-2,5-diphenyl tetrazolium bromide (MTT) for 2 h. After this time, the medium was removed and the resulting crystals were solubilized in DMSO. The absorbance was measured at 570 nm with a reference wavelength of 690 nm. Absorbance readings were subtracted from the value of blank wells, and results were calculated as a percentage of absorbance with respect to control samples. Experiments were performed in duplicate.

2.4. Flow Cytometry Analysis

Cell death and apoptosis were analysed by the annexin-V FITC/propidium iodide (PI) double staining method after 48 h of treatment. Cells were harvested by trypsinization, suspended in the staining buffer (10 mm HEPES ⁄ NaOH, pH 7.4, 140 mm NaCl, 2.5 mm CaCl2), stained with FITC-labeled annexin V and PI for 15 min at RT in the dark and then kept on ice until analysis; 20,000/sample cells were analysed using a MACSQuant 10 instrument. Data were interpreted with MACSQuantify software (Milthenyi Biotech, S.r.l., Bologna, Italy).

2.5. Semi-Quantitative Real-Time Polymerase Chain Reaction (RT-PCR) and Gene Expression Array Cards Analysis

Total RNA was extracted using Aurum Total mini kit (BioRad, Milan Italy). cDNA was synthesized from 1 μg of total RNA using the FirstAid kit (Fermentas, ThermoFisher Scientific, Waltham, MA, USA) and loaded on 384-well microfluidic cards designed to perform probe-based TaqMan real-time PCR on an Applied Biosystems® QuantStudioTM 7 Flex instrument (Thermo Fisher Scientific, Monza, Italy). Cards were configured with selected primers and probe sets to analyze 90 target genes and three housekeeping genes (18s rrna, gliceraldeide-3-fosfato deidrogenasi (GAPDH) and β-actin). Results were evaluated using cloud-based platform software (ThermoFisher Scientific). For semi-quantitative real-time PCR, cDNA were amplified using SsoAdvanced Universal Syber Green Supermix (BioRad) containing 25 pmol of forward and reverse primers using a CXF96 Touch Cycler (BioRad). All samples were tested in triplicate. Amplification of the β-actin transcript from each sample was included as internal control. Sequences of primers (intron spanning) are given separately in the Supporting Information (Table S3).

2.6. Western Blot Analysis

Prepared cells were lysed with RIPA buffer containing proteases and phosphatases inhibitors. Proteins were separated on SDS-polyacrylamide gels, transferred to nitrocellulose membranes and then treated with the following primary antibodies: mouse monoclonal anti-β-catenin (Santa Cruz Biotechnology Inc.), rabbit polyclonal pSer473-AKT and pSer9-GSK3β (Cell Signalling Technology, MA, USA). Anti-β-actin mouse monoclonal antibody (Sigma Aldrich, Merck KGaA, Darmstadt, Germany) was used to normalize protein content. Horseradish peroxide-conjugated goat anti-mouse and goat anti-rabbit antibodies complexes were detected by chemiluminescence (Cell Signalling Technology). Imaging and densitometric analysis were performed with a UVITEC Mini HD9 acquisition system (Alliance UVItec Ltd., Cambridge, UK).

2.7. Immunofluorescence Analysis

Cells on coverslips were fixed with 4% paraformaldehyde for 15 min at room temperature followed by 0.1% Triton X-100 to allow cell permeabilization. Cells were then incubated with anti-α-SMA monoclonal (Sigma Aldrich, Merck Life Science S.r.l. Milan, Italy), or anti-FAP rabbit (Cohesion Bioscience, London, UK) for 1 h. Primary antibodies were visualized using an anti-mouse or anti-rabbit IgG Alexa Fluor 488 (BD Biosciences, Milan, Italy). Fluorescence signals were recorded using a CCD camera (Zeiss, Oberkochen, Germany).

2.8. Cytokines Protein Array

The expression of 20 human cytokines and 41 growth factors were analysed using a commercially available antibody array system (RayBio® C-Series Human Inflammation Array C1 Map, and C-Series Human Growth Factors Array C1 RayBiotech, Inc. Peachtree Corners, GA, USA) that uses membrane-bound cytokine-specific antibodies to capture proteins in biological fluids. The procedure was performed according to the manufacturer’s instructions. Cells were seeded in 10 cm culture dishes and treated (or not) with vismodegib and sonidegib 10 µM for 72 h. After removing the drug, cells (and control untreated proliferating fibroblasts) were maintained in serum-free medium for 48 h before collect conditioned medium. The cytokine array membranes were blocked in 1 × blocking buffer for 30 min at room temperature (RT) and then were incubated with 1 mL of conditioned medium at 4 °C overnight. The medium was then decanted from each container, and the membranes were washed three times with 1 × wash buffer I, followed by two washes with 1 × wash buffer II at RT. Next, the membranes were incubated in biotin-conjugated primary antibodies for 2 h at RT and then washed as described above before incubation in 1:1000-diluted horseradish peroxidase-conjugated streptavidin for 2 h. The membranes were then washed thoroughly and incubated with a chemiluminescent ECL substrate at RT for 5 min. Imaging and densitometric analysis were performed with a UVITEC Mini HD9 acquisition system (Alliance UVItec Ldt, Cambridge, UK).

2.9. Quantification of Vismodegib and Sonidegib Uptake

Cell uptake was determined in three independent experiments performed in three different cell lines. Cells were harvested and counted after 72 h treatments and stored at –80°C until analysis. Cell pellets were extracted after three freezing/thawing cycles in liquid nitrogen and addition of 200 pmoles of d31Cer[NS] 34:0 as the internal standard (iSTD) added in isopropanol containing 0.001% BHT. The cell suspension was extracted twice with 1 mL of ethyl acetate. The upper organic phase was evaporated under nitrogen and the dry extract was dissolved in 200 µL isopropyl alcohol before injection. Calibration curves were prepared with the vismodegib and sonidegib authentic compounds in the concentration range 2–400 nM with the iSTD used for the extraction. Vismodegib and sonidegb were quantified with reversed phase liquid chromatography coupled to a time-of-flight mass spectrometry (LCMS) system (Agilent Technologies, Santa Clara, CA, USA). For the quantitation, the protonated ions [M + H]+ were extracted for vismodegib (m/z 421,0175), sonidegib (m/z 486,1999), and the iSTD (m/z 568,7077). The method was linear (R > 0.9999) in the concentration range between 2–400 nM for both vismodegib and sonidegib.

2.10. Statistical Analysis

Results in the figures are representative of several experiments we performed with at least five cell lines from different donors for each cell type. Quantitative data were reported as mean ± standard deviation (SD). Student’s t-test was used to assess statistical significance with thresholds of * p ≤ 0.05 and ** p ≤ 0.01. Array card results reported in Table 2, Table 3 and Table 4 were obtained with a one-way ANOVA statistical test with thresholds >2.0 and <0.5 fold-change and p value < 0.05.

2.11. Ethics Statement

This study was approved by the Institutional Research Ethics Committee (Istituti Regina Elena e San Gallicano), Rome, Italy, and was performed after informed written consent to collect samples of human material for research.

3. Results

3.1. Patients’ Demographic and Genetic Characterization

A total of 10 NBCCS patients were included in the study. Patients were first diagnosed by clinical examination and family history, and then the genetic characterization confirmed Gorlin Syndrome diagnosis. Consistent with previous evidences demonstrating the absence of hot spot mutation in PTCH1 gene [51,52,53,54], molecular analysis of patients enclosed in this study recorded heterozygotic genetic variation distributed along the entire coding sequence (Table 1). Among these mutations, eight were missense changes resulting in harmful or pathogenic effects based on data from online databases (https://varsome.com, https://clinvarminer.genetics.utah.edu and https://www.lovd.nl, accessed on 20 June 2021), 1 was a deletion resulted in a frameshift of the coding sequence, and one insertion of a surplus CGG triplet was at the site of seven times CGG repeat in a functionally important 5′UTR intron region involved in PTCH1b variant transcription and translation [55,56].
Table
Table 1. Mutations found in nevoid basal cell carcinoma syndrome (NBCCS) cases examined.
Table 1. Mutations found in nevoid basal cell carcinoma syndrome (NBCCS) cases examined.
PatientsAgeSexAffected ExonDNA/Protein ChangeEffect on Sequence
NBCCS-HF163FExon 17c.2635G > T p.Asp879Tyrmissense
NBCCS-HF276MExon 12c.1510T > C (p.Leu503Ser)#missense
NBCCS-HF625FExon 10c.1309G > A (p.Val437Ile)missense
NBCCS-HF757FExon 15c.1510T> C (p.Leu503Ser)#missense
NBCCS-HF864FExon 5c.653delA (p.Gln218Gln fsX219)deletion
NBCCS-HF969MExon 1c.113G > T (p.Gly38Val)missense
NBCCS-HF1078FExon 24c.4172G > A (p.Arg1391Gln)missense
NBCCS-HF1132MExon 10c.1309G > A p.Val437Ilemissense
NBCCS-HF1257FExon 2c.-6_-4dupGGC[1]5′ untranslated region cis-regulatory element
NBCCS-HF1356MExon 19c.3227T > C (p.Ile1076Ser)missense
Mutation analysis of PTCH1 gene was performed by direct sequencing. Data were matched to https://varsome.com, https://clinvarminer.genetics.utah.edu and https://www.lovd.nl (accessed on 20 June 2021) databases to assess the pathogenicity of gene variants. # Indicated patients of the same family.

3.2. Comparative In Vitro Evaluation of Vismodegib and Sonidegib on Fibroblasts Proliferation

To evaluate the cytotoxicity of SMO antagonist, a 72 h MTT assay was performed in NHFs, NBCCS-HFs and CAFs. No obvious cytotoxicity was observed with vismodegib at all the concentration tested. However, with higher concentrations used (20–100 µM), it reduced cell proliferation of every cell types analysed in a dose-dependent manner. By contrast, sonidegib exerted a global stronger effect leading to a progressive decrease of the proliferation rate with lower doses and near to complete dead of cell culture exposed to 50 µM and 100 µM (Figure 1a). Differences in vismodegib and sonidegib toxicity were confirmed by quantification of apoptotic rate by Annexin V/PI staining (Figure 1b). Overall, no significant differences emerged comparing sensibility to both Hh inhibitors in normal, syndromic and BCC-associated fibroblasts. Furthermore, the specificity of SMO-dependent effect was confirmed by small interfering RNA (si-RNA) experiments. In this case, a mean 0.68 ± 0.22 and 0.53 ± 0.16 reduction in SMO mRNA in NHFs and CAFs respectively, impacted on cell proliferation compared to non-specific (NS-RNA). By contrast, syndromic fibroblasts with comparable decrease level of SMO mRNA (0.45 ± 0.02), showed no change in the MTT assay (Supplementary Figure S1a).

3.3. Analysis of Smoothened (SMO) Antagonists on Normal Human Fibroblasts (NHFs), NBCCS-HFs and CAFs Gene Expression Profile

To further evaluate the biological effect of Hh inhibitors on fibroblasts avowing acute cell damage, was selected the dose of 10 µM. We compared the gene expression profile of 9 independent NHFs strains and 9 NBCCS-HFs treated (or not) with vismodegib and sonidegib for 72 h. The study also included 6 BCC-associated CAF cell lines under the same treatment regime. A total of 93 genes of interest related to the Hh pathway, Wnt pathway and to pro-tumoral fibroblast phenotype were evaluated using the gene expression array card system. First, to select genes dissimilarly expressed in untreated cells at basal level, samples were clustered in three different biogroups according to the donor type (NHF, NBCCS-HF and CAF). The level of mRNAs in healthy fibroblasts was used as reference and arbitrarily indicated as 1.0. Among selected genes were identified 24 mRNAs with over two-fold significant difference (>2.0 and p-value < 0.05) in expression levels between NBCCS and NHF samples (19 up-regulated and 5 down-regulated) (Table 2 and Figure 2).
Table
Table 2. Gene expression analysis, 72 h.
Table 2. Gene expression analysis, 72 h.
GeneNBCCS-HFNBCCS-HF NBCCS-HFCAFCAFCAFNHFNHFNHF
VisSonVisSonVisSon
IL1α252.348.123.510.02.16.11.01.01.7
HHIP195.5122.942.74.86.01.91.00.320.54
IL1854.410.911.55.113.520.01.02.72.7
A2M47.121.28.454.257.952.81.01.42.2
CSF213.66.728.40.570.290.621.01.30.52
IGF210.66.611.87.47.68.41.00.870.96
CCL59.34.33.63.43.83.41.01.31.6
IGF17.34.78.06.66.76.81.00.631.3
CCL27.17.26.92.51.71.81.00.580.59
ICAM16.54.44.12.63.22.71.00.950.7
MTSS16.05.53.92.92.32.81.00.830.32
HGF5.12.52.52.00.871.61.00.820.71
MMP14.02.73.30.910.531.11.00.960.61
IL63.73.63.32.52.32.41.01.10.94
IL1β11.53.93.70.550.190.601.00.320.17
SFRP22.97.59.13.95.44.11.00.420.32
VCAM12.410.16.90.252.11.11.00.590.29
CXCL162.44.53.72.92.82.81.00.740.98
CSF31.26.15.50.750.270.121.02.10.34
Gli22.10.71.93.03.31.91.01.51.4
FGF92.02.23.04.02.22.01.01.01.5
PTCH21.11.51.42.11.82.41.01.84.9
WIF10.180.460.750.030.030.021.02.10.08
PTGER30.20.570.450.280.240.221.00.820.7
CES10.411.62.04.33.43.51.01.30.8
Gli30.470.640.840.660.580.641.01.00.91
BDKRB20.470.540.630.670.560.571.00.720.62
SMO1.70.770.872.91.71.41.00.791.4
Wnt7a1.11.63.250.260.280.311.00.540.57
CASP10.590.90.690.450.380.351.01.00.65
PTGIS1.11.71.10.140.160.191.00.680.97
CD401.31.62.70.960.751.01.00.770.59
LTC4S0.870.622.31.31.11.31.00.851.6
Gene expression analysis NBCCS-HFs and CAFs compared to NHFs. Data also report the effect of 72 h treatment with vismodegib and sonidegib (10 µM). Significantly overexpressed mRNA > 2.0 fold-change and p < 0.5 are marked in red, whereas significantly downregulated mRNA < 0.5 fold-change and p < 0.5 are marked in blue. The level of mRNAs in healthy fibroblasts was used as reference and arbitrarily indicated as 1.0.
Six transcripts correspond to genes involved in inflammation (IL1α, IL18, IL1β, IL6, CXCL16, and CCL2), and four mRNAs correspond to growth factors (FGF9, IGF1, IGF2 and HGF). In addition, members of the immunoglobulin superfamily of endothelial adhesion molecules, vascular cell adhesion molecule (VCAM-1), intercellular cell adhesion molecule (ICAM-1), and colony-stimulating factors 2 (CSF2) which stimulate granulocytes and macrophages were increased. Interestingly, Alpha-2-Macroglobulin (A2M) a protease inhibitor and cytokine transporter resulted significantly higher. Due to the capacity of A2M to disrupt inflammatory cascade, its high expression could be explained by the marked pro-inflammatory profile of Gorlin’s patient fibroblasts. Human hedgehog interacting protein (HHIP), a negative regulator of Hh signalling frequently underexpressed in several tumours [57,58,59,60], showed an extremely high level of expression in NBCCS-HFs suggesting a compensatory mechanism acting on Hh pathway. Moreover, our data confirmed previous studies [49,50,51] demonstrating high level of matrix metalloproteinase 1 (MMP1) and of Secreted Frizzled Related Protein 2 (SFRP2), a soluble modulator of Wnt signalling. By contrast, Wnt Inhibitory Factor 1 (WIF1) an inhibitor of the same pathway, Prostaglandin E Receptor 3 (PTGER3), and Bradykinin receptor B2 (BDKRB2) were significantly lower (fold-decrease < 0.5 and p < 0.05) in NBCCS-HFs than in normal fibroblasts (Table 2). Another sonic hedgehog target gene, metastasis suppressor 1 (MTSS1), a cytoskeletal associated protein implicated in cell motility and invasiveness [61,62,63], resulted significantly up-modulated in syndromic fibroblasts. Interestingly, MTSS1 is involved in carcinogenesis not only due to an augmented metastatic capability but also because MTSS1 serves as a co-transcription factor by binding to GLI proteins enhancing Hh target genes transcription [64]. Carboxylesterase1 (CES1), an enzyme involved in lipids hydrolysis, lipoproteins assembly and fatty acyl and cholesterol ester metabolism [54], also implicated in tumour-killing activity of monocytes [65] showed significant lower level in syndromic fibroblasts. A pronounced pro-inflammatory signature of Gorlin syndrome cells was additionally evidenced by the higher mRNA level of several other interleukins and cytokines (CXCL10, IL8, CCL5, CCL19, TNFα and PTGS2) even if this results set did not fully reach statistical significance (>2.0 and p ≥ 0.05 e.g., insignificant). GLI2 and GLI3 displayed opposite regulation (up- and down-regulation, respectively). This is in line with the idea that GLI2 is a positive regulator (activated by SHH) whereas GLI3 is a repressor (reduced by Shh) of Hh targets [6,66]. Overall mRNA profile proven stable Hh over-activation in fibroblasts carrying heterozygotic PTCH1 mutation that persists in vitro. The protein array system used to compare some of the inflammatory proteins released in culture medium by healthy and syndromic cells confirmed differences in IL6 expression (fold-increase 11.3 ± 13.9; p = 0.032) and CXCL8 (fold-increase 3.72 ± 2.8; p = 0.055) (Figure 3a,b) and closely correlated to the mRNA results (Table 2 for IL6 and Table 3 for CXCL8).
Nevertheless, in the case of IL1α and β the corresponding proteins were barely revealed with this detection method. Secretion of TIMP2 protein, a modulator of MMPs, was also significantly higher (2.19 ± 1.1 p = 0.046) in syndromic cells confirming the propensity of NBCCS-HFs to actively remodel the extracellular matrix. Similarly, an array panel for the detection of some mitogenic peptides revealed lessened basic fibroblast growth factor (bFGF) (fold-decrease 0.33 ± 0.30 p = 0.014) and higher insulin growth factor binding protein 6 (IGFBP6; fold-increase 2.1 ± 0.9 p = 0.049) in NBCCS-HFs with respect to normal fibroblasts (Figure 3c,d). Interestingly, CAFs isolated from sporadic BCCs showed an overall expression profile similar to NBCCS-HFs (Table 2). Among the 19 mRNAs up-regulated in syndromic fibroblasts, 15 demonstrated analogous modulation in BCC-associated fibroblasts (IL1α, CCL2, CXCL16, IL6, IL18, ICAM, A2M, SFRP2, FGF9, IGF2, IGF1, HGF, HHIP, MTSS1 and GLI2). As is likely in syndromic patients, CAFs show a moderate increase of additive inflammatory markers (CCL5, and CXCL10 e.g., insignificant). Significant up-modulation of PTCH2 and SMO transcripts confirmed the activation of Hh signalling in BCC-associated CAFs and suggested this anomaly as a possible priming mechanism for Hh signalling activation in absence of genetic aberration. In line with this idea, Walter and collaborators reported SMO overexpression as a mechanism for the activation of Hh signalling in human pancreatic CAFs [48]. NBCCS-HFs and CAFs also share the reduced expression of WIF1 and PTGER3. CAFs additionally presented significant lower level of Wnt7a, caspase-1 and prostacyclin synthase (PTGIS) compared to normal fibroblasts. By contrast with syndromic fibroblasts, CAFs express augmented level of CES1. More data regarding level of mRNA expression uniformly expressed in normal and pathological cells are summarized in Table 3.
Table
Table 3. Gene expression analysis, unmodified genes.
Table 3. Gene expression analysis, unmodified genes.
GeneNBCCS-HFNBCCS-HFNBCCS-HFCAFCAF CAFNHFNHFNHF
VisSonVisSonVisSon
a-SMA0.730.881.00.580.540.631.01.00.86
FAP1.11.231.81.00.860.941.01.20.98
Wnt10.983.10.980.440.470.301.00.840.54
Wnt2b1.20.960.770.870.820.921.01.11.1
Wnt5a1.51.21.50.20.420.471.01.01.0
Wnt63.11.212.30.670.340.951.00.590.76
Wnt9a0.90.400.250.880.880.881.00.510.75
Wnt9b0.010.220.271.00.220.121.00.170.53
Wnt10b0.210.540.210.090.510.1221.00.850.59
Wnt110.531.61.60.640.610.311.00.740.8
Wnt160.381.23.50.440.480.531.01.41.3
TIMP11.81.41.91.40.891.71.01.31.2
TIMP20.841.11.30.931.01.01.01.11.0
TGFβ10.721.10.770.770.771.21.01.00.97
VEGFA1.31.41.91.31.21.21.00.831.0
TNFSF13B1.852.761.861.61.91.51.00.981.2
TNFSF1B0.400.800.640.931.31.71.00.430.37
TNFSF1A0.910.781.20.540.490.491.00.630.75
TNF3.86.111.91.11.41.71.00.880.52
SUFU1.01.01.00.960.990.941.01.11.1
SMAD70.941.81.91.10.971.11.00.130.58
SMAD31.91.31.51.31.31.31.01.11.2
SHH0.110.751.10.291.10.161.00.614.2
SFRP11.10.340.471.91.21.91.00.770.25
PTGS25.45.12.61.81.91.61.00.830.43
PTGS11.41.01.30.540.510.361.00.910.76
PTGIR1.51.22.11.041.11.21.02.11.2
PTGFR0.380.590.890.70.770.591.01.11.2
PTGER20.861.31.70.460.540.321.00.410.30
PTCH10.920.690.921.31.41.31.01.00.95
PLA2G1B0.280.300.380.960.061.51.00.460.30
PDE4B0.480.590.670.500.440.491.00.930.97
PDE4A0.71.21.81.01.11.11.01.00.86
MMP30.861.01.80.740.60.541.01.20.81
MAPK80.810.801.00.800.810.861.01.00.85
MAPK30.860.881.10.870.891.11.00.981.1
MAPK140.760.891.10.690.670.721.00.950.99
MAPK10.890.910.920.640.650.621.00.810.74
IL70.340.763.40.520.961.51.00.951.2
IGFBP61.01.31.90.610.670.901.01.01.6
IGFBP30.890.981.60.640.500.691.00.801.2
Gli11.41.20.321.41.90.941.01.71.2
FGF70.980.931.51.12.41.41.01.21.1
FGF21.21.21.11.41.21.61.00.861.2
FAP1.11.21.81.030.860.941.01.20.98
EGF2.30.650.592.31.70.321.02.70.30
DKK30.870.590.710.540.530.601.00.961.2
DKK10.530.720.911.70.770.591.00.80.93
CXCL819.27.49.51.81.40.821.00.71.4
CXCL103.23.914.119.03.31.31.02.81.4
CTNNB10.941.51.80.821.01.11.01.21.0
CSF10.710.971.10.810.720.791.01.11.0
COL1A20.580.991.11.20.660.621.01.31.6
COL1A11.20.880.931.11.61.41.01.31.4
CCNB11.51.81.20.850.890.531.00.860.65
CCL30.381.23.50.440.480.911.01.41.4
CCL197.10.646.11.30.240.641.01.71.8
CACNA2D10.861.31.22.21.21.61.01.51.0
CACNA1C0.951.10.641.31.31.41.01.21.2
BDKRB10.480.921.10.470.480.371.01.30.77
Summary of data regarding mRNA that resulted similarly (e.g., insignificant or flat) expressed in NHFs, NBCCS-HFs and CAFs at basal level or that are not significantly modified by SMO inhibitors treatment (72 h). The level of mRNAs in healthy fibroblasts was used as reference and arbitrarily indicated as 1.0.
For further statistical analysis, data were pooled in 9 separated biogroups according to cell cultures origin and treatment regime. After 72 h, inhibition of SMO partially reverted the above described “Gorlin syndrome fibroblast expression profile” (Table 2). Nonetheless, most of the deregulated genes did not reached expression level similar to healthy fibroblasts. Notably, only the mRNA level of GLI2 was fully normalized by SMO antagonists. In detail, both vismodegib and sonidegib significantly (two-fold difference and p < 0.05 NBCCS-HFs-treated vs. untreated) counteract deregulation of IL1α, IL1β, IL18, HGF, A2M (by decreasing), CES1 and WIF1 (by increasing) in NBCCS-HFs (Table 2 and Figure 4) that, however, resulted still significantly higher compared to healthy fibroblasts.
Additionally, hyper expression of HHIP, ICAM, MTSS1, MMP1 resulted attenuated by treatments (Table 2). By contrast, the presence of Hh inhibitors exacerbated SFRP2 and VCAM1 overexpression. Also transcripts of Wnt11, CSF3 and leukotriene synthase 4 (LTC4S), comparable at basal level in NBCCS-HFs and healthy fibroblasts, was stimulated by both compounds. Once again, in line with MTT assay results, sonidegib exerted a comprehensive stronger effect compared to vismodegib. For example, focusing on mRNA variation of the four highest up-modulated genes and responsive to the activity of both molecules (IL1α, HHIP, IL1β, and A2M), we observed a mean fold-change of 0.37 ± 0.18 and 0.17 ± 0.05 with vismodegib and sonidegib respectively. Reversion of WIF1, PTGER3, CES1, GLI3 and BDKRB2 hypo-expression in the presence of vismodegib or sonidegib was extremely similar with both drugs (2.36 ± 1.0 and 2.88 ± 1.39). Surprisingly, in CAFs pharmacological inhibition of Hh pathway exerted a moderate effect on the pro-tumorigenic profile (Table 2). IGF1, ICAM, MTSS1, IL6, SFRP2, and CXCL16, were unaffected by treatments and still resulted significantly higher than in normal fibroblasts. Under treatment with both SMO antagonists, IL1α, CCL2, HGF, FGF9 and ICAM1 mRNA tended to expression level measured in healthy fibroblasts, whereas deregulation of IL18 and VCAM were significantly (p < 0.05) aggravated by Hh inhibition (Figure 5).
HHIP and GLI2 gene expression was significantly repressed only in the presence of sonidegib. Notably, in the case of CSF3 we observed divergent regulation in NBCCS-HFs (increase) and CAFs (decrease) (Table 2). The mRNA expression of the major part of inflammatory mediators were not influenced by vismodegib and sonidegib in CAFs, suggesting that in this type of cell this paracrine character is not necessarily driven by Hh signalling. In line with the idea that Hh intracellular signalling is mostly inactive in adult healthy tissue, treatment of normal fibroblasts failed to show any relevant effect with the exception of a IL18 stimulation, a mild but significant decrease of SFRP2 and IL1β with both inhibitors (Figure 6).
Some additional Hh target genes were modified by vismodegib and sonidegib alternatively (Table 2). Finally, the specificity of the effect on mRNA profile observed in the presence of compounds were confirmed by siRNA experiment and following RT-PCR of a sampling of the three most representative mRNA modulated by both molecules (HHIP, IL1β and MTSS1) (Figure S1b and Table 4). Genetic-imposed Hh signalling over-activation represents a quite stable situation that might necessitate a prolonged treatment to achieve prominent repression of SMO activity. Thus, we used a 2-week treatment to simulate the effect of stable Hh inhibition similar to continuative therapeutic regime used in the clinic. Overall, longer treatment impacted firmly on syndromic fibroblast phenotype presenting 10 of the 19 genes responsive to treatments (IL1α, IGF2, IGF1, CCL2, ICAM, MTSS1, MMP1, IL6, SFRP2, and CSF3) with level of expression nearer to normal cells respect to cells treated for 72 h (Table 4 and Supplementary Figure S3).
Table
Table 4. Gene expression analysis, 2 weeks.
Table 4. Gene expression analysis, 2 weeks.
GeneNBCCS-HFNBCCS-HFNBCCS-HFNBCCS-HFNBCCS-HF
Vis 72 hVis 2 weeksSon 72 hSon 2 weeks
IL1α252.348.12.423.54.0
HHIP195.5122.949.042.772.3
IL1854.410.960.711.51.1
A2M47.121.213.08.413.1
CSF213.66.77.828.42.2
IGF210.66.61.911.83.7
CCL59.34.30.283.63.3
IGF17.34.70.88.02.9
CCL27.17.22.06.90.9
ICAM16.54.43.54.12.7
MTSS16.05.51.23.90.28
HGF5.12.54.22.52.5
MMP14.02.71.43.31.8
IL63.73.61.23.30.6
IL1β11.53.98.23.77.5
SFRP22.97.50.69.10.4
VCAM12.410.16.46.911.9
CXCL162.44.52.83.70.49
CSF32.26.10.175.50.28
Gli22.10.521.81.93.1
FGF92.02.21.63.03.6
PTCH21.11.50.21.41.0
WIF10.180.460.120.750.06
PTGER30.20.571.40.451.3
CES10.411.60.82.00.83
Gli30.470.640.650.840.66
BDKRB20.470.540.880.630.62
SMO1.70.771.50.871.7
Wnt7a1.11.60.433.250.7
CASP10.590.90.560.690.29
PTGIS1.11.70.351.10.71
CD401.31.60.822.70.66
LTC4S0.870.621.12.31.6
Gene expression analysis NBCCS-HFs, CAFs compared to NHFs at basal level and following 2 weeks treatment with vismodegib and sonidegib (10 µM). Significantly overexpressed mRNA > 2.0 fold-change and p < 0.5 are marked in red, whereas significantly downregulated mRNA < 0.5 fold-change and p < 0.5 are marked in blue. The level of mRNAs in healthy fibroblasts was used as reference and arbitrarily indicated as 1.0.
In particular, deregulation of IL6, MTSS1, and MMP1 was completely counteracted after 2 weeks of continuous exposure to SMO antagonists. Notably, vismodegib and sonidegib overturned SFRP2 and CSF3 expression only after long treatment whereas, after 72 h treatment, the expression of these two genes resulted extremely high. Long treatment (2 weeks) partially abrogated differences in vismodegib and sonidegib efficacy observed at 72 h. A deeper examination of Table 3 provided additional information on NBCCS-HF and BCC-associated CAF phenotype. Both cell types failed to display a significant different expression of some typical CAF’ markers such as α-smooth muscle actin (α-SMA) and fibroblast activation protein (FAP) at mRNA and protein level suggesting that these genes are not necessarily involved in skin basaloid proliferation promotion by mesenchymal cells (Table 3 and Figure 7a). Similarly, additional evaluation by RT-PCR did not provide evidence of a difference in the level of fibroblast-specific protein 1 (FSP1) and platelet-derived growth factor receptor (PDGFR)-α and β in syndromic fibroblasts (Figure 7b). Consistent with our data, previous transcriptomic analysis of Hh target genes in NBCCS fibroblasts demonstrated no difference in PDGFRα mRNA level compared to healthy cells [49].
Numerous studies have highlighted the importance of molecular crosstalk between the Hh and other signalling cascades, including Wnt pathway. Wnt/β-catenin signalling contributes to the homeostasis of normal skin cells, and its malfunction is implicated in tumour initiation, progression and invasion, but also maintains cancer stem cells which contribute to tumour recurrence [67]. The expression of β-catenin, the key protein regulating Wnt pathway intracellular activity, showed a high heterogenic expression among all fibroblast cell cultures. Altogether, Western blot and densitometric analysis did not evidence statistical difference in the basal level of β-catenin expression in NBCCS-HFs compared to NHFs (data not shown). Vismodegib and sonidegib treatment down-regulates β-catenin at 24 and 72 h in normal fibroblasts whereas NBCCS-HFs resulted less affected especially with the shorter treatment (Figure 7c,d). The mechanism of β-catenin protein reduction did not involve modification of the corresponding mRNA (Table 3), and was associated to an augmented activity of glycogen synthase kinase-3β (GSK3β) since the inactive pSer9-GSK3β was lower in treated cells compared to untreated cells suggesting an accelerated β-catenin ubiquitin-dependent degradation (Figure 7c,d). Correspondingly, pSer473-AKT, an upstream regulator of GSK3β phosphorylation and activity, decreased suggesting the axis pAKT-GSK3β for the control of β-catenin protein turnover and confirming an attenuated effect in the case of syndromic fibroblasts. Again specific interfering of SMO expression confirmed the specificity of down-modulation of the axis pAKT-GSK3β-β-catenin (Figure S1c). Vismodegib and sonidegib treatment exerted a similar effect on CAFs (Figure 7d). Consistent with the synergic regulation of Hh and Wnt pathways, the crosstalk protein SFRP1 carrying a GLI-binding site [68,69] resulted as decreased by SMO antagonists in NBCCS-HF, CAFs and healthy fibroblasts (Table 3). The expression level of β-catenin nuclear partners, the co-transcription factors LEF1, TCF1 and TCF4, did not show any significant variation related to pharmacological treatments (Figure 5). However, a tendency of augmented basal level of expression of LEF/TCF transcription factors was observed in syndromic cells compared to normal cells. The analysis of a wide range of Wnt secreted factors (Wnt5a, Wnt2a, Wnt6, Wnt7a, Wnt9a, Wnt9b, Wnt10b and Wnt11) and members of Dickkopf protein family (DKK1, DKK2 and DKK3) implicated in Wnt/β-catenin signalling miss to identify any specific trait of NBCCS-HFs or variation under Hh inhibitors exposure (Table 2 and Table 3).

3.4. Vismodegib and Sonidegib Uptake Revealed Profound Differences

Vismodegib and sonidegib have been proved to have analogous clinical benefits [36,40]. Herein, disparity in SMO-inhibitory properties observed in vitro, especially the toxic one, represents an unexpected result since these molecules are reported to be similar in terms of IC50 for Hh inhibition in human cell-free assays (3.0 µM and 2.5 µM respectively) [70,71]. In GLI-luciferase assay in human embryonic palatal mesenchyme (HEPM) cells, vismodegib inhibited Hh signalling with an IC50 of 2.8 nM whereas sonidegib shown an IC50 of 12.7 nM [72]. By contrast, in a model using serum-starved human medulloblastoma cells treated with SHH, vismodegib resulted two-fold less potent than sonidegib in restraining GLI1 mRNA production [71]. Moreover, sonidegib displayed similar effectiveness to vismodegib at inhibiting GLI1 and PTCH1 mRNA in in vitro cell cultures and in topically-treated depilated mouse skin [73]. In the same study, SMO-mediated G-protein activation assay and medulloblastoma cells proliferation confirmed similar potency [73]. Thus, overall previous comparative data did not bring out significant inequality in vismodegib and sonidegib effectiveness. Since appreciable off-target effects have been previously excluded at concentrations up to 10 µM [71], we next evaluated the cell internalization of Hh inhibitors by LCMS analysis after 3, 6 and 24 h of drug exposure. At all experimental time points, results evidenced an apparent disparity in the concentration of sonidegib and vismodegib. The sonidegib to vismodegib concentration ratio accounted for 73.3 ± 27.5; 15.1 ± 3.0 and 21.1 ± 1.22 fold-difference at 3, 6 and 24 h respectively (Figure 8). The concentration of vismodegib in the supernatant remained unchanged over time; by contrast, that of sonidegib showed a time-dependent decline opposite to the increasing cell uptake (data not shown).

4. Discussion

Antagonistic activity on Hh signalling represents an important therapeutic opportunity for local advanced, metastatic BCCs and for NBCCS patients presenting a dramatic predisposition to develop multiple BCC during their entire life. Despite the motivating use in therapy, in vitro models are rarely reported to study SMO inhibitors since cell lines are not considered to be predictive in the clinic because Hh signalling is mostly silent in adult tissue and the activation of the Hh pathway is not necessarily conserved in pathological cell cultures. In the case of tumour cells, this could be explained by the absence of the microenvironment, specifically fibroblasts that secrete Hh ligands. Moreover, the fact that multiple intracellular signallings susceptible to culture medium composition are implicated in non-canonical Hh pathway activation might complicate data interpretation. In this study, NBCCS-HFs presenting heterozygous germline mutation of PTCH1 and low passage BCC-CAF cell lines maintained a specific signature ascribable to Hh signalling activity demonstrating the reliability of the model used. Data concerning NBCCS patients confirmed that disease phenotypes may arise from small reductions in PTCH1 activity as in the case of heterozygotic setting presenting one functional copy. All NBCCS-HF cell lines showed a similar Hh signalling signature and comparable sensibility to Hh inhibition (data not shown). The read-out of Hh signalling activation in the oncogenic field is currently in progress. Hence, the definition of its reversion is still being debated. Some groups proposed incremented GLI1 as a reliable marker [74,75,76], whereas others used a panel of multiple genes (two or more), such as GLI1, PTCH1, HIP, cyclinD1, PTCH2, SFRP1 and N-myc to define Hh level of activity [77,78,79]. However, the findings of tumour cells may not be necessarily applicable to normal cells or syndromic fibroblasts carrying a heterozygous loss of function PTCH1. In line with this concept, our results showed that the assessment of Hh pathway activity based on GLI and PTCH1 genes expression is not applicable to mesenchymal cells. GLI transcription factors have overlapping as well as distinctive functions. In humans, GLI2 is the primary activator, GLI3 (lacking the C-terminal transactivation domain) is the primary repressor [80,81], and GLI1 is a target gene that acts as an activator in a positive feedback loop [82]. In this study, GLI1 expression resulted as slightly but not significantly increased in untreated NBCCS-HFs and CAFs compared to healthy fibroblasts. Thus, it is possible that moderate overactivity of Hh signalling in non-tumoral cells does not represent a sufficient input to trigger further amplification of the signalling cascade. Upregulation of GLI2 in CAFs and NBCCS-HFs strengthens the proposed idea that GLI2 regulation has a dominant role on the Hh signalling dynamics [83]. At the same time, a lower abundance of GLI3 in the context of endogenous Hh activation endorses a repressor role for this DNA-binding factor. In this study, the most relevant marker of Hh pathway activation is HHIP overexpression. HHIP is a downstream target gene of Hh signalling responsible for a negative feedback loop due to the capacity to bind and sequester all three Hh ligands away from PTCH receptors [58,84,85,86]. This mechanism probably plays an essential role in moderating the signalling and maintaining a physiological equilibrium in cells of syndromic patients. Accordingly, HHIP expression was forcefully modulated by Hh pathway inhibition and SMO siRNA. Additional target genes, such as CyclinD1, PDGFRα, VEGF, SFRP1, Wnt2 were not differentially expressed in syndromic fibroblasts compared to control fibroblasts confirming the distinctive feature of Hh activity in mesenchymal cells. CAFs feature has been documented in several genodermatoses [87]. Previous studies indicated that NBCCS-fibroblasts carrying constitutive activation of Hh intracellular signalling may promote basaloid proliferation by the production of a supportive tumour microenvironment similar to sporadic BCC [50]. Accordingly, the repertoire of growth factors overproduced by syndromic and BCC-associated fibroblasts fully overlaps (IGF1 and 2, HGF, FGF9). Of note, the production of bFGF is lower in syndromic cells compared to normal fibroblasts. A possible explanation for this observation may reside in the evidence that SMO and GLI1 genes are regulated by bFGF [88]. Thus, contraction of bFGF expression could represent a strategy to avoid detrimental excess in Hh signalling cascade.
A pronounced pro-inflammatory signature is a common character shared by NBCCS-HFs and BCC-CAFs which is consistent with previous studies indicating the Hh pathway as a regulator of cancer immune response [45,89]. However, this aspect seems to be attributable to Hh deregulation only in NBCCS cells, since the expression of the major part of cytokines and interleukins were not perturbed by vismodegib and sonidegib in BCC-stromal fibroblasts. Of note, divergent from other immune response-regulating genes, the expression of CCL2 was rapidly reduced by treatments in CAFs but not in BCC-associated fibroblasts that necessitated prolonged treatment to restore physiological level of this chemokine. Interestingly, the level of CCL2 correlates with the presence of tumour-associated polarized M2 macrophages and poor prognosis in prostate and breast cancer [90,91]. In line with our data, CAFs characterized by the secretion of a high level of pro-inflammatory cytokines are reported to be weakly positive for α-SMA expression [42]. However, this observation seems to be in contrast with a very recent study demonstrating that Hh signalling activity correlates directly to the amount of myofibroblastic CAFs (presenting elevated α-SMA expression) and indirectly to the amount of inflammatory CAFs in pancreatic adenocarcinoma [4].
Another consistent difference between fibroblasts associated with sporadic BCC and syndromic fibroblasts resides in the characteristic over-expression of two Hh signalling components: SMO and PTCH2. Our data are consistent with recent studies in human and mouse models of pancreatic cancer proposing SMO overexpression as a mechanism of activation of this pathway in tumour stroma [48,92]. Of note, even if the global magnitude of Hh target genes deregulation in CAF is lower than in syndromic cells, the efficacy of inhibitors resulted as minor in non-syndromic cells. This is probably related to the fact that targeting CAFs phenotype might be difficult due to the frequent epigenetic permanent modification associated to fibroblast attitude [93,94].
A robust increase of SHH and PTCH2 in normal cells exposed to SMO inhibitors probably aims to restore the minimal basal level of intracellular Hh signalling proving the non-negligible role of minute basic Hh activity in adult cells. However, the overall impact of SMO antagonist treatment was feeble in normal fibroblasts. In line with this idea, data from clinical trials showed that Hh target inhibition in tumours is more pronounced than in the matched normal skin when both tissues were investigated [76]. Opposite deregulated expression of CES1, an enzyme implicated in cellular cholesterol and fatty acids levels [95], might reveal the essential difference of Hh pathway deregulation in NBCCS-HFs and in CAFs, since cholesterol and its oxygenated derivate oxysterols are agonists for SMO [96,97,98] also in absence of Hh ligands [99]. On the contrary, SMO activation is impaired through cholesterol deficiency and sterol depletion [98]. Thus, lipid metabolic abnormalities such as reprogram cholesterol metabolism in a tumor microenvironment may take part in Hh activation. Once again, lipid asset remodelling, in particular lowering the amount of accessible cholesterol (that anti-correlates with PTCH1 activity in normal cells) [100] may act to restrain Hh pathway activity. Additionally, since Hh ligands need to be covalently modified by cholesterol and lipid palmitate to achieve the full activation [101,102], CES-1 modulation might be potentially implicated in limiting SHH bioavailability in NBCCS cells. Vismodegib and sonidegib share essential overlapping Hh inhibitor class-dependent effects. In clinical investigations, these molecules demonstrated similar clinical efficacy, safety, and tolerability profiles [32,103,104]. However, due to the non-coincident study designs, the indirect comparison of pivotal clinical trials used for vismodegib and sonidegib approval (ERIVANCE and BOLT respectively) is enabled and no studies have been reported to evaluate these drugs versus a common comparator [105]. Biological differences between the two Hh inhibitors emerged mostly from pharmacokinetic profiles. First, sonidegib has a longer half-life than vismodegib (29.6 vs. 4–12 days) [106,107]. Both compounds are very highly bound to plasma proteins, but the binding is concentration-dependent for vismodegib [108] and concentration-independent for sonidegib [71,104]. This difference and the major grade of lipophilicity could explain the higher accumulation of sonidegib within the tissue [39,109]. According to clinical data, we observed that the distribution of vismodegib is mainly in the extracellular space whereas sonidegib abundantly penetrates the plasma membrane reaching intracellular concentration about 20-fold higher than vismodegib. The level of intracellular accumulation of SMO inhibitors should in theory be of limited relevance since vismodegib and sonidegib-interacting pocket involves the same SMO-binding cavity that consists in the extracellular stretch of heptahelical membrane domain and a cavity exposed in extracellular space [72,110]. However, our results correlating the intracellular concentration to the efficacy on target genes modulation, strongly suggest that Hh inhibitors may act in the cytoplasm too. A possible explanation is that sonidegib, as reported for two other SMO antagonists, SANT1 and SANT2 [111], blocks the transport of cytoplasmic SMO to the primary cilium of the cell, reducing the level of active SMO. In line with principle, vismodegib been structurally similar to sonidegib might be as well capable to prevent SMO recycling to the membrane but due to limited drug intracellular uptake fails to reach remarkable cytoplasmic concentration. Thus, while inhibiting SMO activity with the well-known canonical mechanism of action, some antagonists might engage cytoplasmic-inactive SMO preventing its primary cilium localization. Therefore, this represents a second mechanism by which inhibitors can attenuate SMO activity potentiating the pharmacological effectiveness. An augmented pharmacological effect for sonidegib, due to the secondary intracellular interaction with SMO, may explain the powerful effect at 72 h of treatment that is partially confirmed after 2 weeks of drug exposure. This augmented efficacy might be useful in therapy to reconsider active doses, reduce side effects, and enlarge the possible spectrum of target patients. Although vismodegib and sonidegib have only been approved for cases of advanced BCCs, the question has been raised whether they would be useful for the greater majority of BCCs that present the typical of-fice setting (i.e., non-advanced). Moreover, Hh inhibitors might be used as neoadjuvant therapy to decrease primary tumour size before excision reducing the complexity and extent of closure [112].

5. Conclusions

In conclusion, in this study we provided evidence that targeting the Hh pathway impacting on fibroblast tumour supportive functions might be considered a therapeutic option for BCC independent of the mutation status of Hh components in tumour cells. The advantage of targeting stromal cells in addition to neoplastic cells resides in a more durable response since normal cells have a relatively stable genetic constitution compared to cancer cells that easily develop secondary resistance to target therapies. However, further investigations will address the question of how CAF targeting could be improved. Also the evidence that a cancer-prone phenotype of NBCCS fibroblasts was significantly repressed by SMO inhibitors strengthens the use of this class of compound in the management of syndromic patients.

Supplementary Materials

The following are available online: https://www.mdpi.com/article/10.3390/cancers13225858/s1, Figures S1: (a) MTT assay of NBCCS-HFs, CAFs and NHFs transfected with specific SMO siRNA for 72 h. Data represent the average of experiments performed with 2 different SMO siRNA. (b) Semi-quantitative RT-PCR was performed to assess the level of SMO, IL1β, MTSS1 and HHIP level of expression in cells transfected with SMO siRNA (72 h). Histograms report mean ± SD of 2 set of experiments performed with 2 different SMO siRNA. (c) One representative western blot detecting β-catenin, pSer9-GSK3β and pSer473-AKT after 72 h treatment with siRNA (NS siRNA or SMO siRNA); Figure S2: Volcano plot represents gene expression analysis comparing NBCCS HFs treated with vismodegib (a) and sonidegib (b) for 2 weeks (n = 9). (b) Analysis was performed using three different housekeeping genes (β actin, GAPDH and 18S). One way ANOVA statistical test with thresholds >2.0 fold change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold change and p < 0.05 reported in green; any fold difference with p≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black; Table S1: Healthy donors characteristics; Table S2: BCC patients characteristics; Table S3: Primers for Real Time Quantitative Polymerase Chain Reaction.

Author Contributions

B.B. designed the study, conducted part of the experiments, analyzed data, and wrote the manuscript; L.E. patients selection, clinical data curation and critically revision of the paper. S.C., C.A. and E.C. performed part of the experiments. E.M., selected some BCC patients and collected all healthy tissue sample; V.S. and A.C. selected NBCCS and some BCC patients and collected tissue sample. C.C. examined all biopsies and separated perilesional and lesional area of skin cancers. P.F. coordinated the clinical unit. All authors have read and agreed to the published version of the manuscript.

Funding

Research was supported by the grant from the Italian Ministry of Health.

Institutional Review Board Statement

The Institutional Research Ethics Committee (Istituti Regina Elena e San Gallicano) approved all research activities involving human subjects (protocol number: Prot. CE/286/06).

Informed Consent Statement

The Declaration of Helsinki Principles was followed and patients gave written informed consent to collect samples of human material for research.

Data Availability Statement

The datasets used or analysed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We are grateful to Alessia Lupino for editing the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bath-Hextall, F.; Leonardi-Bee, J.; Smith, C.; Meal, A.; Hubbard, R. Trends in Incidence of Skin Basal Cell Carcinoma. Additional Evidence from a UK Primary Care Database Study. Int. J. Cancer 2007, 121, 2105–2108. [Google Scholar] [CrossRef] [PubMed]
  2. Telfer, N.R.; Colver, G.B.; Morton, C.A.; British Association of Dermatologists. Guidelines for the Management of Basal Cell Carcinoma. Br. J. Dermatol. 2008, 159, 35–48. [Google Scholar] [CrossRef] [PubMed]
  3. Mathias, S.D.; Chren, M.M.; Colwell, H.H.; Yim, Y.M.; Reyes, C.; Chen, D.M.; Fosko, S.W. Assessing Health-Related Quality of Life for Advanced Basal Cell Carcinoma and Basal Cell Carcinoma Nevus Syndrome: Development of the First Disease-Specific Patient-Reported Outcome Questionnaires. JAMA Dermatol. 2014, 150, 169–176. [Google Scholar] [CrossRef] [PubMed]
  4. Steele, N.G.; Biffi, G.; Kemp, S.B.; Zhang, Y.; Drouillard, D.; Syu, L.; Hao, Y.; Oni, T.E.; Brosnan, E.; Elyada, E.; et al. Inhibition of Hedgehog Signaling Alters Fibroblast Composition in Pancreatic Cancer. Clin. Cancer Res. 2021, 27, 2023–2037. [Google Scholar] [CrossRef]
  5. Wu, F.; Zhang, Y.; Sun, B.; McMahon, A.P.; Wang, Y. Hedgehog Signaling: From Basic Biology to Cancer Therapy. Cell. Chem. Biol. 2017, 24, 252–280. [Google Scholar] [CrossRef] [Green Version]
  6. Matise, M.P.; Joyner, A.L. Gli Genes in Development and Cancer. Oncogene 1999, 18, 7852–7859. [Google Scholar] [CrossRef] [Green Version]
  7. Regl, G.; Kasper, M.; Schnidar, H.; Eichberger, T.; Neill, G.W.; Ikram, M.S.; Quinn, A.G.; Philpott, M.P.; Frischauf, A.M.; Aberger, F. The Zinc-Finger Transcription Factor GLI2 Antagonizes Contact Inhibition and Differentiation of Human Epidermal Cells. Oncogene 2004, 23, 1263–1274. [Google Scholar] [CrossRef] [Green Version]
  8. Kasper, M.; Regl, G.; Frischauf, A.M.; Aberger, F. GLI Transcription Factors: Mediators of Oncogenic Hedgehog Signalling. Eur. J. Cancer 2006, 42, 437–445. [Google Scholar] [CrossRef]
  9. Bigelow, R.L.; Chari, N.S.; Unden, A.B.; Spurgers, K.B.; Lee, S.; Roop, D.R.; Toftgard, R.; McDonnell, T.J. Transcriptional Regulation of Bcl-2 Mediated by the Sonic Hedgehog Signaling Pathway through Gli-1. J. Biol. Chem. 2004, 279, 1197–1205. [Google Scholar] [CrossRef] [Green Version]
  10. Pelullo, M.; Zema, S.; Nardozza, F.; Checquolo, S.; Screpanti, I.; Bellavia, D. Wnt, Notch, and TGF-β Pathways Impinge on Hedgehog Signaling Complexity: An Open Window on Cancer. Front. Genet. 2019, 10, 711. [Google Scholar] [CrossRef] [Green Version]
  11. Bijlsma, M.F.; Roelink, H. Skin-Derived Vitamin D(3) Protects Against Basal Cell Carcinoma. J. Invest. Dermatol. 2017, 137, 2469–2471. [Google Scholar] [CrossRef] [Green Version]
  12. Cai, K.; Na, W.; Guo, M.; Xu, R.; Wang, X.; Qin, Y.; Wu, Y.; Jiang, J.; Huang, H. Targeting the Cross-Talk between the Hedgehog and NF-κB Signaling Pathways in Multiple Myeloma. Leuk. Lymphoma 2019, 60, 772–781. [Google Scholar] [CrossRef]
  13. Ng, J.M.; Curran, T. The Hedgehog’s Tale: Developing Strategies for Targeting Cancer. Nat. Rev. Cancer. 2011, 11, 493–501. [Google Scholar] [CrossRef] [Green Version]
  14. Epstein, E.H. Basal cell carcinomas: Attack of the hedgehog. Nat. Rev. Cancer. 2008, 10, 743–754. [Google Scholar] [CrossRef]
  15. Lam, C.W.; Xie, J.; To, K.F.; Ng, H.K.; Lee, K.C.; Yuen, N.W.; Lim, P.L.; Chan, L.Y.; Tong, S.F.; McCormick, F. A Frequent Activated Smoothened Mutation in Sporadic Basal Cell Carcinomas. Oncogene 1999, 18, 833–836. [Google Scholar] [CrossRef] [Green Version]
  16. Bonilla, X.; Parmentier, L.; King, B.; Bezrukov, F.; Kaya, G.; Zoete, V.; Seplyarskiy, V.B.; Sharpe, H.J.; McKee, T.; Letourneau, A.; et al. Genomic Analysis Identifies New Drivers and Progression Pathways in Skin Basal Cell Carcinoma. Nat. Genet. 2016, 48, 398–406. [Google Scholar] [CrossRef]
  17. Crowson, A.N. Basal Cell Carcinoma: Biology, Morphology and Clinical Implications. Mod. Pathol. 2006, 19 (Suppl. 2), S127–S147. [Google Scholar] [CrossRef]
  18. Heitzer, E.; Lassacher, A.; Quehenberger, F.; Kerl, H.; Wolf, P. UV Fingerprints Predominate in the PTCH Mutation Spectra of Basal Cell Carcinomas Independent of Clinical Phenotype. J. Investig. Dermatol. 2007, 127, 2872–2881. [Google Scholar] [CrossRef] [Green Version]
  19. D’Amato, C.; Rosa, R.; Marciano, R.; D’Amato, V.; Formisano, L.; Nappi, L.; Raimondo, L.; Di Mauro, C.; Servetto, A.; Fulciniti, F.; et al. Inhibition of Hedgehog Signalling by NVP-LDE225 (Erismodegib) Interferes with Growth and Invasion of Human Renal Cell Carcinoma Cells. Br. J. Cancer 2014, 111, 1168–1179. [Google Scholar] [CrossRef] [Green Version]
  20. Freitas, R.D.; Dias, R.B.; Vidal, M.T.A.; Valverde, L.F.; Gomes Alves Costa, R.; Damasceno, A.K.A.; Sales, C.B.S.; Siquara da Rocha, L.O.; Dos Reis, M.G.; Soares, M.B.P.; et al. Inhibition of CAL27 Oral Squamous Carcinoma Cell by Targeting Hedgehog Pathway with Vismodegib or Itraconazole. Front. Oncol. 2020, 10, 563838. [Google Scholar] [CrossRef]
  21. Benvenuto, M.; Masuelli, L.; De Smaele, E.; Fantini, M.; Mattera, R.; Cucchi, D.; Bonanno, E.; Di Stefano, E.; Frajese, G.V.; Orlandi, A.; et al. In Vitro and in Vivo Inhibition of Breast Cancer Cell Growth by Targeting the Hedgehog/GLI Pathway with SMO (GDC-0449) Or GLI (GANT-61) Inhibitors. Oncotarget 2016, 7, 9250–9270. [Google Scholar] [CrossRef] [Green Version]
  22. Williams, J.A. Hedgehog Signaling Pathway as a Target for Therapeutic Intervention in Basal Cell Carcinoma. Drug News. Perspect. 2003, 16, 657–662. [Google Scholar] [CrossRef]
  23. Maas, S.M.; Lombardi, M.P.; van Essen, A.J.; Wakeling, E.L.; Castle, B.; Temple, I.K.; Kumar, V.K.; Writzl, K.; Hennekam, R.C. Phenotype and Genotype in 17 Patients with Goltz-Gorlin Syndrome. J. Med. Genet. 2009, 46, 716–720. [Google Scholar] [CrossRef] [Green Version]
  24. Onodera, S.; Saito, A.; Hasegawa, D.; Morita, N.; Watanabe, K.; Nomura, T.; Shibahara, T.; Ohba, S.; Yamaguchi, A.; Azuma, T. Multi-Layered Mutation in Hedgehog-Related Genes in Gorlin Syndrome may Affect the Phenotype. PLoS ONE 2017, 12, e0184702. [Google Scholar] [CrossRef]
  25. Johnson, R.L.; Rothman, A.L.; Xie, J.; Goodrich, L.V.; Bare, J.W.; Bonifas, J.M.; Quinn, A.G.; Myers, R.M.; Cox, D.R.; Epstein, E.H., Jr.; et al. Human Homolog of Patched, a Candidate Gene for the Basal Cell Nevus Syndrome. Science 1996, 272, 1668–1671. [Google Scholar] [CrossRef] [Green Version]
  26. Evans, D.G.; Farndon, P.A.; Burnell, L.D.; Gattamaneni, H.R.; Birch, J.M. The Incidence of Gorlin Syndrome in 173 Consecutive Cases of Medulloblastoma. Br. J. Cancer 1991, 64, 959–961. [Google Scholar] [CrossRef] [Green Version]
  27. Kimonis, V.E.; Goldstein, A.M.; Pastakia, B.; Yang, M.L.; Kase, R.; DiGiovanna, J.J.; Bale, A.E.; Bale, S.J. Clinical Manifestations in 105 Persons with Nevoid Basal Cell Carcinoma Syndrome. Am. J. Med. Genet. 1997, 69, 299–308. [Google Scholar] [CrossRef] [Green Version]
  28. Shanley, S.; Ratcliffe, J.; Hockey, A.; Haan, E.; Oley, C.; Ravine, D.; Martin, N.; Wicking, C.; Chenevix-Trench, G. Nevoid Basal Cell Carcinoma Syndrome: Review of 118 Affected Individuals. Am. J. Med. Genet. 1994, 50, 282–290. [Google Scholar] [CrossRef]
  29. Lo Muzio, L. Nevoid Basal Cell Carcinoma Syndrome (Gorlin Syndrome). Orphanet J. Rare Dis. 2008, 3, 32. [Google Scholar] [CrossRef] [Green Version]
  30. Quinn, A.G.; Campbell, C.; Healy, E.; Rees, J.L. Chromosome 9 Allele Loss Occurs in both Basal and Squamous Cell Carcinomas of the Skin. J. Investig. Dermatol. 1994, 102, 300–303. [Google Scholar] [CrossRef] [Green Version]
  31. Tang, J.Y.; Ally, M.S.; Chanana, A.M.; Mackay-Wiggan, J.M.; Aszterbaum, M.; Lindgren, J.A.; Ulerio, G.; Rezaee, M.R.; Gildengorin, G.; Marji, J.; et al. Inhibition of the Hedgehog Pathway in Patients with Basal-Cell Nevus Syndrome: Final Results from the Multicentre, Randomised, Double-Blind, Placebo-Controlled, Phase 2 Trial. Lancet Oncol. 2016, 17, 1720–1731. [Google Scholar] [CrossRef]
  32. Skvara, H.; Kalthoff, F.; Meingassner, J.G.; Wolff-Winiski, B.; Aschauer, H.; Kelleher, J.F.; Wu, X.; Pan, S.; Mickel, L.; Schuster, C.; et al. Topical Treatment of Basal Cell Carcinomas in Nevoid Basal Cell Carcinoma Syndrome with a Smoothened Inhibitor. J. Invest. Dermatol. 2011, 131, 1735–1744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Sekulic, A.; Migden, M.R.; Lewis, K.; Hainsworth, J.D.; Solomon, J.A.; Yoo, S.; Arron, S.T.; Friedlander, P.A.; Marmur, E.; Rudin, C.M.; et al. Pivotal ERIVANCE Basal Cell Carcinoma (BCC) Study: 12-Month Update of Efficacy and Safety of Vismodegib in Advanced BCC. J. Am. Acad. Dermatol. 2015, 72, 1021–1026. [Google Scholar] [CrossRef] [PubMed]
  34. Chiang, A.; Jaju, P.D.; Batra, P.; Rezaee, M.; Epstein, E.H., Jr.; Tang, J.Y.; Sarin, K.Y. Genomic Stability in Syndromic Basal Cell Carcinoma. J. Invest. Dermatol. 2018, 138, 1044–1051. [Google Scholar] [CrossRef] [Green Version]
  35. Sharpe, H.J.; Pau, G.; Dijkgraaf, G.J.; Basset-Seguin, N.; Modrusan, Z.; Januario, T.; Tsui, V.; Durham, A.B.; Dlugosz, A.A.; Haverty, P.M.; et al. Genomic Analysis of Smoothened Inhibitor Resistance in Basal Cell Carcinoma. Cancer. Cell. 2015, 27, 327–341. [Google Scholar] [CrossRef] [Green Version]
  36. Nguyen, T.T.L.; Atwood, S.X. Illuminating Alternative Strategies to Treat Targeted Chemotherapy-Resistant Sporadic Basal Cell Carcinoma. J. Investig. Dermatol. 2018, 138, 1017–1019. [Google Scholar] [CrossRef]
  37. Wahid, M.; Jawed, A.; Mandal, R.K.; Dar, S.A.; Khan, S.; Akhter, N.; Haque, S. Vismodegib, Itraconazole and Sonidegib as Hedgehog Pathway Inhibitors and their Relative Competencies in the Treatment of Basal Cell Carcinomas. Crit. Rev. Oncol. Hematol. 2016, 98, 235–241. [Google Scholar] [CrossRef]
  38. Atwood, S.X.; Sarin, K.Y.; Whitson, R.J.; Li, J.R.; Kim, G.; Rezaee, M.; Ally, M.S.; Kim, J.; Yao, C.; Chang, A.L.; et al. Smoothened Variants Explain the Majority of Drug Resistance in Basal Cell Carcinoma. Cancer Cell. 2015, 27, 342–353. [Google Scholar] [CrossRef] [Green Version]
  39. Danial, C.; Sarin, K.Y.; Oro, A.E.; Chang, A.L. An Investigator-Initiated Open-Label Trial of Sonidegib in Advanced Basal Cell Carcinoma Patients Resistant to Vismodegib. Clin. Cancer Res. 2016, 22, 1325–1329. [Google Scholar] [CrossRef] [Green Version]
  40. Brancaccio, G.; Pea, F.; Moscarella, E.; Argenziano, G. Sonidegib for the Treatment of Advanced Basal Cell Carcinoma. Front. Oncol. 2020, 10, 582866. [Google Scholar] [CrossRef]
  41. Dummer, R.; Ascierto, P.A.; Basset-Seguin, N.; Dréno, B.; Garbe, C.; Gutzmer, R.; Hauschild, A.; Krattinger, R.; Lear, J.T.; Malvehy, J.; et al. Sonidegib and Vismodegib in the Treatment of Patients with Locally Advanced Basal Cell Carcinoma: A Joint Expert Opinion. J. Eur. Acad. Dermatol. Venereol. 2020, 34, 1944–1956. [Google Scholar] [CrossRef]
  42. Zollinger, M.; Lozac’h, F.; Hurh, E.; Emotte, C.; Bauly, H.; Swart, P. Absorption, Distribution, Metabolism, and Excretion (ADME) of ¹⁴C-Sonidegib (LDE225) in Healthy Volunteers. Cancer Chemother. Pharmacol. 2014, 74, 63–75. [Google Scholar] [CrossRef]
  43. Yoshida, G.J. Regulation of Heterogeneous Cancer-Associated Fibroblasts: The Molecular Pathology of Activated Signaling Pathways. J. Exp. Clin. Cancer Res. 2020, 39, 112. [Google Scholar] [CrossRef]
  44. Bellei, B.; Migliano, E.; Picardo, M. A Framework of Major Tumor-Promoting Signal Transduction Pathways Implicated in Melanoma-Fibroblast Dialogue. Cancers 2020, 12, 3400. [Google Scholar] [CrossRef]
  45. Kalluri, R. The Biology and Function of Fibroblasts in Cancer. Nat. Rev. Cancer. 2016, 16, 582–598. [Google Scholar] [CrossRef]
  46. 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]
  47. Takabatake, K.; Shimo, T.; Murakami, J.; Anqi, C.; Kawai, H.; Yoshida, S.; Wathone Oo, M.; Haruka, O.; Sukegawa, S.; Tsujigiwa, H.; et al. The Role of Sonic Hedgehog Signaling in the Tumor Microenvironment of Oral Squamous Cell Carcinoma. Int. J. Mol. Sci. 2019, 20, 5779. [Google Scholar] [CrossRef] [Green Version]
  48. Shan, T.; Chen, S.; Chen, X.; Lin, W.R.; Li, W.; Ma, J.; Wu, T.; Ji, H.; Li, Y.; Cui, X.; et al. Prometastatic Mechanisms of CAF-Mediated EMT Regulation in Pancreatic Cancer Cells. Int. J. Oncol. 2017, 50, 121–128. [Google Scholar] [CrossRef] [Green Version]
  49. Walter, K.; Omura, N.; Hong, S.M.; Griffith, M.; Vincent, A.; Borges, M.; Goggins, M. Overexpression of Smoothened Activates the Sonic Hedgehog Signaling Pathway in Pancreatic Cancer-Associated Fibroblasts. Clin. Cancer Res. 2010, 16, 1781–1789. [Google Scholar] [CrossRef] [Green Version]
  50. Valin, A.; Barnay-Verdier, S.; Robert, T.; Ripoche, H.; Brellier, F.; Chevallier-Lagente, O.; Avril, M.F.; Magnaldo, T. PTCH1 +/− Dermal Fibroblasts Isolated from Healthy Skin of Gorlin Syndrome Patients Exhibit Features of Carcinoma Associated Fibroblasts. PLoS ONE 2009, 4, e4818. [Google Scholar] [CrossRef] [Green Version]
  51. Gache, Y.; Brellier, F.; Rouanet, S.; Al-Qaraghuli, S.; Goncalves-Maia, M.; Burty-Valin, E.; Barnay, S.; Scarzello, S.; Ruat, M.; Sevenet, N.; et al. Basal Cell Carcinoma in Gorlin’s Patients: A Matter of Fibroblasts-Led Protumoral Microenvironment? PLoS ONE 2015, 10, e0145369. [Google Scholar] [CrossRef]
  52. Ponti, G.; Manfredini, M.; Pastorino, L.; Maccaferri, M.; Tomasi, A.; Pellacani, G. PTCH1 Germline Mutations and the Basaloid Follicular Hamartoma Values in the Tumor Spectrum of Basal Cell Carcinoma Syndrome (NBCCS). Anticancer Res. 2018, 38, 471–476. [Google Scholar]
  53. Klein, R.D.; Dykas, D.J.; Bale, A.E. Clinical Testing for the Nevoid Basal Cell Carcinoma Syndrome in a DNA Diagnostic Laboratory. Genet. Med. 2005, 7, 611–619. [Google Scholar] [CrossRef] [Green Version]
  54. Evans, T.; Boonchai, W.; Shanley, S.; Smyth, I.; Gillies, S.; Georgas, K.; Wainwright, B.; Chenevix-Trench, G.; Wicking, C. The Spectrum of Patched Mutations in a Collection of Australian Basal Cell Carcinomas. Hum. Mutat. 2000, 16, 43–48. [Google Scholar] [CrossRef]
  55. Lian, J.; Bahitham, W.; Panigrahi, R.; Nelson, R.; Li, L.; Watts, R.; Thiesen, A.; Lemieux, M.J.; Lehner, R. Genetic Variation in Human Carboxylesterase CES1 Confers Resistance to Hepatic Steatosis. Biochim. Biophys. Acta Mol. Cell. Biol. Lipids 2018, 1863, 688–699. [Google Scholar] [CrossRef]
  56. Tietze, J.K.; Pfob, M.; Eggert, M.; von Preußen, A.; Mehraein, Y.; Ruzicka, T.; Herzinger, T. A Non-Coding Mutation in the 5′ Untranslated Region of Patched Homologue 1 Predisposes to Basal Cell Carcinoma. Exp. Dermatol. 2013, 22, 834–835. [Google Scholar] [CrossRef] [Green Version]
  57. Ozretić, P.; Bisio, A.; Musani, V.; Trnski, D.; Sabol, M.; Levanat, S.; Inga, A. Regulation of Human PTCH1b Expression by Different 5′ Untranslated Region Cis-Regulatory Elements. RNA Biol. 2015, 12, 290–304. [Google Scholar] [CrossRef] [Green Version]
  58. Song, Y.; Tian, Y.; Zuo, Y.; Tu, J.C.; Feng, Y.F.; Qu, C.J. Altered Expression of PTCH and HHIP in Gastric Cancer through their Gene Promoter Methylation: Novel Targets for Gastric Cancer. Mol. Med. Rep. 2013, 7, 1159–1168. [Google Scholar] [CrossRef] [Green Version]
  59. Chuang, P.T.; McMahon, A.P. Vertebrate Hedgehog Signalling Modulated by Induction of a Hedgehog-Binding Protein. Nature 1999, 397, 617–621. [Google Scholar] [CrossRef]
  60. Li, P.; Zhang, X.; Murphy, A.J.; Costa, M.; Zhao, X.; Sun, H. Downregulation of Hedgehog-Interacting Protein (HHIP) Contributes to Hexavalent Chromium-Induced Malignant Transformation of Human Bronchial Epithelial Cells. Carcinogenesis 2021, 42, 136–147. [Google Scholar] [CrossRef]
  61. Wang, X.; Ma, W.; Yin, J.; Chen, M.; Jin, H. HHIP Gene Overexpression Inhibits the Growth, Migration and Invasion of Human Liver Cancer Cells. J. BUON 2020, 25, 2424–2429. [Google Scholar] [PubMed]
  62. Kayser, G.; Csanadi, A.; Kakanou, S.; Prasse, A.; Kassem, A.; Stickeler, E.; Passlick, B.; Zur Hausen, A. Downregulation of MTSS1 Expression is an Independent Prognosticator in Squamous Cell Carcinoma of the Lung. Br. J. Cancer 2015, 112, 866–873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Mertz, K.D.; Pathria, G.; Wagner, C.; Saarikangas, J.; Sboner, A.; Romanov, J.; Gschaider, M.; Lenz, F.; Neumann, F.; Schreiner, W.; et al. MTSS1 is a Metastasis Driver in a Subset of Human Melanomas. Nat. Commun. 2014, 5, 3465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Xie, F.; Ye, L.; Chen, J.; Wu, N.; Zhang, Z.; Yang, Y.; Zhang, L.; Jiang, W.G. The Impact of Metastasis Suppressor-1, MTSS1, on Oesophageal Squamous Cell Carcinoma and its Clinical Significance. J. Transl. Med. 2011, 9, 95. [Google Scholar] [CrossRef] [Green Version]
  65. Callahan, C.A.; Ofstad, T.; Horng, L.; Wang, J.K.; Zhen, H.H.; Coulombe, P.A.; Oro, A.E. MIM/BEG4, a Sonic Hedgehog-Responsive Gene that Potentiates Gli-Dependent Transcription. Genes Dev. 2004, 18, 2724–2729. [Google Scholar] [CrossRef] [Green Version]
  66. Markey, G.M. Carboxylesterase 1 (Ces1): From Monocyte Marker to Major Player. J. Clin. Pathol. 2011, 64, 107–109. [Google Scholar] [CrossRef]
  67. Lee, J.; Platt, K.A.; Censullo, P.; Ruiz i Altaba, A. Gli1 is a Target of Sonic Hedgehog that Induces Ventral Neural Tube Development. Development 1997, 124, 2537–2552. [Google Scholar] [CrossRef]
  68. Yang, K.; Wang, X.; Zhang, H.; Wang, Z.; Nan, G.; Li, Y.; Zhang, F.; Mohammed, M.K.; Haydon, R.C.; Luu, H.H.; et al. The Evolving Roles of Canonical WNT Signaling in Stem Cells and Tumorigenesis: Implications in Targeted Cancer Therapies. Lab. Invest. 2016, 96, 116–136. [Google Scholar] [CrossRef] [Green Version]
  69. Katoh, Y.; Katoh, M. Hedgehog Signaling Pathway and Gastrointestinal Stem Cell Signaling Network (Review). Int. J. Mol. Med. 2006, 18, 1019–1023. [Google Scholar] [CrossRef] [Green Version]
  70. Bovolenta, P.; Esteve, P.; Ruiz, J.M.; Cisneros, E.; Lopez-Rios, J. Beyond Wnt Inhibition: New Functions of Secreted Frizzled-Related Proteins in Development and Disease. J. Cell. Sci. 2008, 121, 737–746. [Google Scholar] [CrossRef] [Green Version]
  71. Scales, S.J.; de Sauvage, F.J. Mechanisms of Hedgehog Pathway Activation in Cancer and Implications for Therapy. Trends Pharmacol. Sci. 2009, 30, 303–312. [Google Scholar] [CrossRef]
  72. Pan, S.; Wu, X.; Jiang, J.; Gao, W.; Wan, Y.; Cheng, D.; Han, D.; Liu, J.; Englund, N.P.; Wang, Y.; et al. Discovery of NVP-LDE225, a Potent and Selective Smoothened Antagonist. ACS Med. Chem. Lett. 2010, 1, 130–134. [Google Scholar] [CrossRef] [Green Version]
  73. Bariwal, J.; Kumar, V.; Dong, Y.; Mahato, R.I. Design of Hedgehog Pathway Inhibitors for Cancer Treatment. Med. Res. Rev. 2019, 39, 1137–1204. [Google Scholar] [CrossRef]
  74. Lauressergues, E.; Heusler, P.; Lestienne, F.; Troulier, D.; Rauly-Lestienne, I.; Tourette, A.; Ailhaud, M.C.; Cathala, C.; Tardif, S.; Denais-Laliève, D.; et al. Pharmacological Evaluation of a Series of Smoothened Antagonists in Signaling Pathways and After Topical Application in a Depilated Mouse Model. Pharmacol. Res. Perspect. 2016, 4, e00214. [Google Scholar] [CrossRef]
  75. Stecca, B.; Mas, C.; Clement, V.; Zbinden, M.; Correa, R.; Piguet, V.; Beermann, F.; Ruiz i Altaba, A. Melanomas Require HEDGEHOG-GLI Signaling Regulated by Interactions between GLI1 and the RAS-MEK/AKT Pathways. Proc. Natl. Acad. Sci. USA 2007, 104, 5895–5900. [Google Scholar] [CrossRef] [Green Version]
  76. Bhattacharya, R.; Kwon, J.; Ali, B.; Wang, E.; Patra, S.; Shridhar, V.; Mukherjee, P. Role of Hedgehog Signaling in Ovarian Cancer. Clin. Cancer Res. 2008, 14, 7659–7666. [Google Scholar] [CrossRef] [Green Version]
  77. Rodon, J.; Tawbi, H.A.; Thomas, A.L.; Stoller, R.G.; Turtschi, C.P.; Baselga, J.; Sarantopoulos, J.; Mahalingam, D.; Shou, Y.; Moles, M.A.; et al. A Phase I, Multicenter, Open-Label, First-in-Human, Dose-Escalation Study of the Oral Smoothened Inhibitor Sonidegib (LDE225) in Patients with Advanced Solid Tumors. Clin. Cancer Res. 2014, 20, 1900–1909. [Google Scholar] [CrossRef] [Green Version]
  78. Karhadkar, S.S.; Bova, G.S.; Abdallah, N.; Dhara, S.; Gardner, D.; Maitra, A.; Isaacs, J.T.; Berman, D.M.; Beachy, P.A. Hedgehog Signalling in Prostate Regeneration, Neoplasia and Metastasis. Nature 2004, 431, 707–712. [Google Scholar] [CrossRef]
  79. Thayer, S.P.; di Magliano, M.P.; Heiser, P.W.; Nielsen, C.M.; Roberts, D.J.; Lauwers, G.Y.; Qi, Y.P.; Gysin, S.; Fernández-del Castillo, C.; Yajnik, V.; et al. Hedgehog is an Early and Late Mediator of Pancreatic Cancer Tumorigenesis. Nature 2003, 425, 851–856. [Google Scholar] [CrossRef] [Green Version]
  80. Lee, Y.; Kawagoe, R.; Sasai, K.; Li, Y.; Russell, H.R.; Curran, T.; McKinnon, P.J. Loss of Suppressor-of-Fused Function Promotes Tumorigenesis. Oncogene 2007, 26, 6442–6447. [Google Scholar] [CrossRef] [Green Version]
  81. Tsanev, R.; Tiigimägi, P.; Michelson, P.; Metsis, M.; Østerlund, T.; Kogerman, P. Identification of the Gene Transcription Repressor Domain of Gli3. FEBS Lett. 2009, 583, 224–228. [Google Scholar] [CrossRef] [Green Version]
  82. Sasaki, H.; Nishizaki, Y.; Hui, C.; Nakafuku, M.; Kondoh, H. Regulation of Gli2 and Gli3 Activities by an Amino-Terminal Repression Domain: Implication of Gli2 and Gli3 as Primary Mediators of Shh Signaling. Development 1999, 126, 3915–3924. [Google Scholar] [CrossRef]
  83. Smelkinson, M.G. The Hedgehog Signaling Pathway Emerges as a Pathogenic Target. J. Dev. Biol. 2017, 5, 14. [Google Scholar] [CrossRef] [Green Version]
  84. Cohen, M.; Kicheva, A.; Ribeiro, A.; Blassberg, R.; Page, K.M.; Barnes, C.P.; Briscoe, J. Ptch1 and Gli Regulate Shh Signalling Dynamics Via Multiple Mechanisms. Nat. Commun. 2015, 6, 6709. [Google Scholar] [CrossRef]
  85. Bosanac, I.; Maun, H.R.; Scales, S.J.; Wen, X.; Lingel, A.; Bazan, J.F.; de Sauvage, F.J.; Hymowitz, S.G.; Lazarus, R.A. The Structure of SHH in Complex with HHIP Reveals a Recognition Role for the Shh Pseudo Active Site in Signaling. Nat. Struct. Mol. Biol. 2009, 16, 691–697. [Google Scholar] [CrossRef]
  86. Holtz, A.M.; Peterson, K.A.; Nishi, Y.; Morin, S.; Song, J.Y.; Charron, F.; McMahon, A.P.; Allen, B.L. Essential Role for Ligand-Dependent Feedback Antagonism of Vertebrate Hedgehog Signaling by PTCH1, PTCH2 and HHIP1 during Neural Patterning. Development 2013, 140, 3423–3434. [Google Scholar] [CrossRef] [Green Version]
  87. Chuang, P.T.; Kawcak, T.; McMahon, A.P. Feedback Control of Mammalian Hedgehog Signaling by the Hedgehog-Binding Protein, Hip1, Modulates Fgf Signaling during Branching Morphogenesis of the Lung. Genes Dev. 2003, 17, 342–347. [Google Scholar] [CrossRef] [Green Version]
  88. Zamarrón, A.; García, M.; Río, M.D.; Larcher, F.; Juarranz, Á. Effects of Photodynamic Therapy on Dermal Fibroblasts from Xeroderma Pigmentosum and Gorlin-Goltz Syndrome Patients. Oncotarget 2017, 8, 77385–77399. [Google Scholar] [CrossRef] [Green Version]
  89. Zhu, Z.X.; Sun, C.C.; Ting Zhu, Y.; Wang, Y.; Wang, T.; Chi, L.S.; Cai, W.H.; Zheng, J.Y.; Zhou, X.; Cong, W.T.; et al. Hedgehog Signaling Contributes to Basic Fibroblast Growth Factor-Regulated Fibroblast Migration. Exp. Cell Res. 2017, 355, 83–94. [Google Scholar] [CrossRef]
  90. Jain, S.; Song, R.; Xie, J. Sonidegib: Mechanism of Action, Pharmacology, and Clinical Utility for Advanced Basal Cell Carcinomas. OncoTargets Ther. 2017, 10, 1645–1653. [Google Scholar] [CrossRef] [Green Version]
  91. Ueno, T.; Toi, M.; Saji, H.; Muta, M.; Bando, H.; Kuroi, K.; Koike, M.; Inadera, H.; Matsushima, K. Significance of Macrophage Chemoattractant Protein-1 in Macrophage Recruitment, Angiogenesis, and Survival in Human Breast Cancer. Clin. Cancer Res. 2000, 6, 3282–3289. [Google Scholar] [PubMed]
  92. Qian, B.Z.; Li, J.; Zhang, H.; Kitamura, T.; Zhang, J.; Campion, L.R.; Kaiser, E.A.; Snyder, L.A.; Pollard, J.W. CCL2 Recruits Inflammatory Monocytes to Facilitate Breast-Tumour Metastasis. Nature 2011, 475, 222–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Tian, H.; Callahan, C.A.; DuPree, K.J.; Darbonne, W.C.; Ahn, C.P.; Scales, S.J.; de Sauvage, F.J. Hedgehog Signaling is Restricted to the Stromal Compartment during Pancreatic Carcinogenesis. Proc. Natl. Acad. Sci. USA 2009, 106, 4254–4259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Chen, X.; Song, E. Turning Foes to Friends: Targeting Cancer-Associated Fibroblasts. Nat. Rev. Drug Discov. 2019, 18, 99–115. [Google Scholar] [CrossRef]
  95. Albrengues, J.; Bertero, T.; Grasset, E.; Bonan, S.; Maiel, M.; Bourget, I.; Philippe, C.; Herraiz Serrano, C.; Benamar, S.; Croce, O.; et al. Epigenetic Switch Drives the Conversion of Fibroblasts into Proinvasive Cancer-Associated Fibroblasts. Nat. Commun. 2015, 6, 10204. [Google Scholar] [CrossRef] [Green Version]
  96. Xu, J.; Xu, Y.; Xu, Y.; Yin, L.; Zhang, Y. Global Inactivation of Carboxylesterase 1 (Ces1/Ces1g) Protects Against Atherosclerosis in Ldlr (-/-) Mice. Sci. Rep. 2017, 7, 17845. [Google Scholar] [CrossRef] [Green Version]
  97. Luchetti, G.; Sircar, R.; Kong, J.H.; Nachtergaele, S.; Sagner, A.; Byrne, E.F.; Covey, D.F.; Siebold, C.; Rohatgi, R. Cholesterol Activates the G-Protein Coupled Receptor Smoothened to Promote Hedgehog Signaling. Elife 2016, 5, e20304. [Google Scholar] [CrossRef]
  98. Kinnebrew, M.; Iverson, E.J.; Patel, B.B.; Pusapati, G.V.; Kong, J.H.; Johnson, K.A.; Luchetti, G.; Eckert, K.M.; McDonald, J.G.; Covey, D.F.; et al. Cholesterol Accessibility at the Ciliary Membrane Controls Hedgehog Signaling. Elife 2019, 8, e50051. [Google Scholar] [CrossRef]
  99. Raleigh, D.R.; Sever, N.; Choksi, P.K.; Sigg, M.A.; Hines, K.M.; Thompson, B.M.; Elnatan, D.; Jaishankar, P.; Bisignano, P.; Garcia-Gonzalo, F.R.; et al. Cilia-Associated Oxysterols Activate Smoothened. Mol. Cell 2018, 72, 316–327.e5. [Google Scholar] [CrossRef] [Green Version]
  100. Blassberg, R.; Jacob, J. Lipid Metabolism Fattens Up Hedgehog Signaling. BMC Biol. 2017, 15, 95. [Google Scholar] [CrossRef] [Green Version]
  101. Radhakrishnan, A.; Rohatgi, R.; Siebold, C. Cholesterol Access in Cellular Membranes Controls Hedgehog Signaling. Nat. Chem. Biol. 2020, 16, 1303–1313. [Google Scholar] [CrossRef]
  102. Taylor, F.R.; Wen, D.; Garber, E.A.; Carmillo, A.N.; Baker, D.P.; Arduini, R.M.; Williams, K.P.; Weinreb, P.H.; Rayhorn, P.; Hronowski, X.; et al. Enhanced Potency of Human Sonic Hedgehog by Hydrophobic Modification. Biochemistry 2001, 40, 4359–4371. [Google Scholar] [CrossRef]
  103. Pepinsky, R.B.; Zeng, C.; Wen, D.; Rayhorn, P.; Baker, D.P.; Williams, K.P.; Bixler, S.A.; Ambrose, C.M.; Garber, E.A.; Miatkowski, K.; et al. Identification of a Palmitic Acid-Modified Form of Human Sonic Hedgehog. J. Biol. Chem. 1998, 273, 14037–14045. [Google Scholar] [CrossRef] [Green Version]
  104. Sekulic, A.; Migden, M.R.; Basset-Seguin, N.; Garbe, C.; Gesierich, A.; Lao, C.D.; Miller, C.; Mortier, L.; Murrell, D.F.; Hamid, O.; et al. Long-Term Safety and Efficacy of Vismodegib in Patients with Advanced Basal Cell Carcinoma: Final Update of the Pivotal ERIVANCE BCC Study. BMC Cancer 2017, 17, 332. [Google Scholar] [CrossRef]
  105. Lear, J.T.; Migden, M.R.; Lewis, K.D.; Chang, A.L.S.; Guminski, A.; Gutzmer, R.; Dirix, L.; Combemale, P.; Stratigos, A.; Plummer, R.; et al. Long-Term Efficacy and Safety of Sonidegib in Patients with Locally Advanced and Metastatic Basal Cell Carcinoma: 30-Month Analysis of the Randomized Phase 2 BOLT Study. J. Eur. Acad. Dermatol. Venereol. 2018, 32, 372–381. [Google Scholar] [CrossRef]
  106. Migden, M.; Farberg, A.S.; Dummer, R.; Squittieri, N.; Hanke, C.W. A Review of Hedgehog Inhibitors Sonidegib and Vismodegib for Treatment of Advanced Basal Cell Carcinoma. J. Drugs Dermatol. 2021, 20, 156–165. [Google Scholar] [CrossRef]
  107. Dummer, R.; Guminksi, A.; Gutzmer, R.; Lear, J.T.; Lewis, K.D.; Chang, A.L.S.; Combemale, P.; Dirix, L.; Kaatz, M.; Kudchadkar, R.; et al. Long-Term Efficacy and Safety of Sonidegib in Patients with Advanced Basal Cell Carcinoma: 42-Month Analysis of the Phase II Randomized, Double-Blind BOLT Study. Br. J. Dermatol. 2020, 182, 1369–1378. [Google Scholar] [CrossRef]
  108. Passarelli, A.; Galdo, G.; Aieta, M.; Fabrizio, T.; Villonio, A.; Conca, R. A Vismodegib Experience in Elderly Patients with Basal Cell Carcinoma: Case Reports and Review of the Literature. Int. J. Mol. Sci. 2020, 21, 8596. [Google Scholar] [CrossRef]
  109. Lorusso, P.M.; Jimeno, A.; Dy, G.; Adjei, A.; Berlin, J.; Leichman, L.; Low, J.A.; Colburn, D.; Chang, I.; Cheeti, S.; et al. Pharmacokinetic Dose-Scheduling Study of Hedgehog Pathway Inhibitor Vismodegib (GDC-0449) in Patients with Locally Advanced Or Metastatic Solid Tumors. Clin. Cancer Res. 2011, 17, 5774–5782. [Google Scholar] [CrossRef] [Green Version]
  110. Graham, R.A.; Hop, C.E.; Borin, M.T.; Lum, B.L.; Colburn, D.; Chang, I.; Shin, Y.G.; Malhi, V.; Low, J.A.; Dresser, M.J. Single and Multiple Dose Intravenous and Oral Pharmacokinetics of the Hedgehog Pathway Inhibitor Vismodegib in Healthy Female Subjects. Br. J. Clin. Pharmacol. 2012, 74, 788–796. [Google Scholar] [CrossRef] [Green Version]
  111. Pietrobono, S.; Stecca, B. Targeting the Oncoprotein Smoothened by Small Molecules: Focus on Novel Acylguanidine Derivatives as Potent Smoothened Inhibitors. Cells 2018, 7, 272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Rohatgi, R.; Milenkovic, L.; Corcoran, R.B.; Scott, M.P. Hedgehog Signal Transduction by Smoothened: Pharmacologic Evidence for a 2-Step Activation Process. Proc. Natl. Acad. Sci. USA 2009, 106, 3196–3201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 13 05858 g001 550
Figure 1. (a) MTT assay of NHFs (n = 11), NBCCS-HFs (n = 9) and CAFs (n = 6) after 72 h treatment with vismodegib and sonidegib (range 5–100 μM). Data are reported as mean ± SD. Experiments were performed in duplicate. (b) AnnexinV/iodide propidium staining evidenced apoptosis cell death only with sonidegib treatment. Dot plots show one representative experiment performed 48 h after Hh inhibitors treatment. * p ≤ 0.05 and ** p ≤ 0.01.
Figure 1. (a) MTT assay of NHFs (n = 11), NBCCS-HFs (n = 9) and CAFs (n = 6) after 72 h treatment with vismodegib and sonidegib (range 5–100 μM). Data are reported as mean ± SD. Experiments were performed in duplicate. (b) AnnexinV/iodide propidium staining evidenced apoptosis cell death only with sonidegib treatment. Dot plots show one representative experiment performed 48 h after Hh inhibitors treatment. * p ≤ 0.05 and ** p ≤ 0.01.
Cancers 13 05858 g001
Cancers 13 05858 g002 550
Figure 2. (a) Volcano plot represents basal gene expression analysis comparing NBCCS-HFs (n = 9) mRNA to NHFs (n = 9) mRNA level of expression. (b) Volcano plot represents basal gene expression analysis comparing CAFs (n = 6) mRNA to NHFs (n = 9) mRNA level of expression. Analysis was performed combining three different housekeeping genes (β-actin, GAPDH and 18 s). One-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Figure 2. (a) Volcano plot represents basal gene expression analysis comparing NBCCS-HFs (n = 9) mRNA to NHFs (n = 9) mRNA level of expression. (b) Volcano plot represents basal gene expression analysis comparing CAFs (n = 6) mRNA to NHFs (n = 9) mRNA level of expression. Analysis was performed combining three different housekeeping genes (β-actin, GAPDH and 18 s). One-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Cancers 13 05858 g002
Cancers 13 05858 g003 550
Figure 3. (a) One representative protein array assay comparing inflammatory factors released by NBCCS-HFs and NHFs. (b) Histograms report densitometric analysis of four independent experiments mean ± SD. (c) One representative protein array assay comparing growth factors released by NBCCS-HFs and NHFs. (d) Histograms report densitometric analysis of three independent experiments mean ± SD. * p ≤ 0.05.
Figure 3. (a) One representative protein array assay comparing inflammatory factors released by NBCCS-HFs and NHFs. (b) Histograms report densitometric analysis of four independent experiments mean ± SD. (c) One representative protein array assay comparing growth factors released by NBCCS-HFs and NHFs. (d) Histograms report densitometric analysis of three independent experiments mean ± SD. * p ≤ 0.05.
Cancers 13 05858 g003
Cancers 13 05858 g004 550
Figure 4. Volcano plot represents gene expression analysis comparing NBCCS-HFs treated with vismodegib (a) and sonidegib (b) for 72 h (n = 9) to NBCCS-HFs untreated (n = 9). (b) Analysis was performed using three different housekeeping genes (β-actin, GAPDH and 18S). A one-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Figure 4. Volcano plot represents gene expression analysis comparing NBCCS-HFs treated with vismodegib (a) and sonidegib (b) for 72 h (n = 9) to NBCCS-HFs untreated (n = 9). (b) Analysis was performed using three different housekeeping genes (β-actin, GAPDH and 18S). A one-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Cancers 13 05858 g004
Cancers 13 05858 g005 550
Figure 5. Volcano plot represents gene expression analysis comparing CAFs treated with vismodegib (a) and sonidegib (b) for 72 h (n = 6) to CAFs untreated (n = 6). (b) Analysis was performed using three different housekeeping genes (β-actin, GAPDH and 18S). One-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Figure 5. Volcano plot represents gene expression analysis comparing CAFs treated with vismodegib (a) and sonidegib (b) for 72 h (n = 6) to CAFs untreated (n = 6). (b) Analysis was performed using three different housekeeping genes (β-actin, GAPDH and 18S). One-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Cancers 13 05858 g005
Cancers 13 05858 g006 550
Figure 6. Volcano plot represents gene expression analysis comparing NHFs treated with vismodegib (a) and sonidegib (b) for 72 h (n = 9) to NHFs untreated (n = 9). (b) Analysis was performed using three different housekeeping genes (β-actin, GAPDH and 18S). A one-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Figure 6. Volcano plot represents gene expression analysis comparing NHFs treated with vismodegib (a) and sonidegib (b) for 72 h (n = 9) to NHFs untreated (n = 9). (b) Analysis was performed using three different housekeeping genes (β-actin, GAPDH and 18S). A one-way ANOVA statistical test with thresholds >2.0 fold-change and p < 0.05 defined significant increase are reported in red; significant decrease <0.5 fold-change and p < 0.05 reported in green; any fold-difference with p ≥ 0.5, e.g., insignificant are reported in grey; ≤2.0 or ≥0.5 difference e.g., flat reported in black.
Cancers 13 05858 g006
Cancers 13 05858 g007 550
Figure 7. (a) Immunofluorescence analysis of α-SMA and FAP protein expression in NBCCS and control cells. Original magnification 20×. (b) Semi-quantitative RT-PCR was performed to assess the level of PDGFRα, β and FSP1 mRNA. Histograms report mean ± SD from 9 NBCCS and 9 NHF cell lines. Value from healthy fibroblasts was arbitrarily indicated as 1.0. (c) One representative Western blot detecting β-catenin, pSer9-GSK3β and pSer473-AKT of NHFs and NBCCS-HFs. Histograms show densitometric analysis mean ± SD of nine independent experiments. * p < 0.05; ** p < 0.01. (d) Semi-quantitative RT-PCR was performed to assess the mRNA level of expression Wnt-β-catenin pathway transcription factors LEF1, TCF1 and TCF4. Histograms report mean ± SD from 9 NBCCS and 9 NHF cell lines. Value from healthy fibroblasts was arbitrarily indicated as 1.0. (e) One representative Western blot detecting β-catenin, pSer9-GSK3β and pSer473-AKT of CAFs. Histograms show densitometric analysis mean ± SD of five independent experiments. * p < 0.05; ** p < 0.01.
Figure 7. (a) Immunofluorescence analysis of α-SMA and FAP protein expression in NBCCS and control cells. Original magnification 20×. (b) Semi-quantitative RT-PCR was performed to assess the level of PDGFRα, β and FSP1 mRNA. Histograms report mean ± SD from 9 NBCCS and 9 NHF cell lines. Value from healthy fibroblasts was arbitrarily indicated as 1.0. (c) One representative Western blot detecting β-catenin, pSer9-GSK3β and pSer473-AKT of NHFs and NBCCS-HFs. Histograms show densitometric analysis mean ± SD of nine independent experiments. * p < 0.05; ** p < 0.01. (d) Semi-quantitative RT-PCR was performed to assess the mRNA level of expression Wnt-β-catenin pathway transcription factors LEF1, TCF1 and TCF4. Histograms report mean ± SD from 9 NBCCS and 9 NHF cell lines. Value from healthy fibroblasts was arbitrarily indicated as 1.0. (e) One representative Western blot detecting β-catenin, pSer9-GSK3β and pSer473-AKT of CAFs. Histograms show densitometric analysis mean ± SD of five independent experiments. * p < 0.05; ** p < 0.01.
Cancers 13 05858 g007
Cancers 13 05858 g008 550
Figure 8. Quantification of vismodegib and sonidegib uptake. Intracellular concentration of vismodegib and sonidegib was evaluated after 3, 6 and 24 h of drugs exposure. Histograms show mean ± SD of three independent experiments.
Figure 8. Quantification of vismodegib and sonidegib uptake. Intracellular concentration of vismodegib and sonidegib was evaluated after 3, 6 and 24 h of drugs exposure. Histograms show mean ± SD of three independent experiments.
Cancers 13 05858 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eibenschutz, L.; Caputo, S.; Camera, E.; Carbone, A.; Silipo, V.; Migliano, E.; Aurizi, C.; Cota, C.; Frascione, P.; Bellei, B. Evaluation of Hedgehog Pathway Inhibition on Nevoid Basal Cell Carcinoma Syndrome Fibroblasts and Basal Cell Carcinoma-Associated Fibroblasts: Are Vismodegib and Sonidegib Useful to Target Cancer-Prone Fibroblasts? Cancers 2021, 13, 5858. https://doi.org/10.3390/cancers13225858

AMA Style

Eibenschutz L, Caputo S, Camera E, Carbone A, Silipo V, Migliano E, Aurizi C, Cota C, Frascione P, Bellei B. Evaluation of Hedgehog Pathway Inhibition on Nevoid Basal Cell Carcinoma Syndrome Fibroblasts and Basal Cell Carcinoma-Associated Fibroblasts: Are Vismodegib and Sonidegib Useful to Target Cancer-Prone Fibroblasts? Cancers. 2021; 13(22):5858. https://doi.org/10.3390/cancers13225858

Chicago/Turabian Style

Eibenschutz, Laura, Silvia Caputo, Emanuela Camera, Anna Carbone, Vitaliano Silipo, Emilia Migliano, Caterina Aurizi, Carlo Cota, Pasquale Frascione, and Barbara Bellei. 2021. "Evaluation of Hedgehog Pathway Inhibition on Nevoid Basal Cell Carcinoma Syndrome Fibroblasts and Basal Cell Carcinoma-Associated Fibroblasts: Are Vismodegib and Sonidegib Useful to Target Cancer-Prone Fibroblasts?" Cancers 13, no. 22: 5858. https://doi.org/10.3390/cancers13225858

APA Style

Eibenschutz, L., Caputo, S., Camera, E., Carbone, A., Silipo, V., Migliano, E., Aurizi, C., Cota, C., Frascione, P., & Bellei, B. (2021). Evaluation of Hedgehog Pathway Inhibition on Nevoid Basal Cell Carcinoma Syndrome Fibroblasts and Basal Cell Carcinoma-Associated Fibroblasts: Are Vismodegib and Sonidegib Useful to Target Cancer-Prone Fibroblasts? Cancers, 13(22), 5858. https://doi.org/10.3390/cancers13225858

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