Abstract
Transmembrane protease serine 4 (TMPRSS4) is upregulated in various kinds of human cancers, including pancreatic cancer. However, its biological function in pancreatic ductal adenocarcinoma (PDAC) remains unclear. In the current study, real-time qPCR, immunohistochemical staining, Western blotting, and database (Cancer Genome Atlas and Gene Expression) analysis revealed remarkable overexpression of TMPRSS4 in PDAC tissue as compared to non-tumor tissue. The TMPRSS4 overexpression was associated with poor prognosis of PDAC patients. Moreover, multivariate analysis revealed that TMPRSS4 serves as an independent risk factor in PDAC. We performed gain-and loss-of-function analysis and found that TMPRSS4 promotes cellular proliferation and inhibits apoptosis of PDAC cells both in vitro and in vivo. Furthermore, we showed that TMPRSS4 might promote cell proliferation and inhibit apoptosis through activating ERK1/2 signaling pathway in pancreatic cancer cells. These findings were validated by using ERK1/2 phosphorylation inhibitor SCH772984 both in vitro and in vivo. Taken together, this study suggests that TMPRSS4 is a proto-oncogene, which promotes initiation and progression of PDAC by controlling cell proliferation and apoptosis. Our findings indicate that TMPRSS4 could be a promising prognostic biomarker and a therapeutic target for the treatment of pancreatic cancer.
Introduction
Pancreatic ductal adenocarcinoma (PDAC) is the most common type of pancreatic cancer (PC) and the fourth most lethal malignancy. It has a poor prognosis with the five-year survival rate of only 10% (). There has not been much progress in the early prognosis and treatment in pancreatic cancer as compared to other malignancies, such as breast and colorectal cancers (). Clinical management of pancreatic cancer poses several challenges such as prevention, early diagnosis at a curable stage, and effective biomarkers for its precise detection (). Therefore, there is an unmet need for the development of sensitive, as well as specific therapeutic targets and an in-depth understanding of the regulatory mechanisms that mediate the development of pancreatic cancer.
Transmembrane protease serine 4 (TMPRSS4), previously referred to as TMPRSS3, is a protein-coding gene located in the human chromosome 11q23.3 (). TMPRSS4 belongs to the type II transmembrane serine protease family. This family of proteins is anchored to the plasma membrane by the transmembrane domain. They play crucial roles in the execution of diverse cellular functions (, ). Cancer cells have distinct characteristics of abnormal cell proliferation, apoptosis inhibition, active migration, and invasion. Recent studies suggest that TMPRSS4 facilitates cancerous tumor growth and metastasis in multiple malignancies. A positive correlation between the overexpression of TMPRSS4 and poor prognosis was demonstrated in a variety of human cancers (–). In lung cancer, TMPRSS4 promotes cell growth, and imparts drug resistance to chemotherapy (). In pancreatic cancer TMPRSS4 is overexpressed and involved in metastasis formation and tumor invasion ().
Additionally, TMPRSS4 enhances cell invasion by activating NF-kB/MMP-9 signaling pathway in human gastric cancer cells (). Also, TMPRSS4 mediates an epithelial to mesenchymal transition phenotype by the activation of distinct signaling pathways such as FAK, ERK1/2, Akt, Src, and Rac1 (). Previous studies indicate that TMPRSS4 might play a crucial role in the development and progression of tumors. However, the underlying molecular mechanisms of TMPRSS4 in pancreatic cancer development and progression remain unclear.
In this study, we unraveled that overexpression of TMPRSS4 is associated with poor prognosis in pancreatic cancer. We also found that TMPRSS4 mediates cell proliferation and inhibits apoptosis in pancreatic cancer cells. These findings suggest that TMPRSS4 is involved in the development and progression of pancreatic cancer. Further in-depth investigation of TMPRSS4 might lead to its development as an efficient prognostic biomarker and therapeutic target for pancreatic cancer.
Materials and Methods
GEO and TCGA Data Mining and Analysis
We probed the GEO database (https://www.ncbi.nlm.nih.gov/geo/) to explore the expression of TMPRSS4 in PDAC and non-cancerous tissues. We found overexpression of TMPRSS4 in the datasets of GSE62165 (118 PDAC samples and 13 non-tumor samples) (), GSE15471 (36 paired PDAC tissues) (, ), GSE62452 (69 pancreatic tumors and 61 adjacent non-tumors) (), and GSE16515 (36 pancreatic tumor samples and 16 non-tumor samples) (). We validated the aberrant upregulation of TMPRSS4 by integrating analysis of the TCGA (The Cancer Genome Atlas) data and GTEx data (https://xenabrowser.net/datapages/) (). The association between the TMPRSS4 expression and the prognosis of PDAC patients was investigated in the 149 (all 177 pancreatic cancer) and 36 PDAC cases from TCGA and GSE16515, respectively.
Tissue Specimens
The current study is approved by the Southwest Hospital Ethics Committee. Tissue specimens were procured from the Southwest Hospital, Third Military Medical University in Chongqing, China. We performed immunohistochemical analysis in a tissue microarray containing 97 cases of histopathologically diagnosed PDAC and 12 normal pancreatic tissues. Additionally, 14 fresh-frozen PDAC tissue samples and matched adjacent non-tumor tissues were taken for RT-qPCR and Western blot analysis.
Cell Lines and Cell Culture
Seven pancreatic cancer cell lines, AsPC-1, BxPC-3, Capan-1, CFPAC-1, PANC-1, SW1990 (Shanghai Institute of Biochemistry and Cell Biology, China), and Hs766t (ATCC, Manassas, VA, USA), and an immortalized human pancreatic ductal epithelial cell lines: HPDE6-C7 (BeNa Culture Collection, China), were used in this study. These cell lines were maintained in the RPMI 1640 medium (Gibco Invitrogen, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (FBS) (Gibco Invitrogen, Grand Island, NY, USA) and 1% antibiotics, in a humidified atmosphere with 5% carbon dioxide (CO2) and 95% air at 37°C.
RNA Extraction and Real-Time PCR
The total RNA was isolated from the PDAC tissue and PDAC cells with the Trizol reagent (Invitrogen, Camarillo, CA, USA) as per the manufacturer's instruction. Furthermore, RNA concentration was evaluated, and stored at −80°C refrigerator to avoid degradation. Prime Script RT reagent Kit (TaKaRa) and SYBR Premix Ex Taq II (TaKaRa) was used to detect the TMPRSS4 mRNA expression and β-actin was used as an internal control. The 2−ΔΔCT method was used to calculate the relative expression of TMPRSS4. The primers used for the TMPRSS4 was (): Forward: 5′-CCGATGTGTTCAACTGGAAG−3′ and Reverse: 5′-CCCATCCAATGATCCAGAGT-3′.
Western Blot
Western blot was performed as previously described () using antibodies included those targeting TMPRSS4 (1:800; Proteintech, USA), ERK1/2 (1:1,000; Cell Signaling Technology, USA), phospho-ERK1/2 (p-ERK1/2) (1:1,000; Cell Signaling Technology, USA), Bax (1:1,000; Cell Signaling Technology, USA), Bcl-2 (1:1,000; Cell Signaling Technology, USA), cleaved caspase-3 (1:1,000; Cell Signaling Technology, USA), and β-actin (1:5,000; Cell Signaling Technology, USA). The horseradish peroxidase-conjugated antibody was used as the secondary antibody (anti-mouse 1:5,000; Cell Signaling Technology, anti-rabbit; 1:5,000; Cell Signaling Technology). Protein levels were normalized against the endogenous control, β-actin.
Immunohistochemistry
Immunohistochemistry was performed as per the manufacturer's instruction (Maixin, Fuzhou, China). The staining intensity and percentage of positively stained cells were used as two essential criteria for IHC scoring. The staining intensity was calculated as follows: 0 (no staining), 1 (weak staining), 2 (moderate staining), and 3 (strong staining). The percentage of the positively stained area was calculated as follows: 0 (no staining), 1 (1–10% staining), 2 (10–50% staining), and 3 (more than 50% staining). The staining index (SI) was calculated by multiplying the values of staining intensity and the percentage of positively stained cells. This evaluation method was employed to categorize the tumor samples into low expression (SI ≤ 4) and high expression (SI ≥ 6) groups.
Establishing the Knockdown and Overexpression of TMPRSS4 Cell Lines
TMPRSS4 knockdown and overexpression AsPC-1 and PANC-1 cell lines were established by lentivirus-mediated transfection. Human TMPRSS4 siRNA and scrambled siRNA were purchased from RiboBio Co. (Guangzhou, China). The transient transfection was performed using Lipofectamine 3000 (Invitrogen) as per the manufacturer's instructions. We cloned the TMPRSS4-targeting short hairpin RNA (shRNA) and TMPRSS4 coding sequence fragments into the pLKD-CMV-G&PR-U6 vector and pRLenti-EF1a-mCherry-P2A-Puro-CMV-3Flag vector (OBiO, Shanghai, China), to generate shTMPRSS4 lentivirus and TMPRSS4 overexpressing lentivirus, respectively. The PDAC cells were infected with concentrated lentivirus as per the manufacturer's instructions. The eGFP-positive or mCherry positive virus-infected cells were selected using puromycin.
Inhibition of the ERK1/2 Pathway
The AsPC-1 and PANC-1 cells that stably overexpressed TMPRSS4 and the control cells were treated with the specific ERK1/2 phosphorylation inhibitor-SCH772984 (2 mM; CSNpharm, USA) for 24 h. Control cells were incubated with the same volume of dimethyl sulfoxide (DMSO). The efficacy of SCH772984 inhibition was analyzed by Western blot. These treated cells were used for further analysis.
Cell Proliferation Analysis
The proliferative ability of PDAC cells was determined by the Cell Counting Kit 8 (CCK-8) assay and the 5-Ethynyl-2′-deoxyuridine (EdU) immunofluorescent assay. The CCK-8 and EdU assay were performed as per the manufacturer's instructions (Dojinodo, Shanghai, China) (RiboBio, Guangzhou, China). The CCK-8 assay was performed by seeding the 5,000 stably transfected cells into 96-well plates and adding the 10 μl of Cell Counting Kit-8 solution to each well-followed by incubation for 2 h at 37°C. The OD value was detected at 24, 48, 72, and 96 h after the cells were seeded. The EdU assay was performed by seeding the 5,000 stably transfected cells into 96-well plates, followed by 48 h of incubation with EdU and immunofluorescence staining.
Apoptosis Assay
Cellular apoptosis was quantified by flow cytometry using an Annexin V-FITC/PI Staining kit as per the manufacturer's instruction (KeyGEN BioTECH, Jiangsu, China). Briefly, 5 × 105 cells were harvested and suspended in 500 μl of binding buffer containing 5 μl annexin V-FITC and 5 μl PI. These cells were incubated for 15 min at room temperature in the dark. The apoptotic cells were detected by the FACSCalibur (BD Biosciences).
RNA Sequencing Analysis
The total RNA was isolated from the AsPC-1 cells that were stably transfected with the shTMPRSS4 lentivirus using TRIzol. Subsequently, RNA-seq analysis was performed (Genedenovo Co., Guangzhou, China) in three shTMPRSS4-treated cells along with the appropriate control samples using the Illumina HiSeq 2500 (Illumina, USA). The experimental procedure is as follows: Total RNA was extracted using Trizol reagent kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's protocol. RNA quality was assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and checked using RNase free agarose gel electrophoresis. After total RNA was extracted, eukaryotic mRNA was enriched by Oligo(dT) beads, while prokaryotic mRNA was enriched by removing rRNA by Ribo-ZeroTM Magnetic Kit (Epicentre, Madison, WI, USA). Then the enriched mRNA was fragmented into short fragments using fragmentation buffer and reverse transcripted into cDNA with random primers. Second-strand cDNA were synthesized by DNA polymerase I, RNase H, dNTP and buffer. Then the cDNA fragments were purified with QiaQuick PCR extraction kit (Qiagen, Venlo, The Netherlands), end repaired, poly (A) added, and ligated to Illumina sequencing adapters. The ligation products were size selected by agarose gel electrophoresis, PCR amplified, and sequenced using Illumina HiSeq2500 by Gene Denovo Biotechnology Co. (Guangzhou, China).
Animal Experiments
Four to 6 weeks old female BALB/c nude mice were purchased from the Peking University Animal Center (Beijing, China). TMPRSS4 knockdown or overexpressing mouse xenograft tumor models were established. For TMPRSS4 knockdown tumor model, the mice were divided into two groups (n = 8 per group). They were subcutaneously inoculated with 5 × 106 AsPC-1/negative control or AsPC-1/TMPRSS4-shRNA cells in the left oxter. The tumor length (L) and width (W) were monitored every week using a digital Vernier caliper. The tumor volume was determined using the equation (L*W2)/2. Mice were sacrificed 4 weeks post injection, and tumors were excised and weighed. The mice were sacrificed 4 weeks post injection.
The TMPRSS4 overexpressing tumor model were established to explore the association between the TMPRSS4 and ERK1/2 signaling pathways in vivo. Two groups of mice (n = 8/group) were subcutaneously inoculated with 5 × 106 AsPC-1/ TMPRSS4 or AsPC-1/Vector cells in the left dorsal flank. One week after injection, mice were randomized into four subgroups (n = 4 in each group, with mean tumor volumes of 8.74 ± 2.84, 7.48 ± 1.63, 7.38 ± 0.86, and 7.49 ± 3.35 mm3 for the TMPRSS4+DMSO, Vector+DMSO, TMPRSS4+SCH772984, and Vector+SCH772984 group, respectively) and treated with (i, ii) DMSO control, (iii, iv) 25 mg/kg SCH772984 per day by intraperitoneal injection for 3 weeks. The tumor growth was monitored, and tumor excision was performed as mentioned previously. This animal study was approved by the Southwest Hospital Ethics Committee.
Statistical Analysis
The statistical analyses were performed using SPSS statistical software (version 22.0, Chicago, IL, USA). Data were presented as the mean ± standard deviation (SD) of three independent experiments and analyzed with Student's t-test. Statistical tests included in this study for data analysis were the Chi-square test, the log-rank test, and the Student's t-test. The Cox proportional hazards regression model was employed to distinguish independent prognostic factors. All tests were two-sided, and statistical significance was defined as p < 0.05. Figures were generated using the GraphPad Prism software version 8.0.
Results
TMPRSS4 Overexpression in PDAC Is Associated With Poor Prognosis
The mRNA expression profiles of the GSE62165, GSE15471, GSE62452, and GSE16515 datasets from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) were analyzed and showed a significant overexpression of TMPRSS4 mRNA in pancreatic cancer tissue as compared to the non-tumor tissue (Figures 1A–D). These findings were validated by data from the Cancer Genome Atlas dataset (TCGA) and the Genotype-Tissue Expression (GTEx) (Figure 1E). Additionally, we performed RT-qPCR to evaluate the TMPRSS4 expression in 14 paired PDAC tissues collected from our hospital. The result demonstrated that TMPRSS4 was upregulated in PDAC tissue, which further confirmed the above data (Figure 1F).
Figure 1
Moreover, analysis of the GSE62452 dataset revealed that TMPRSS4 overexpression is correlated with poor prognosis in pancreatic cancer (Figure 1G). However, in patients with 149 PDAC pathologic type of TCGA dataset pancreatic cancer, Kaplan-Meier analysis showed no significant difference between TMPRSS4 expression and overall survival (Figure 1H). We then further explored the relationship between TMPRSS4 expression and prognosis in all 177 pancreatic cancer patients in the TCGA dataset, and found that they were negatively correlated and TCGA dataset.
In line with these findings, Western blot analysis demonstrated the upregulation of TMPRSS4 in seven pancreatic cancer cell lines (AsPC-1, BxPC-3, Capan-1, CFPAC-1, Hs766t, PANC-1 and SW1990) and in the PDAC tissue in contrast to the primary normal human pancreatic duct epithelial cell (HPDE6-C7) or adjacent non-tumor tissue (Figures 1I,J). IHC staining was followed to investigate the clinical significance of TMPRSS4 expression in PDAC in a tissue microarray containing 97 PDAC tissue and 12 normal pancreatic tissue. TMPRSS4 expression was markedly upregulated in pancreatic cancer tissue as compared to the normal pancreatic tissue (Figures 1K,L). The Kaplan–Meier survival curve and log-rank test demonstrated that TMPRSS4 overexpression was significantly correlated with the poor overall survival in PDAC patients (Figure 1M). We statistically analyzed the IHC score and found that the TMPRSS4 level was strongly correlated with age, tumor size and differentiation (Table 1), We further evaluated the prognostic significance of TMPRSS4 expression by employing univariable and multivariable Cox proportional hazards analysis (Table 2). In univariable analysis, patients with the TMPRSS4 overexpression exhibited the lowest OS (HR, 2.392, p = 0.003) (Table 2). We selected the variables whose p-value was <0.1 in Cox univariate analysis and put them into Cox multivariate regression model for further analysis. The outcome of our analysis is demonstrated in Table 2, which was in line with our previous analysis (HR, 1.898, p = 0.033).
Table 1
| Characteristics | Low group (n = 47) | High group (n = 50) | P-value |
|---|---|---|---|
| Age(years) | |||
| <60 years | 23 (48.9) | 13 (26.0) | 0.033 |
| ≥ 60 years | 24 (51.1) | 37 (74.0) | |
| Gender | |||
| female | 22 (46.8) | 16 (32.0) | 0.199 |
| male | 25 (53.2) | 34 (68.0) | |
| Tumor size (cm) | |||
| 0–2 | 7 (14.9) | 7 (14.0) | 0.021 |
| ≥2, <5 | 38 (80.9) | 31 (62.0) | |
| ≥5 | 2 (4.3) | 12 (24.0) | |
| Differentiation | |||
| Moderate-Well | 40 (85.1) | 31 (62.0) | 0.019 |
| Poor | 7 (14.9) | 19 (38.0) | |
| Clinical stage | |||
| I-II | 44 (93.6) | 42 (84.0) | 0.241 |
| III-IV | 3 (6.4) | 8 (16.0) | |
| Perineuronal invasion | |||
| No | 36 (76.6) | 32 (64.0) | 0.257 |
| Yes | 11 (23.4) | 18 (36.0) | |
Clinical characteristics of pancreatic cancer stratified with TMPRSS4 expression.
Table 2
| Clinical factor | Univariable analysis | Multivariable analysis | ||||
|---|---|---|---|---|---|---|
| HR | 95%CI | P-value | HR | 95%CI | P-value | |
| Age (≥60 vs. <60 years) | 1.786 | 0.984–3.241 | 0.056 | |||
| Gender (male vs. female) | 1.649 | 0.922–2.948 | 0.092 | |||
| Differentiation (Poor vs. Moderate+Well) | 2.155 | 1.163–3.994 | 0.015 | 2.353 | 1.227–4.513 | 0.01 |
| Tumor size (vs. 0–2 cm) | ||||||
| ≥ 2, <5 | 2.161 | 0.664–7.035 | 0.201 | - | - | - |
| ≥ 5 | 5.133 | 1.393–18.914 | 0.014 | |||
| Perineuronal invasion (Positive vs. Negative) | 1.062 | 0.565–1.995 | 0.851 | |||
| Clinical stage (III-IV vs. I-II) | 6.938 | 3.051–15.779 | <0.001 | 7.590 | 3.199–18.009 | <0.001 |
| TMPRSS4 expression (High vs. Low) | 2.392 | 1.357–4.215 | 0.003 | 1.898 | 1.052–3.426 | 0.033 |
Cox univariable and multivariable analysis of clinicopathological variables and TMPRSS4 expression in relation to OS in pancreatic cancer patients.
In total TCGA pancreatic cancer, TMPRSS4 expression is associated with prognosis. Probably due to ethnic differences, the expression level of TMPRSS4 was associated with overall survival in this study, which was not in 149 PDAC pathologic type of TCGA. These findings suggest that TMPRSS4 is overexpressed in pancreatic cancer tissue, and it can be an independent prognostic factor in PDAC patients.
TMPRSS4 Promotes Cell Proliferation and Inhibits Cell Apoptosis in Pancreatic Cancer Both in vitro and in vivo
The overexpression of TMPRSS4 in multiple malignancies suggests that it could be an oncogenic promoter in PDAC. To explore the biological function of TMPRSS4 in PDAC development and progression, the AsPC-1 and PANC-1 cell lines with TMPRSS4 knockdown and overexpression were established by lentivirus-mediated transfection.
The EdU and CCK-8 assays were performed and the results indicated a significant decrease in cell proliferation in the pancreatic cancer cells with reduced TMPRSS4 expression (Figures 2A,D). In contrast, cells with induced TMPRSS4 overexpression promoted cellular proliferation (Figures 2B,E). Subsequently, we performed flow cytometry assay to quantify the cellular apoptosis (Figures 2F–H). The percentage of apoptotic cells was significantly higher in siTMPRSS4 pancreatic cancer cells than in the control cells (Figure 2F) but lower in the TMPRSS4 overexpressing cells (Figure 2G).
Figure 2
Furthermore, TMPRSS4 downregulation was found to be associated with the reducing tumor volume or weight of the subcutaneous xenografts in nude mice (Figures 3A–C). IHC staining of xenograft tissues showed reduced expression of proliferating cell nuclear antigen (PCNA), a biomarker for cell proliferation (Figure 3D).
Figure 3
We also examined the expression of apoptosis-related proteins in the AsPC-1-shTMPRSS4, PANC-1-shTMPRSS4, and control cells (Figure 4B). We found that pro-apoptotic proteins such as Bax, cleaved caspase 3 were increased, and anti-apoptosis protein Bcl-2 was decreased in the TMPRSS4 silencing cells. Conversely, overexpression of TMPRSS4 led to the reduced pro-apoptotic proteins and increased anti-apoptotic proteins. These data indicate that overexpression of TMPRSS4 inhibited apoptosis in the PDAC cells.
Figure 4
Furthermore, TMPRSS4 downregulation was found to be associated with the reducing tumor volume or weight of the subcutaneous xenografts in nude mice (Figures 3A–C). IHC staining of xenograft tissues showed reduced expression of proliferating cell nuclear antigen (PCNA), a biomarker for cell proliferation (Figure 3D). Conversely, upregulation of TMPRSS4 demonstrated opposite effects (Figures 6B–E, SS4+DMSO vs. vector+DMSO). These results suggest that TMPRSS4 plays critical roles in regulating cellular proliferation and apoptosis in PDAC cells.
TMPRSS4 Activates ERK1/2 Signaling Pathway in Pancreatic Cancer
We performed RNA-seq to investigate which signaling pathways are associated with the TMPRSS4 expression level. Interestingly, the Gene Set Enrichment Analysis (GSEA) revealed that the TMPRSS4 overexpression is strongly correlated with genes associated with Mitogen-activated protein kinases (MAPKs) signaling pathways (Figure 4A). Moreover, previous studies have validated the fact that the activation of ERK1/2 axis contributes to the tumorigenic effect of TMPRSS4 (, ). Therefore, we probed the correlation between the TMPRSS4 and the three classical MAPKs signaling pathways, that is, ERK1/2, JNK, and p38 MAPKs signaling pathways, and found that only the ERK1/2 pathway had a positive correlation with the TMPRSS4 expression (Figure 4A). The phosphorylation of ERK1/2 was decreased in TMPRSS4 knockdown cells; whereas it was elevated in TMPRSS4 overexpressing cells (Figure 4B).
The ERK1/2 Signaling Pathway Is Essential for the TMPRSS4-Induced Proliferative and Apoptotic Behavior in PDAC
We treated the TMPRSS4-overexpressing cells with the ERK1/2 phosphorylation inhibitor, that is, SCH772984, to evaluate the effect of ERK1/2 on the TMPRSS4-induced cell proliferation and apoptosis inhibition. As depicted in Figure 5, the TMPRSS4-induced proliferative and anti-apoptotic behavior was reversed by the SCH772984 treatment. After blocking the ERK1/2 phosphorylation, pro-apoptotic proteins Bax, and cleaved caspase three levels were increased whereas anti-apoptotic protein Bcl-2 level was diminished in the TMPRSS4-overexpressing PDAC cells (Figure 6A). In human PDAC xenograft AsPC-1 models, oncogenic effect induced by the TMPRSS4 overexpression was remarkably inhibited by SCH772984 treatment. In line with the EdU and CCK-8 assays results in vitro, IHC staining of the xenograft tumor tissue demonstrated that the PCNA expression was significantly decreased in the SCH772984 treated xenografts (Figure 6E, SS4+DMSO vs. SS4+SCH). Taken together, these results further validate the fact that the activation of the ERK1/2 signaling pathway affects the functional outcome of TMPRSS4, that is, promotes cell proliferation and inhibits apoptosis in pancreatic cancer.
Figure 5
Figure 6
Discussion
In the present study, we found that TMPRSS4 is significantly overexpressed in PDAC tissues and cell lines. Moreover, overexpression of TMPRSS4 is positively correlated with tumor size and differentiation, also a reduced overall survival rate. Additionally, we performed gain-and loss-of-function analysis and found that TMPRSS4 promotes cellular proliferation and inhibits apoptosis in PDAC. Importantly, the outcome of this study revealed the critical role of TMPRSS4 in mediating cell proliferation and apoptosis in pancreatic cancer cells via the ERK1/2 signaling pathway.
TMPRSS4 was found to be overexpressed in multiple human malignancies such as thyroid, lung, breast, pancreatic, gastric, colon, and other cancers. As an oncogenic protein, TMPRSS4 significantly contributes to the development of cancerous tumors. A recent study suggests that TMPRSS4 regulates both proliferation and invasion through Slug and cyclin D1 in prostate cancer cells (). Similarly, TMPRSS4 facilitates cellular proliferation via CREB phosphorylation in thyroid cancer (). Additionally, TMPRSS4 plays an influential role in radiosensitivity and chemosensitivity by disrupting cell cycle and apoptosis in lung cancer and triple-negative breast cancer (, ). These results are in line with our findings that TMPRSS4 accelerates cell proliferation and inhibits apoptosis in PDAC.
Cancer cells are characterized by sustained cell proliferation and diminished apoptosis (). In the present study, we investigated the role of TMPRSS4 in cell proliferation and apoptosis in PDAC cells. Intriguingly, TMPRSS4 might activate the MAPKs signaling pathways to promote cellular proliferation and suppress apoptosis of PDAC cells.
MAPKs signaling pathways play crucial roles in multiple biological processes by converting the extracellular stimuli into cellular responses. These pathways were reported to be overactivated in human cancers. The MAPKs signaling kinases are chiefly divided into three families: extracellular signal-regulated kinase (ERK), Jun N-terminal kinase (JNK), and p38 MAPKs (, ). Previous reports suggest that TMPRSS4 might assist the cellular proliferation, progression, and invasion in lung and colon cancer cells by activating the ERK1/2 and p38 MAPK signaling pathways (). Similarly, Yunhee Lee et al., reported that the TMPRSS4 induced AP-1 activation is mediated by Axl overexpression and subsequent ERK1/2 and JNK signaling pathway activation, which in turn promotes prostate cancer cell proliferation and invasion (). Interestingly, in line with these findings, our GSEA results also showed that a higher TMPRSS4 level correlated with gene signatures associated with MAPKs signaling pathways. Further in-depth analysis revealed that among the three classical MAPKs signaling pathways, the ERK1/2 signaling pathway was positively correlated with TMPRSS4 expression, but not JNK and p38 MAPKs signaling pathways.
In the results of genes differentially expressed by RNA-seq, we found MAPKs signaling pathway-related genes such as ACNG4, JUND, TGFA, PTK2B, etc (The annotated informations are in Supplementary Material 3). The bioinformatics analysis revealed that the functions of TMPRSS4 are mediated via the ERK1/2 signaling pathway. Our results demonstrated that TMPRSS4 knocked down leads to decreased p-ERK1/2 expression, while TMPRSS4 overexpression leads to increased p-ERK1/2 expression in AsPC-1 and PANC-1 cells. We found that the apoptosis-related proteins Bax, and cleaved caspase 3 were overexpressed while the apoptosis-inhibited protein Bcl2 was downregulated in TMPRSS4 knocked down AsPC-1 and PANC-1 cells. Moreover, we found a decreased expression of Bax and cleaved caspase3 and an increased expression of Bcl2 in TMPRSS4 overexpressing cells.
Caspase 3, a crucial component of the caspase cascade, is activated in both the extrinsic and intrinsic death signal pathways (). Bcl2 family proteins are key regulators of apoptosis and encompass both proapoptotic proteins (Bid, Bak, Bax, and Bim) and anti-apoptotic proteins (Bcl-2, Bcl-XL, and Mcl-1) ().
The TMPRSS4 overexpression activates the ERK1/2 signaling pathway and facilitated the subsequent cell survival by upregulating anti-apoptotic protein, Bcl2 and downregulating apoptotic proteins, Bax and cleaved caspase 3. Previous studies have shown that ERK1/2 plays a crucial role in cell apoptosis. Adrian Achuthan et al. reported that suppression of ERK1/2 activity resulted in p90 ribosomal-S6 kinase deactivation and the Bad dephosphorylation, which in turn promoted caspase-3 activity and the subsequent cell apoptosis (39). Another study stated that ACY-1215, an HDAC6-selective inhibitor, inhibited cellular proliferation, and promoted apoptosis via PI3K/AKT/mTOR and ERK pathways (40). In this study, we validated that the ERK1/2 signaling pathway is responsible for the TMPRSS4-induced cell proliferative and apoptotic behaviors, that were reversed by ERK inhibitor SCH772984 treatment.
In conclusion, TMPRSS4 is overexpressed in PDAC tissues and cell lines. TMPRSS4, as an oncogene in PDAC, promotes cellular proliferation and inhibits apoptosis by activating the ERK1/2 pathway. Our results suggest that TMPRSS4 protein could serve as a potential biomarker in PDAC prognosis and a potential therapeutic target for PDAC treatment.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: https://www.ncbi.nlm.nih.gov/geo/(GEO database); https://xenabrowser.net/datapages/(TCGA data); https://xenabrowser.net/datapages/(GTEx data).
Ethics statement
The studies involving human participants were reviewed and approved by Southwest Hospital Ethics Committee. The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by Southwest Hospital Ethics Committee.
Author contributions
JG and WH contributed to the conception, design, acquisition of data and data analysis, obtaining the findings, and drafting the manuscript. JZ, XW, TT, LY, YZ, SL, JY, and LZ contributed to acquisition of data and technical support. YF and HW offered conception, design and critical revision of the manuscript for important intellectual content. All authors read and approved the final manuscript.
Funding
This work was supported by the Key-Area Research and Development Program of Guangdong Province, China (2020B010165004) for YF and the National Key R&D Program of China (No. 2017YFC1308600) and the National Natural Science Foundation of China (No. 81672382) for HW.
Acknowledgments
Chuanming Xie, Yujun Zhang, Jiejuan Lai, Li Shuai of the Institute of Hepatopancreatobiliary Surgery, Southwest Hospital, Army Medical University in Chongqing received our sincere thanks for technical guidance and experimental suggestions in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.628353/full#supplementary-material
References
1.
HowladerNNooneAMKrapchoMMillerDBrestAYuMet al. (editors). SEER Cancer Statistics Review, 1975-2017. Bethesda, MD: National Cancer Institute (2019). Available online at: https://seer.cancer.gov/csr/1975_2017/
2.
Nevala-PlagemannCHidalgoMGarrido-LagunaI. From state-of-the-art treatments to novel therapies for advanced-stage pancreatic cancer. Nature Rev Clin Oncol. (2019) 17:108–23. 10.1038/s41571-019-0281-6
3.
KleeffJKorcMApteMLa VecchiaCJohnsonCDBiankinAVet al. Pancreatic cancer. Nat Rev Dis Primers. (2016) 2:16022. 10.1038/nrdp.2016.22
4.
WallrappCHahnelSMuller-PillaschFBurghardtBIwamuraTRuthenburgerMet al. A novel transmembrane serine protease (TMPRSS3) overexpressed in pancreatic cancer. Cancer Res. (2000) 60:2602–6.
5.
de AberasturiALCalvoA. TMPRSS4: an emerging potential therapeutic target in cancer. Br J Cancer. (2015) 112:4–8. 10.1038/bjc.2014.403
6.
MartinCEListK. Cell surface-anchored serine proteases in cancer progression and metastasis. Cancer Metastasis Rev. (2019) 38:357–87. 10.1007/s10555-019-09811-7
7.
LarzabalLNguewaPAPioRBlancoDSanchezBRodriguezMJet al. Overexpression of TMPRSS4 in non-small cell lung cancer is associated with poor prognosis in patients with squamous histology. Br J Cancer. (2011) 105:1608–14. 10.1038/bjc.2011.432
8.
ChikaishiYUramotoHKoyanagiYYamadaSYanoSTanakaF. TMPRSS4 Expression as a Marker of Recurrence in Patients with Lung Cancer. Anticancer Res. (2016) 36:121–7.
9.
VillalbaMDiaz-LagaresARedradoMde AberasturiALSeguraVBodegasMEet al. Epigenetic alterations leading to TMPRSS4 promoter hypomethylation and protein overexpression predict poor prognosis in squamous lung cancer patients. Oncotarget. (2016) 7:22752–69. 10.18632/oncotarget.8045
10.
FanXLiangYLiuYBaiYYangCXuS. The upregulation of TMPRSS4, partly ascribed to the downregulation of miR-125a-5p, promotes the growth of human lung adenocarcinoma via the NFkappaB signaling pathway. Int J Oncol. (2018) 53:148–58. 10.3892/ijo.2018.4396
11.
VillalbaMExpositoFPajaresMJSainzCRedradoMRemirezAet al. TMPRSS4: A Novel Tumor Prognostic Indicator for the Stratification of Stage IA Tumors and a Liquid Biopsy Biomarker for NSCLC Patients. J Clin Med. (2019) 8:2134. 10.3390/jcm8122134
12.
ChengDKongHLiY. TMPRSS4 as a poor prognostic factor for triple-negative breast cancer. Int J Mol Sci. (2013) 14:14659–68. 10.3390/ijms140714659
13.
LiangBWuMBuYZhaoAXieF. Prognostic value of TMPRSS4 expression in patients with breast cancer. Med Oncol. (2013) 30:497. 10.1007/s12032-013-0497-8
14.
LiXMLiuWLChenXWangYWShiDBZhangHet al. Overexpression of TMPRSS4 promotes tumor proliferation and aggressiveness in breast cancer. Int J Mol Med. (2017) 39:927–35. 10.3892/ijmm.2017.2893
15.
ChengDLiangBLiY. High TMPRSS4 expression is a predictor of poor prognosis in cervical squamous cell carcinoma. Cancer Epidemiol. (2013) 37:993–7. 10.1016/j.canep.2013.08.009
16.
DaiWZhouQXuZZhangE. Expression of TMPRSS4 in patients with salivary adenoid cystic carcinoma: correlation with clinicopathological features and prognosis. Med Oncol. (2013) 30:749. 10.1007/s12032-013-0749-7
17.
HuangAZhouHZhaoHQuanYFengBZhengM. High expression level of TMPRSS4 predicts adverse outcomes of colorectal cancer patients. Med Oncol. (2013) 30:712. 10.1007/s12032-013-0712-7
18.
LuoZYWangYYZhaoZSLiBChenJF. The expression of TMPRSS4 and Erk1 correlates with metastasis and poor prognosis in Chinese patients with gastric cancer. PLoS ONE. (2013) 8:e70311. 10.1371/journal.pone.0070311
19.
ShengHShenWZengJXiLDengL. Prognostic significance of TMPRSS4 in gastric cancer. Neoplasma. (2014) 61:213–7. 10.4149/neo_2014_027
20.
ChengYWangKGengLSunJXuWLiuDet al. Identification of candidate diagnostic and prognostic biomarkers for pancreatic carcinoma. EBioMedicine. (2019) 40:382–93. 10.1016/j.ebiom.2019.01.003
21.
ExpositoFVillalbaMRedradoMde AberasturiALCirauquiCRedinEet al. Targeting of TMPRSS4 sensitizes lung cancer cells to chemotherapy by impairing the proliferation machinery. Cancer Lett. (2019) 453:21–33. 10.1016/j.canlet.2019.03.013
22.
JinJShenXChenLBaoLWZhuLM. TMPRSS4 promotes invasiveness of human gastric cancer cells through activation of NF- kB/MMP-9 signaling. Biomed Pharmacother. (2016) 77:30–6. 2016/01/23. 10.1016/j.biopha.2015.11.002
23.
JankyRBindaMMAllemeerschJVan den BroeckAGovaereOSwinnenJVet al. Prognostic relevance of molecular subtypes and master regulators in pancreatic ductal adenocarcinoma. BMC Cancer. (2016) 16:632. 10.1186/s12885-016-2540-6
24.
BadeaLHerleaVDimaSODumitrascuTPopescuI. Combined gene expression analysis of whole-tissue and microdissected pancreatic ductal adenocarcinoma identifies genes specifically overexpressed in tumor epithelia. Hepato Gastroenterol. (2008) 55:2016–27.
25.
IdichiTSekiNKuraharaHYonemoriKOsakoYAraiTet al. Regulation of actin-binding protein ANLN by antitumor miR-217 inhibits cancer cell aggressiveness in pancreatic ductal adenocarcinoma. Oncotarget. (2017) 8:53180–93. 10.18632/oncotarget.18261
26.
YangSHePWangJSchetterATangWFunamizuNet al. A Novel MIF Signaling Pathway Drives the Malignant Character of Pancreatic Cancer by Targeting NR3C2. Cancer Res. (2016) 76:3838–50. 10.1158/0008-5472.can-15-2841
27.
LiLZhangJWJenkinsGXieFCarlsonEEFridleyBLet al. Genetic variations associated with gemcitabine treatment outcome in pancreatic cancer. Pharmacogenet Genomics. (2016) 26:527–37. 10.1097/fpc.0000000000000241
28.
GoldmanMJCraftBHastieMRepečkaKMcDadeFKamathAet al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. (2020) 38:675–8. 10.1038/s41587-020-0546-8
29.
JungHLeeKPParkSJParkJHJangYSChoiSYet al. TMPRSS4 promotes invasion, migration and metastasis of human tumor cells by facilitating an epithelial-mesenchymal transition. Oncogene. (2008) 27:2635–47. 10.1038/sj.onc.1210914
30.
LiJWuHLiWYinLGuoSXuXet al. Downregulated miR-506 expression facilitates pancreatic cancer progression and chemoresistance via SPHK1/Akt/NF-kappaB signaling. Oncogene. (2016) 35:5501–14. 10.1038/onc.2016.90
31.
LeeYKoDMinHJKimSBAhnHMLeeYet al. TMPRSS4 induces invasion and proliferation of prostate cancer cells through induction of Slug and cyclin D1. Oncotarget. (2016) 7:50315–32. 10.18632/oncotarget.10382
32.
GuanHLiangWLiuJWeiGLiHXiuLet al. Transmembrane protease serine 4 promotes thyroid cancer proliferation via CREB phosphorylation. Thyroid. (2015) 25:85–94. 10.1089/thy.2014.0155
33.
AssaniGSegboJYuXYessoufouAXiongYZhouFet al. Downregulation of TMPRSS4 Enhances Triple-Negative Breast Cancer Cell Radiosensitivity Through Cell Cycle and Cell Apoptosis Process Impairment. Asian Pacific J Cancer Prev. (2019) 20:3679–87. 10.31557/apjcp.2019.20.12.3679
34.
HanahanDWeinbergRA. Hallmarks of cancer: the next generation. Cell. (2011) 144:646–74. 10.1016/j.cell.2011.02.013
35.
McKinnonBDKocbekVNirgianakisKBersingerNAMuellerMD. Kinase signalling pathways in endometriosis: potential targets for non-hormonal therapeutics. Hum Reprod Update. (2016) 22:382–403. 10.1093/humupd/dmv060
36.
WagnerEFNebredaAR. Signal integration by JNK and p38 MAPK pathways in cancer development. Nat Rev Cancer. (2009) 9:537–49. 10.1038/nrc2694
37.
NagataS. Apoptosis and Clearance of Apoptotic Cells. Annu Rev Immunol. (2018) 36:489–517. 10.1146/annurev-immunol-042617-053010
38.
ChipukJEGreenDR. How do BCL-2 proteins induce mitochondrial outer membrane permeabilization?Trends Cell Biol. (2008) 18:157–64. 10.1016/j.tcb.2008.01.007
39.
AchuthanAAslamASMNguyenQLamPYFleetwoodAJFryeATet al. Glucocorticoids promote apoptosis of proinflammatory monocytes by inhibiting ERK activity. Cell Death Dis. (2018) 9:267. 10.1038/s41419-018-0332-4
40.
CaoJLvWWangLXuJYuanPHuangSet al. Ricolinostat (ACY-1215) suppresses proliferation and promotes apoptosis in esophageal squamous cell carcinoma via miR-30d/PI3K/AKT/mTOR and ERK pathways. Cell Death Dis. (2018) 9:817. 10.1038/s41419-018-0788-2
Summary
Keywords
TMPRSS4, proliferation, apoptosis, ERK1/2, pancreatic ductal adenocarcinoma
Citation
Gu J, Huang W, Zhang J, Wang X, Tao T, Yang L, Zheng Y, Liu S, Yang J, Zhu L, Wang H and Fan Y (2021) TMPRSS4 Promotes Cell Proliferation and Inhibits Apoptosis in Pancreatic Ductal Adenocarcinoma by Activating ERK1/2 Signaling Pathway. Front. Oncol. 11:628353. doi: 10.3389/fonc.2021.628353
Received
11 November 2020
Accepted
23 February 2021
Published
18 March 2021
Volume
11 - 2021
Edited by
Laura Pellegrini, Institute of Oncology Research (IOR), Switzerland
Reviewed by
Tao Li, Hubei University of Medicine, China; Didier Dreau, University of North Carolina at Charlotte, United States; Swati Venkat, Roswell Park Comprehensive Cancer Center, United States
Updates
Copyright
© 2021 Gu, Huang, Zhang, Wang, Tao, Yang, Zheng, Liu, Yang, Zhu, Wang and Fan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huaizhi Wang whuaizhi@gmail.comYingfang Fan fanxifan@smu.edu.cn
This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Oncology
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.