- Open Access
Epiregulin reprograms cancer-associated fibroblasts and facilitates oral squamous cell carcinoma invasion via JAK2-STAT3 pathway
Journal of Experimental & Clinical Cancer Researchvolume 38, Article number: 274 (2019)
Local resident normal fibroblasts (NFs) are the major source of cancer-associated fibroblasts (CAFs), which are distinguishable from NFs by their tumor-supportive properties. However, the mechanism and the effects underlying the transition of NFs to CAFs in oral squamous cell carcinoma (OSCC) remain unclear.
Five pairs of matching primary NFs and CAFs derived from OSCC patients were sent for RNA sequencing. Epiregulin (EREG) expression was analyzed by IHC in fibroblasts from OSCC patients. The role of EREG in the NF-CAF transition and the consequential effects on OSCC progression were examined by upregulation/downregulation of EREG in NFs/CAFs both in vitro and in vivo.
Here, we identified epiregulin (EREG) as the most remarkably upregulated gene in CAFs. High EREG expression in CAFs correlated with higher T stage, deeper invasion and inferior worst pattern of invasion (WPOI) in OSCC patients and predicted shorter overall survival. Overexpression of EREG in NFs activated the CAF phenotype. Mechanistically, the JAK2/STAT3 pathway was enhanced by EREG in parallel with increased IL-6 expression, which could be inhibited by the JAK2 inhibitor AG490. Recombinant IL-6 upregulated the JAK2/STAT3/EREG pathway in a feedback loop. Moreover, EREG-induced CAF activation promoted the epithelial-mesenchymal transition (EMT) necessary for migration and invasion, which was dependent on JAK2/STAT3 signaling and IL-6. In vivo, EREG expression in stroma fibroblasts promoted tumor growth with high stromal α-SMA, phospho-JAK2/STAT3, and IL-6 expression and upregulated EMT in HSC3 cells.
EREG is essential for the NF-CAF transformation needed to induce EMT of tumor cells in a JAK2-STAT3- and IL-6-dependent manner in OSCC.
It is widely accepted that tumor development depends on a complex interaction between malignant cells and their microenvironment . Cancer-associated fibroblasts (CAFs), as a major cellular component of the tumor stroma, coevolve with cancer cells and contribute to cancer cell growth, survival, and metastasis [2, 3]. Often marked by smooth muscle actin (SMA), vimentin and N-cadherin, CAFs acquire a phenotype similar to that of myofibroblasts, which are activated in wound healing and fibrosis, with a morphology and function that differ from those of normal fibroblasts (NFs) . In oral squamous cell carcinoma (OSCC), CAFs have been recognized as a poor prognostic factor, promoting invasion, metastasis, and recurrence of OSCC [5,6,7,8]. There are 5 cellular origins of CAFs, among which NFs are a major source [9, 10]. Despite the overwhelming consensus on the tumor-supportive role of CAFs, the molecular mechanisms underlying the transition from NFs to CAFs remain largely unclear. Further clarification of the mechanism underlying the transition from NFs to CAFs is important for cancer diagnosis and therapy.
Epiregulin (EREG), which is mainly secreted by fibroblasts, is a member of the epidermal growth factor (EGF) family of peptide growth factors and promotes tumor development, migration and invasion [11,12,13,14]. It has been reported that EREG expression is increased in healing wounds [13, 15]. However, whether EREG is involved in the NF-CAF transition remains unknown. As CAFs function similarly to myofibroblasts in wound healing, we hypothesized that EREG might participate in the transition of NFs to the CAF phenotype from the very start of NF activation, thus promoting tumor progression.
This study aimed to determine the role of EREG in the NF-CAF transition in the progression of OSCC and to identify the underlying mechanism. By understanding the NF-CAF transition, we can improve our knowledge of the supportive role of the tumor stroma, thus providing new insights into cancer treatment.
Material & Methods
Primary NFs/CAFs separation, cultivation and RNA-sequence
OSCCs, and matching normal oral biopsies were collected following ethical approval (2017NL-013(KS)) and written informed consent. 11 matching CAFs, and NFs were isolated according to a method previously described  and forwarded for morphologic examination. After appropriate NFs/CAFs were isolated and chosen, 5 pairs of matched NFs and CAFs were sent for library construction, and RNA sequencing was performed.
We performed the RNA-seq with the help of Novel Bioinformatics Co., Ltd. (Shanghai, China). Firstly, total RNA was extracted by Trizol reagent (Invitrogen) separately. The RNA quality was checked by Bioanalyzer 2200 (Aligent) and kept at − 80 °C.The RNA with RIN > 8.0 is right for cDNA library construction. Secondly, the complementary DNA (cDNA) libraries for single-end sequencing were prepared using Ion Total RNA-Seq Kit v2.0 (Life Technologies). The cDNA libraries were then processed for the Proton Sequencing process according to the commercially available protocols. Thirdly, Mapping of single-end reads. The clean reads were then aligned to human genome (version: GRCh38.p1) using the MapSplice program (v2.2.0). Finally, pathway analysis was performed and we also applied EBseq algorithm to filter the differentially expressed genes, the significant analysis and FDR analysis were performed. Besides, we presented gene co-expression Networks to find the relations among different mRNA and LncRNA. Results of LncRNA have been reported in a previous publication . Results of mRNAs with the most significant changes are listed in Additional file 1: Table S1. The data in the list are the mean value of the 5 pairs of samples.
Patients and tissue samples
For the analysis of expression of EREG, 104 consecutive unselected T1-4N0M0 OSCC patients treated in the Department of Oral and Maxillofacial Surgery at Nanjing Stomatology Hospital from 1 January, 2009 to 31 December, 2010 were included into the study. Including criteria: (1) diagnosed with primary OSCC by hematoxylin and eosin staining, and staged as T1-4N0M0 according to the 8th version of American Joint Committee on Cancer (AJCC) TNM staging system  by experienced pathologists from the Department of Pathology at Nanjing Stomatology Hospital; (2) underwent primary surgical treatments in the Department of Oral and Maxillofacial Surgery at Nanjing Stomatology Hospital; (4) paraffin-fixed sample available; (3) with complete clinical-pathological information and follow-up data. Excluding criteria: (1) diagnosed with autoimmune or other malignant diseases and pregnant; (2) had adjuvant treatment before or after surgery. The ethical approval for this study was obtained from the Research Ethics Committee of Nanjing Stomatology Hospital (approval number: 2017NL-013(KS)). The clinical pathologic characteristics of the 104 OSCC patients are listed in Table 1.
Immunohistochemistry staining and evaluation
Immunofluorescence of cryosections or immunohistochemistry of paraffin-embedded tissue was done as previously described . Primary antibodies were rabbit anti-human EREG (Abcam), phospho-JAK2 (Proteintech), phospho-STAT3 (CST), SMA (Abcam).
EREG expression in CAFs was identified on the basis of SMA staining. CAFs were defined as SMA positive, spindled-shaped cells in the stroma (data of SMA staining not shown). EREG staining was assessed only in fibroblasts in tumor stroma by the evaluation of staining intensity and the percentage of EREG positive fibroblasts according to criteria modified from Zhang J et al. (Fig. 2a). Briefly, EREG expression was scored by a combination of staining intensity (score: 0 = no staining, 1 = weak staining, 2 = medium staining, 3 = strong staining) and percentage of EREG positive cells (score 0 = 0%~ 5% positive cells, 1 = 6%~ 33% positive cells, 2 = 34%~ 66% positive cells and 3 = 67%~ 100% positive cells) by a method modified from Zhang J, et al. (Fig. 2a) . EREG staining index was taken as the product of staining intensity and percentage. The staining index was further divided into low and high: score 0–4 (negative to medium) was defined as EREG low, score 6 and 9 (strong) was defined as high EREG expression. Two observers examined the images independently without the knowledge of patients’ information.
Construction of plasmid and siRNA
EREG over expression plasmid was purchased from Genechem Co.,Ltd. (Shanhai, China). The plasmid served as vector was pcDNA3.1+ (Amp). Amplified DNA fragments were inserted into the EcoR I/Xhol sites of pcDNA3.1+ vector. The amplification sequences were: EREG-F: GAATTCatgaccgcggggaggaggatggag; EREG-R: CTCGAGtcagacttgcggcaactctggatc. Empty plasmid was used as control.
The siRNAs were synthesized by Ribobio Co. Ltd. (Guangzhou, China) with following sequences: siNC: 5′-UUCUCCGAACGUGUCACGUdTdT-3′, siEREG-1: 5′-CCAGGAGAGUCCAGUGAUAdTdT-3′, siEREG-2: 5′-CCACCAACCUUUAAGCAAAdTdT-3′, siEREG-3: 5′-UACACUUUGUUAUUGACACUUdTdT-3′.
Plasmid-mediated overexpression of EREG in NFs, and siRNA-mediated interference of EREG in CAFs, were carried out using Lipofectamine® 2000 Transfection Reagent (Thermo Fisher Scientific) according to manufacturer’s instructions. Transient transfection was applied in all in vitro experiments.
Quantitative reverse transcription PCR
Cells were lysed with TRIzol® reagent (Life Technologies), total RNA was extracted following the manufacturer’s instructions, and cDNA synthesis was conducted using High-Capacity cDNA Archive Kit system (Applied Biosystems). Quantitative PCR (q-PCR) was then conducted using iTaq™ Universal SYBR® Green Supermix (BIO-RAD). Comparative 2-△Ct or 2-△△Ct method was used to quantify the relative mRNA expression as indicated.
Western blots were performed using an SDS–PAGE electrophoresis system as described previously , employing rabbit anti-human EREG (Abcam), rabbit anti-human SMA (Abcam), rabbit anti-human E-cadherin (CST), rabbit anti-human N-cadherin (CST), rabbit anti-human vimentin (Proteintech), rabbit anti-human p-JAK2 (CST), rabbit anti-human JAK2 (CST), rabbit anti-human p-STAT3 (CST), mouse anti-human STAT3 (CST), rabbit anti-human p-NF-kB (CST), rabbit anti-human NF-kB (CST), rabbit anti-human p-p38 (CST), rabbit anti-human p38 (CST), rabbit anti-human p-akt (CST), rabbit anti-human akt (CST), rabbit anti-human p-erk1/2 (CST), rabbit anti-human erk1/2 (CST), rabbit anti-human snail (CST), antibodies. The blots were re-probed for rabbit anti-human GAPDH (Bioworld) to control for protein loading and transfer.
Immunofluorescence analyses were carried out as previously described : primary antibodies, goat anti-human EREG (R&D), rabbit anti-human SMA (Abcam) and mouse anti-human CK (Abcam); secondary antibodies, AlexaFluor 647 donkey anti-goat IgG and AlexaFluor 488 goat anti-rabbit 594 goat anti-mouse IgG (Life Technologies).
Conditioned medium (CM) collection
NFs and CAFs were seeded into 12-well plates at identical density. NFs were transfected with control plasmid (pcDNA3.1) or EREG-overexpression plasmid (pcDNA3.1 EREG), and CAFs were transfected with control siRNA (siNC) or EREG siRNA (siEREG). The culture medium were collected 48 h after transfection, centrifuged at 3000RPM for 10 min to remove suspended cells, and stored in − 80 °C until Use.
Cells were transfected with over-expression plasmid or siRNA for EREG, followed with or without treatment of AG490 or rIL6 as indicated. After 48 h, the supernatants were collected for the measurement of IL6 with an ELISA kit (Dakewei, China) according to the manufacturer’s instructions.
Primary NFs and CAFs were transfected with pcDNA3.1-control, pcDNA3.1-EREG siNC or siEREG as indicated. The treated cells were either seeded in 96-well, or their culture medium was used to treat HSC3 cells seeded in 96-well. Cell viabilities were determined at 0, 24, 48, 72 and 96 h after transfection (in case of primary NFs/CAFs) or CM treatment (in case of HSC3 cells).
Colony formation assay
30/ml HSC3 cells were seeded in 12 well plates. Conditioned-medium from NFs + pcDNA3.1/pcDNA3.1-EREG, or CAFs+siNC/siEREG 24 h after transfection were used to treat HSC3 cells. The CM was changed every 3 days. After 1 to 2 weeks of cultivation, when visible colonies were formed, the medium was removed. Cells were 5% formalin-fixed after washed with PBS for 3 times, then stained with crystal violet.
Wound healing assays
HSC3 cells were seeded in 12-well plates until 100% fusion. After overnight starvation with serum-free DF12 medium, a scratch was made using a micropipette tip and cells were washed to remove detached cells and debris. Photographs of the same area of the wound were taken at 0 and 18 h for measuring the closure of the wound after the treatment of the indicated serum-free CM from NFs/CAFs.
Invasion and migration assays
The invasion and migration ability of HSC3 cells were tested in Transwell (boyden chamber) with an 8-um pore size of the polycarbonate membranes. Primary NFs and CAFs were placed in the lower chamber at a concentration of 3 × 10^4 cells/ml (0.6 ml/well). The NFs were then transfected with or without EREG-overexpression plasmids, and CAFs transfected with or without siEREG RNA. AG490 and recombinant human IL6 (Abcam) were added as indicated. After 5 h of incubation, the transfection mix was replaced by 0.5% FBS. Simultaneously, HSC3 cells, which had been previously starved overnight, were re-suspended with serum-free DF12 culture medium at a concentration of 10^6 cells/ml and then added to the upper chamber (200ul/well). For invasion assay, 40ul matrigel mix (matrigel: DMEM = 1: 2) was seeded at the bottom of Transwell chamber 1 h prior to the re-suspension of HSC3 cells to let the gel solidify. The number of migrated and invaded cells were qualified by blinded counting of migrated/invaded cells on the lower surface of the membrane, with five fields per chamber, at 24 h and 48 h, respectively.
3D culture invasion model and evaluation of carcinoma cell invasion
3D culture invasion model was constructed and based on a method previously described by Igarashi et al. [21, 22]. Matrices containing fibroblasts were prepared by the addition of 1.4 parts of DF12 and 0.6 parts of FBS to 1 part of collagen type I solution (Thermo Fisher, 3 mg/ml), and neutralized with 1 M NaOH solution. Fibroblasts were then embedded into this solution to the final cell density of 5*10^5/ml, and seeded on top with HSC3 at the density of 4.5*10^5/cm2 in triplicates. After 5 days of cultivation, 3D constructs were harvested, formalin-fixed, and paraffin-embedded, and sent for slide generation and HE staining.
To investigate the role of EREG in fibroblasts in vivo, stable EREG over-expression in NFs and EREG knockdown in CAFs were realized by lentivirus vector as previously decribed [17, 23]. Then 106 conditional NFs/CAFs were subcutaneously co-injected with 106 HSC3 in rear flank of four-week old male BALB/c-nu/nu T cells-deficient mice (Cavens, Changzhou, China). Tumor volumes (mm3) were measured at various time points until termination of the experiments. Tumors were collected for further experiments. Sixteen mice were used to construct xenografts with HSC3 and fibroblasts (four in each group). At the end of the study, mice were sacrificed and xenografted tumors were harvested. The experiments were executed in compliance with institutional guidelines and regulations.
All images represent at least three independent experiments. The results are presented as the mean ± SEM of at least three independent experiments. Statistical analysis was performed with SPSS® v 19 (SPSS, Chicago, IL, USA). Survival of patients was plotted using the Kaplan–Meier method and analyzed using the log-rank test, with patients censored at last follow-up. Kruskal–Wallis and Mann– Whitney U- and t-tests were used to compare groups, as appropriate. P < 0.05 was considered statistically significant for all tests (*p < 0.05, **p < 0.01, **p < 0.001, ***p < 0.0001).
EREG expression is higher in CAFs than in matched NFs
As an abnormal activated form of NFs, CAFs differ from NFs with distinct phenotypes and functions. Exploration of the regulatory molecules that drive the transition from NFs to CAFs is urgently needed. Therefore, we carried out genomic analysis to identify differentially expressed genes between matching NFs and CAFs obtained from 5 OSCC patients as previously reported . A series of differentially expressed genes were identified, among which Ereg was one of the genes with the most significant differential expression (Fig. 1a, Additional file 1: Table S1. The data are the mean value of the 5 matched samples).
To validate the results of the genomic studies, we then examined Ereg expression at the gene level in another 5 matching NFs and CAFs by quantitative PCR and found that Ereg was significantly higher in CAFs than the matching NFs (Fig. 1b). Furthermore, immunoblotting assays and immunofluorescence results also showed that EREG was significantly higher in CAFs than in matched NFs at the protein level, indicating that EREG may be involved in the transition from NFs to CAFs (Fig. 1c). Furthermore, the EREG expression profiles in normal mucosa tissue and tumor tissue were assessed by immunofluorescence staining in frozen sections. The results demonstrated that EREG expression was also significantly higher in fibroblasts from the cancer stroma (CAFs) than in fibroblasts from the stroma around normal mucosa (NFs; Fig. 1e). These results all proved that EREG expression was much higher in CAFs than in matched NFs.
High EREG expression in fibroblasts is related to poor prognosis in OSCC patients
It is well recognized that CAFs can promote the progression of tumors. However, the clinical outcomes of EREG expression in fibroblasts remain unclear. Considering the high level of EREG in CAFs, we first analyzed the correlation between EREG expression and the clinicopathological characteristics of patients with OSCC (n = 104). Chi-square analysis revealed that high EREG expression in CAFs was closely related to higher T staging (p = 0.0041), deeper depth of invasion (DOI, p < 0.0001) and inferior worst pattern of invasion (WPOI, p < 0.0001) (Table 1). Then, the prognostic value of EREG expression in fibroblasts was assessed. Kaplan–Meier analysis showed that high EREG expression in CAFs was significantly correlated with inferior overall survival (OS), and disease-free survival (DFS) (Fig. 2b, p = 0.0398, and p = 0.0314). The results demonstrated that high EREG expression in fibroblasts correlated with poor prognosis in OSCC patients.
EREG promotes the transition from NFs to CAFs
In addition, we ascertained the role of EREG during the NF-CAF transition. First, we modified EREG expression in NFs and CAFs separately and evaluated changes in the cell phenotypes. High N-cadherin, vimentin, and SMA expression was used to mark the CAF phenotype. Three siRNAs (siEREG-1, siEREG-2, and siEREG-3) were designed, and their interference rates were 64.9, 58.3, and 48.9%, respectively. Then, siRNA-1, which had the highest interference rate, was chosen for subsequent experiments (hereafter referred to as siEREG; Additional file 2: Figure S1B). Plasmid overexpression of EREG in NFs was also successful (Additional file 2: Figure S1A). In addition to the confirmation of EREG expression at the gene level, we also confirmed the overexpression and interference of EREG by ELISA at the protein level (Additional file 2: Figure S1C). Overexpression of EREG in NFs significantly augmented the expression of N-cadherin, vimentin, and SMA, while siRNA modulation of EREG levels in CAFs resulted in reduced N-cadherin, vimentin, and SMA expression (Fig. 3a&b), indicating that EREG was involved in the transition from NFs to CAFs. However, the specific mechanism underlying EREG participation in the NF-CAF transition remains to be further elucidated.
EREG regulates the NF-CAF transition via the JAK2-STAT3 pathway
To ascertain the pathway of the EREG-regulated NF-CAF transition, major pathways, including the JAK2-STAT3 pathway, NF-kB pathway, p38-MAPK pathway, PI3K/AKT pathway and erk1/2 pathway, were tested in unmodified NFs and EREG-overexpressing NFs. The results showed that both the levels of phospho-JAK2 and phospho-STAT3 and the ratios of phospho-JAK2/JAK2 and phospho-STAT3/STAT3 were significantly increased by the overexpression of EREG in NFs, indicating that the JAK2-STAT3 pathway was the most activated pathway after EREG overexpression (Fig. 3c-e). Although a close relationship between the erk1/2 pathway and EREG has been reported , our results demonstrated that erk1/2 did not play a significant role in the EREG-mediated NF-CAF transition (Fig. 3c-e).
To further confirm the role of the JAK2-STAT3 pathway in the EREG-mediated NF-CAF transition, we treated EREG-overexpressing NFs with AG490, an inhibitor of the JAK2-STAT3 pathway, and found that AG490 could reverse EREG-mediated CAF marker augmentation (Fig. 3f-h). Moreover, AG490 reduced EREG-mediated CAF activation (Additional file 3: Figure S2A). IL-6 is an important downstream factor of the JAK2-STAT3 pathway . To further confirm the role of the JAK2-STAT3 pathway, we added recombinant human IL-6 (rIL-6) into the culture medium to construct a rescue assay. The results showed that while siEREG significantly lowered the expression of N-cadherin, vimentin and SMA, rIL-6 upregulated their expression (Fig. 3f-h). IL-6 levels were also assessed by ELISA. The results demonstrated that the IL-6 level was increased after EREG overexpression and was downregulated after EREG interference (Fig. 3i). In clinical OSCC samples, we also found that high EREG expression in CAFs was correlated with higher phospho-JAK2 and phospho-STAT3 expression (Fig. 3j-k). These results all demonstrate that EREG activates NFs in a JAK2-STAT3-dependent way. In addition to the change in protein markers, in OSCC, it has been reported that CAFs promote tumor migration  and invasion [7, 26] and are related to poor prognosis . However, the effects of the EREG-mediated NF-CAF transition during OSCC progression are still unknown.
CAFs-derived EREG promotes the invasion and migration of OSCC cells
Then we investigated whether EREG expression in fibroblasts affected their tumor supportive functions. Ereg-expression level was assessed in 4 human OSCC cell lines and one normal epithelium cell line HACAT by qPCR. All 4 OSCC cell lines showed higher Ereg expression level compared with HACAT, and HSC3 was chosen for further experiments as it expressed the lowest level of Ereg among OSCC cell lines (Additional file 2: Figure S1D). Then, CM obtained from NFs/CAFs were added in to HSC3, and we found that compared with CM obtained from NFs, CM obtained from CAFs significantly improved the cell viabilities and promoted proliferation of HSC3 (Fig. 4a&b). However, modulation of EREG-expression in NFs/CAFs (transfected with plasmid or siRNA as indicated) did not significantly affect their supporting role of HSC3 viabilities and proliferation (Fig. 4a&b).
In our retrospective study of immunohistochemistry EREG staining in CAFs of OSCC, we already found that high EREG expression in CAFs of OSCC was correlated to a deeper invasion depth and inferior WPOI (Table 1). To confirm these results in vitro and to study the EREG-mediated NF activation in terms of their invasion-supportive abilities, in transwell migration assays, we found that CAFs with EREG interfered had an attenuated effect in inducing HSC3 migration, while up-regulation of EREG in NFs augmented their migration inducing ability (Fig. 4c&e). Moreover, the inhibition of JAK2-STAT3 pathway by AG490 could attenuate the migration-supportive role of EREG-high NFs, while rIL6 could rescue the migration-supportive role of EREG-low CAFs. The invasion abilities of HSC3 were assessed by transwell invasion assays. EREG-overexpressed NF increased HSC3 cell invasion more than two-folds than NF-control, while AG490 inhibited this invasion-supportive ability. Invasion-supportive ability was attenuated in CAF-siEREG, but could be rescued by rIL6 (Fig. 4c&f). Again, AG490 treatment attenuated the EREG-mediated migration and invasion-supportive abilities of CAFs (Additional file 3: Figure S2B).
To further confirm the invasion-inducing ability of EREG, we constructed a 3D invasion model with conditional NFs/CAFs, and found that augmented EREG expression in NFs led to deeper invasion. CAFs supported HSC3 invasion most significantly, with deeper invasion and more diffused spreading. This ability was attenuated after EREG interference, but could be restored by rIL6 (Fig. 4d&g). Meanwhile, overexpression of EREG in NFs significantly increased their invasion-supportive ability, but could be interfered after inhibition of JAK2-STAT3 pathway through AG490 (Fig. 4d&g).
To further reveal the mechanism underlying the EREG-mediated invasion-supportive role of fibroblasts, we analyzed the epithelial-mesenchymal transition (EMT) markers of HSC3 cells after treatment with CM from conditional NFs/CAFs. The results showed that compared with NFs, CAFs could significantly promote EMT in HSC3 cells in an EREG- and JAK2-STAT3-dependent way (Fig. 4h&i).
Together, these results indicate that the NF-CAF transition, which is induced by EREG, affects the migration and invasion of OSCC cells through EMT in a JAK2-STAT3-dependent way.
EREG expression in fibroblasts promotes OSCC progression with stronger invasive character in vivo
To further examine the tumorigenic role of EREG in vivo, we established a xenograft model in which 1*106 HSC3 cells were subcutaneously coinjected with 106 conditioned fibroblasts into the flank of nude mice, which were allowed to develop measurable tumors. One of the most important functional characteristics of CAFs is their tumor-supportive nature in vivo. Indeed, HSC3 cells coinjected with CAFs demonstrated rapid in vivo tumor growth, but this tumor-supportive effect could be attenuated by EREG knockdown (Fig. 5a-c). On the other hand, NFs demonstrated an in vivo tumor-supportive role after EREG overexpression (Fig. 5a-c). EREG knockdown and EREG overexpression were confirmed by qPCR and IHC (Fig. 5d & h). The EREG-mediated NF-CAF transition was confirmed in vivo, with IHC showing a positive correlation between EREG expression and SMA expression in fibroblasts (Fig. 5h). Moreover, EREG modulated the NF-CAF transition through JAK2-STAT3 pathway activation (Fig. 5e-h). In addition, WB confirmed that EREG regulated the EMT-promoting abilities of fibroblasts in vivo (Fig. 5i&j). Thus, these data provide novel insights into the involvement of EREG in the EREG-JAK2-STAT3-IL-6 axis in fibroblasts, which supports OSCC invasion through induction of EMT in the tumor microenvironment.
Available evidence suggests that CAFs, as an irreversible active form of NFs, promote the progression of OSCC . However, the molecular mechanisms underlying the NF-CAF transition remain largely unclear.
The activation of NFs occurs through a multistep process. During wound healing, spindle-shaped quiescent resting fibroblasts (NFs) are reversibly activated to facilitate repair and regeneration. Normal activated fibroblasts (NAFs) gain smooth muscle actin (SMA) and vimentin and became stellate in shape. Activated fibroblasts can be further activated until there is “no way back” and they finally transit to CAFs. Such CAFs gain enhanced proliferative properties and are a functionally diverse population, adding to the dynamic complexity of the tumor microenvironment milieu .
To identify key participants in NF-CAF reprogramming, 5 pairs of matching NFs and CAFs were sent for genomic analysis. We reviewed genes with the most significant changes and identified epiregulin (EREG) as a strong candidate. From its initial discovery, EREG has been closely linked with fibroblasts . It has been reported as the most significantly altered gene between NFs and CAFs in other studies [6, 13]. It is overexpressed not only in CAFs but also in the stroma during wound healing [13, 15], indicating that it could take part in NF activation from the very start, activating even the most quiescent NFs.
In a cohort consisting of 104 OSCC patients, we found that high levels of EREG expression in fibroblasts were significantly correlated with poorer OS and DFS, deeper DOI and higher levels of WPOI in OSCC patients. DOI and WPOI are prominent indicators for tumor invasion and are closely related to OSCC prognosis [30, 31]. These results indicate that EREG secreted by CAFs is related to inferior OSCC prognosis and promotes OSCC invasion.
Next, we identified the causal relationship between EREG expression in fibroblasts and their CAF-like phenotype. Overexpression of EREG led to the activation of NFs, resulting in augmented expression of CAF markers (SMA, N-cadherin, and vimentin) and upregulated cell viability. Moreover, CAFs, a naturally irreversible activated form of NFs, were deprived of their CAF phenotype by EREG knockdown.
Through years of study, the tumor-promoting functions of the JAK2-STAT3 pathway have been well recognized in various malignancies . It has been reported that JAK-STAT can promote the proliferation  and activation of fibroblasts . IL-1 induces the JAK-STAT pathway to shape CAF heterogeneity in pancreatic ductal adenocarcinoma . CAFs promote bladder cancer EMT via paracrine IL-6 . However, further evidence is needed regarding JAK-STAT-mediated fibroblast activation in oral cancer. Our study revealed that, in the abovementioned EREG-mediated NF activation, the JAK2-STAT3 pathway is the most significantly activated pathway. EREG could activate fibroblasts in a JAK2-STAT3-dependent way.
In addition to protein marker alteration, a widely recognized functional characteristic of CAFs is their supportive role in tumor progression. Existing evidence supports the idea that EREG is involved in carcinogenesis [12, 13] and metastasis  and correlates with poor prognosis [38, 39]. Fibroblast-derived EREG promotes epithelial cell proliferation in cholesteatoma . In our study, even though the EREG level was positively correlated with fibroblast proliferation ability, colony formation and CCK8 assays revealed that modulating EREG in NFs/CAFs did not significantly alter their effects on oral cancer cell proliferation. On the other hand, wound-healing assays, transwell assays and 3D invasion models demonstrated that altered EREG expression in fibroblasts significantly influences their supportive role in cancer cell migration and invasion, with more significant effects on invasion. High EREG levels in fibroblasts could promote their invasion-supportive abilities by encouraging OSCC EMT in a JAK2-STAT3 pathway-dependent way.
Last, we confirmed our results in vivo. In a xenograft tumor model, EREG expression in fibroblasts had a tumor-supportive function in an EREG-related pattern. Phospho-JAK2 and phospho-STAT3 levels were related to EREG expression.
In recent years, CAFs have become a novel target in cancer treatment [40,41,42]. The identification of EREG as a pivotal participant in the NF-CAF transition highlights critical pathways that could be targeted for novel therapeutic interventions.
In this study, we showed that EREG is an important protein that could reprogram NFs into CAFs via the JAK2-STAT3 pathway. The EREG expression level was related to the invasion-supportive role of fibroblasts and predicted a poor clinical outcome of OSCC patients. Thus, EREG overexpression in fibroblasts is implicated in the NF-CAF transition and in OSCC progression (Fig. 5k).
American Joint Committee on Cancer
depth of invasion
epidermal growth factor
epidermal growth factor receptor
oral squamous cell carcinoma
worst pattern of invasion
Mueller MM, Fusenig NE. Friends or foes - bipolar effects of the tumour stroma in cancer. Nat Rev Cancer. 2004;4(11):839–49.
Hanahan D, Coussens LM. Accessories to the crime: functions of cells recruited to the tumor microenvironment. Cancer Cell. 2012;21(3):309–22.
Xu J, Lu Y, Qiu S, Chen ZN, Fan Z. A novel role of EMMPRIN/CD147 in transformation of quiescent fibroblasts to cancer-associated fibroblasts by breast cancer cells. Cancer Lett. 2013;335(2):380–6.
Marsh T, Pietras K, McAllister SS. Fibroblasts as architects of cancer pathogenesis. Biochim Biophys Acta. 2013;1832(7):1070–8.
Li H, Zhang J, Chen SW, et al. Cancer-associated fibroblasts provide a suitable microenvironment for tumor development and progression in oral tongue squamous cancer. J Transl Med. 2015;13:198.
Costea DE, Hills A, Osman AH, et al. Identification of two distinct carcinoma-associated fibroblast subtypes with differential tumor-promoting abilities in oral squamous cell carcinoma. Cancer Res. 2013;73(13):3888–901.
Elmusrati AA, Pilborough AE, Khurram SA, Lambert DW. Cancer-associated fibroblasts promote bone invasion in oral squamous cell carcinoma. Br J Cancer. 2017;117(6):867–75.
Kanemaru A, Yamamoto K, Kawaguchi M, et al. Deregulated matriptase activity in oral squamous cell carcinoma promotes the infiltration of cancer-associated fibroblasts by paracrine activation of protease-activated receptor 2. International journal of cancer Journal international du cancer. 2017;140(1):130–41.
Polanska UM, Acar A, Orimo A. Experimental generation of carcinoma-associated fibroblasts (CAFs) from human mammary fibroblasts. Journal of visualized experiments : JoVE. 2011(56):e3201.
Mitra AK, Zillhardt M, Hua Y, et al. MicroRNAs reprogram normal fibroblasts into cancer-associated fibroblasts in ovarian cancer. Cancer discovery. 2012;2(12):1100–8.
Yoshikawa M, Kojima H, Yaguchi Y, Okada N, Saito H, Moriyama H. Cholesteatoma fibroblasts promote epithelial cell proliferation through overexpression of epiregulin. PLoS One. 2013;8(6):e66725.
Bauer AK, Velmurugan K, Xiong KN, Alexander CM, Xiong J, Brooks R. Epiregulin is required for lung tumor promotion in a murine two-stage carcinogenesis model. Mol Carcinog. 2016.
Neufert C, Becker C, Tureci O, et al. Tumor fibroblast-derived epiregulin promotes growth of colitis-associated neoplasms through ERK. J Clin Invest. 2013;123(4):1428–43.
Riese DJ 2nd, Cullum RL. Epiregulin: roles in normal physiology and cancer. Semin Cell Dev Biol. 2014;28:49–56.
Roy S, Khanna S, Rink C, Biswas S, Sen CK. Characterization of the acute temporal changes in excisional murine cutaneous wound inflammation by screening of the wound-edge transcriptome. Physiol Genomics. 2008;34(2):162–84.
Liu Y, Hu T, Shen J, et al. Separation, cultivation and biological characteristics of oral carcinoma-associated fibroblasts. Oral Dis. 2006;12(4):375–80.
Ding L, Ren J, Zhang D, et al. A novel stromal lncRNA signature reprograms fibroblasts to promote the growth of oral squamous cell carcinoma via LncRNA-CAF/interleukin-33. Carcinogenesis. 2018.
Tirelli G, Gatto A, Boscolo Nata F, et al. Prognosis of oral cancer: a comparison of the staging systems given in the 7th and 8th editions of the American joint committee on Cancer staging manual. Br J Oral Maxillofac Surg. 2018;56(1):8–13.
Liang D, Xiao-Feng H, Guan-Jun D, et al. Activated STING enhances Tregs infiltration in the HPV-related carcinogenesis of tongue squamous cells via the c-Jun/CCL22 signal. Biochim Biophys Acta. 2015;1852(11):2494–503.
Underwood TJ, Hayden AL, Derouet M, et al. Cancer-associated fibroblasts predict poor outcome and promote periostin-dependent invasion in oesophageal adenocarcinoma. J Pathol. 2015;235(3):466–77.
Igarashi M, Irwin CR, Locke M, Mackenzie IC. Construction of large area organotypical cultures of oral mucosa and skin. Journal of oral pathology & medicine : official publication of the International Association of Oral Pathologists and the American Academy of Oral Pathology. 2003;32(7):422–30.
Yamada R, Kitajima K, Arai K, Igarashi M. Cytokeratin expression of engrafted three-dimensional culture tissues using epithelial cells derived from porcine periodontal ligaments. Journal of oral pathology & medicine : official publication of the International Association of Oral Pathologists and the American Academy of Oral Pathology. 2014;43(8):637–45.
Wang J, Guan X, Guo F, et al. miR-30e reciprocally regulates the differentiation of adipocytes and osteoblasts by directly targeting low-density lipoprotein receptor-related protein 6. Cell Death Dis. 2013;4:e845.
Rokavec M, Oner MG, Li H, et al. IL-6R/STAT3/miR-34a feedback loop promotes EMT-mediated colorectal cancer invasion and metastasis. J Clin Invest. 2014;124(4):1853–67.
Zhou B, Chen WL, Wang YY, et al. A role for cancer-associated fibroblasts in inducing the epithelial-to-mesenchymal transition in human tongue squamous cell carcinoma. Journal of oral pathology & medicine : official publication of the International Association of Oral Pathologists and the American Academy of Oral Pathology. 2014;43(8):585–92.
Hassona Y, Cirillo N, Heesom K, Parkinson EK, Prime SS. Senescent cancer-associated fibroblasts secrete active MMP-2 that promotes keratinocyte dis-cohesion and invasion. Br J Cancer. 2014;111(6):1230–7.
Ding L, Zhang Z, Shang D, et al. Alpha-smooth muscle actin-positive myofibroblasts, in association with epithelial-mesenchymal transition and lymphogenesis, is a critical prognostic parameter in patients with oral tongue squamous cell carcinoma. Journal of oral pathology & medicine : official publication of the International Association of Oral Pathologists and the American Academy of Oral Pathology. 2014;43(5):335–43.
Kalluri R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 2016;16(9):582–98.
Toyoda H, Komurasaki T, Uchida D, et al. Epiregulin. A novel epidermal growth factor with mitogenic activity for rat primary hepatocytes. J Biol Chem. 1995;270(13):7495–500.
Almangush A, Bello IO, Keski-Santti H, et al. Depth of invasion, tumor budding, and worst pattern of invasion: prognostic indicators in early-stage oral tongue cancer. Head Neck. 2014;36(6):811–8.
Almangush A, Bello IO, Coletta RD, et al. For early-stage oral tongue cancer, depth of invasion and worst pattern of invasion are the strongest pathological predictors for locoregional recurrence and mortality. Virchows Arch. 2015;467(1):39–46.
Sansone P, Bromberg J. Targeting the interleukin-6/Jak/stat pathway in human malignancies. Journal of clinical oncology : official journal of the American Society of Clinical Oncology. 2012;30(9):1005–14.
Liu LD, Dong CH, Shi HJ, et al. A novel type II membrane receptor up-regulated by IFN-alpha in fibroblasts functions in cell proliferation through the JAK-STAT signalling pathway. Cell Prolif. 2006;39(2):93–103.
Hendrayani SF, Al-Khalaf HH, Aboussekhra A. The cytokine IL-6 reactivates breast stromal fibroblasts through transcription factor STAT3-dependent up-regulation of the RNA-binding protein AUF1. J Biol Chem. 2014;289(45):30962–76.
Biffi G, Oni TE, Spielman B, et al. IL1-induced JAK/STAT signaling is antagonized by TGFbeta to shape CAF heterogeneity in pancreatic ductal adenocarcinoma. Cancer discovery. 2019;9(2):282–301.
Goulet CR, Champagne A, Bernard G, et al. Cancer-associated fibroblasts induce epithelial-mesenchymal transition of bladder cancer cells through paracrine IL-6 signalling. BMC Cancer. 2019;19(1):137.
Liu S, Ye D, Xu D, et al. Autocrine epiregulin activates EGFR pathway for lung metastasis via EMT in salivary adenoid cystic carcinoma. Oncotarget. 2016;7(18):25251–63.
Jacobs B, De Roock W, Piessevaux H, et al. Amphiregulin and epiregulin mRNA expression in primary tumors predicts outcome in metastatic colorectal cancer treated with cetuximab. Journal of clinical oncology : official journal of the American Society of Clinical Oncology. 2009;27(30):5068–74.
Stahler A, Heinemann V, Giessen-Jung C, et al. Influence of mRNA expression of epiregulin and amphiregulin on outcome of patients with metastatic colorectal cancer treated with 5-FU/LV plus irinotecan or irinotecan plus oxaliplatin as first-line treatment (FIRE 1-trial). International journal of cancer Journal international du cancer. 2016;138(3):739–46.
Mertens JC, Fingas CD, Christensen JD, et al. Therapeutic effects of deleting cancer-associated fibroblasts in cholangiocarcinoma. Cancer Res. 2013;73(2):897–907.
Chan JSK, Sng MK, Teo ZQ, Chong HC, Twang JS, Tan NS. Targeting nuclear receptors in cancer-associated fibroblasts as concurrent therapy to inhibit development of chemoresistant tumors. Oncogene. 2017.
Zhang J, Liu J. Tumor stroma as targets for cancer therapy. Pharmacol Ther. 2013;137(2):200–15.
We thank Dr. Rie Yamada, Dr. Kayoko Kitajima, Dr. Kyoko Arai and Professor Masaru Igarashi (Department of Endodontics, The Nippon Dental University School of Life Dentistry at Niigata, Japan) for their generous help and precious advice in the construction of 3D invasion model.
Availability of data and supporting materials
Please contact author for data requests.
This work was supported by the National Natural Science Foundation of China (No. 81772880, 81702680), Nanjing Medical Science and Technique Development Foundation(No.YKK16164, QRX17083), Jiangsu Provincial Key Medical Discipline(since 2017), Nanjing Municipal Key Medical Laboratory Constructional Project Funding (Since 2016), Center of Nanjing Clinical Medicine of tumor project (Since 2014).
Ethics approval and consent to participate
The ethical approval for this study was obtained from the Research Ethics Committee of Nanjing Stomatology Hospital, number: 2017NL-013(KS).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. The levels and fold-changes of differentially expressed mRNAs. Table S2. Primer sequences. (DOC 133 kb)
Figure S1. A, Relative Ereg mRNA levels in NFs after transfection with the control plasmid pcDNA3.1 or the EREG-overexpressing plasmid pcDNA3.1-EREG, showing successful EREG overexpression in NFs. B, Relative EREG mRNA levels in NFs after transfection with control siRNA (siNC) or EREG interference siRNA. EREG siRNA-1 (hereby referred to as siEREG) was chosen for further experiments. C, ELISA showing increased EREG expression of NFs after EREG overexpression and decreased EREG expression of CAFs after EREG interference. D, Ereg mRNA expression was analyzed in 1 normal epithelial cell line (HACAT) and 4 different oral squamous cell carcinoma cell lines (OSCC3, SCC4, HSC3, and Cal27) by quantitative RT-PCR. HSC3 cells, having the lowest Ereg mRNA level, were chosen for further research. (TIF 3286 kb)
Figure S2. AG490 reduces EREG-mediated CAF activation. AG490 was added to the CM of CAFs transfected with pcDNA3.1-EREG, and treatment was carried out for 48 h. Then, CAFs were sent for protein extraction for WB analysis (A) or seeded for transwell assays (B). The results revealed that AG490 reduced EREG-mediated CAF activation (A) and pro-migration and pro-invasion abilities (B). (TIF 1444 kb)