Abstract
Lei Wang1,2,*, Joy Q. Jin2,*, Yong Zhou2,3, Ziqiang Tian1, David M. Jablons2 and Biao He2
1Department of Thoracic Surgery, Fourth Hospital of Hebei Medical University, Shijiazhuang, Hebei 050011, China
2Thoracic Oncology Program, Department of Surgery, Helen Diller Family Comprehensive Cancer Center, University of California, San Francisco, CA 94115, USA
3Department of Respiratory Diseases, Sir Run Run Shaw Hospital, Zhejiang University, Hangzhou, Zhejiang 310016, China
*These authors have contributed equally to this work
Correspondence to:
Biao He, email: [email protected]
David M. Jablons, email: [email protected]
Keywords: sonic hedgehog; Gli; epithelial-mesenchymal transition; esophageal adenocarcinoma
Received: September 19, 2017 Accepted: November 09, 2017 Published: December 01, 2017
ABSTRACT
Esophageal adenocarcinoma (EAC) accounts for the most esophageal cancer cases in the US, and is notoriously aggressive. This study examines the role of Sonic Hedgehog (SHh)/Gli signaling in the regulation of epithelial-mesenchymal transition (EMT), a process tied to invasion and metastasis, in EAC.
Gli/EMT protein expression levels were examined by western blot in paired EAC patient tissues (n = 24) and cell lines (OE19, OE33). Functional analyses were performed (siRNA, treatment with Gli-inhibitor, AKT-inhibitor, and N-Shh recombinant proteins) to investigate SHh/Gli signaling and EMT, cell cycle, and prognostic markers in EAC cell lines. MTS, luciferase reporter, qRT-PCR, western blot, wound healing, and transwell assays were executed to analyze pathway activity, cell migration, and invasion.
Aberrant Gli1/2 expression was found in EAC patient tissues, and was significantly associated with increased EMT and AKT pathway activity. Stimulation of SHh/Gli resulted in EMT signaling, including expression of E-cadherin, N-cadherin, Vimentin, β-catenin, Snail, and Slug, as well as cell cycle progression at mRNA and protein levels in EAC cell lines. Gli inhibition via small molecule administration and siRNA significantly reduced EMT, decreasing cell mobility and invasion. Both Gli and AKT inhibition rescued E-cadherin expression and suppressed AKT phosphorylation.
This study provides evidence for a strong association between aberrant Gli1/2 expression and AKT/EMT markers in EAC; activated SHh/Gli signaling may be a critical component in promoting cell survival, metastases, and resistance to chemotherapy, and represents a promising avenue to target tumor proliferation and mobility.
INTRODUCTION
Esophageal carcinoma is one of the leading causes of cancer-related deaths worldwide, with low 5-year survival rates at approximately 12-20% [1, 2]. Esophageal adenocarcinoma (EAC) is one of the most common subtypes of esophageal cancer [1, 3]. Patients suffer from poor prognoses and high incidence of recurrence and metastasis, even following positive response to traditional therapeutic strategies including surgery, chemotherapy, and radiotherapy [4, 5]. Thus, more meaningful treatment strategies are desperately needed to target the mechanisms leading to tumor invasion and metastasis in EAC.
The epithelial-mesenchymal transition (EMT) is a process in which cells lose their cell-cell adhesive properties and gain migratory and invasive potential; it is essential for events in embryonic development, wound healing, fibrosis, cancer progression, and metastasis [6–8]. EMT involves multiple complex changes in the distribution and function of proteins, including E-cadherin, an adhesive protein inactivated in numerous cancers [8]. This process is regulated by a number of converging signaling cascades, including the Sonic Hedgehog (SHh) pathway, and confers critical traits required for seeding metastasis and developing stem cell properties that allow new cancer cell colonies to be launched [9]. Acquisition of EMT features has been associated with poor prognosis and chemotherapeutic resistance—accordingly, further knowledge can improve our understanding of tumor recurrence and metastasis to identify potential therapeutic targets in EAC [10–13].
Activation of SHh signaling has been implicated in the tumorigenesis and metastasis of various cancers [14–26]. The canonical cascade is initiated by Patched (Ptch) and Smoothened (Smo); Sonic Hedgehog ligand (Shh) binding to Ptch allows release of Smo, causing active full-length Gli to enter the nucleus and activate transcription of target genes, including Gli1, Ptch1, and Cyclin D1, in a context- and cell-type dependent manner [27]. Among the Gli family of transcription factors (TFs), Gli1 and Gli2 are considered activators, while Gli3 serves as a repressor [27]. Non-canonical Gli activation independent of Shh activation has also been noted in many cancer cells types, owing to stimulation by other oncogenic signaling pathways such as transforming growth factor β (TGF-β), epidermal growth factor receptor (EGFR), RAS, and AKT/PI3K [28–32]. As Gli TFs constitute the final effectors of the SHh pathway, and are implicated in multiple other oncogenic signaling cascades, they represent an important downstream target for potential cancer therapeutics. Efforts have been made to develop inhibitors of the SHh pathway, including Vismodegib, a Smo inhibitor, as well as GANT-61, a Gli inhibitor that regulates Gli-dependent transcription, to promote anti-cancer activity [33, 34]. In our lab, we have also developed a novel Gli inhibitor (Gli-i) that specifically blocks Gli1 and Gli2 transcriptional activity with significant efficacy [35].
However, the relationship between EMT and SHh/Gli activation has not previously been studied in EAC, and existing data from other solid tumor types are controversial. While some studies have shown that SHh/Gli inhibition block EMT, the exact mechanisms have not been elucidated. In melanoma and pancreatic cancers, results suggest the role of Gli in facilitating cancer migration and invasion by regulation of E-cadherin [36, 37]. On the other hand, another conflicting report proposes the inhibition of Gli in promoting the same EMT characteristics in pancreatic cancer [38]. In lung squamous cell cancer (SCC), Gli expression is inversely correlated to that of E-cadherin [39]. Studies in melanoma and hepatocellular carcinoma have linked Gli1 to vascular/capsular invasion, advanced tumor stage, and upregulation of matrix metalloprotease (MMP)-2 and MMP-9, while siRNA silencing of Gli1 successfully reduced invasion and increased E-cadherin expression [40, 41]. Our lab recently studied upregulated signaling in lung cancer, investigating Gli1’s inverse correlation with E-cadherin; inhibition of the SHh pathway upregulates E-cadherin expression and suppresses lung cancer cell migration [39].
In this study, aberrant Gli activation was studied in both EAC tissue samples and cell lines. Gli and EMT markers were found to be inversely correlated, and inhibition of the former minimized migration and invasion of EAC cells. Gli suppression induced upregulated E-cadherin expression and downregulated phosphorylated AKT, suggesting Gli may be critical for the metastasis and recurrence of esophageal adenocarcinomas.
RESULTS
Gli is upregulated in EAC tissue samples, and correlates with EMT and AKT pathway markers
Previous studies have suggested presence of increased Gli signaling components in EAC, but have not investigated the relationship between aberrant Gli activation and EMT. Protein expression of Gli1, Gli2, and key EMT pathway markers was first examined in 24 matched EAC and normal patient tissue samples from the Thoracic Oncology Program at University of California, San Francisco (UCSF). In western blot analysis, 91.7% (22/24) and 87.5% (21/24) of tumor samples showed higher Gli1 and Gli2 expression than in paired normal tissue, respectively (Figure 1A, Supplementary Figure 1A). Overexpression of N-cadherin, a biomarker associated with increased EMT [42], was observed in 83.3% (20/24) in tumor samples, while a loss of E-cadherin was seen in 95.8% (23/24) of tumor tissues. Activation of β-catenin and vimentin, both correlated with EMT invasive and proliferative characteristics [43, 44], was observed in 83.3% (20/24) of the patient sample cohort. AKT pathway proteins were subsequently analyzed via western blot (Figure 1B, Supplementary Figure 1B); activation of m-TOR, p-AKT, and p-S6K1 was observed in 87.5% (21/24), 75% (18/24), and 83.3% (20/24) of the tumor tissue samples, respectively. M-TOR activation has previously been shown to induce overexpression of downstream effectors, including p-S6K1 (activated S6K1); both factors have been linked to crosstalk with Gli signaling, and are associated with enhanced aggression [45, 46]. On the other hand, p-AKT plays an important role in tumorigenesis and resistance to therapy in EAC [47]. Thus, our results demonstrate Gli1 and Gli2 hyperactivation in EAC patient tissue samples, with strong correlations to 7 pathway markers implicated in upregulated EMT signaling.
Figure 1: Gli1/2, EMT, and AKT pathway markers are upregulated in EAC tissue samples. (A) Western blots of Gli1 and Gli2 protein expression in 24 matched pairs of esophageal tumor (T) and normal (N) tissues. GAPDH served as a loading control; results of 6 numbered, representative sample pairs are shown here. EMT markers (E-cadherin, N-cadherin, Vimentin, and β-catenin) were also examined. (B) Western blots of AKT pathway protein (p-MTOR, p-AKT, and p-S6K1) expression in 24 matched pairs of esophageal T and N tissues. GAPDH served as a loading control; results of 6 representative sample pairs are shown.
Gli inhibition and siRNA knockdown reduce EMT and cell cycle activity in EAC cell lines
To further understand the role of Gli in affecting EMT and cell cycle activity, two EAC cell lines, OE19 and OE33, were subjected to cell viability, reporter, and qRT-PCR analyses. Inhibition of Gli activity was achieved by treating cells with a Gli inhibitor (Gli-i) developed by our lab, which has been shown to specifically target Gli1 and Gli2 transcriptional activity [35]. Stimulation of the SHh/Gli pathway was achieved by applying recombinant Shh proteins (N-Shh) to the cultured cells.
MTS assays using a serial dilution of Gli-i (with DMSO as vehicle control) in both OE19 and OE33 EAC cell lines yielded IC50 values of 363.2 nM and 464.3 nM, respectively (Figure 2A). Results suggest Gli inhibition may drastically reduce EAC cell viability; luciferase reporter assays were then performed to measure Gli-mediated transcriptional activity in both cell lines (Figure 2B). As expected, N-Shh stimulation of cells significantly increased reporter activity in OE19 and OE33 (p < 0.05), while Gli-i treatment of transfected cells showed over 30% less relative reporter activity in both cell lines. When subjected to both Gli-i and N-Shh treatment, OE19 and OE33 continued to demonstrate decreased Gli transcriptional activity. qRT-PCR analyses of EMT markers in cells treated in the same conditions (DMSO, N-Shh, Gli-i, or Gli-i + N-Shh) produced consistent results (Figure 2C). While N-Shh stimulation increased relative mRNA expression of CTNNB1 (gene coding for β-catenin), CDH2 (gene coding for N-cadherin), Vimentin, Snail, Slug, and Zeb1, administration of Gli-i in OE19 and OE33 cells significantly reduced mRNA expression levels of all markers tested (p < 0.05). All six genes (CTNNB1, CDH2, Vimentin, Snail, Slug, and Zeb1) are positively associated with EMT activity; Snail, Slug, and Zeb1 function to repress E-cadherin, which typically adheres adjacent epithelial cells in carcinomas, leading to increased metastasis and EMT [48, 49].
Figure 2: Gli inhibition reduces EMT activity in EAC cell lines. (A) MTS cell viability assay in (i) OE19, and (ii) OE33 EAC cell lines. Cells were subjected to a serial dilution of Gli-i with DMSO control over a 3-day period, yielding IC50 values of 363.2 nM and 464.3 nM, respectively. (B) Luciferase reporter assay in OE19 and OE33, treated with DMSO (control), N-Shh (0.5 mg/mL), Gli-i (500 nmol/L), or Gli-i (500 nmol/L) and N-Shh (0.5 mg/mL) stimulation, for 24 hours. Results are expressed as fold activity, i.e. the ratio of luciferase activity induced in Gli-transfected cells relative to basal luciferase activity in control transfected OE19/33 cells (p-values of < 0.05, < 0.01, or < 0.001 were indicated as *, **, or ***, respectively). (C) qRT-PCR expression levels of EMT markers in (i) OE19 (CTNNB1, CDH2, Vimentin), and (ii) OE33 (Vimentin, CDH2, CTNNB1, Snail, Slug, Zeb1). Cells were treated with DMSO, N-Shh, Gli-i, or Gli-i + N-Shh prior to RNA extraction according to the conditions described above (p-values of < 0.05, < 0.01, or < 0.001 were indicated as *, **, or ***, respectively).
To determine protein expression patterns of EMT and related markers following Gli-i treatment, western blots were then conducted in OE19 and OE33 (Figure 3). Gli inhibition is associated with decreased p-AKT expression (normalized against total AKT protein levels), as well as positive correlation with p21, and negative correlation with Snail and Slug (Figure 3A). The relationship of Gli-i treatment with increased p21 is especially promising given that p21, a Cyclin-dependent kinase (Cdk) inhibitor, promotes cell cycle arrest and negatively regulates cell cycle activities [50]. Additionally, Cyclin D1 showed a decrease in expression level following Gli-i administration, suggesting a decrease in TGF-β mediated cell migration (Figure 3B) [51]. Further examination of p-ERK revealed decreased expression (normalized against total ERK protein levels) following Gli-i incubation, even with N-Shh stimulation (Figure 3B). Activation of ERK is required for TGF-β1 induced EMT [52]. EMT markers, including E-cadherin, N-cadherin, and β-catenin, showed results consistent with Figure 1; Gli-i induced a decrease in E-cadherin and increase in expression of the latter two markers.
Figure 3: Gli inhibition and siRNA knockdown reduce EMT and cell cycle activity in EAC cell lines. (A) Western blots of OE19 and OE33 cells, pretreated with either DMSO (control) or Gli-i (500 nmol/L). Expression of AKT, p-AKT, Snail, Slug, and p21 was examined, with GAPDH serving as a loading control. (B) Western blots of EMT and cell cycle signaling markers in (i) OE19 (Gli1, p-S6K1, ERK, p-ERK, E-cadherin, N-cadherin, Cyclin D1), and (ii) OE33 (Gli1, Gli2, m-TOR, E-cadherin, N-cadherin, β-catenin, Cyclin D1), with GAPDH as a control. Cells were pretreated with DMSO, N-Shh (0.5 mg/mL), Gli-i (500 nmol/L), or Gli-i (500 nmol/L) with N-Shh (0.5 mg/mL) stimulation prior to protein extraction. (C) Western blots of OE19 and OE33 cells, subjected to siRNA treatment. Cells were given control siRNA (100 nM) or Gli1 + Gli2 siRNA (100 nM each) for 48 hours prior to total protein extraction. Expression of Gli1, Gli2, E-cadherin, N-cadherin, and β-catenin was examined, with GAPDH serving as an internal control.
In addition, functional analyses with siRNA knockdown of both Gli1 and Gli2 yielded western blot data consistent with earlier parts of the figure, showing decreased Gli1, Gli2, and EMT activity (as measured by E-cadherin, N-cadherin, and β-catenin) in both OE19 and OE33 cell lines (Figure 3C). Taken together, these results strongly suggest that inhibition of Gli signaling may suppress EMT and cell cycle progression in EAC.
Gli signaling promotes cell migration and invasion in EAC
To assess the effects of Gli signaling on cell migration and invasion in EAC, wound healing and transwell assays were conducted in OE19 and OE33 cell lines, subjected to treatment with DMSO (vehicle control), N-Shh, Gli-i, and Gli-i + N-Shh in the same manner as above. Over the course of 24 hours, OE33 cells exhibited significantly slower migration when treated with Gli-i (Figure 4A). Cells stimulated by N-Shh demonstrated greater mobility during the same period, although such cells still migrated slower than the DMSO control when Gli-i was administered (Figure 4B). Furthermore, in transwell invasion assays, both OE19 and OE33 demonstrated decreased invasion by approximately 50% following Gli-i treatment (Figure 5A, 5B). Consistent with data in earlier figures, these results indicate that Gli signaling may promote both cell migration and invasion through loss of E-cadherin in EAC.
Figure 4: Gli signaling promotes cell migration in EAC. (A) Wound healing assays of OE33 EAC cells, treated with DMSO (control), N-Shh (0.5 mg/mL), Gli-i (500 nmol/L), or Gli-i (500 nmol/L) with N-Shh (0.5 mg/mL) stimulation. Representative pictures taken at 0 and 24 hours under a light microscope are shown (100x). (B) Quantification of wound healing assays, expressed as fold change in comparison to the control group (p-values of < 0.05, < 0.01, or < 0.001 were indicated as *, **, or ***, respectively).
Figure 5: Gli signaling promotes cell invasion in EAC. (A) Transwell invasion assay of EAC OE19 and OE33 cells, treated with DMSO (control), N-Shh (0.5 mg/mL), Gli-i (500 nmol/L), or Gli-i (500 nmol/L) with N-Shh (0.5 mg/mL) stimulation. Representative images captured by a light microscope are shown (100x). (B) Quantification of transwell invasion assays, expressed as fold change in comparison to the control group (p-values of < 0.05, < 0.01, or < 0.001 were indicated as *, **, or ***, respectively).
AKT inhibition reduces EMT and cell cycle activity in EAC cell lines
Further examination on AKT signaling was conducted in EAC cell lines, following initial investigation (Figures 1-3). OE19 and OE33 EAC cells were administered MK2206, an AKT inhibitor (AKT-i), and analyzed by western blot to determine protein levels of relevant pathway and EMT markers after 30 hours of incubation. In OE19 (Figure 6A, Supplementary Figure 2A), AKT-i treated cells exhibited decreased m-TOR, p-AKT, and p-S6K1, as well as inhibited EMT activity (increased E-cadherin, decreased N-cadherin and Vimentin); Cyclin D1 levels were lower in both AKT-i and AKT-i + N-Shh stimulated cells. Consistent with earlier results, AKT inhibition yielded decreased m-TOR, N-cadherin, and Vimentin, alongside increased E-cadherin levels in OE33 (Figure 6B, Supplementary Figure 2B). N-Shh has previously been shown to promote EMT, cell migration, and invasion [53]. As with Gli-i, N-Shh used in combination with AKT-i rescues the effects induced by the latter, suggesting that SHh/Gli signaling may regulate EMT via the AKT pathway in EAC. Thus, inhibition of AKT and Gli signaling may serve as a promising avenue for decreased EMT and cell cycle activity, both of which are linked to greater aggression and viability in EAC.
Figure 6: AKT inhibition reduces EMT and cell cycle activity in EAC cell lines. Western blots of EMT and cell cycle signaling markers in (A) OE19 (m-TOR, p-AKT, p-S6K1, E-cadherin, N-cadherin, Vimentin, Cyclin D1), and (B) OE33 (m-TOR, E-cadherin, N-cadherin, Vimentin), with GAPDH as a loading control. Cells were pretreated with DMSO, N-Shh (0.5 mg/mL), AKT-i (1 μmol/L), or AKT-i (1 μmol/L) with N-Shh (0.5 mg/mL) stimulation prior to protein extraction.
DISCUSSION
SHh/Gli signaling plays an important role in solid malignancies; downstream pathway effectors are tied to EMT, and many serve as potential therapeutic targets in the prevention of tumor aggression and metastasis. This study provides strong preliminary evidence that in vitro Gli inhibition impairs the migratory, invasive, and proliferative properties of EAC cells by reducing EMT. EAC tissue specimens from patients demonstrated inverse correlations between Gli1/2 and key EMT markers, including cell adhesion protein E-cadherin, in addition to a positive correlation with N-cadherin, Vimentin, and β-catenin, all three of which are linked to increased EMT [42-44, 54]. Data also indicated strong relationships to increased m-TOR, p-S6K1, and phosphorylated AKT, which play important roles in tumorigenesis, resistance to therapy, and cancer aggression [45–47]. After demonstrating that high levels of SHh/Gli activity were linked to the EMT pathway and its regulators, functional study showed the important effects of Gli-inhibition on reversing these expression patterns. Administration of a novel Gli-i small molecule, in addition to siRNA knockdown and AKT-i treatment, decreased expression of both proteins involved in cell cycle progression and thus proliferation, as well as markers responsible for the characteristic properties of EMT. Furthermore, when used in combination with Gli-i and AKT-i, N-Shh can rescue the effects of such inhibitors, suggesting that SHh/Gli signaling may regulate EMT via the AKT pathway in EAC. These findings were corroborated by results in both wound healing and transwell invasion assays following Gli-i treatment. Investigation via western blot and qRT-PCR paint a more complete picture of potential mechanisms responsible for declining migratory capabilities following SHh signal suppression.
EMT regulation is a sophisticated process; the SHh/Gli pathway may interact at various levels, including direct transcription of Snail, an EMT marker also investigated in this study [53]. Crosstalk between SHh and TGF-β signaling may play an important role in EMT regulation, while at the protein level, interactions between Gli and β-catenin at cell-cell junctions may directly contribute to cell mobility [53]. Thus, the molecular mechanism behind SHh pathway regulation of EMT in EAC requires further exploration.
Upstream of Gli, Vismodegib, a Smo inhibitor preventing activation of SHh target genes, has already been approved for basal cell carcinoma treatment, and is looking to expand its use to other cancers [55, 56]. While Vismodegib shows promise for reducing tumor size and lengthening progression free survival times, two mutations (D477G, E518K) conferring resistance mechanisms have already been identified in relapsed patients [57, 58]. Second generation Smo inhibitors, including Itraconazole and arsenic trioxide, both inhibit the growth of medulloblastoma and basal cell carcinoma, prolonging survival of mice with intracranial drug-resistant Smo [59, 60]. In addition to development of upstream SHh inhibitors, we propose that downstream Gli inhibition may prove even more effective than current Smo targeting that is currently in development. Gli suppression has the potential to circumvent upstream resistance mechanisms that have posed challenges thus far, given evidence of (1) non-canonical Gli activity independent of SHh in cancer, and (2) Gli as a target of other converging oncogenic signaling pathways [61–66]. While further research is still required, including additional patient tissue studies, Gli inhibition is a promising area for further investigation and potential drug development, owing to its close connection to EMT and regulation of tumor mobility in EAC.
This study provides evidence for aberrant upregulation of the Gli signaling pathway and a strong association between expression of Gli versus AKT, EMT, and cell cycle checkpoint markers in EAC. These findings suggest that activated SHh/Gli signaling may be a critical component in promoting cell survival, metastases, and resistance to chemotherapy. Inhibition of Gli and AKT pathway activity may thus serve as a potential therapeutic strategy for the treatment of human esophageal adenocarcinoma.
MATERIALS AND METHODS
Tissue specimens
Tissue specimens were collected from 24 patients who underwent surgical resection for esophageal adenocarcinoma (EAC) at the Thoracic Oncology Program at the University of California, San Francisco (UCSF). The study was approved by the Committee on Human Research (CHR approval number: H8714-11647-10) at UCSF; written, informed consent was obtained from each patient before specimen collection. These tissue samples were snap-frozen in liquid nitrogen immediately after resection and kept at -170°C before use.
Cell culture and drug treatment
Human esophageal adenocarcinoma (EAC) cell lines OE19 and OE33 were purchased from American Type Culture Collection (ATCC, Manassas, VA, USA) and cultured in RPMI 1640 (Life Technologies, Carlsbad, CA, USA), supplemented with 10% fetal bovine serum (FBS) and antibiotics. Cells were seeded one day before treatment with Gli inhibitor (Gli-i) or AKT inhibitor MK-2206 2HCI (AKT-i) (Selleck Chemicals, Houston, TX, USA), with or without N-Shh recombinant proteins (eBioscience, San Diego, CA, USA), and incubated for 30 hours, with DMSO vehicle as a control. Cells were then cultured for various periods of time prior to harvesting for subsequent analyses, including western blot, qRT-PCR, proliferation, migration, and invasion assays.
Protein extraction and western blot
OE19 and OE33 cells were pretreated with DMSO (control), N-Shh (0.5 mg/mL), Gli-i (500 nmol/L), or AKT-i (1 μmol/L), with or without subsequent stimulation by N-Shh (0.5mg/mL) for 24 hours. Total protein was then extracted from cell lines using M-PER Mammalian Protein Extraction Reagent (ThermoFisher Scientific, Waltham, MA, USA) and Complete Protease Inhibitor Cocktails (Roche, Lewes, UK), according to manufacturers’ protocols. Protein concentrations were determined using the Pierce BCA Protein Assay Kit (ThermoFisher Scientific), with 10 μg per sample run on 4-20% gradient SDS-polyacrylamide gels before transfer to Immobilon-P membranes (Millipore, Billerica, MA, USA). Nitrocellulose membranes were blocked in 5% nonfat milk and probed with the primary antibodies overnight at 4°C. The membranes were incubated with appropriate secondary antibodies, followed by detection using an ECL blotting analysis system. Antibodies applied to detect protein expression are as follows: Gli1 (1:100), Cyclin D1 (1:200), Snail (1:100), and Slug (1:150) (Santa Cruz Biotechnology, Dallas, TX, USA); Gli2 (1:150) and phosphorylated m-TOR (1:200) (Abcam, Cambridge, UK); N-Cadherin (1:200), E-cadherin (1:100), and β-catenin (1:200) (EMD Millipore, Hayward, CA, USA); Vimentin (1:200), p-AKT (1:200), AKT (1:200), p-ERK (1:200), and ERK (1:200) (Cell Signaling, Danvers, MA, USA); and GAPDH (1:500) (Sigma-Aldrich, St. Louis, MO, USA).
Luciferase reporter assay
To measure Gli-mediated Shh transcriptional activity, luciferase reporter constructs, 8x wild-type Gli binding site (8x Gliwt Luc) or 8x mutant Gli binding site (8x Glimut Luc) plasmids, and human Gli1 and Gli2 expression vectors (pcDNA3.1-Gli1/2) were co-transfected into EAC cell lines OE19 and OE33 in 24-well plates. The Renilla luciferase pRL-TK plasmid (Promega, Madison, WI, USA) was co-transfected to normalize for transfection efficiency according to manufacturer’s protocol. All transfections were performed using Lipofectamine2000 (Invitrogen, Carlsbad, CA, USA) in accordance with the manufacturer’s instructions. OE19 and OE33 cells were pretreated with DMSO, N-Shh (0.5 mg/mL), and/or Gli-i (500 nmol/L) for 30 minutes. After 24 hours, cells were lysed and subjected to the luciferase assays described above. Results are expressed as fold induction, which is the ratio of luciferase activity induced in Gli-transfected cells relative to basal luciferase activity in control transfected cells.
RNA extraction and qRT-PCR
Total RNA was isolated from OE19 and OE33 cultured cells using an RNeasy kit (Qiagen, Hilden, Germany). Genomic DNA contamination was eliminated by DNase I treatment, and reverse transcription was conducted with 500 ng RNA using the iScript cDNA synthesis kit (Bio-Rad, Hercules, CA, USA). The resulting cDNA were analyzed with real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) using Gene Expression Assays in a 7900 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) for 40 cycles (96°C for 15 seconds and 60°C for 1 minute). Expression of Vimentin, CDH2, CTNNB1, Snail, Slug, Zeb1, and GAPDH was quantified using commercially available primer and probe sequences (Applied Biosystems), and Relative Quantification Software (Applied Biosystems). Gene expression was normalized to GAPDH expression; experiments were performed in triplicate and repeated three times.
siRNA transfection
OE19 and OE33 cells were plated in six-well plates with fresh media without antibiotics for 24 hours before transfection. Transfection was performed using Lipofectamine2000 (Life Technologies, Carlsbad, CA) according to manufacturer’s protocol, with siRNAs targeting Gli1 (Assay ID 107670) and Gli2 (Assay ID 109640) at a total concentration of 50 nmol/L (Life Technologies). The efficiency of siRNA transfection was evaluated by western blot.
MTS survival, wound healing, and transwell invasion assays
In MTS survival assays, treated cells were plated in 96-well plates at a density of 500-1000 cells/well, with medium changed daily. Logarithmically growing cells were treated with increasing doses of Gli-i or AKT-i, and a DMSO vehicle control for 3 days. Cells were subsequently assessed for cell viability using CellTiter-Glo Luminescent Cell Viability Assay reagent (Promega, Madison, WI, USA) according to manufacturer’s instructions. Luminescence was measured using a GloMax-96 Microplate Luminometer (Promega), and percent cell survival was calculated based on the reading of untreated cells as 100% using GraphPad Prism 6.0 software to generate dose-response curves and IC50 values. All samples were conducted in triplicate, and data were standardized to that of DMSO treated cells.
In migration assays, four wounds were made in each treatment condition, and cell migration was determined by the average differences in distance between 0 and 24 hours following treatment. Transwell invasion assays (ThermoFisher Scientific) were performed according to manufacturer’s protocols. All experiments were conducted more than three times, and representative results were included in the text.
Statistical analysis
Two-sided student t-tests were performed for wound healing (migration) and transwell (invasion) assay analyses, while ANOVA and Scheffe’s tests were used to determine significance in luciferase reporter and qRT-PCR data. P-values of < 0.05, < 0.01, or < 0.001 were indicated as *, **, or ***, respectively, in figures.
Abbreviations
AKT-i: AKT inhibitor; Cdk: Cyclin-dependent kinase; EAC: esophageal adenocarcinoma; EGFR: epidermal growth factor receptor; EMT: epithelial-mesenchymal transition; Gli-i: Gli inhibitor; MMP: matrix metalloprotease; N: normal; p-: phosphorylated; Ptch: Patched; q-RTPCR: real-time reverse transcription polymerase chain reaction; SCC: squamous cell cancer; SHh: Sonic Hedgehog (pathway); Shh: Sonic Hedgehog (ligand); Smo: Smoothened; T: tumor; TF: transcription factor; TGF-β: transforming growth factor β; UCSF: University of California, San Francisco.
Author contributions
LW and JJ performed cell culture, western blot, cell viability, reporter, siRNA, migration and transwell assays, data analysis, and drafted the manuscript. YZ and ZT carried out qRT-PCR and conducted data and statistical analyses. DMJ participated in the study design. BH conceived of the study, participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
ACKNOWLEDGMENTS
The authors acknowledge the important contributions of the other members of the University of California San Francisco Thoracic Oncology Program and Biao He’s research team, including Fleur Leguay, BS; Ngoc Hoang, BS; Luis Acevedo, BS.
CONFLICTS OF INTEREST
The authors declare that they have no competing interests.
FUNDING
This work was supported by the Eileen D. Ludwig Endowed for Thoracic Oncology Research (to B. He); The Kazan, McClain, Abrams, Fernandez, Lyons, Greenwood, Harley & Oberman Foundation, the Ziegelmam Family Foundation, and the Barbara Isackson Lung Cancer Research Fund (to D. M. Jablons); the Science & Technology Support Project of Hebei Province (132077127D) (to L. Wang).
REFERENCES
1. Cancer Genome Atlas Research Network, Analysis Working Group: Asan University, BC Cancer Agency, Brigham and Women’s Hospital, Broad Institute, Brown University, Case Western Reserve University, Dana-Farber Cancer Institute, Duke University, Greater Poland Cancer Centre, Harvard Medical School, Institute for Systems Biology, KU Leuven, et al. Integrated genomic characterization of oesophageal carcinoma. Nature. 2017; 541:169–75. https://doi.org/10.1038/nature20805.
2. Shaheen N, Ransohoff DF. Gastroesophageal reflux, barrett esophagus, and esophageal cancer: Clinical applications. JAMA. 2002; 287:1982–6.
3. Napier KJ, Scheerer M, Misra S. Esophageal cancer: A review of epidemiology, pathogenesis, staging workup and treatment modalities. World J Gastrointest Oncol. 2014; 6:112–20. https://doi.org/10.4251/wjgo.v6.i5.112.
4. Chen X, Yang CS. Esophageal adenocarcinoma: A review and perspectives on the mechanism of carcinogenesis and chemoprevention. Carcinogenesis. 2001; 22:1119–29.
5. Rubenstein JH, Shaheen NJ. Epidemiology, diagnosis, and management of esophageal adenocarcinoma. Gastroenterology. 2015; 149:302–17 e1. https://doi.org/10.1053/j.gastro.2015.04.053.
6. De Craene B, Berx G. Regulatory networks defining emt during cancer initiation and progression. Nature Reviews Cancer. 2013; 13:97–110. https://doi.org/10.1038/nrc3447.
7. Hong KO, Kim JH, Hong JS, Yoon HJ, Lee JI, Hong SP, Hong SD. Inhibition of akt activity induces the mesenchymal-to-epithelial reverting transition with restoring e-cadherin expression in kb and koscc-25b oral squamous cell carcinoma cells. J Exp Clin Cancer Res. 2009; 28:28. https://doi.org/10.1186/1756-9966-28-28.
8. Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. J Clin Invest. 2009; 119:1420–8. https://doi.org/10.1172/JCI39104.
9. Larue L, Bellacosa A. Epithelial-mesenchymal transition in development and cancer: Role of phosphatidylinositol 3’ kinase/akt pathways. Oncogene. 2005; 24:7443–54. https://doi.org/10.1038/sj.onc.1209091.
10. Soltermann A, Tischler V, Arbogast S, Braun J, Probst-Hensch N, Weder W, Moch H, Kristiansen G. Prognostic significance of epithelial-mesenchymal and mesenchymal-epithelial transition protein expression in non-small cell lung cancer. Clinical Cancer Research. 2008; 14:7430–7. https://doi.org/10.1158/1078-0432.ccr-08-0935.
11. Yang MH, Wu MZ, Chiou SH, Chen PM, Chang SY, Liu CJ, Teng SC, Wu KJ. Direct regulation of twist by hif-1 alpha promotes metastasis. Nature Cell Biology. 2008; 10:295–305. https://doi.org/10.1038/ncb1691.
12. Vendrell JA, Thollet A, Nguyen NT, Ghayad SE, Vinot S, Bieche I, Grisard E, Josserand V, Coll JL, Roux P, Corbo L, Treilleux I, Rimokh R, et al. Znf217 is a marker of poor prognosis in breast cancer that drives epithelial-mesenchymal transition and invasion. Cancer Research. 2012; 72:3593–606. https://doi.org/10.1158/0008-5472.can-11-3095.
13. Van den Broeck A, Vankelecom H, Van Eijsden R, Govaere O, Topal B. Molecular markers associated with outcome and metastasis in human pancreatic cancer. Journal of Experimental & Clinical Cancer Research. 2012; 31:68. https://doi.org/10.1186/1756-9966-31-68.
14. Lauth M, Bergstrom A, Shimokawa T, Toftgard R. Inhibition of gli-mediated transcription and tumor cell growth by small-molecule antagonists. Proc Natl Acad Sci USA. 2007; 104:8455–60.
15. Lauth M, Toftgard R. Non-canonical activation of gli transcription factors. Cell Cycle. 2007; 6:2458–63.
16. Lauth M, Toftgard R. The hedgehog pathway as a drug target in cancer therapy. Curr Opin Investig Drugs. 2007; 8:457–61.
17. Mimeault M, Batra SK. Frequent deregulations in the hedgehog signaling network and cross-talks with the epidermal growth factor receptor pathway involved in cancer progression and targeted therapies. Pharmacological Reviews. 62:497–524. https://doi.org/10.1124/pr.109.002329.
18. Stanton BZ, Peng LF. Small-molecule modulators of the sonic hedgehog signaling pathway. Molecular Biosystems. 6:44–54. https://doi.org/10.1039/b910196a.
19. Tostar U, Malm CJ, Meis-Kindblom JM, Kindblom LG, Toftgard R, Unden AB. Deregulation of the hedgehog signalling pathway: A possible role for the ptch and sufu genes in human rhabdomyoma and rhabdomyosarcoma development. J Pathol. 2006; 208:17–25. https://doi.org/10.1002/path.1882.
20. Kinzler KW, Bigner SH, Bigner DD, Trent JM, Law ML, Obrien SJ, Wong AJ, Vogelstein B. Identification of an amplified, highly expressed gene in a human glioma. Cytogenetics and Cell Genetics. 1987; 46:639.
21. Chi SM, Huang SH, Li CX, Zhang XL, He NG, Bhutani MS, Jones D, Castro CY, Logrono R, Haque A, Zwischenberger J, Tyring SK, Zhang HW, et al. Activation of the hedgehog pathway in a subset of lung cancers. Cancer Letters. 2006; 244:53–60. https://doi.org/10.1016/j.canlet.2005.11.036.
22. Thompson MC, Fuller C, Hogg TL, Dalton J, Finkelstein D, Lau CC, Chintagumpala M, Adesina A, Ashley DM, Kellie SJ, Michael DT, Curran T, Gajjar A, et al. Genornics identifies medulloblastoma subgroups that are enriched for specific genetic alterations. Journal of Clinical Oncology. 2006; 24:1924–31. https://doi.org/10.1200/jco.2005.04.4974.
23. Thayer SP, di Magliano MP, Heiser PW, Nielsen CM, Roberts DJ, Lauwers GY, Qi YP, Gysin S, Fernandez-del Castillo CF, Yajnik V, Antoniu B, McMahon M, Warshaw AL, et al. Hedgehog is an early and late mediator of pancreatic cancer tumorigenesis. Nature. 2003; 425:851–6. https://doi.org/10.1038/nature02009.
24. Taylor MD, Liu L, Raffel C, Hui CC, Mainprize TG, Zhang XY, Agatep R, Chiappa S, Gao LZ, Lowrance A, Hao AH, Goldstein AM, Stavrou T, et al. Mutations in sufu predispose to medulloblastoma. Nature Genetics. 2002; 31:306–10. https://doi.org/10.1038/ng916.
25. Li H, Lui N, Cheng T, Tseng HHK, Yue DS, Giroux-Leprieur E, Do HT, Sheng Q, Jin JQ, Luh TW, Jablons DM, He BA. Gli as a novel therapeutic target in malignant pleural mesothelioma. PLoS One. 2013; 8. https://doi.org/10.1371/journal.pone.0057346.
26. Long B, Zhu HL, Zhu CX, Liu T, Meng WT. Activation of the hedgehog pathway in chronic myelogeneous leukemia patients. Journal of Experimental & Clinical Cancer Research. 2011; 30:8. https://doi.org/10.1186/1756-9966-30-8.
27. Riobo NA, Manning DR. Pathways of signal transduction employed by vertebrate hedgehogs. Biochem J. 2007; 403:369–79. https://doi.org/10.1042/BJ20061723.
28. Brennan D, Chen X, Cheng L, Mahoney M, Riobo NA. Noncanonical hedgehog signaling. Vitam Horm. 2012; 88:55–72. https://doi.org/10.1016/B978-0-12-394622-5.00003-1.
29. Javelaud D, Pierrat MJ, Mauviel A. Crosstalk between tgf-beta and hedgehog signaling in cancer. FEBS Lett. 2012; 586:2016–25. https://doi.org/10.1016/j.febslet.2012.05.011.
30. Perrot CY, Javelaud D, Mauviel A. Overlapping activities of tgf-beta and hedgehog signaling in cancer: Therapeutic targets for cancer treatment. Pharmacol Ther. 2013; 137:183–99. https://doi.org/10.1016/j.pharmthera.2012.10.002.
31. Mangelberger D, Kern D, Loipetzberger A, Eberl M, Aberger F. Cooperative hedgehog-egfr signaling. Front Biosci (Landmark Ed). 2012; 17:90–9.
32. Riobo NA, Lu K, Ai X, Haines GM, Emerson CP Jr. Phosphoinositide 3-kinase and akt are essential for sonic hedgehog signaling. Proc Natl Acad Sci U S A. 2006; 103:4505–10. https://doi.org/10.1073/pnas.0504337103.
33. Sandhiya S, Melvin G, Kumar SS, Dkhar SA. The dawn of hedgehog inhibitors: Vismodegib. J Pharmacol Pharmacother. 2013; 4:4–7. https://doi.org/10.4103/0976-500X.107628.
34. Agyeman A, Jha BK, Mazumdar T, Houghton JA. Mode and specificity of binding of the small molecule gant61 to gli determines inhibition of gli-dna binding. Oncotarget. 2014; 5:4492–503. https://doi.org/10.18632/oncotarget.2046.
35. Bosco-Clement G, Zhang F, Chen Z, Zhou HM, Li H, Mikami I, Hirata T, Yagui-Beltran A, Lui N, Do HT, Cheng T, Tseng HH, Choi H, et al. Targeting gli transcription activation by small molecule suppresses tumor growth. Oncogene. 2014; 33:2087–97. https://doi.org/10.1038/onc.2013.164.
36. Alexaki VI, Javelaud D, Van Kempen LCL, Mohammad KS, Dennler S, Luciani F, Hoek KS, Juarez P, Goydos JS, Fournier PJ, Sibon C, Bertolotto C, Verrecchia F, et al. Gli2-mediated melanoma invasion and metastasis. Journal of the National Cancer Institute. 2010; 102:1148–59. https://doi.org/10.1093/jnci/djq257.
37. Inaguma S, Kasai K, Hashimoto M, Ikeda H. Gli1 modulates emt in pancreatic cancer-letter. Cancer Research. 2012; 72:3702–3. https://doi.org/10.1158/0008-5472.can-12-0379.
38. Joost S, Almada LL, Rohnalter V, Holz PS, Vrabel AM, Fernandez-Barrena MG, McWilliams RR, Krause M, Fernandez-Zapico ME, Lauth M. Gli1 inhibition promotes epithelial-to-mesenchymal transition in pancreatic cancer cells. Cancer Research. 2012; 72:88–99. https://doi.org/10.1158/0008-5472.can-10-4621.
39. Yue D, Li H, Che J, Zhang Y, Tseng HH, Jin JQ, Luh TM, Giroux-Leprieur E, Mo M, Zheng Q, Shi H, Zhang H, Hao X, et al. Hedgehog/gli promotes epithelial-mesenchymal transition in lung squamous cell carcinomas. J Exp Clin Cancer Res. 2014; 33:34. https://doi.org/10.1186/1756-9966-33-34.
40. Faiao-Flores F, Alves-Fernandes DK, Pennacchi PC, Sandri S, Vicente AL, Scapulatempo-Neto C, Vazquez VL, Reis RM, Chauhan J, Goding CR, Smalley KS, Maria-Engler SS. Targeting the hedgehog transcription factors gli1 and gli2 restores sensitivity to vemurafenib-resistant human melanoma cells. Oncogene. 2017; 36:1849–61. https://doi.org/10.1038/onc.2016.348.
41. Chen JS, Li HS, Huang JQ, Zhang LJ, Chen XL, Wang Q, Lei J, Feng JT, Liu Q, Huang XH. Down-regulation of gli-1 inhibits hepatocellular carcinoma cell migration and invasion. Mol Cell Biochem. 2014; 393:283–91. https://doi.org/10.1007/s11010-014-2071-x.
42. Mosnier JF, Kandel C, Cazals-Hatem D, Bou-Hanna C, Gournay J, Jarry A, Laboisse CL. N-cadherin serves as diagnostic biomarker in intrahepatic and perihilar cholangiocarcinomas. Mod Pathol. 2009; 22:182–90. https://doi.org/10.1038/modpathol.2008.123.
43. Mendez MG, Kojima S, Goldman RD. Vimentin induces changes in cell shape, motility, and adhesion during the epithelial to mesenchymal transition. FASEB J. 2010; 24:1838–51. https://doi.org/10.1096/fj.09-151639.
44. Jiang YG, Luo Y, He DL, Li X, Zhang LL, Peng T, Li MC, Lin YH. Role of wnt/beta-catenin signaling pathway in epithelial-mesenchymal transition of human prostate cancer induced by hypoxia-inducible factor-1alpha. Int J Urol. 2007; 14:1034–9. https://doi.org/10.1111/j.1442-2042.2007.01866.x.
45. Karlsson E, Magic I, Bostner J, Dyrager C, Lysholm F, Hallbeck AL, Stal O, Lundstrom P. Revealing different roles of the mtor-targets s6k1 and s6k2 in breast cancer by expression profiling and structural analysis. PLoS One. 2015; 10:e0145013. https://doi.org/10.1371/journal.pone.0145013.
46. Wang Y, Ding Q, Yen CJ, Xia W, Izzo JG, Lang JY, Li CW, Hsu JL, Miller SA, Wang X, Lee DF, Hsu JM, Huo L, et al. The crosstalk of mtor/s6k1 and hedgehog pathways. Cancer Cell. 2012; 21:374–87. https://doi.org/10.1016/j.ccr.2011.12.028.
47. Saeed N, Shridhar R, Hoffe S, Almhanna K, Meredith KL. Akt expression is associated with degree of pathologic response in adenocarcinoma of the esophagus treated with neoadjuvant therapy. J Gastrointest Oncol. 2016; 7:158–65. https://doi.org/10.3978/j.issn.2078-6891.2015.067.
48. Tomizawa Y, Wu TT, Wang KK. Epithelial mesenchymal transition and cancer stem cells in esophageal adenocarcinoma originating from barrett’s esophagus. Oncol Lett. 2012; 3:1059–63. https://doi.org/10.3892/ol.2012.632.
49. Liu Y, El-Naggar S, Darling DS, Higashi Y, Dean DC. Zeb1 links epithelial-mesenchymal transition and cellular senescence. Development. 2008; 135:579–88. https://doi.org/10.1242/dev.007047.
50. Abbas T, Dutta A. P21 in cancer: Intricate networks and multiple activities. Nat Rev Cancer. 2009; 9:400–14. https://doi.org/10.1038/nrc2657.
51. Dai M, Al-Odaini AA, Fils-Aime N, Villatoro MA, Guo J, Arakelian A, Rabbani SA, Ali S, Lebrun JJ. Cyclin d1 cooperates with p21 to regulate tgfbeta-mediated breast cancer cell migration and tumor local invasion. Breast Cancer Res. 2013; 15:R49. https://doi.org/10.1186/bcr3441.
52. Xie L, Law BK, Chytil AM, Brown KA, Aakre ME, Moses HL. Activation of the erk pathway is required for tgf-beta1-induced emt in vitro. Neoplasia. 2004; 6:603–10. https://doi.org/10.1593/neo.04241.
53. Zhang J, Tian XJ, Xing J. Signal transduction pathways of emt induced by tgf-beta, shh, and wnt and their crosstalks. J Clin Med. 2016; 5. https://doi.org/10.3390/jcm5040041.
54. Nakajima S, Doi R, Toyoda E, Tsuji S, Wada M, Koizumi M, Tulachan SS, Ito D, Kami K, Mori T, Kawaguchi Y, Fujimoto K, Hosotani R, et al. N-cadherin expression and epithelial-mesenchymal transition in pancreatic carcinoma. Clin Cancer Res. 2004; 10:4125–33. https://doi.org/10.1158/1078-0432.CCR-0578-03.
55. Azmi AS. Unveiling the role of nuclear transport in epithelial-to-mesenchymal transition. Curr Cancer Drug Targets. 2013; 13:906–14.
56. Dlugosz A, Agrawal S, Kirkpatrick P. Vismodegib. Nat Rev Drug Discov. 2012; 11:437–8. https://doi.org/10.1038/nrd3753.
57. Amakye D, Jagani Z, Dorsch M. Unraveling the therapeutic potential of the hedgehog pathway in cancer. Nat Med. 2013; 19:1410–22. https://doi.org/10.1038/nm.3389.
58. Macha MA, Batra SK, Ganti AK. Profile of vismodegib and its potential in the treatment of advanced basal cell carcinoma. Cancer Manag Res. 2013; 5:197–203. https://doi.org/10.2147/CMAR.S45976.
59. Kim J, Tang JY, Gong R, Kim J, Lee JJ, Clemons KV, Chong CR, Chang KS, Fereshteh M, Gardner D, Reya T, Liu JO, Epstein EH, et al. Itraconazole, a commonly used antifungal that inhibits hedgehog pathway activity and cancer growth. Cancer Cell. 2010; 17:388–99. https://doi.org/10.1016/j.ccr.2010.02.027.
60. Ally MS, Ransohoff K, Sarin K, Atwood SX, Rezaee M, Bailey-Healy I, Kim J, Beachy PA, Chang AL, Oro A, Tang JY, Colevas AD. Effects of combined treatment with arsenic trioxide and itraconazole in patients with refractory metastatic basal cell carcinoma. JAMA Dermatol. 2016; 152:452–6. https://doi.org/10.1001/jamadermatol.2015.5473.
61. Lauth M, Toftgard R. Non-canonical activation of gli transcription factors: implications for targeted anti-cancer therapy. Cell Cycle. 2007; 6:2458–63.
62. Mimeault M, Batra SK. Frequent deregulations in the hedgehog signaling network and cross-talks with the epidermal growth factor receptor pathway involved in cancer progression and targeted therapies. Pharmacol Rev. 2010; 62:497–524. https://doi.org/10.1124/pr.109.002329.
63. Guo X, Wang XF. Signaling cross-talk between tgf-beta/bmp and other pathways. Cell Res. 2009; 19:71–88. https://doi.org/10.1038/cr.2008.302.
64. Schnidar H, Eberl M, Klingler S, Mangelberger D, Kasper M, Hauser-Kronberger C, Regl G, Kroismayr R, Moriggl R, Sibilia M, Aberger F. Epidermal growth factor receptor signaling synergizes with hedgehog/gli in oncogenic transformation via activation of the mek/erk/jun pathway. Cancer Res. 2009; 69:1284–92. https://doi.org/10.1158/0008-5472.can-08-2331.
65. Pasca di Magliano M, Sekine S, Ermilov A, Ferris J, Dlugosz AA, Hebrok M. Hedgehog/ras interactions regulate early stages of pancreatic cancer. Genes Dev. 2006; 20:3161–73. https://doi.org/10.1101/gad.1470806.
66. Stecca B, Mas C, Clement V, Zbinden M, Correa R, Piguet V, Beermann F, Ruiz IAA. Melanomas require hedgehog-gli signaling regulated by interactions between gli1 and the ras-mek/akt pathways. Proc Natl Acad Sci U S A. 2007; 104:5895–900. https://doi.org/10.1073/pnas.0700776104.