1. Introduction
Pancreatic cancer has been identified as one of the most aggressive malignant tumors in the human alimentary system and the 5-year survival rate of patients with pancreatic cancer has been found to be approximately 11% [
1]. Fewer than 20% of the patients with pancreatic cancer are candidates for curative resections of tumors, whereas most of the patients are found to have unresectable disease or distant metastasis, of which liver metastasis account for 50% [
2]. Even among patients who undergo curative operations, 80% of cases experience recurrence within 2 years and the liver is the most frequently involved organ postoperatively [
3,
4]. Unfortunately, applicable therapies such as systemic chemotherapy, stereotactic body radiation therapy, transarterial chemoembolization, radiofrequency ablation, cryoablation and ethanol injection, have achieved unsatisfactory effects on liver metastasis in current clinical practice [
5]. There is an imperative need to respond this challenge by developing innovative targets in treatments of liver metastasis in pancreatic cancer.
In previous studies, multiple novel molecular targets with respect to liver metastasis had been found in human malignancies [
6,
7,
8,
9,
10], which improved patients’ outcomes and highlighted the roles of molecular targets in liver metastasis [
1]. Clinical trials combining different molecular inhibitors with chemotherapies in patients with local advanced or metastatic disease of pancreatic cancer are ongoing, although initial results are discouraging [
11]. Due to transitions in the tumor microenvironment (TME) from primary tumors to liver metastasis and corresponding variations in the molecular profiling, distinctive molecular targets found in liver metastasis may provide active areas of research and clinical benefits. Several studies have made some breakthrough in these fields [
12,
13], but no doubt more attempts will be expected.
The crosstalk between tumors cells and the TME, has been posited to induce dynamic variations in molecular profiling, and therefore we attempted to identify differentially expressing genes (DEGs) in liver metastasis compared to primary tumors of pancreatic cancer. Out of 29 candidate genes, proprotein convertase subtilisin/kexin type 6 (PCSK6), was found to be an up-regulated DEG in liver metastasis in this study, and was closely related to tumor progression [
14]. Current evidence has shown that PCSK6 participates in proliferation, migration and apoptosis of tumor cells in different human malignancies [
15,
16], but its role in the liver metastasis of pancreatic cancer remains unclear.
In this study, we used CRISPR/Cas9 technology to knock out PCSK6 in pancreatic cancer in vitro and in vivo, aiming to explore the functions of PCSK6 and the corresponding mechanism involved in liver metastasis, to update the knowledge regarding the role of PCSK6-targeted therapeutics in the treatments of the liver metastasis in pancreatic cancer.
2. Materials and Methods
2.1. Bioinformatics Analysis of Liver-Metastasis-Related Genes
Transcriptome data were collected from Gene Expression Omnibus (GEO,
www.ncbi.nlm.nih.gov/geo/, accessed on 15 August 2021). The searching criteria were defined as pancreatic cancer OR pancreatic tumor OR pancreatic ductal adenocarcinoma (keyword), Homo sapiens (organisms), and expression profiling by array (study type).
DEGs were generated between liver metastasis and primary tumor in each data series, respectively using the online tool GEO2R (
www.ncbi.nlm.nih.gov/geo/geo2r, accessed on 15 August 2021), with parameters including adjusted
p-values and the Benjamini Hochberg (false discovery rate). Log transformation was applied to the data. DEGs with Log2 FC ≥ 1 or ≤−1 and adjusted
p-values < 0.05 were considered statistically significant. Overlapping DEGs from different data series were obtained using Venn diagrams (bioinformatics.psb.ugent.be/webtools/Venn, accessed on 29 August 2021).
Overlapping DEGs were validated in The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) database using an online tool GEPIA2 (
gepia2.cancer-pku.cn/#index, accessed on 10 September 2021). Adjusted
p-values < 0.01 and Log2 FC ≥ 1 or ≤−1 were defined as cutoff values.
To exclude the influence of mesenchymal components, DEGs obtained in the previous step were verified in cell lines of pancreatic cancer using Cancer Cell line Encyclopedia (CCLE) through an online tool DepMap Portal (depmap.org/portal/, accessed on 12 September 2021). The mRNA expression profiles of 13 cell lines from CCLE (ASPC1, BXPC3, CAPAN1, CAPAN2, CFPAC1, HPAC, HPAFII, HS766T, MIAPACA2, PANC1, SU8686, SUIT2 and SW1990) were applied. The mRNA expression of each gene was presented as transcripts per kilobase million (TPM). Log2(TPM + 1) < 2 was defined as negative expression. The DEGs in the previous step with Log2(TPM + 1) < 2 in all 13 cell lines were excluded. The verified DEGs were represented by heat plots using an online tool Heatmapper (
www.heatmapper.ca/, accessed on 13 September 2021). Lasso Cox regression (R package glmnet, R version 4.2.1) was performed, followed by univariate and multivariate Cox regression (R package survival, R version 4.2.1), to find prognostic genes among the candidate genes generated above. For overall survival analysis, cohorts were grouped into high and low expressions of DEGs using the median expressions of candidate genes in the TCGA database. Kaplan–Meier survival plots were drawn, and Log-rank tests were employed (R package survminer, R version 4.2.1). The cutoff of
p-values was defined as <0.05.
2.2. Cell Culture
The human pancreatic cancer cell line SUIT-2 was obtained from the Japanese Collection of Research Bioresources Cell Bank (JCRB1094, RRID: CVCL_3172) and cultured in minimum essential medium (Gibco, 31095-029) supplemented with 10% fetal bovine serum (Gibco, A3160801). Another human pancreatic cell line MIA PaCa-2 was obtained from ATCC (CRM-CRL-1420, RRID: CVCL_0428) and cultured in Dulbecco’s modified Eagle’s medium (Gibco, 31966-021) supplemented with 10% fetal bovine serum and 2.5% horse serum (Gibco, 16050-130). Both cell lines were incubated in a humidified atmosphere with 5% CO2 at 37 °C.
2.3. CRISPR/Cas9 Gene Editing
In this study, PCSK6 was knocked out in the cell lines mentioned above using the CRISPR/Cas9 gene editing system as described previously [
17]. Briefly, forward and reverse strands of oligonucleotides for the designed sgRNAs targeting PCSK6 were synthesized by Eurofins Genomics (
Table S1). The sgRNA oligos were cloned into the pSpCas9(BB)-2APuro (PX459) V2.0 vector, which was kindly provided by Feng Zhang (Addgene #62988,
http://n2t.net/addgene:62988, RRID: Addgene_62988, accessed on 14 September 2021). The vectors with inserted sgRNA oligos were transformed into Endura DUOs Electrocompetent Cells (BioSEARCH Technologies, 60242-2, Middleton, WI, USA). The CRISPR plasmids were isolated and validated via DNA sequencing. Cell lines were transfected with the validated CRISPR plasmids using a Lipofectamine 3000 Transfection Kit (Invitrogen, L3000-008, Waltham, MA, USA) according to the manufacturer’s instructions. Puromycin (Invivogen, ant-pr-1, 10 μg/mL, Waltham, MA, USA) was used to eliminate un-transfected cells 48 h after transfection and the generated pool cells were seeded in a 96-well plate to establish single clones. The finally generated cells included SUIT-WT (SUIT-2 wild type), SUIT-NC (negative control transfected with nonsense sgRNA), SUIT-SC1 (single clone with knockout of PCSK6 by sgRNA1), SUIT-SC2 (single clone with knockout of PCSK6 by sgRNA2), MIA-WT (MIA PaCa-2 wild type), MIA-NC (negative control transfected with nonsense sgRNA), MIA-SC1 (single clone with knockout of PCSK6 by sgRNA1) and MIA-SC2 (single clone with knockout of PCSK6 by sgRNA2). RT-PCR and genomic DNA sequencing were both performed to validate the single clones.
2.4. RT-PCR
Quantitative analysis of PCSK6 mRNA was performed by means of RT-PCR. Briefly, total RNAs were extracted from cells using a NucleoSpin RNA Plus Kit (MACHEREY-NAGEL, 740984.250, Dueren, Germany). The concentration and purity of the obtained RNAs were determined using the NanoDrop 2000 Spectrophotometer. RNAs were converted to cDNAs using a High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, 4368814, Baltics, Vilnius, Lithuania). Real time PCR was performed via the Bio-Rad CFX Connect Real-Time System, using a PowerSYBR Green PCR Master Mix kit (Applied Biosystems, 4367659, Woolston Warrington, UK) following the manufacturer’s instructions. 2
−∆∆ct values were calculated for the relative quantification of mRNAs. Forward and reverse primers for real-time PCR were synthesized by Eurofins Genomics (Ebersberg, Germany,
Table S1).
2.5. Construction of Plasmid for Genomic DNA Sequencing
Genomic DNA was isolated using a NucleoSpin Kit (MACHEREY-NAGEL, 740952.250, Dueren, Germany). Amplifications of sgRNAs targeting regions in the obtained genomic DNA were performed with Q5
® Hot Start High-Fidelity 2X Master Mix (New England Biolabs, M0494L, Frankfurt, Germany) according to the manufacturer’s instructions and primers for sgRNAs targeting regions were synthesized by Eurofins Genomics (
Table S1). Products of DNA amplifications were purified with Wizard
® SV Gel and a PCR Clean-Up System (Promega, A9281, Madison, WI, USA), followed by ligation of the inserted DNA and vector using NEB PCR Cloning Kits (New England Biolabs, E1203S, Frankfurt, Germany). The generated vectors in above step were transformed into Competent E. coli (New England Biolabs, C3019, Frankfurt, Germany) and isolated using GeneJET Plasmid miniprep Kit (Thermo Scientific, K0503, Baltics UAB, Vilnius, Lithuania). DNA Sequencing of the obtained plasmids was performed by Eurofins Genomics.
2.6. Proliferation Assay
Cells were seeded on a 96-well plate (1000 cells/well) and stained with DAPI (Invitrogen, H3570, Eugene, OR, USA) in quintuplicate at 6 h, 24 h, 48 h, 72 h, 96 h and 120 h after seeding, respectively. Images were obtained using an Evos FL Auto 2 imaging system (Invitrogen by Thermo Fisher Scientific, Bothell, WA, USA). Counts of cells were determined using Thermo Scientific HCS Studio Cellomics Scan software (Version 6.60-1.00x, build 8153) and ratios of cell proliferation were calculated using the following formula: mean of cell counts (N hours)/mean of cell counts (6 h).
2.7. Colony Formation Assay
Cells were seeded on 6-well plates (100 cells/well) in triplicate and incubated for 8 days. Then, the medium was removed, and cells were fixed with 4% formaldehyde solution for 20 min. Cells were stained with 0.1% crystal violet for 15 min. Colonies consisting of more than 50 cells were counted visually and mean counts of colonies were calculated.
2.8. Cell Cycle Analysis
A total of 106 cells in the exponential phase were collected and fixed in 70% cooled ethanol for 30 min. Cells were incubated with Propidium Iodide Solution (Biolegend, 421301, München, Germany) and RNse A (MACHEREY-NAGEL, R2045S, Dueren, Germany) for 15 min at room temperature, followed by adding cooled PBS to terminate the reaction. Stained cell samples were examined by means of flow cytometry (BD Biosciences LSRII, San Jose, CA, USA) and data were processed using FlowJo v10.8 Software (BD Life Sciences, Ashland, OR, USA).
2.9. Migration Assay
After starving for 24 h, 5 × 104 cells were seeded in superior chamber of 8.0 μm Corning FluoroBlok Cell Culture Insert (Corning, 351152, Durham, NC, USA) adding medium without fetal bovine serum, whereas medium with 10% fetal bovine serum was added to the inferior chamber. Cells on the inferior surface of the superior chamber were stained with DAPI (Invitrogen, H3570) 20 h after seeding and images were obtained using an Evos FL Auto 2 imaging system (Invitrogen by Thermo Fisher Scientific, Bothell, WA, USA). Counts of migrated cells in each field (30 fields in total) were determined using Thermo Scientific HCS Studio Cellomics Scan software (Version 6.60-1.00x, Build 8153) and mean migrated cells were calculated according to the following formula: cell counts per field = total counts/30 fields.
2.10. Wound-Healing Assay
In the wound-healing assay, 106 cells were seeded in 6-well plates and starved in medium without fetal bovine serum before drawing a straight line in each well to establish the wound. Images were obtained at 6 different regions of the wound in each well, at baseline and 20 h after establishing the wound, using an Evos FL Auto 2 imaging system (Invitrogen by Thermo Fisher Scientific, Bothell, WA, USA). The areas of the wounds in the 6 above mentioned regions of each well were determined using imageJ (1.53 k, National Institute of Health, Bethesda, MD, USA). The ratios of wound change (healing ratio) were calculated using the following formula: (area of region baseline–area of region 20 h)/area of regions baseline.
2.11. Immunofluorescence
In this step, 105 Cells were seeded on cell culture slides (SPL Life Sciences, 30104, Gyeonggi-do, Korea) and fixed with 4% formaldehyde solution, followed by permeabilization with 0.1% Triton X-100. Non-specific antigen was blocked by 10% Goat Serum. All the steps above were performed at room temperature. Primary E-cadherin antibodies (CST, 3195, RRID: AB_2291471) were incubated overnight at 4 °C and Isotype Normal Rabbit IgG (RD, AB-105-C) was utilized as a control. Secondary Goat anti Rabbit IgG antibody with Alexa Fluor 488 (Invitrogen, A11034, RRID: AB_2576217, Eugene, Oregon, USA) was incubated for 1 h at room temperature. F-actin was stained with Phalloidin FITC Reagent (Abcam, ab235137, Berlin, Germany) for 1 h. The nucleus was stained with DAPI (Thermo Scientific, V13069211, Eugene, OR, USA) for 10 min. Images were acquired and processed using an Evos FL Auto 2 imaging system (Invitrogen by Thermo Fisher Scientific, Bothell, WA, USA). Quantitative analyses of E-cadherin staining at cell–cell junctions were completed using imageJ.
2.12. Gene Set Enrichment Analysis (GSEA) and Gene Ontology (GO) Analysis
The expression profiles of pancreatic cancer cell lines were downloaded from CCLE through the DepMap Portal. The cell lines were ranked in decreasing order according to the median expression of PCSK6. The three cell lines ranked at the top (SUIT-2, SU8686 and ASPC1) and three cell lines at the bottom (SW1990, HS766T and CFPAC1) were assigned to PCSK6-high group and PCSK6-low group, respectively. In total, 16,383 genes in PCSK6-high group and PCSK6-low group were included in GSEA. Hallmark gene sets from the Molecular Signatures Database (MSigDB v7.4
www.gsea-msigdb.org/gsea/msigdb, accessed on 15 September 2021) were used to perform GSEA (GSEA v4.2.3). NES was generated and normalized
p-values of the gene set <0.05 were considered statistically significant.
The R package clusterProfilter was adopted to perform GO analysis. Candidate pathways generated by GSEA were sorted into biological process (BP), molecular function (MF) and cell component (CC) categories. A network plot was established based on the target BP categories to provide an overview of the relationships among these categories and to identify core genes.
Potential proteolysis sites or motifs within target proteins were analyzed via the ELM database (the Eukaryotic Linear Motif resource for Functional Sites in Proteins), which provide a repository of annotated motif data and is a valuable tool for motif prediction (
http://elm.eu.org/, accessed on 25 September 2021).
2.13. Western Blot
The total protein was lysed using the PierceTM RIPA Buffer (Thermo Scientific, 89900, Rockford, IL, USA) mixed with HaltTM Protease and Phosphatase Inhibitor Cocktail (Thermo Scientific, 78442, Rockford, IL, USA). The quantitation of proteins was carried out using a PierceTM BCA Protein Assay Kit (Thermo scientific, 23227, Rockford, IL, USA). Protein samples were prepared by mixing an equal amount of protein with LDS Sample Buffer (NovexTM by life technologies, B0008, Carlsbad, CA, USA) and Sample Reducing Agent (NovexTM by life technologies, B0009) before denaturation. Gel electrophoresis was performed using a Mini Gel Tank chamber system (Invitrogen) with BoltTM 4–12% Bis-Tris gel (Invitrogen, NW04122BOX, Carlsbad, CA, USA), followed by transferring the protein to the Nitrocellulose Blotting membrane (AmershamTM, 10600003, Germany). 5% BSA or nonfat milk was used to block nonspecific antigens, and primary antibodies were incubated overnight at 4 °C, including phospho-p44/42 ERK1/2 (CST, 4370, RRID: AB_2315112), p44/42 ERK (CST, 4695, RRID: AB_390779) and GAPDH (CST, 5174, RRID: AB_10622025). HRP-linked Anti Rabbit IgG (CST, 7074, RRID: AB_2099233) was used as secondary antibody. Examination of the signal was performed using Amersham Imager 600 (Pittsburgh, PA, USA) with SignalFire™ ECL Reagent (CST, 6883S, Frankfurt, Germany).
2.14. Tumor Xenograft Model of Intrasplenic Injection
Intrasplenic injection was performed to establish liver metastasis. Briefly, 6-week-old NSG mice (NOD.Cg-Prkdcscid/l2rgtm1Wjl/Szj, The Jackson Laboratory, RRID: IMSR_JAX: 005557) were prepared via inhalation anesthesia with 2% Isoflurane, and intraperitoneal analgesia with 1% ketamine chloride and 0.1% xylazine chloride. A left flank incision was made to expose and move the spleen out of the abdominal cavity. Then, 4 × 104 tumor cells resuspended in 20 μL PBS were injected into the lateral end of the spleen with a 30 G insulin syringe, followed by pressing the injection site with a cotton swab for 1 min to prevent leakage and facilitate hemostasis. Splenic vessels and gastrosplenic ligament were ligated 5 min after injection, and the spleen was removed to avoid tumor formation in the spleen. The abdominal wall was closed with a 5-0 suture (ETHICON, 1666H, San Lorenzo, Puerto Rico, USA). Three weeks after establishing the model, mice were sacrificed, and an exploration was carried out to rule out the planting lesions of tumors within the abdominal cavity. The whole liver was resected for pathological analysis, including macroscopic examinations and microscopic investigations of histological sections with hematoxylin and eosin (HE) staining. Four NSG mice were applied in each group (SUIT-SC and SUIT-NC). Three cross-sections of each liver sample were subjected to HE staining, and each section was scanned and transformed into a digital image. A total of 12 digital images for each group were utilized to calculate the counts of liver metastatic niches. The sizes of metastatic niches were determined using NDP.view 2.7.39 software. All the animal experiments were approved by the governments of Mittelfranken or Unterfranken, Germany. All the procedures mentioned above, animal housing, and animal husbandry were performed in accordance with the ethical regulations of the governments of Mittelfranken or Unterfranken, Germany.
2.15. Statistical Analysis
The statistical methods used in bioinformatics analysis have been mentioned above. Moreover, the variables in this study were continuous variables and are presented as means ± standard deviation (SD). The two-tailed unpaired Student’s t-test was utilized for comparisons between two groups. One-way ANOVA was utilized for comparisons among three or more groups, followed by multiple comparisons between groups using the Scheffe method. p-values < 0.05 were considered statistically significant. Different levels of p-values were defined as * <0.05, ** <0.01 and *** <0.001. All the statistical analyses were performed using the SPSS software 26.0 (SPSS Inc, Chicago, IL, USA).
4. Discussion
The TME is comprised of tumor cells and the mesenchyme, including cellular components and non-cellular components. The interactions between tumor cells and the mesenchyme were found to induce variations in molecular profiles of tumor cells and remodel the TME, finally promoting tumor progression [
18]. The biological behaviors of tumor cells are governed by alterations of driver genes, however, epigenetic, post-transcriptional, and post-translational modifications coupled with the TME have emerged as equally critical determinants during the course of tumor development [
19]. Therefore, the TME of liver metastatic niches may provide us with new clues to find genes related to the liver metastasis of pancreatic cancer. We carried out unique bioinformatics analysis of the mRNA expression profiles in both liver metastasis and primary tumors, and identified 29 candidate genes. One of the obtained genes, PCSK6, was significantly prognostic and this was chosen for further study.
PCSK6 locates at 15q26.3 and is the coding gene of PACE4, a member of the proprotein convertase family and a kind of calcium-dependent serine proteinase which participates in post-translational modifications [
20]. The proprotein convertase family consists of nine currently known members, including PC1, PC2, furin, PC4, PC5, PACE4, PC7, SKI-1, and PCSK9. With the exception of SKI-1 and PCSK9, the remaining seven members of the proprotein convertase family have been found to activate various protein precursors through proteolysis at specific single or paired basic amino acids within the motif (R/K)Xn(R/K)↓, however, there are variations in cleavage sites for each proprotein convertase [
14]. PACE4 is synthesized as an inactive precursor and then it is converted into a mature enzyme via both autoactivation in the endoplasmic reticulum and activation by other proteinases, which are prerequisites for processing its substrates [
21]. Up-regulation of PCSK6 in the liver metastasis of pancreatic cancer was identified in this study. The detailed functions of PCSK6 in tumor cells and its impact on the liver metastasis of pancreatic cancer remain to be determined. CRISPR/Cas9, a cutting-edge gene editing technology, was utilized to inactivate PCSK6 in this study. Through the accurate binding of sgRNA guided Cas9 nuclease and genome DNA, exons of PCSK6 were edited at designated sites, which were proven via genome DNA sequencing. Variations in base sequences among exons resulted in the decreased efficiency of mRNA transcription and significant inhibitions of PCSK6 mRNA expression, which were confirmed via RT-PCR.
The results that single clones of both SUIT-2 and MIA PaCa-2 showed lower proliferation rates and fewer colony formation counts compared to controls, indicated that the inactivation of PCSK6 could inhibit the growth of tumor cells in pancreatic cancer. G0/G1 cell cycle arrest with decreased ratios and decreased trends in the S/G2 phase, were demonstrated in single clones of SUIT-2 and MIA PaCa-2, respectively, at least in part validating the above results. Similar findings were also obtained in skin cell tumors [
15,
22] and prostate cancer [
23,
24]. PACE4 acts through calcium-dependent serine proteases and shows a preference for the cleavage sites of (R/K)Xn(R/K), accordingly its substrates include various growth factors, metalloproteinases, etc. [
14]. We postulated that the inactivation of PCSK6 attenuated the interactions between PACE4 and its substrates, leading to immature growth factors and decreased downstream signaling, mediating the growth of tumor cells.
In addition, the inactivation of PCSK6 in single clones of both cell lines could also suppress cell migration, according to the results of the migration assay and wound-healing assay in this study. Most previous researchers focused on the effect of PCSK6 on cell growth or apoptosis [
23,
24,
25], whereas one study noted that the elevated expression of PCSK6 enhanced the migratory capacity of skin tumor cells [
22], which was consistent with our findings in relation to pancreatic cancer. Molecular reprogramming in tumor cells may induce modifications of cell–cell junctions or the cell skeleton, which play a critical role in cell movement, deformation, etc. [
26]. E-cadherin is a component of the stable adherens junctions connecting the actin skeletons of adjacent cells [
27]. Proprotein convertases including PC7, furin, and PACE4 have been proven to regulate cell–cell junctions by processing E-cadherin during blastocyst formation [
28]. Moreover, we found specific motifs or cleavage sites for proprotein convertases within the amino acids sequence of E-cadherin (
Table S5). Therefore, we examined the expression of E-cadherin, and it was worth noting that, instead of decreasing the expression of E-cadherin, the inactivation of PCSK6 enhanced the expression of E-cadherin at sites of cell–cell junctions in single clones of SUIT-2, although not in MIA PaCa-2, which showed an extremely low expression of E-cadherin [
29,
30]. Previous studies have shown that, when one of the proprotein convertases is inactivated, other members can process immature E-cadherin efficiently [
31]. Accordingly, a possible explanation for our observations was that the maturation of E-cadherin was not dependent on PACE4, and other proprotein convertases probably played a compensatory role in processing premature E-cadherin. Alternatively, other unknown mechanisms might be responsible for the increased expression of E-cadherin after the inactivation of PCSK6 in SUIT-2 cells. Repression of the Raf-MEK-ERK signaling pathway would lead to restoration of E-cadherin [
32,
33] and this might be one of the reasons for the up-regulation of E-cadherin observed in this study. Furthermore, F-actin, acting in coordination with adherens junctions to direct cell morphology and cell movement, was found to be up-regulated in single clones of both cell lines. Specifically, transformation of cell morphology from spindle-shaped to cubic or round, as well as regaining of a glandule-like structure and apicobasal polarity, were features highlighted after the inactivation of PCSK6. Variation in cell–cell junctions and the cell skeleton enable tumor cells to adapt to different phases of tumor progression and to various mesenchymal components in the TME [
27]. In summary, the inactivation of PCSK6 led to the decreased migratory capacity of tumor cells in pancreatic cancer, and reprogramming of cell–cell junctions or the cell skeleton might be one of the causes for this effect.
The final step in the metastatic process is the generation of colonies in the liver, which requires tumor cells to overcome serial obstacles to outgrow into a macroscopic metastasis, including getting access to the liver through blood vessels, gaining a nutritional supply or growth stimuli from surrounding tissue, remodeling or establishing the TME, and escaping from immune elimination [
19]. In this study, tumor xenograft model of intrasplenic injection was used to simulate the colonization of tumor cells in the liver. Prominent macroscopic metastatic lesions were found in the NC group of SUIT-2 but not in the SC group. Microscopic metastatic niches were found in both groups of SUIT-2, whereas fewer counts and smaller-sized metastatic niches were observed in the SC group compared to the NC group. Considering that PACE4 can process and activate growth stimuli and adhesion molecules in the extracellular matrix (ECM) in the manner of a secretory protein [
14], the results showing that the inactivation of PCSK6 made tumor cells less competent to establish metastatic niches in the liver and finally retarded tumor growth in vivo, might in part be attributed to the presence of fewer extravasating cells from the blood wall, fewer anchored cells in liver tissue, or less active growth stimuli surrounding tumor niches. Small molecule inhibitors targeting PACE4 have been investigated in relation to arthritis disease [
34], human skin tumors [
15] and human prostate cancer [
35]. Accordingly, our preclinical animal model will provide evidence and options for future evaluations of small-molecule inhibitors of PACE4 or their combined use with chemotherapy in the treatment of liver metastasis in pancreatic cancer.
KRAS, CDKN2A, TP53 and SMAD4 are the most frequently observed driver genes involved with genetic mutations in pancreatic cancer, and KRAS mutations occur in 90% of the patients [
36]. Through GSEA, we found that the high expression of PCSK6 was accompanied by the up-regulation of Kras signaling pathway compared to the group with a relatively low expression of PCSK6, in the cell lines of pancreatic cancer. The Raf-MEK1/2-ERK1/2 axis is identified as one of the downstream signaling pathways of Kras. Previously studies have linked furin, another member of the proprotein convertase family, to the activation of the Raf-MEK1/2-ERK1/2 pathway via endoproteolytic cleavage in KRAS or BRAF mutant cells [
37]. As a ubiquitously distributed proprotein convertase in human tissues, PACE localizes subcellularly to the trans-Golgi network or the cell surface, and its substrates include hormones, tyrosine phosphatases, growth factors, metalloproteinases, and adhesion molecules in the extracellular matrix (ECM) or on the cell surface [
38]. One possible mechanism is that PACE4-processing growth factors or other substrates on the cell surface or in the ECM might exert an effect on the Raf-MEK1/2-ERK1/2 axis. We evaluated the impact of PCSK6 inactivation on the Raf-MEK1/2-ERK1/2 axis in this study. Strikingly, compared to controls, single clones of both cell lines showed extremely low activities of phosphorylated ERK1/2 in this study, which is a critical modulator regulating the transcriptions of multiple genes correlated with the cell cycle, apoptosis, cell migration, cell invasion, etc. [
39]. Since SUIT-2 and MIA PaCa-2 in this study are both cell lines with KRAS mutations and they activate the downstream Raf-MEK1/2-ERK1/2 axis, regardless of whether the upstream components such as RTK are muted due to immature ligands, there may be another reasonable explanation for the notable suppression of the Raf-MEK1/2-ERK1/2 axis after PCSK6 inactivation. Therefore, another potential mechanism is that Kras, Raf, MEK1/2, and ERK1/2 are potential substrates of proprotein convertases according to our analysis based on the predictions of the ELM database, which provide theoretical evidence for interactions between PACE4 and Kras or the Raf-MEK1/2-ERK1/2 axis, however, the current literature does not provide sufficient evidence for the interaction of PACE4 with Kras or the Raf-MEK1/2-ERK1/2 axis intracellularly [
39,
40]. More investigations will be necessary to confirm this.
Previous researchers found that PCSK6 regulated cell apoptosis via endoplasmic reticulum stress and the mitochondria signaling pathway [
25]. We have advanced the knowledge by demonstrating that the inactivation of PCSK6 can significantly inhibit the Raf-MEK1/2-ERK1/2 axis, having a repressive effect on cell growth, cell migration and the generation of metastatic niches in the liver in pancreatic cancer. The epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor has demonstrated limited survival benefits for patients with pancreatic cancer due to mutated KRAS [
41], however, the inactivation of PCSK6 may be a novel target of treatment, bypassing mutated KRAS and retarding tumor progression in pancreatic cancer, providing that a direct interaction between PACE4 and the Raf-MEK1/2-ERK1/2 axis can be proven in the future.