- Open Access
Exosomal DNAJB11 promotes the development of pancreatic cancer by modulating the EGFR/MAPK pathway
Cellular & Molecular Biology Letters volume 27, Article number: 87 (2022)
Pancreatic ductal adenocarcinoma (PDAC) is a malignant tumor with invasive and metastatic characteristics and poor prognosis. Intracellular protein homeostasis is associated with invasion and metastasis of pancreatic cancer, but the specific molecular mechanism remains unclear. Our previous studies have revealed that DNAJB11, a key protein in protein homeostasis, is secreted by exosomes in the supernatant of dissociated pancreatic cancer cells with high metastasis. The results from transcriptome sequencing and co-immunoprecipitation (Co-IP)-based liquid chromatography with tandem mass spectrometry (LC–MS/MS) showed that depletion of DNAJB11 levels could increase HSPA5 expression and induce endoplasmic reticulum stress through the PRKR-like endoplasmic reticulum kinase signaling pathway in pancreatic cancer cells. Furthermore, exosomal DNAJB11 promoted cell development of PC cells in vitro and in vivo. In addition, exosomal DNAJB11 could regulate the expression of EGFR and activate the downstream MAPK signaling pathway. Clinical blood samples were collected to evaluate the potential of exosome DNAJB11 as a diagnostic biomarker and therapeutic target for the treatment of pancreatic cancer. This study could provide a new theoretical basis and potential molecular targets for the treatment of pancreatic cancer.
Pancreatic cancer (PC) is a lethal malignancy with low resection and high recurrence rates. In addition, PC is dormant at the early stage, with an overall 5-year survival rate of < 5% . However, since 2000, the PC incidence has increased by approximately 1% per year. Furthermore, by 2030, PC is expected to be the second leading cause of cancer deaths worldwide [2, 3]. Current treatment options, including surgical resection, conventional chemotherapy, radiotherapy, and combination therapies, remain suboptimal , probably owing to the tumor microenvironment (TME) of PC . Therefore, it is important to discuss the cancer progression mechanism and identify novel biomarkers for the early prevention and detection of PC.
Exosomes are cup-like extracellular vesicles 40–150 nm in diameter that transmit functional cargoes (proteins, mRNAs, and miRNAs) to recipient cells and regularize their functioning or pathological progression . Accumulating evidence supports that exosomes could affect the TME and play an important role in tumor progression and metastasis [7, 8]. Furthermore, the expression level of pancreatic cancer-derived exosomes is different from patients without PC . The endoplasmic reticulum (ER) is an important organelle in cells that coordinates protein folding and is the point of origin for protein secretion. The ER, like a senior engineer, maintains protein homeostasis through quality control . The unfolded protein response (UPR) is activated when the dynamic process is broken (ER stress). Cancer cells are exposed to multiply factors and trigger the UPR. Then, the UPR stimulates the transcription of genes required for tumor survival and growth [11, 12]. Heat Shock Protein Family A (Hsp70) Member 5 (HSPA5, also called BiP or GRP78) is a major regulator of ER homeostasis and controls the activation of transmembrane endoplasmic reticulum stress (ERS) sensors (e.g., IRE1, ATF6, and PERK) . HSPA5 is responsible for maintaining the “stemness” properties in PC attributed to destructive properties like metastasis . In our previous study, the DnaJ Heat Shock Protein Family (Hsp40) Member B11 (DNAJB11), a co-chaperone for HSPA5, recruited HSPA5 and other chaperones to the substrate for ERAD (ER-associated degradation) to activate UPR , was secreted to the extracellular region through exosome method. There was a significant difference in expression levels between the conditioned medium (CM) of the invasive and the less invasive pancreatic cancer cell lines . Hence, DNAJB11 may contribute to the progression and metastasis of PC. However, the molecular mechanisms underlying the efficacy of exosomal DNAJB11 in PC have not been characterized.
In this study, we revealed that exosomal DNAJB11 can enhance the invasive ability of poorly invasive PC cells, and further explained the effect of UPR in pancreatic cancer. Next, exosomal DNAJB11 contributed to cancer progression, and exosomal DNAJB11 can act as a diagnostic biomarker for PC.
Materials and methods
Ethics statement and human plasma
Ethics approval was obtained from the committee of the Ethical Committee of the Affiliated Shengjing Hospital of China Medical University with certificate number 2019PS481K for the animal studies and 2017PS46K for the blood specimens. In total, 90 plasma specimens were collected in the study. Before collecting blood specimens, no patients received other treatment like chemotherapy, radiotherapy, or surgery. All the patients and healthy individuals provided informed written consent. All animal experiments were conducted following the standards of the Ethical Committee of Shengjing Hospital.
Cell lines and cell culture
The culture conditions of the cells were as previously described . AsPC-1 (cat. no. TCHu 8), Capan-2 (cat. no. SCSP-568), BxPC-3 (cat. no. TCHu 12), and MIA Paca-2 (cat. no. SCSP-568) were purchased from the Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences (Shanghai, China). HPDE6-C7 (H6c7, cat. no. ECA001), which is a normal pancreatic duct epithelial cell line, was purchased from Kerafast (Boston, MA, USA). The cells were cultured in RPMI-1640 (Gibco-BRL, Grand Island, NY, USA) supplemented with 10% fetal bovine serum + 1% streptomycin. Cells were maintained at 37 °C in an atmosphere containing 5% CO2.
Exosomes extraction from plasma and medium, TEM, and NTA analysis
Plasma exosomes were isolated from PC, benign pancreatic disease, and healthy individuals by ExoQuick (SBI, CA, USA) following the manufacturer’s instructions . Cell culture medium from each cell line was used for exosome isolation using the differential ultracentrifugation method. The exosome pellets were resuspended in 1 ml of PBS. Purified exosomes (10 μg) from each cell resuspended in 500 μl PBS (Gibco-BRL, NY, USA) were used to treat correspondence cells according to the assay. Exosomes were prepared for identification by western blotting, transmission electron microscopy (TEM), and nanoparticle tracking analysis (NTA). Exosomes were visualized using TEM (JEM1230, JEOL, Akishima, Japan). The concentration of the samples (particles/mL) and size distribution (in nanometers) of the exosome samples were determined using ZetaView (Particle Metrix, Meerbusch, Germany) and the corresponding ZetaView 8.04.02. software.
The tissue microarray (TMA) containing 90 PC tissue samples and 56 matched adjacent paracancerous tissue samples with survival times was purchased from Outdo Biotech (Shanghai, China). The samples were collected from November 2004 to December 2008. The immunohistochemistry assay and the analysis of clinicopathological features were based on tissue microarray. Detailed clinical and pathologic information of patients is presented in Additional file 2: Table S2. The intensity of staining was divided into four scores: 0 (none), 1 (weak), 2 (moderate), and 3 (strong). The immunohistochemistry results were expressed as the percentage of positive cells and staining intensity. The cutoff score was 2.0, and a score of ≥ 2.0 was defined as positive.
Construction of DNAJB11-knockdown cells in AsPC-1 cells
Lentiviruses containing shRNA targeting human DNAJB11 (Hanheng Biotechnology, Shanghai, China) were transfected into AsPC-1 to knock down DNAJB11. Stable transfection with a lentiviral vector was performed following the manufacturer’s protocols. Knockdown efficacy was determined by quantitative polymerase chain reaction (qPCR) and western blot.
Sequencing was performed by Amogene (Xiamen, China). First, total RNA was quantified and purified. Three micrograms of RNA of each sample was used for RNA-seq library preparation. Second, RNA-seq libraries were conducted using the NEBNext Ultra II RNA Library Prep Kit (NEB, E7760, Illumina, San Diego, CA, USA) according to the manufacturer’s protocol. In brief, poly-T oligo-attached magnetic beads were used to purify mRNA. First- and second-strand cDNA synthesis was performed. The library fragments were analyzed by AMPure XP system (Beckman Coulter, Brea, CA, USA). Next, PCR was conducted with universal PCR primers, index (X) primer, and Phusion High-Fidelity DNA polymerase. Finally, PCR products were purified, and RNA integrity was assessed on an Agilent 2100 Bioanalyser (Agilent, Santa Clara, CA, USA).
For the quantification of gene expression level and differential expression analysis, the read numbers mapped to each gene were counted using Cufflinks V2.2.1.(http://cole-trapnell-lab.github.io/cufflinks/) . Next, the fragments per kilobase of transcript per million mapped reads of each gene was calculated on the basis of gene length and its mapped read counts. Finally, differentially expressed genes (DEGs) were analyzed using Cufflinks. The Benjamini–Hochberg correction was performed to adjust the P-values. Corrected P-value < 0.05 and absolute fold change > 1.5 were selected as the threshold of differential expression.
Co-immunoprecipitation (Co-IP) was performed following the instructions on the cross-linking CoIP kit (Thermo Fisher, Waltham, MA, USA). First, the cells were lysed with 500 µL RIPA buffer containing protease and phosphatase inhibitor on ice for 30 min. Next, the cell lysate supernatant was collected by centrifugation and incubated with anti-DNAJB11 antibody overnight at 4 °C. Next, protein A agarose was added, and cell lysates were further mixed and incubated for 4 h, after being centrifuged at 3,000g for 3 min and washed thrice with lysis buffer. Finally, the bound proteins were eluted by boiling for 5 min in 2× loading buffer and separated by SDS-PAGE followed by liquid chromatography and tandem mass spectrometry (LC–MS/MS). LC–MS/MS was performed as previously described with minor modifications . Protein digestion was performed following the standard protocol, and 5 μL of each sample was loaded for LC–MS/MS. The Eksigent nanoLC 415 high-performance liquid chromatograph (Eksigent, Redwood, CA, USA) with the 5600 + TripleTOF (AB Sciex, Framingham, CA, USA) was used to analyze protein fractions. Mobile phases consisted of A (98% H2O and 2% acetonitrile with 0.1% formic acid) and B (2% H2O and 98% acetonitrile with 0.1% formic acid). Peptides were loaded onto a C18 trap column (100 μm i.d. × 2 cm) and separated by a capillary separation column (75 μm i.d. × 15 cm) at 0.3 μL/min with a 90 min gradient. The MS data were acquired using Data Independent Acquisition (DIA) mode, 30 data-dependent MS/MS scans per full scan, rolling collision energy, dynamic exclusion (exclusion duration 15 s), MS/MS scan range of 100–1500 m/z, and scan time of 50 ms. The ProteinPilot software (Version 4.5, ABSciex) was used for data analysis. The Unused Protein Score was used to measure protein confidence.
Measurement of ECAR
Glycolytic ability was measured using a Seahorse XF96 analyzer (Seahorse Biosciences, Agilent) according to the manufacturer’s instructions. Twenty-thousand cells per well were plated in 96-well XF microplates and cultured for 6 h. All measurements were recorded at set time intervals and normalized to total protein content. Extracellular acidification rate (ECAR) after oligomycin treatment indicates glycolytic capacity. ECAR was calculated using the Wave software.
Western blot, immunohistochemical assay
Western blotting was performed as described previously . Samples of equivalent total protein (20 μg) were loaded. Primary antibodies against CD63, CD9 (Santa Cruz Biotechnology, Santa Cruz, CA, USA), c-MYC, p-EGFR, p-MEK, p-ERK(Cell Signaling Technology, Danvers, MA, USA), GAPDH, DNAJB11, HSPA5, ATF6, IRE1, XBP1, PERK, ATF4, EGFR, Raf-1, MEK, and ERK1/2 (ProteinTech Group, Rosemont, IL, USA) were used. Immunohistochemical (IHC) assay was performed following a previously described procedure .
Colony formation, invasion and migration assays, cell apoptosis
A colony formation assay was performed to assess cell proliferation. Transwell and wound healing assays evaluated cell invasion and migration. The assays were performed as previously described . Apoptosis assay was performed using an apoptosis detection kit according to the manufacturer’s protocol. The cells were cultured in six-well plates at a density of 1 × 105 cells per well. Following digestion with trypsin and washing twice with PBS, the cells pellet were incubated with 5 µl Annexin V-FITC and 5 µl PI for 15 min at room temperature in the dark. The cell suspensions were analyzed by flow cytometry (FACScan; BD Biosciences) and FlowJo software (version 10.0; FlowJo LLC). The total amount of apoptosis was the results of early apoptosis plus late apoptosis.
The subcutaneous tumor model was used to evaluate cell proliferation in vivo according to a previously described method . First, 18 4-week-old female BALB/c nude mice were obtained from Huafukang Biotechnology Co (Beijing, China). Next, 1 × 106 cells from each group in 0.1 mL PBS were injected into the right axillary region of each mouse to generate 60 mm3 tumors. The tumor diameter was measured twice a week. Volume was calculated as volume = 0.5 × L × W2, where L is the long axis of the tumor and W is the short axis of the tumor.
Enzyme-linked immunosorbent assay (ELISA)
The DNAJB11 protein level of plasma-derived exosomes was determined using the human DNAJB11 ELISA kit (Cloud-Clone, Wuhan, China) following the manufacturer’s protocol.
The Gene Expression Profiling Interactive Analysis (GEPIA) database (http://gepia.cancer-pku.cn) provided the DEGs based on the integrated RNA-seq dataset from Genotype-Tissue Expression (GTEx) databases and The Cancer Genome Atlas (TCGA) . UALCAN (http://ualcan.path.uab.edu/index.html)  was used to assist DNAJB11 expression across TCGA cancers. The Kaplan–Meier (KM) plotter database included the genes that were associated with overall survival (OS) and relapse-free survival (RFS) of patients’ information . The R packages “clusterProfiler” and “ggplot2” were used to run and visualize Gene Ontology (GO) analysis, KEGG pathways, and Gene Set Enrichment Analysis (GSEA). GO terms with corrected P-value less than 0.05 were significantly enriched by differentially expressed genes. Multivariate Cox analysis was performed to identify the independent prognostic factor on the basis of the common clinicopathological data.
The R version 3.6.1 and RStudio software were used for statistical analyses and figure outputs. Other statistical analyses and graphics were performed and generated using GraphPad Prism 9.0. The data were analyzed using two-tailed Student’s t-test and one-way analysis of variance (ANOVA). P-value < 0.05 was considered statistically significant in this study.
Exosomal DNAJB11 from plasma obtained from patients with PC was correlated with pancreatic cancer progression
Plasma exosomes from patients with malignant pancreatic cancer (MP-exo), patients with benign pancreatic disease (BP-exo), and healthy people (N-exo) were isolated using an exosome isolation kit. The morphology and size of the exosomes were assessed using nanoparticle tracking analysis (NTA) and transmission electron microscopy (TEM) (Fig. 1A). In this study, patients with PC had higher protein concentrations than healthy individuals (Fig. 1B). Furthermore, the plasma exosomal markers CD63, CD9, and exosomal DNAJB11 were observed using western blotting (Fig. 1C). Moreover, the sizes of different groups showed no significant differences (Fig. 1D). The corresponding clinical information of the research group is provided in Additional file 1: Table S1. We collected plasma samples from patients with malignant pancreatic cancer (MP, n = 31), benign pancreatic disease (BP, n = 28; biliary pancreatitis ; acute pancreatitis ; chronic pancreatitis ; pancreatic cystadenoma ; and pancreatic islet cell tumor), and normal controls (N, n = 31). The samples were used to assess whether exosomal DNAJB11 could be a diagnostic biomarker for pancreatic cancer. The result showed that exosomal DNAJB11 was significantly higher in the MP group than the BP group (P < 0.01) and N (P < 0.0001) group (Fig. 1E). Furthermore, the receiver operating characteristic curve was obtained to assess the diagnostic role of exosomal DNAJB11 for pancreatic cancer. The areas under the curve (AUC) indicated the diagnostic ability between the MP and N groups (AUC 0.8044), the MP and BP groups (AUC 0.735), and the BP and N groups (AUC 0.7615) (Fig. 1F–H). Thus, exosomal DNAJB11 is a potential diagnostic biomarker of PC.
Survival analysis and clinical features of DNAJB11 expression in PC tissue
Multiple bioinformatics tools were performed for analyzing the clinical data on DNAJB11 from the GTEx database and TCGA. DNAJB11 was markedly expressed in most cancers, including PC by UALCAN and GEPIA (Fig. 2A, B). KM plotter analysis was used to analyze the relationship between DNAJB11 expression and survival; 74 patients with a high expression of DNAJB11 indicated poor overall survival (OS) compared with those with low expression (P = 0.026, Fig. 2C), but relapse-free survival (RFS) was not significantly different (P = 0.13) (Fig. 2D). To characterize DNAJB11 expression in PC, we purchased a PC tissue microarray containing 90 PC tissue samples and 56 matched adjacent paracancerous tissue samples. Integral images of immunohistochemical staining of the TMA with an anti-DNAJB11 were captured. Higher levels of DNAJB11 were observed in PC tissue, compared with the matched adjacent tissues. Corresponding staining of pancreatic cancer and adjacent tissues are shown in Fig. 2E. The immunohistochemical (IHC) score of every sample in the TMA was determined (Additional file 2: Table S2). DNAJB11 expression was significantly higher in PC compared with the adjacent tissues (Table 1). Subsequently, the effects of DNAJB11 expression on the pathological characteristics and prognosis of patients with PC were analyzed. The survival outcome of the patients was related to the DNAJB11 levels (Fig. 2F). Patients with high expression of DNAJB11 showed poorer OS than those with low expression of DNAJB11 (P < 0.009). To study the prognosis-related factors, we used multivariate Cox regression. The results showed that the expression of DNAJB11 was an independent risk factor for PC (Fig. 2G). The high DNAJB11 expression levels correlate with tumor size (Table 2). The results showed that DNAJB11 was overexpressed in PC and closely associated with the survival of patients.
Exosomes containing DNAJB11 promote PC cell proliferation, migration, and invasion
To characterize the functional role of exosomal DNAJB11 during PC progression, first, DNAJB11 expression was measured in four human PC cell lines (AsPC-1, BxPC-3, MIA paca-2, and Capan-2) and a normal pancreatic cell line (HPDE6-C7) (Fig. 3A). Compared with the other cell lines, the expression of DNAJB11 was higher in the AsPC-1. Therefore, AsPC-1 cell lines were chosen for further assays, and the efficacy of DNAJB11 knockdown was determined by western blot analysis. We extracted the AsPC-1, shDNAJB AsPC-1, and Capan-2 cell exosomes, respectively. Furthermore, the shDNAJB AsPC-1 and Capan-2 cells were treated with AsPC-1-derived exosomes (50 µg/ml). The presence of exosomal markers and the corresponding expression of exosomal DNAJB11 in the different groups were determined by western blot analysis (Fig. 3B). The results indicated that DNAJB11 could be delivered into other cells through the exosomes. Further, we verified the effects of exosomal DNAJB11 in PC progression by comparing the functions of wild-type (WT) and DNAJB11-knockdown exosomes from AsPC-1 cells in vitro. The shDNAJB11 group showed poorer proliferation, migration, and invasion abilities than the control group. In addition, the abilities were recovered after treatment with AsPC-1-exo. Furthermore, exosomes containing DNAJB11 significantly stimulated proliferation, migration, and invasion compared with shDNAJB11-exo (Fig. 3C–E). shDNAJB11 group exhibited a significant increase in cellular apoptosis rate compared with AsPC-1. However, exosomal DNAJB11 significantly decreased the cellular apoptosis rates compared with Capan-2 (Fig. 3F).
DNAJB11 mediates transcriptomic changes in PC cells
To characterize DNAJB11-mediated transcriptomic changes in PC cells, transcriptomic sequencing (RNA-seq) was conducted after DNAJB11 knockdown. On analysis, a total of 24,760 expressed genes were detected, 837 genes were upregulated, and 1499 genes were downregulated (Additional file 3: Table S3). Gene Ontology enrichment analysis was performed on all the DEGs, and detailed information is provided in Additional file 2: Table S2. Volcano and clustering heat maps were performed on all the detected genes (Additional file 8: Fig S1B). GO and KEGG analysis was used to analyze the biological processes (BP), cellular components (CC), molecular functions (MF), and enriched pathways (Fig. 4A–F). The top ten upregulated and downregulated enriched terms of DEGs in the BP, CC, and MF terms are shown in Fig. 4C, D.
Furthermore, 32 genes were enriched in ERS in the upregulated group, whereas 59 were enriched in the small-molecule catabolic process. In addition, KEGG analysis of DEGs enriched the pathway involved in cell metabolism (Fig. 4E, F). The results illustrate that DNAJB11 has an important role in cell metabolism. We performed the GSEA analyses based on the DEGs to further investigate the underlying molecular mechanisms. The DEGs were enriched for 50 hallmark gene sets (Additional file 4: Table S4), correlating with proliferation (HALLMARK_E2F_TARGETS) and interferon response (HALLMARK_ INTERFERON_GAMMA). The result indicated that the proliferation signature was linked with the gain of the inflammatory interferon response in the PC tissue . Notably, the glycolysis gene set was ranked higher in the upregulated group (Fig. 4G). The Seahorse array was performed to analyze metabolism. Consistent with the GSEA results, DNAJB11 was involved in the energy metabolism process. Thus, inhibiting the DNAJB11 expression could reduce glycolysis (Fig. 4H, I, Additional file 5: Table S5). In addition, the TNFA-NFKB signaling pathway and RAS signaling pathway are presented in upregulated and downregulated groups. The observations suggested that DNAJB11 could regulate the progression of PC by increasing the RAS signaling pathway.
CoIP-base proteomics analysis of DNAJB11
To understand the DNAJB11 molecular mechanisms in PC cells, we performed a Co-IP-based proteomics assay to discover its interacting proteins. IgG served as a negative control for Co-IP. Two protein samples obtained from previous experiments were qualitatively analyzed by LC–MS/MS. On the basis of the search results, a Venn diagram was drawn (Fig. 5A). There were 157 proteins in the DNAJB11-IP group, and 32 in the IgG group (Table 3). Furthermore, 24 proteins were identified in both groups (Fig. 5A, Additional file 6: Table S6). In addition, 133 proteins in the DNAJB11-IP group were analyzed by five online bioinformatics databases (Fig. 5B–F, Additional file 7: Table S7). WoLF PSORT (https://wolfpsort.hgc.jp/)  was used to predict the localization of proteins (Fig. 5B). The results showed that DNAJB11-IP mostly exited in the cytosol and nucleus. The GO and KEGG analysis showed proteins enriched in translation, transport, and catabolism (Fig. 5C) as well as cellular process and binding, and structural molecule activity (Fig. 5D). Therefore, we ran InterProScan  to analyze domain architectures. Results showed that there were 14 proteins that include P-loop-containing nucleoside triphosphate hydrolase (Fig. 5E). The genes were placed into clusters of orthologous groups (COG) categories according to their functional annotation. Thirty-four proteins were enriched in the COG category for translation, ribosomal structure, and biogenesis (COG J), and 14 proteins were in the posttranslational modification, protein turnover, and chaperones (COG O) (Fig. 5F).
Blocking DNAJB11 induced ER stress and suppressed the EGFR signaling pathway
DNAJB11 plays an important role in the UPR. Our data showed that blocking DNAJB11 could induce ER stress by PERK-mediated UPR (Fig. 6A). The functional network integration was performed using GeneMANIA, and DNAJB11 was high compared with HSPA5 (Fig. 6B) . The results were validated by western blot analysis (Fig. 6C). In the shDNAJB11 group, HSPA5, PERK, and ATF4 expression levels increased, but not ATF6, XBP1, and IRE1. Thus, blocking DNAJB11 activates ER stress through PERK-ATF4 of the UPR. Therefore, the mechanism by which exosomal DNAJB11 promotes cancer cell growth was verified. However, compared with AsPC-1 cells, we did not find a significant difference in the DNAJB11/PERK signaling pathway in shDNAJB11-AsPC-1 and Capan-2 cells after incubation with exosomal DNAJB11. This is probably because overexpression was more complex. Therefore, the UPR was regulated by another gene concurrently and not by DNAJB11 only. In addition, GSEA analyses revealed that DNAJB11 depletion could increase EGFR expression (Fig. 6D). Furthermore, our Co-IP-based MS data showed that DNAJB11 interacts with EGFR. UALCAN was also performed to analyze the correlation between DNAJB11 and EGFR (P < 0.001) (Fig. 6E). Combining the RNA-seq results, we hypothesized that DNAJB11 induced EGFR/RAF/MAPK signaling to promote cancer cell growth. Indeed, blocking DNAJB11 suppresses the RAF/MPAK signaling by reducing the EGFR expression level. Moreover, EGFR and p-EGFR protein levels were elevated in shDNAJB11-AsPC-1 and Capan-2 cells after the uptake of DNAJB11-exo (Fig. 6F). Significant differences were also seen in p-MEK and p-ERK protein expression in each experimental group compared with the corresponding control groups. Our data reveal that DNAJB11 induces the activation of the EGFR signaling pathway and suppresses the PERK-ATF4 ER stress signaling pathways.
Exosomal DNAJB11 promotes pancreatic cancer growth in a nude mouse model
To further investigate the effect of exosomal DNAJB11 in vivo, subcutaneous tumors were generated in nude mice with 1 × 107 Capan-2 cells. The nude mice treated with normal AsPC-1-exo showed rapid and larger tumor growth than those treated with shDNAJB11-exo (Fig. 7A–D). Furthermore, immunohistochemical assays showed that DNAJB11 expression was higher than in those injected with shDNAJB11-exo, indicating that exosomes enriched in DNAJB11 increased tumor growth (Fig. 7E).
Pancreatic cancer is an aggressive malignancy worldwide. Patients diagnosed with early-stage pancreatic cancer should receive therapy for a good clinical outcome [29, 30]. Early detection of PC is important in the prognosis and survival of patients with PC. In this study, we revealed that not only was the expression of DNAJB11 in the PC tissue high, but also exosomal DNAJB11 was higher in plasma exosomes from patients with PC than healthy individuals. Exosomal DNAJB11 and DNAJB11 may play an important role in PC development.
DNAJB11, a co-chaperone to HSPA5, is a soluble resident ER glycoprotein . The N terminus of DNAJB11 was expressed in the cytosol and the nucleus, and the C terminus was targeted to the ER lumen and membrane-anchored [32, 33]. Thus, the conformation is vital in the process of unfolded or misfolded protein transfer. Besides, the N terminus includes the J domain, which provides a dynamic surface to form functional interactions with the HSPA5 . Furthermore, the binding of DNAJB11 and HSPA5 could quickly protect the unfolded protein from aggregation, and DNAJB11 and HSPA5 are more frequently dispersed than other chaperones (e.g., ERdj4, ERdj5). Previous studies have reported that DNAJB11 is highly upregulated by unfolded secretory protein stress in the ER and could regulate a partial HSPA5’s gating activity to stimulate ATPase activity up to fivefold [36,37,38]. In addition, HSPA5 activates the UPR by increasing the probability that ERS sensors recognize unfolded proteins . The UPR plays an important role in cancer progression [40, 41]. Previous studies established that DNAJB11 promotes cancer growth in hepatocellular carcinoma . However, the interaction between DNAJB11 and UPR on PC has not been studied.
To elucidate further molecular mechanisms of DNAJB11 in patients with PC, IHC and TCGA database analysis of patients with PC were performed. DNAJB11 expression was positively correlated with advanced status and poor survival outcomes. Furthermore, DNAJB11 was associated with EGFR. As a result, DNAJB11 probably upregulated EGFR expression and directly promoted angiogenesis in PC tissues. In addition, there was a significant increase in DNAJB11 in the plasma exosomes from patients with PC. Hence, DNAJB11 overexpression contributed to the development and progression of PC.
In our study, we found that DNAJB11 reduced the expression of HSPA5 in pancreatic cancer cells. In shDNAJB11 cells, there was a corresponding increase in HSPA5 expression. However, in the Co-IP MS experiment, the direct interaction between DNAJB11 and HSPA5 was confirmed, as previously reported by other studies . Furthermore, DNAJB11 expression enhanced ATP hydrolysis, thus inducing a conformational change in occupied HSPA5, which locks the unfolded protein, and DNAJB11 release from HSPA5 complexes well before folding is complete [37, 43]. DNAJB11 has “handcuff”-like function with respect to HSPA5. When DNAJB11 was overexpressed, occupied HSPA5 could not bind to other proteins. However, the binding of HSPA5 and DNAJB11 was more resistant to mutational disruption than previously identified for other chaperones. In shDNAJB11 cells, UPR activity was increased in our study. The depletion of DNAJB11 increased the possibility of HSPA5 interacting with the unfolded protein and inducing UPR to protect against apoptosis from cancer cells.
It is known that the ER is the quality controller of proteins, and long-term ERS produces large amounts of unfolded protein. The unfolded proteins can be refolded to generate functional proteins or, when they accumulate improperly folded proteins, can exceed the threshold and be degraded by proteasome or trigger cell death signals by activating UPR, a phenomenon also called imbalance of ER homeostasis . Over decades, the potential role of DNAJB11 in human chronic diseases has been investigated [45, 46]. Cell survival under chronic ERS that develops hepatic carcinomas is a symptom of aggressive cancer . Recently, it has been reported that DNAJB11 upregulation decreases the potential antitumor activity of celecoxib . ER-associated degradation (ERAD) enables cancer cells to tolerate toxic glycoprotein stress, thereby helping them to survive . Pancreatic cancer cells maintain homeostasis by inducing the repeated use of abnormally folded proteins in cells, contributing to the progression of tumor tissues. We demonstrated that DNAJB11 is involved in ERAD by increasing ATF4 expression. ATF4 promotes pancreatic cancer progression and the mechanism of gemcitabine resistance . Thus, silencing DNAJB11 expression triggers cellular apoptosis.
In addition, we compared the metabolism of shDNAJB11 and that of the parent cell line. The knockdown group showed low glycolysis. In addition, hypoxia was another cellular response by ERS noted. In this study, a comprehensive analysis of our RNA-seq and Co-IP-based MS data, seven genes (ALPP, CD55, CGN, GSN, ITGAV, TFRC, and VIL1) were found in both groups. Transferrin receptor (TFRC) is important for iron uptake and is related to ferroptosis, a non-apoptotic form of cell death. The complex relationship between DNAJB11 and TFRC in pancreatic cancer needs to be elucidated in future studies.
Therefore, we provided new insight into DNAJB11, identifying it as an exosome with important extracellular biological function. DNAJB11 in exosomes was protected from degradation. Like a sponge, it took up extracellular unfolded proteins and performed its function. However, the mechanisms of exosome–recipient cell interaction are unclear. In this study, exosomal protein altered the ability of invasion and metastasis in receipt cells. Previous studies report a similar molecular mechanism, that DNAJB11 is directly co-secreted with unfolded proteins to the extracellular environment and maintains proteostasis [51, 52]. Exosomal DNAJB11 activated the EGFR/MAPK signaling pathway in targeted PC cells. DNAJB11 is expressed in the upregulated differential expression group both intracellularly and extracellularly. Therefore, it is conceivable that interference with DNAJB11 can improve the tumor cells and their surrounding microenvironment simultaneously. This study also has certain limitations. For example, there is no description of how exosomes are transported into another cell. The different types of exosome were used to treat corresponding cells, and differences in function were assessed. The differences in function observed may be due to differential expression of DNAJB11. Hence, exosomal DNAJB11 is a potential biomarker for PC development.
In summary, exosomal DNAJB11 is crucial in maintaining extracellular signaling, altered metabolic behavior, and adhesion properties of pancreatic cancer cells. Exosomal DNAJB11 may be a potential biomarker in PC development and progression. However, DNAJB11 is a crucial regulator upon ER stress and negatively controls the UPR signaling pathway. DNAJB11 promoted cancer development through the EGFR/MAPK signaling pathway, providing new insight into the development of pancreatic cancer. DNAJB11 inhibitors have been proposed as targets to impair the survival of cancer cells.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
Unfolded protein response
DnaJ Heat Shock Protein Family (Hsp40) Member B11
Heat Shock Protein Family A (Hsp70) Member 5
Gene Set Enrichment Analysis
The Cancer Genome Atlas
Transmission electron microscopy
Differentially expressed genes
Liquid chromatography with tandem mass spectrometry
Extracellular acidification rates
Gene expression profiling interactive analysis
Nanoparticle tracking analysis
Cancer Genome Atlas Research Network. Integrated genomic characterization of pancreatic ductal adenocarcinoma. Cancer Cell. 2017;32(2):185-203.e13.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71(1):7–33.
Rahib L, Smith BD, Aizenberg R, Rosenzweig AB, Fleshman JM, Matrisian LM. Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the United States. Cancer Res. 2014;74(11):2913–21.
Hidalgo M. Pancreatic cancer. N Engl J Med. 2010;362(17):1605–17.
Truong LH, Pauklin S. Pancreatic cancer microenvironment and cellular composition: current understandings and therapeutic approaches. Cancers (Basel). 2021. https://doi.org/10.3390/cancers13195028.
Liao W, Du Y, Zhang C, Pan F, Yao Y, Zhang T, et al. Exosomes: the next generation of endogenous nanomaterials for advanced drug delivery and therapy. Acta Biomater. 2019;86:1–14.
Chang CH, Pauklin S. Extracellular vesicles in pancreatic cancer progression and therapies. Cell Death Dis. 2021;12(11):973.
Weidle UH, Birzele F, Kollmorgen G, Ruger R. The multiple roles of exosomes in metastasis. Cancer Genomics Proteomics. 2017;14(1):1–15.
Ariston Gabriel AN, Wang F, Jiao Q, Yvette U, Yang X, Al-Ameri SA, et al. The involvement of exosomes in the diagnosis and treatment of pancreatic cancer. Mol Cancer. 2020;19(1):132.
Needham PG, Guerriero CJ, Brodsky JL. Chaperoning endoplasmic reticulum-associated degradation (ERAD) and protein conformational diseases. Cold Spring Harb Perspect Biol. 2019. https://doi.org/10.1101/cshperspect.a033928.
Koumenis C, Wouters BG. “Translating” tumor hypoxia: unfolded protein response (UPR)-dependent and UPR-independent pathways. Mol Cancer Res. 2006;4(7):423–36.
Faitova J, Krekac D, Hrstka R, Vojtesek B. Endoplasmic reticulum stress and apoptosis. Cell Mol Biol Lett. 2006;11(4):488–505.
Huang W, Liao CC, Han Y, Lv J, Lei M, Li Y, et al. Co-activation of Akt, Nrf2, and NF-κB signals under UPRER in torpid Myotis ricketti bats for survival. Commun Biol. 2020;3(1):658.
Dauer P, Sharma NS, Gupta VK, Durden B, Hadad R, Banerjee S, et al. ER stress sensor, glucose regulatory protein 78 (GRP78) regulates redox status in pancreatic cancer thereby maintaining “stemness.” Cell Death Dis. 2019;10(2):132.
Shen Y, Meunier L, Hendershot LM. Identification and characterization of a novel endoplasmic reticulum (ER) DnaJ homologue, which stimulates ATPase activity of BiP in vitro and is induced by ER stress. J Biol Chem. 2002;277(18):15947–56.
Liu P, Weng Y, Sui Z, Wu Y, Meng X, Wu M, et al. Quantitative secretomic analysis of pancreatic cancer cells in serum-containing conditioned medium. Sci Rep. 2016;6:37606.
Liu P, Kong L, Jin H, Wu Y, Tan X, Song B. Differential secretome of pancreatic cancer cells in serum-containing conditioned medium reveals CCT8 as a new biomarker of pancreatic cancer invasion and metastasis. Cancer Cell Int. 2019;19:262.
Wu M, Tan X, Liu P, Yang Y, Huang Y, Liu X, et al. Role of exosomal microRNA-125b-5p in conferring the metastatic phenotype among pancreatic cancer cells with different potential of metastasis. Life Sci. 2020;255: 117857.
Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013;31(1):46–53.
Jin H, Liu P, Wu Y, Meng X, Wu M, Han J, et al. Exosomal zinc transporter ZIP4 promotes cancer growth and is a novel diagnostic biomarker for pancreatic cancer. Cancer Sci. 2018;109(9):2946–56.
Kong L, Liu P, Fei X, Wu T, Wang Z, Zhang B, et al. A prognostic prediction model developed based on four CpG sites and weighted correlation network analysis identified DNAJB1 as a novel biomarker for pancreatic cancer. Front Oncol. 2020;10:1716.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.
Nagy A, Lanczky A, Menyhart O, Gyorffy B. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci Rep. 2018;8(1):9227.
Ligorio M, Sil S, Malagon-Lopez J, Nieman LT, Misale S, Di Pilato M, et al. Stromal microenvironment shapes the intratumoral architecture of pancreatic cancer. Cell. 2019;178(1):160-75.e27.
Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, et al. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35:W585–7.
Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, et al. InterProScan: protein domains identifier. Nucleic Acids Res. 2005;33:W116–20.
Mostafavi S, Ray D, Warde-Farley D, Grouios C, Morris Q. GeneMANIA: a real-time multiple association network integration algorithm for predicting gene function. Genome Biol. 2008;9(Suppl 1):S4.
Poruk KE, Firpo MA, Adler DG, Mulvihill SJ. Screening for pancreatic cancer: why, how, and who? Ann Surg. 2013;257(1):17–26.
Pereira SP, Oldfield L, Ney A, Hart PA, Keane MG, Pandol SJ, et al. Early detection of pancreatic cancer. Lancet Gastroenterol Hepatol. 2020;5(7):698–710.
Jin Y, Zhuang M, Hendershot LM. ERdj3, a luminal ER DnaJ homologue, binds directly to unfolded proteins in the mammalian ER: identification of critical residues. Biochemistry. 2009;48(1):41–9.
Lau PP, Villanueva H, Kobayashi K, Nakamuta M, Chang BH, Chan L. A DnaJ protein, apobec-1-binding protein-2, modulates apolipoprotein B mRNA editing. J Biol Chem. 2001;276(49):46445–52.
Yu M, Haslam RH, Haslam DB. HEDJ, an Hsp40 co-chaperone localized to the endoplasmic reticulum of human cells. J Biol Chem. 2000;275(32):24984–92.
Chen KC, Qu S, Chowdhury S, Noxon IC, Schonhoft JD, Plate L, et al. The endoplasmic reticulum HSP40 co-chaperone ERdj3/DNAJB11 assembles and functions as a tetramer. EMBO J. 2017;36(15):2296–309.
Behnke J, Mann MJ, Scruggs FL, Feige MJ, Hendershot LM. Members of the Hsp70 family recognize distinct types of sequences to execute ER quality control. Mol Cell. 2016;63(5):739–52.
Awad W, Estrada I, Shen Y, Hendershot LM. BiP mutants that are unable to interact with endoplasmic reticulum DnaJ proteins provide insights into interdomain interactions in BiP. Proc Natl Acad Sci USA. 2008;105(4):1164–9.
Shen Y, Hendershot LM. ERdj3, a stress-inducible endoplasmic reticulum DnaJ homologue, serves as a cofactor for BiP’s interactions with unfolded substrates. Mol Biol Cell. 2005;16(1):40–50.
Guo F, Snapp EL. ERdj3 regulates BiP occupancy in living cells. J Cell Sci. 2013;126(Pt 6):1429–39.
Onn A, Ron D. Modeling the endoplasmic reticulum unfolded protein response. Nat Struct Mol Biol. 2010;17(8):924–5.
Urra H, Dufey E, Avril T, Chevet E, Hetz C. Endoplasmic reticulum stress and the hallmarks of cancer. Trends Cancer. 2016;2(5):252–62.
Bartoszewska S, Collawn JF. Unfolded protein response (UPR) integrated signaling networks determine cell fate during hypoxia. Cell Mol Biol Lett. 2020;25:18.
Pan J, Cao D, Gong J. The endoplasmic reticulum co-chaperone ERdj3/DNAJB11 promotes hepatocellular carcinoma progression through suppressing AATZ degradation. Future Oncol. 2018;14(29):3001–13.
Jin Y, Awad W, Petrova K, Hendershot LM. Regulated release of ERdj3 from unfolded proteins by BiP. EMBO J. 2008;27(21):2873–82.
Chen X, Iliopoulos D, Zhang Q, Tang Q, Greenblatt MB, Hatziapostolou M, et al. XBP1 promotes triple-negative breast cancer by controlling the HIF1α pathway. Nature. 2014;508(7494):103–7.
Jordan P, Arrondel C, Bessieres B, Tessier A, Attie-Bitach T, Guterman S, et al. Bi-allelic pathogenic variations in DNAJB11 cause Ivemark II syndrome, a renal–hepatic–pancreatic dysplasia. Kidney Int. 2021;99(2):405–9.
Huynh VT, Audrezet MP, Sayer JA, Ong AC, Lefevre S, Le Brun V, et al. Clinical spectrum, prognosis and estimated prevalence of DNAJB11-kidney disease. Kidney Int. 2020;98(2):476–87.
Huang J, Viswakarma N, Yu S, Jia Y, Bai L, Vluggens A, et al. Progressive endoplasmic reticulum stress contributes to hepatocarcinogenesis in fatty acyl-CoA oxidase 1-deficient mice. Am J Pathol. 2011;179(2):703–13.
Tsutsumi S, Namba T, Tanaka KI, Arai Y, Ishihara T, Aburaya M, et al. Celecoxib upregulates endoplasmic reticulum chaperones that inhibit celecoxib-induced apoptosis in human gastric cells. Oncogene. 2006;25(7):1018–29.
Tax G, Lia A, Santino A, Roversi P. Modulation of ERQC and ERAD: a broad-spectrum spanner in the works of cancer cells? J Oncol. 2019;2019:8384913.
Wei L, Lin Q, Lu Y, Li G, Huang L, Fu Z, et al. Cancer-associated fibroblasts-mediated ATF4 expression promotes malignancy and gemcitabine resistance in pancreatic cancer via the TGF-β1/SMAD2/3 pathway and ABCC1 transactivation. Cell Death Dis. 2021;12(4):334.
Buck TM, Brodsky JL. Escaping the endoplasmic reticulum: why does a molecular chaperone leave home for greener pastures? EMBO J. 2015;34(1):1–3.
Genereux JC, Qu S, Zhou M, Ryno LM, Wang S, Shoulders MD, et al. Unfolded protein response-induced ERdj3 secretion links ER stress to extracellular proteostasis. EMBO J. 2015;34(1):4–19.
The present study was funded by the National Natural Science Foundation of China (grant number 81902953), Liaoning Provincial Natural Science Foundation Project (grant number 180530068), the Outstanding Young Doctor Fund of China Medical University (grant number QGZD2018050), and 345 Talent Project of Shengjing Hospital of China Medical University.
Ethics approval and consent to participate
This research protocol was approved by the Ethics Committee of the Affiliated Shengjing Hospital of China Medical University with certificate number 2019PS481K for the animal studies and 2017PS46K for the human studies..
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.
Additional file 1:
Table S1. Clinical information of the 90 patients run through the ELISA assay.
Additional file 2: Table S2.
Clinical and pathologic information of patients in TMA.
Additional file 3: Table S3.
A complete list of all genes identified in the transcriptome sequencing.
Additional file 4: Table S4.
Enrichment analysis of DEGs.
Additional file 5:
Table S5. The results of the Seahorse assay.
Additional file 6:
Table S6. List of proteins identified in CoIP-based LC–MS/MS.
Additional file 7:
Table S7. Enrichment analysis of proteins identified in CoIP-seq.
Additional file 8:
Figure S1. The flow chart of transcriptome sequencing.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Liu, P., Zu, F., Chen, H. et al. Exosomal DNAJB11 promotes the development of pancreatic cancer by modulating the EGFR/MAPK pathway. Cell Mol Biol Lett 27, 87 (2022). https://doi.org/10.1186/s11658-022-00390-0
- Pancreatic cancer
- Exosomal protein
- Signal transduction
- Early diagnosis