Orchestrating a biomarker panel with lncRNAs and mRNAs for predicting survival in pancreatic ductal adenocarcinoma

The low survival of patients with pancreatic ductal adenocarcinoma (PDAC) makes the treatment of this disease one of the most challenging task in modern medicine. Here, by mining a large‐scale cancer genome atlas data set of pancreatic cancer tissues, we identified 21 long noncoding RNAs (lncRNAs) that significantly associated with overall survival in patients with PDAC (P < .01). Further analysis revealed that 8 lncRNAs turned out to be independently correlated with patients’ overall survival, and the risk score could be calculated based on their expression. To obtain a better predicting power, we integrated lncRNA data with a total of 410 differently expressed messenger RNAs (mRNAs) screened from PDAC and normal tissues in gene expression omnibus (GEO) database. The integration resulted in a much better panel including 8 lncRNAs (RP3.470B24.5, CTA.941F9.9, RP11.557H15.3, LINC00960, AP000479.1, LINC00635, LINC00636, and AC073133.1) and 8 mRNAs (DHRS9, ONECUT1, OR8D4, MT1M, TCN1, MMP9, DPYSL3, and TTN) to predict prognosis. A functional evaluation showed that these lncRNAs might play roles in pancreatic secretion, cell adhesion, and proteolysis. Using normal and pancreatic cancer cell lines, we confirmed that a majority of identified lncRNAs and mRNAs showed altered expressions in pancreatic cancer cells. Especially, LINC01589, LINC00960, TCN1, and MT1M showed a profoundly increased expression in pancreatic cancer cells, which suggests their potentially important role in pancreatic cancer. The results of our work indicate that lncRNAs have vital roles in PADC and provide new insights to integrate multiple kinds of markers in clinical practices.     

Prognostic significance of microRNA-141 expression and its tumor suppressor function in human pancreatic ductal adenocarcinoma

Increasing evidence shows that dysregulation of microRNAs is correlated with tumor development. This study was performed to determine the expression of miR-141 and investigate its clinical significance in pancreatic ductal adenocarcinoma (PDAC). Taqman quantitative RT-PCR was used to detect miR-141 expressions in 94 PDAC tissues and 16 nontumorous pancreatic tissues. Correlations between miR-141 expression and clinicopathologic features and prognosis of patients were statistically analyzed. The effects of miR-141 expression on growth and apoptosis of PDAC cell line (PANC-1) were determined by MTT, colony formation, and flow cytometry assays. Potential target genes were identified by luciferase reporter and Western blot assays. The expression level of miR-141 in PDAC tissues was significantly lower than that in corresponding nontumorous tissues. Downregulation of miR-141 correlated with poorer pT and pN status, advanced clinical stage, and lymphatic invasion. Also, low miR-141 expression in PDAC tissues was significantly correlated with shorter overall survival, and multivariate analysis showed that miR-141 was an independent prognostic factor for PDAC patients. Further, functional researches suggested that miR-141 inhibits growth and colony formation, and enhances caspase-3-dependent apoptosis in PANC-1 cells by targeting Yes-associated protein-1 (YAP1). Therefore, miR-141 is an independent prognostic factor for PDAC patients, and functions as a tumor suppressor gene by targeting YAP1.     

Enhancer of zeste homologue 2 (EZH2) down-regulates RUNX3 by increasing histone H3 methylation

Overexpression of enhancer of zeste homologue 2 (EZH2) occurs in various malignancies and is associated with a poor prognosis, especially because of increased cancer cell proliferation. In this study we found an inverse correlation between EZH2 and RUNX3 gene expression in five cancer cell lines, i.e. gastric, breast, prostate, colon, and pancreatic cancer cell lines. Chromatin immunoprecipitation assay showed an association between EZH2 bound to the RUNX3 gene promoter, and trimethylated histone H3 at lysine 27, and HDAC1 (histone deacetylase 1) bound to the RUNX3 gene promoter in cancer cells. RNA interference-mediated knockdown of EZH2 resulted in a decrease in H3K27 trimethylation and unbound HDAC1 and an increase in expression of the RUNX3 gene. Restoration of RUNX3 expression was not associated with any change in DNA methylation status in the RUNX3 promoter region. RUNX3 was repressed by histone deacetylation and hypermethylation of a CpG island in the promoter region and restored by trichostatin A or/and 5-aza-2'-deoxycytidine. Immunofluorescence staining confirmed restoration of expression of the RUNX3 protein after knockdown of EZH2 and its restoration resulted in decreased cell proliferation. In vivo, an inverse relationship between expression of the EZH2 and RUNX3 proteins was observed at the individual cell level in gastric cancer patients in the absence of DNA methylation in the RUNX3 promoter region. The results showed that RUNX3 is a target for repression by EZH2 and indicated an underlying mechanism of the functional role of EZH2 overexpression on cancer cell proliferation.     

LINC01133 inhibits breast cancer invasion and metastasis by negatively regulating SOX4 expression through EZH2

Mounting evidence highlights long non‐coding RNAs (lncRNAs) as crucial regulators in multiple types of biological processes and contributing to tumourigenesis. LINC01133, located in chromosome 1q23.2, was a recently identified novel lncRNA with a length of 1154nt. It was involved in the development of colorectal cancer and non‐small cell lung cancer. However, its clinical relevance, biological functions and potential molecular mechanism in breast cancer are still unclear. In this study, we found that the LINC01133 expression was significantly down‐regulated in breast cancer samples and was associated with progression and poor prognosis of breast cancer. Further experiments demonstrated that overexpression of LINC01133 inhibited invasion and metastasis in breast cancer both in vitro and in vivo. Mechanistic investigations revealed that LINC01133 repressed SOX4 expression by recruiting EZH2 to SOX4 promoter. Moreover, rescue experiments further confirmed that LINC01133 functional acted as an anti‐oncogene, at least partly, via repressing SOX4 in breast cancer. Taken together, these findings imply that LINC01133 could serve as a novel prognostic biomarker and potential therapeutic target for breast cancer.     

RASAL2 mediated the enhancement of YAP1/TIAM1 signaling promotes malignant phenotypes of pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is characterized by a high incidence of metastasis and dismal prognosis. As a member of Gas-Gap gene, RASAL2 is involved in the hydrolysis of RAS-GTP to RAS-GDP and abnormal expression in human cancers. Here we firstly described the function of RASAL2 on PDAC to enrich the knowledge of RAS family.We interestingly observed that RASAL2 expression was upregulated in PDAC at both mRNA and protein levels, and high expression of RASAL2 predicted a poor prognosis in PDAC patients. Additionally, RASAL2 promoted malignant behaviors of PDAC in vitro and in vivo. To determine the mechanistic roles of RASAL2 signaling and its potential as a therapeutic target in PDAC, we clarified that RASAL2 could accumulate the TIAM1 expression in different level through inhibiting YAP1 phosphorylation, increased TIAM1 mRNA expression and suppressed ubiquitination of TIAM1 protein. In conclusion, RASAL2 enhances YAP1/TIAM1 signaling and promotes PDAC development and progression.     

Five key lncRNAs considered as prognostic targets for predicting pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis, and the 5‐year survival rate was only 7.7%. To improve prognosis, a screening biomarker for early diagnosis of pancreatic cancer is in urgent need. Long non‐coding RNA (lncRNA) expression profiles as potential cancer prognostic biomarkers play critical roles in development of tumorigenesis and metastasis of cancer. However, lncRNA signatures in predicting the survival of a patient with PDAC remain unknown. In the current study, we try to identify potential lncRNA biomarkers and their prognostic values in PDAC. LncRNAs expression profiles and corresponding clinical information for 182 cases with PDAC were acquired from The Cancer Genome Atlas (TCGA). A total of 14 470 lncRNA were identified in the cohort, and 175 PDAC patients had clinical variables. We obtained 108 differential expressed lncRNA via R packages. Univariate and multivariate Cox proportional hazards regression, lasso regression was performed to screen the potential prognostic lncRNA. Five lncRNAs have been recognized to significantly correlate with OS. We established a linear prognostic model of five lncRNA (C9orf139, MIR600HG, RP5‐965G21.4, RP11‐436K8.1, and CTC‐327F10.4) and divided patients into high‐ and low‐risk group according to the prognostic index. The five lncRNAs played independent prognostic biomarkers of OS of PDAC patients and the AUC of the ROC curve for the five lncRNAs signatures prediction 5‐year survival was 0.742. In addition, targeted genes of MIR600HG, C9orf139, and CTC‐327F10.4 were explored and functional enrichment was also conducted. These results suggested that this five‐lncRNAs signature could act as potential prognostic biomarkers in the prediction of PDAC patient's survival.     

Deregulation of the MiR-193b-KRAS Axis Contributes to Impaired Cell Growth in Pancreatic Cancer

Modulation of KRAS activity by upstream signals has revealed a promising new approach for pancreatic cancer therapy; however, it is not clear whether microRNA-associated KRAS axis is involved in the carcinogenesis of pancreatic cancer. Here, we identified miR-193b as a tumor-suppressive miRNA in pancreatic ductal adenocarcinoma (PDAC). Expression analyses revealed that miR-193b was downregulated in (10/11) PDAC specimens and cell lines. Moreover, we found that miR-193b functioned as a cell-cycle brake in PDAC cells by inducing G1-phase arrest and reducing the fraction of cells in S phase, thereby leading to dampened cell proliferation. miR-193b also modulated the malignant transformation phenotype of PDAC cells by suppressing anchorage-independent growth. Mechanistically, KRAS was verified as a direct effector of miR-193b, through which the AKT and ERK pathways were modulated and cell growth of PDAC cells was suppressed. Taken together, our findings indicate that miR-193b-mediated deregulation of the KRAS axis is involved in pancreatic carcinogenesis, and suggest that miR-193b could be a potentially effective target for PDAC therapy.     

Long non‐coding RNA LINC01535 promotes cervical cancer progression via targeting the miR‐214/EZH2 feedback loop

Long non‐coding RNAs (lncRNAs) have shown critical roles in multiple cancers via competitively binding common microRNAs. miR‐214 has been proved to play tumour suppressive roles in various cancers, including cervical cancer. In this study, we identified that lncRNA LINC01535 physically binds miR‐214, relieves the repressive roles of miR‐214 on its target EZH2, and therefore up‐regulates EZH2 protein expression. Intriguingly, we also found that EZH2 directly represses the expression of miR‐214. Thus, miR‐214 and EZH2 form double negative regulatory loop. Through up‐regulating EZH2, LINC01535 further represses miR‐214 expression. Functional experiments showed that enhanced expression of LINC01535 promotes cervical cancer cell growth, migration and invasion in vitro and cervical cancer xenograft growth in vivo. Reciprocally, LINC01535 knockdown suppresses cervical cancer cell growth, migration and invasion. Activation of the miR‐214/EZH2 regulatory loop by overexpression of miR‐214 or silencing of EZH2 reverses the roles of LINC01535 in promoting cervical canc`er cell growth, migration and invasion in vitro and cervical cancer xenograft growth in vivo. Clinically, LINC01535 is significantly up‐regulated in cervical cancer tissues and correlated with advanced clinical stage and poor prognosis. Moreover, the expression of LINC01535 is reversely associated with the expression of miR‐214 and positively associated with the expression of EZH2 in cervical cancer tissues. In conclusion, this study reveals that LINC01535 promotes cervical cancer progression via repressing the miR‐214/EZH2 regulatory loop.     

Identification of multi-SH3 domain-containing protein interactome in pancreatic cancer: a yeast two-hybrid approach

Pancreatic ductal adenocarcinoma (PDAC) is a fatal disease that shows minimal response to chemotherapy. Genetic changes involved in the progression of PDAC concern genes that encode proteins related to signal transduction networks. This fact reveals the importance in identifying the role and the relations between multiple signaling cascades in PDAC. One of the major factors that modulate signaling events is multidomain scaffold proteins that function by binding several proteins simultaneously, inducing their proximity and influencing the outcome of signaling. A particular group among them, containing multiple Src homology 3 (SH3) domains that can bind proteins containing proline-rich motifs, was associated to different aspects of cancer cell homeostasis. In this work, using a microarray-based analysis, we have shown that 13 multiple SH3 domain containing scaffold proteins are expressed in PDAC cells. Using a yeast two-hybrid approach, we have identified proteins that interact with these adaptor proteins. Among them we have found several molecules that modulate cell proliferation and survival (CIZ1, BIRC6, RBBP6), signaling (LTBP4, Notch2, TOM1L1, STK24) and membrane dynamics (PLSCR1, DDEF2, VCP). Our results indicate that interactions mediated by multi-SH3 domain-containing proteins could lead to the formation of dynamic protein complexes that function in pancreatic cancer cell signaling. The identification of such protein complexes is of paramount importance in deciphering pancreatic cancer biology and designing novel therapeutic approaches.     

S100 Calcium-Binding Protein A6 Promotes Epithelial-Mesenchymal Transition through β-Catenin in Pancreatic Cancer Cell Line

The pathogenesis of pancreatic ductal adenocarcinoma (PDAC) remains poorly understood. S100 calcium-binding protein A6 (S100A6) has been associated with PDAC; however, the effect of S100A6 on PDAC migration and invasion has not yet been explored. In this study, Panc-1 cells were transfected with a plasmid to induce overexpression of S100A6, and β-catenin was knocked down using a specific short hairpin RNA (shRNA). The wound-healing and Transwell assays demonstrated that S100A6 promoted PDAC cell migration and invasion. Furthermore, β-catenin shRNA inhibited the migration and invasion of PDAC cells. We confirmed that S100A6 induces PDAC cell migration and invasion via activation of β-catenin in vitro. Assessment of mRNA and protein levels revealed that S100A6 induces increased expression of β-catenin, N-cadherin and vimentin, and decreased expression of E-cadherin in PDAC cells. β-catenin shRNA also altered the expression of epithelial-mesenchymal transition (EMT)-related markers in PDAC cells. Specifically, expression of E-cadherin was increased, whereas expression of N-cadherin and vimentin was decreased. Finally, we demonstrated that S100A6 alters the expression of EMT-related markers via β-catenin activation. In conclusion, S100A6 induces EMT and promotes cell migration and invasion in a β-catenin-dependent manner. S100A6 may therefore represent a novel potential therapeutic target for the treatment of pancreatic cancer.     

Prognostic and Functional Significance of MAP4K5 in Pancreatic Cancer

Objectives

MAP4K5 plays an important role in regulating a range of cellular responses and is involved in Wnt signaling in hematopoietic cells. However, its functions in human malignancies have not been studied. The major objectives of this study are to examine the expression, functions and clinical significance of MAP4K5 in pancreatic ductal adenocarcinoma (PDAC).

Materials and Methods

The expression levels of MAP4K5, E-cadherin, vimentin, and carboxylesterase 2 (CES2) were examined by immunohistochemistry in 105 PDAC and matched non-neoplastic pancreas samples from our institution. The RNA sequencing data of 112 PDAC patients were downloaded from the TCGA data portal. Immunoblotting and RNA sequencing analysis were used to examine the expression of MAP4K5 and E-cadherin in pancreatic cancer cell lines. The effect of knockdown MAP4K5 using siRNA on the expression of CDH1 and vimentin were examined by Real-time RT-PCR in Panc-1 and AsPC-1 cells. Statistical analyses were performed using IBM SPSS Statistics.

Results

MAP4K5 protein is expressed at high levels specifically in the pancreatic ductal cells of 100% non-neoplastic pancreas samples, but is decreased or lost in 77.1% (81/105) of PDAC samples. MAP4K5-low correlated with the loss of E-cadherin (P = 0.001) and reduced CES2 expression (P = 0.002) in our patient populations. The expression levels of MAP4K5 mRNA directly correlated with the expression levels of CDH1 mRNA (R = 0.2490, P = 0.008) in the second cohort of 112 PDAC patients from The Cancer Genome Atlas (TCGA) RNA-seq dataset. Similar correlations between the expression of MAP4K5 and E-cadherin were observed both at protein and mRNA levels in multiple pancreatic cancer cell lines. Knockdown MAP4K5 led to decreased CDH1 mRNA expression in Panc-1 and AsPC-1 cells. MAP4K5-low correlated significantly with reduced overall survival and was an independent prognosticator in patients with stage II PDAC.

Conclusions

MAP4K5 expression is decreased or lost in majority of PDACs. The strong associations between low MAP4K5 expression and loss of E-cadherin, reduced CES2 expression and decreased overall survival may suggest an important role of MAP4K5 in epithelial-to-mesenchymal transition, chemotherapy resistance and tumor progression in pancreatic cancer. Targeting impaired MAP4K5 signaling may represent a new therapeutic strategy for pancreatic cancer.