1. Introduction
Pancreatic cancer is one of the most common causes of cancer-related deaths worldwide. Its two major histological subtypes are pancreatic ductal adenocarcinoma (PDAC), accounting for 90% of all cases, and pancreatic neuroendocrine neoplasm (panNEN), which makes up 3-5% of all cases. PanNEN is classified into well-differentiated pancreatic NE tumor and poorly-differentiated pancreatic NE carcinoma (panNEC). Although PDAC and panNEN are commonly thought to be different diseases with distinct biology, cell of origin, and genomic abnormalities, the idea that PDAC and panNEC share common cells of origin has been gaining support. This is supported by molecular profiling data suggesting that panNEC is genetically and phenotypically related to PDAC [
1]. PanNENs represent a heterogeneous group of epithelial tumors with NE differentiation (NED) that are classified into well-differentiated pancreatic NE tumors (panNETs), including G1, G2, and G3 tumors, and poorly differentiated panNECs [
2]. PanNETs can be regarded as a unique category, where G1-G2 tumors may progress to G3 tumors mainly driven by DAXX/ATRX mutations [
2]. Conversely, panNECs display totally different histomolecular features more closely related to PDAC, including TP53 and Rb alterations [
2] was therefore not surprising that two PDAC-derived tumor cell lines of the quasi-mesenchymal subtype, PANC-1 and MIA PaCa-2, were demonstrated to harbor a NED phenotype [
3]. PANC-1 expresses CK5.6, MNF-116, vimentin (VIM), chromogranin A (CHGA), neural cell adhesion molecule (NCAM/CD56) and somatostatin receptor-2 (SSTR2) but not E-cadherin (ECAD), synaptophysin (SYP) or neurotrophin receptor-1 (NTR1), while MIA PaCa-2 expresses CK5.6, AE1/AE3, ECAD, VIM, CHGA, SYP, SSTR2 and NTR1 but not NCAM [
3]. In addition to the NED markers, we [
4] and others [
5] demonstrated that PANC-1 and MIA PaCa-2 cells express genes associated with endocrine/neuroendocrine differentiation such as
MAFA,
NEURODI,
PDX1, and
NEUROG3. Of note, the PANC-1 cell line has been shown to exhibit epithelial-mesenchymal plasticity (EMP). Parental cultures comprise several clonal subpopulations with different EM transdifferentiation (EMT) phenotypes of which some are more epithelial and others more mesenchymal in nature [
6]. The importance of EMT is also considered in other gastrointestinal (GI) cancers such as colorectal NEC (coloNEC) and a distinct subtype of ovarian cancer (OC), e.g., the heterogeneous small cell carcinoma of the ovary hypercalcemic type (SCCOHT). The SCCOHT represents a rare form of an aggressive ovarian tumor, which is predominantly observed in young women. This malignant neoplasia is associated with paraneoplastic hypercalcemia and affected patients often have a lethal outcome already within a few months after diagnosis [
7,
8]. Cellular models for this tumor entity were established and characterized as the SCCOHT-1 [
9] and BIN-67 cell lines [
10,
11]. Indeed, the characterization of SCCOHT-1 cells and chemotherapeutic responses has revealed a constitutive expression of NED markers such as NCAM in the original patient tumor and derived mouse xenograft tumors [
9,
12]. Thus, the heterogeneity and plasticity of SCCOHT-1 cells may also involve EMP and maturation along a NED phenotype [
13].
A series of studies suggests that NED is closely associated with EMT in different tumor types, mainly prostate, lung and colon/pancreas. In prostate carcinoma (PCa), androgen deprivation has been shown to activate both EMT and NE transdifferentiation programs [
14]. Many factors have been associated with the onset and progression of NED in clinically typical prostate adenocarcinomas including loss of androgen receptor (AR) expression and/or signaling, conventional therapy, and dysregulated cytokine function. The AR is a critical driver of tumor progression as well as therapeutic response in patients with metastatic castrate-resistant PCa (CRPC) [
15]. Emerging evidence suggests that the acquisition of EMT and a cancer stem cell (CSC) phenotype are associated with the development of NED in PCa [
16]. EMT and NED may also be induced by androgen-targeted therapy [
14] and is considered a resistance mechanism to treatments in PCa [
17]. The resulting NE PCa (NEPC) is highly aggressive exhibiting reactivation of developmental programs associated with EMT induction and stem cell-like characteristics.
Interestingly, an inverse relationship between NED and EMT has also been described in some tumors, i.e., small cell lung carcinoma (SCLC), Merkel cell carcinoma and gastroenteropancreatic (GEP)-NET [
18]. In SCLC, an association was revealed between the loss of NED and EMT induction [
19] as inferred from the observation that the low NED subtype had undergone EMT and had activated - amongst others - the TGF-β pathway. However, differential effects of TGF-β on both programs were observed in SCLC in that TGF-β seems to be required for promoting EMT but not NED. In the panNET cell lines, BON-1 (BON), QGP-1 and NT-3 previous work has shown that the ECAD and VIM expression profiles are indicative of a well-differentiated epithelial phenotype [
20]. In BON cells, TGF-β1 has been shown to control proliferation and NED through the SST/SSTR system [
21,
22]. Of note, disrupting either the TGF-β or SST signaling pathway resulted in NED-mesenchymal transition, which is characterized by the loss of NED markers, decreased ECAD and elevated VIM expression. This inverse correlation of TGF-β sigaling activity and EMT was surprising since TGF-β is known as one of the most potent inducers of EMT. In PCa, aberrant TGF-β signaling accelerates progression in a transgenic mouse model via effects on EMT and NED driving tumor progression to CRPC [
23].
Given the concept of a positive association between EMT and NED mainly arising from the PCa model, an important issue remains whether the reverse process, mesenchymal-epithelial transition (MET) affects NED. We have recently shown that PDAC-derived tumor cells can be forced to undergo MET in response to a transdifferentiation culture through exposure to a combination of three cytokines, IL-1β, IFN-γ and TNF-α (TDC-IIT) [6}. Moreover, a member of the TGF-β superfamily of growth and differentiation factors, bone morphogenetic protein-7 (BMP-7), has been reported to be able to induce MET in adult renal fibroblasts of the injured kidney [
24], hepatic stellate cells [
25] and melanoma cells [
26], generating functional epithelial cells [
27]. Interestingly, BMP-7 has been shown to induce MET through downregulation of the EMT-related transcription factor SNAIL [
28], SNAIL-induced α-smooth muscle actin, and concomitant upregulation of ECAD. BMP-7 also acts as an inhibitor of fibrotic progression in many organs through activation of the SMAD1/5 arm of TGF-β signaling and inhibition of TGF-β–mediated EMT [
29,
30] via suppression of canonical TGF-β/SMAD2/3 signaling [
31].
Given the above-mentioned findings, particularly the inverse relationship between EMT and NED in panNET, we decided to study the effects of EMT and MET inducers on EMP and NED in pancreatic cells (normal duct cells and PDAC-derived tumor cells of the epithelial and quasi-mesenchymal subtype), and in cell lines of SCCOHT. Specifically, we asked whether forced conversion from epithelial to mesenchymal (via treatment with TGF-β1) or vice versa (via stimulation with TDC-IIT or BMP-7) will affect NED markers in the same way as EMT markers.
4. Discussion
The relationship between EMT and NED, and the role of TGF-β signaling in controlling these differentiation programs appears to vary among different tumor entities, and even between cancers affecting the same organ, such as the pancreas. While these associations have been studied in detail in a panNET cell line [
21], only little information was available for PDAC in this respect. For other tumor entities such as SCCOHT no such data, whatsoever, were available. We thus initially focussed on the pancreatic model comparing poorly differentiated tumor cells of the quasi-mesenchymal subtype (PANC-1, MIA PaCa-2) with moderately differentated ones of the epithelial subtype (BxPC-3) and with the presumed progenitor cells, non-transformed pancreatic ductal epithelial cells (HPDE6c7). Based on our assumption that NED is associated with EMT and a mesenchymal phenotype, we hypothesised that HPDE6c7 and BxPC-3 cells should not exhibit a NED phenotype. Indeed, HPDE6c7 expressed much lower levels of NED-associated markers than PANC-1 and MIA PaCa-2. The same was true for BxPC-3 with respect to CHGA, SYP and NSE although, unlike HPDE6c7, not for SSTR2 and SSTR5. Hence, we have shown in PANC-1 and MIA PaCa-2 cells endogenous expression NED markers, some of which had previously been detected by flow cytometry [
3].
PANC-1 cells exhibit EMP, meaning that parental cultures of this cell line consist of a mixture of subclones each displaying a different EMT phenotype despite being genetically identical. A previous histomorphological subtyping with a panel of epithelial and mesenchymal markers has shown that these clones can be grossly classified as epithelial (E-type), mesenchymal (M-type), or mixed [
6]. We, therefore, considered it appropriate to test whether NED markers are enriched in either the E or the M-type clones. Monitoring seven (three E-type and four M-type) single cell-derived clones for expression of SYP, CHGA, NCAM, NSE, GLUT2, SSTR2, and SSTR5 showed that M-type clones (P1C3, P3D10, P4B9, P2E8) present with higher levels of these NED markers, except for GLUT2, than E-type clones (P4B11, P3D2, P1G7). This led us to conclude that NED is preferentially associated with a mesenchymal phenotype.
Next, we sought to know if TGF-β1, a powerful promoter of EMT, and BMP-7, another member of the TGF-β superfamily of growth and differentiation factors and promoter of MET, impact EMT and NED-associated gene expression in pancreatic tumor cells. While PANC-1 cells are highly sensitive to this growth factor, MIA PaCa-2 cells are refractory due to a defective type II receptor. These cells could thus only be employed to study the effects of BMP-7. Treatment of PANC-1 cells with TGF-β1 downregulated ECAD and RAC1b and upregulated VIM and SNAIL. However, both PANC-1 and MIA PaCa-2 cells responded to treatment with BMP-7 with induction of SNAIL and VIM, while RAC1b was only suppressed in PANC-1 and RAC1 remained unaltered in both cell lines. The effect of TGF-β1 on VIM and SNAIL was more potent than that of BMP-7. Cotreatment with both growth factors acted in either an additive or synergistic manner to suppress ECAD (only in PANC-1 since MIA PaCa-2 are ECAD-null) and RAC1b, and to enhance SNAIL and VIM expression at both the RNA and protein level. We thus concluded that TGF-β1 and BMP-7 in PANC1, and BMP-7 in MIA PaCa-2 cells, can induce EMT, and that BMP-7 can synergize with TGF-β1 in EMT induction. The observation that BMP-7
promoted EMT was surprising since this growth factor has been identified in other cellular models as either an inhibtor of EMT or even a promoter of MET [
24,
25,
26,
27,
28].
Both TGF-β1 and BMP-7 were capable of inducing SYP protein in PANC-1, however, only TGF-β1 appears to accomplish this by a transcriptional mechanism. Intriguingly, a very strong inductive effect of TGF-β1 or BMP-7 was observed on SSTR2 in PANC-1 cells that could be further enhanced synergistically by combined treatment. In MIA PaCa-2 cells, however, SSTR2 mRNA levels remained unchanged in response to TGF-β1 or BMP-7, while the related SSTR5 was strongly downregulated by TGF-β1, but upregulated by BMP-7 in both cell lines. In addition, TGF-β1 treatment of PANC-1 induced CHGA, NCAM, and NSE, while BMP-7 treatment only upregulated NSE, but downregulated CHGA. These results clearly indicate that at least in control of SSTR5 and CHGA, BMP-7 can also act in an antagonistic fashion to TGF-β1.
SSTR2 and 5 represent not only established markers of NED but also possess tumor-relevant functions. SSTR2 is an inhibitory G protein-coupled receptor, the expression of which is lost in most human pancreatic cancers [
35]. Of note, murine Sstr2 has been identified as a transcriptional target of TGF-β [
35], which suggested the possibility that loss of SMAD4 accounts for the loss of SSTR2 expression in human PDAC. This event may contribute to a growth advantage of tumor cells [
35] and is consistent with findings that SSTR2 exhibits anti-tumor properties. Here, we have confirmed - for the first time - in human PDAC-derived tumor cells a strong positive regulation of SSTR2 by TGF-β1. Also, for the first time,
SSTR5 was identified here as a negative transcriptional target gene of TGF-β1 as evidenced by the dramatic downregulation of its mRNA. This mode of regulation suggests the possibility that SSTR5 normally antagonises TGF-β-dependent EMT or even other cellular responses to TGF-β, such as growth arrest. Hence, while
SSTR2 qualifies as a gene involved in growth arrest in accordance with the proposed anti-tumor function, the reverse may be true for
SSTR5. Moreover, while SSTR2 appears to be involved in mesenchymal conversion, SSTR5 may have a role in promoting MET or an epithelial phenotype. An interesting finding in this context came from a study with a highly invasive paclitaxel-resistant OC cell line. This cell line expresses CD105/endoglin, a stem cell marker and TGF-β co-receptor that may promote EMT and metastasis of OC by inhibiting expression of ECAD. Of note, after CD105 knockdown, the expression of both SSTR5 and ECAD (amongst others) was markedly upregulated [
36]. Conversely, coactivation of SSTR2 in PDAC cells led to increased expression of mesenchymal markers and decreased expression of an epithelial marker [
37]. Moreover, the expression of SSTR2 (along with those of SNAIL, SLUG and VIM) was associated with invasive non-functioning NETs of the pituitary [
38].
Treatment of PANC-1 or MIA PaCa-2 cells with BMP-7 downregulated ECAD and upregulated VIM and SNAIL. This pro-EMT effect was quite surprising as BMP-7 has been identified previously as a MET-inducing (and thus anti-EMT) factor in a range of different cell types such as alveolar type II cells, adult renal epithelial tubular cells and fibroblasts, hepatic stellate cells, and melanoma cells [
24,
25,
26,
27,
28]. Moreover, BMP-7 stimulation of PANC-1 or MIA PaCa-2 cells also altered the expression of some NED genes,
i.e.,
SSTR2, the same way as the EMT-associated genes. The SSTR2 mRNA was induced by BMP-7, although the extent of induction was not as great as that with TGF-β1. In contrast to the suppressive effect of TGF-β1, BMP-7 upregulated SSTR5 mRNA, which would be consistent with a pro-epithelial effect based on the above proposed function for SSTR5. However, the observation that BMP-7 induced only a few NED markers, while others were either inhibited or remained unaffected questions its role as a general promoter of NED. Surprisingly, upon combined stimulation of PANC-1 cells with TGF-β1 and BMP-7, additive effects on induction of
SSTR2 were noted. Together, this clearly suggests an association of EMT and NED through TGF-β signaling, while BMP-7 only partially shares this ability in common with TGF-β.
Prompted by the newly discovered positive association of EMT and TGF-β signaling with NED, we evaluated in another set of experiments the possibility that, conversely, MET in PANC-1 cells is associated with a loss of NED. In agreement with this hypothesis, we observed that during IIT-induced MET most NED markers, except SSTR2, were downregulated along with VIM, RAC1 and SNAIL at the mRNA (
Figure 5) and protein [
6] level, while concomitantly, the epithelial markers ECAD, CLDN4, GRHL2, OVOL2, CK19 and RAC1b were all upregulated [
6]. This clearly shows a simultaneous loss of NED and mesenchymal markers during MET and the acquisition of an epithelial phenotype. Mechanistically, this may – at least in part - be mediated through inhibition of the SMAD2/3 arm of TGF-β signaling, since all markers that are responsive to TGF-β1 treatment were also affected by TDC-IIT but in an antagonistic fashion.
As control for cells with a strong NED phenotype, we employed the panNET cell lines, BON and NT-3, which are both epithelial in nature [
19]. This is in sharp contrast to PDAC cells with NED, which are poorly differentiated/quasi-mesenchymal. In this study, we have, therefore, revealed fundamental differences between two major types of pancreatic cancer, PDAC and panNET, with respect to the association of NED with EMT and TGF-β signaling. In panNET, NED is associated with a well-differentiated epithelial phenotype and functional TGF-β and SST signaling, and defective TGF-β/SST signaling causes loss of NED with mesenchymal conversion [
21]. However, in the present study the reverse situation is operating; NED occurs in poorly differentiated mesenchymal cells and can still be further enhanced by activation of TGF-β or BMP-7 signaling provided the cells have retained responsiveness to these growth factors. Thus, the role of SST/SSTR in TGF-β signaling may differ between panNET (BON, NT-3) and PDAC, which is also supported by the antagonistic regulation of
SSTR5 by TGF-β1. We have thus identified a newly distinguishing feature between panNET and PDAC. Consequently, it will be highly intriguing to test these properties in panNEC, a pancreatic cancer entity that combines features of both panNET and PDAC. Intriguingly, we have shown recently that not only treatment with TGF-β but also the stimulation with SST, or the SST analogs octreotide and lanreotide, was able to regulate a set of NED genes and alter the NED state [
22].
While the phenotypic association between EMT and NED seems to be well established in some cancers, this is not the case for the underlying molecular mechanism(s). Initial insights came from the PCa model [
39] with the identification of microRNA-147b as an inducer of NED through targeting the ribosomal protein PRS15A [
40]. More recently, activation of NFκB-STAT3 signaling by tumor protein D52, isoform 3 (TPD52) has been found to induce distinct NED features (as measured by CHGA and NSE) through EMT under androgen-depleted conditions [
39]. Moreover, the authors were able to show that TPD52 also positively regulates EMT of PCa cells towards NED (as revealed by induction of N-cadherin, VIM and ZEB1, another EMT-associated transcription factor) via activation of NFκB-STAT3. These changes were orchestrated by SNAIL, since silencing of
SNAI1 in TPD52-positive cells blocked the progression of NED [
39]. SNAIL may thus promote tumor aggressiveness in PCa cells through multiple processes; induction of EMT to promote migration, while, in turn, induction of NED promotes tumor proliferation through a paracrine mechanism [
41]. In addition, ZEB1 has been shown to promote NED in PCa [
42]. Liu and colleagues investigated the molecular mechanisms by which androgen deprivation therapy (ADT) induces NED in advanced Pca and found transmembrane protein 1 (MCTP1) to be abundantly expressed in samples from patients with advanced PCa. Of note, after ADT, MCTP1 through SNAIL promoted EMT, NED and cell migration of PC3 and C4-2 PCa cells [
43].
Therapeutic targeting of SNAIL or ZEB1 may thus prove beneficial in abrogating not only EMT but also NED [
41]. Apart from NFκB-STAT3 [
39], other signaling pathways are involved in acquiring NED characteristic features of PCa cells, such as AMPK/SIRT1-p38MAPK-IL6 [
44] and ERK [
45]. Of note, MEK-ERK signaling is also activated by TGF-β and is critically involved in driving TGF-β-dependent EMT, migration, invasion and metastasis [
46,
47,
48], and TGF-β1-induced downregulation of
CDH1 and upregulation of
SNAI1 [
33]. We are currently carrying out ERK immunoblot analysis of PANC-1 subclones to reveal whether the extent of ERK1/2 activation corresponds with NED marker expression. Interestingly, treatment of MIA PaCa-2 and PANC-1 cells with grape seed proanthocyanidins (GSPs) resulted in decreased phosphorylation of ERK1/2, inactivation of NFκB, reversal of EMT, upregulation of ECAD, downregulation of NCAD and VIM, and reduced cell migration [
49]. It is thus conceivable that GSP-induced inhibition of ERK activation also reduces NED in these cells.
The importance of EMT is also considered in other GI cancers such as coloNEC. Current efforts are therefore underway to reveal an association between EMT and NED in model cell lines of this disease, e.g., the coloNEC-derived cell lines, SS-2 and LCC-18. Here, we have studied a non-GI cancer, namely small cell hypercalcemic ovarian cancer, represented by the cell lines SCCOHT-1 [
9] and BIN-67 [
10,
11]. SCCOHT-1 cells constitutively express NED markers such as NCAM in the original patient tumor and derived mouse xenograft tumors [
9,
12]. Moreover, both cell lines developed biallelic deleterious SMARC A4 gene mutations whereby phenotypic and genetic similarities were observed between SCCOHT and highly malignant childhood-onset atypical teratoid/rhabdoid tumors (AT/RTs) of the central nervous system [
50,
51]. This heterogeneity and plasticity of SCCOHT-1 indicates the potential for transdifferentiation involving EMT and maturation along a NED phenotype [
13]. We found both cell lines to be refractory to stimulation with TGF-β1 for any of the above-mentioned EMT/NED genes but have partially retained sensitivity to BMP-7 as evidenced by induction of ECAD in BIN-67 and SNAIL, VIM, CHGA and SSTR2 in SCCOHT-1 cells.
Since both PANC-1 and MIA PaCa-2 cells are highly invasive and metastatic [
52], the NED phenotype besides the mesenchymal subtype may contribute to this property. In addition, the TGF-β1 effects on various NED markers may have therapeutic significance as it was shown that PANC-1 and MIA PaCa-2 cell lines when subjected to (fractionated) radiation upregulate not only the expression of NED genes [
5] but also induce the synthesis and secretion of TGF-β [
53]. Both radiotherapy and chemotherapy induce TGF-β activity, possibly promoting metastatic progression, and high levels of TGF-β are associated with resistance to anticancer treatments [
54]. Therefore, irradiation-induced secretion of TGF-β1 by tumor cells may account for changes in NED marker expression. Given that the resulting induction of EMT and NED may enhance tumor invasion and metastasis in PDAC, the concomitant application of TGF-β inhibitors between radiation cycles should be considered to prevent an unwanted increase in tumor aggressiveness [
54].