USP10 is associated with poor prognosis in PDAC
In order to investigate the expression of USP10 in PDAC tissues, we initially performed bioinformatics analysis of the TCGA and GTEx PDAC transcriptome databases. The results validated that the expression of USP10 mRNA in PDAC was significantly higher than in normal pancreatic tissues (Fig. 1A). Further Kaplan-Meier survival analysis indicated a relationship between elevated USP10 mRNA expression and worse overall survival (OS) in PDAC patients (Fig. 1B).
To further clarify the clinical significance of USP10 in PDAC, we collected two PDAC tissue microarrays from two independent medical centers (Group 1: 90 cases and Group 2: 46 cases) (Table 1). Subsequently, we performed IHC to identify the protein expression of USP10 in PDAC tissues (Fig. 1C-D). Subsequent survival analysis revealed that patients in both Group 1 and 2 with high USP10 expression had a shorter overall survival (MST: 15.0m vs. 21.8m, P = 0.007; 14.0m vs. 20.0m, P = 0.047) (Fig. 1E-F). Furthermore, there was a significant correlation between elevated USP10 expression and the occurrence of lymph node metastasis and distant metastasis, as indicated in Table 2. Overall, these results indicate that the elevated USP10 expression in PDAC tissues is associated with invasion, metastasis and poor survival.
Table 1
Baseline characteristics of patients with PDAC
Characteristics | Group 1 (n = 90) | Group 2 (n = 46) |
Age | | |
≥ 65 | 22 (32.4%) | 22 (47.8%) |
< 65 | 68 (67.6%) | 24 (52.8%) |
Gender | | |
Male | 53 (58.9%) | 22 (47.8%) |
Female | 37 (41.1%) | 24 (52.8%) |
CA19-9 | | |
≥ 27U/ml | 72 (80.0%) | 40 (87.0%) |
< 27U/ml | 18 (20.0%) | 6 (13.0%) |
Tumor size | | |
≥ 3cm | 74 (82.2%) | 36 (78.3%) |
< 3cm | 16 (17.8%) | 10 (21.7%) |
Stage | | |
I | 24 (26.7%) | 15 (32.6%) |
II | 28 (31.1%) | 21 (45.7%) |
III | 10 (11.1%) | 3 (6.5%) |
IV | 28 (31.1%) | 7 (15.2%) |
Table 2
Correlation between USP10 expression and clinical characteristics in patients with PDAC
Characteristics | Group 1 (n = 88) | Group 2 (n = 46) |
| USP10 High (n = 45) | USP10 Low (n = 45) | p | USP10 High (n = 21) | USP10 Low (n = 25) | p |
Age | | | | | | |
≥ 65 | 10 | 12 | 0.624 | 9 | 13 | 0.536 |
< 65 | 35 | 33 | 12 | 12 |
Gender | | | | | | |
Male | 25 | 28 | 0.643 | 9 | 13 | 0.536 |
Female | 20 | 27 | 12 | 12 |
CA19-9 | | | | | | |
≥ 27U/ml | 32 | 40 | 0.035 | 16 | 24 | 0.079 |
< 27U/ml | 13 | 5 | 5 | 1 |
Tumor size | | | | | | |
≥ 3cm | 37 | 37 | 1.000 | 19 | 17 | 0.084 |
< 3cm | 8 | 8 | 2 | 8 |
Node Metastasis | | | | | | |
Yes | 35 | 23 | 0.008 | 15 | 10 | 0.042 |
No | 10 | 22 | 6 | 15 |
Metastasis | | | | | | |
Yes | 8 | 20 | 0.012 | 5 | 2 | 0.220 |
No | 37 | 25 | 16 | 23 |
USP10 promotes a malignant phenotype in pancreatic cancer cells
To clarify the role of USP10 in the proliferation, invasion, and metastasis of PDAC cells, we stably knocked down USP10 in BxPC-3 and AsPC-1 cell lines and stably overexpressed USP10 in PANC-1 cell lines (Fig. 2A). Firstly, the results of colony formation and subcutaneous tumor formation assays indicated that overexpression of USP10 significantly promoted colony formation and tumorigenesis of PDAC cells, while USP10 knockdown had a significant inhibitory effect (Fig. 2B-E). Cell invasion and migration experiments further revealed that overexpression of USP10 significantly upregulated the invasion and migration capacity of PDAC cells, while knockdown had an inhibitory effect as well (Fig. 2F-G). In summary, these results suggest that USP10 can significantly promote tumorigenesis, invasion, and migration of PDAC cells.
USP10 was involved in the FAK signaling pathway
To clarify the specific mechanisms by which USP10 promotes the malignant phenotype of PDAC cells, we initially performed GSEA analysis using the PAAD-TCGA transcriptome database. The results revealed a significant association between USP10 and the FAK signaling pathway (Fig. 3A). To further explore the above bioinformatics analysis results, we subsequently analyzed FAK pathway-related protein expression in stable cell lines. The results showed that overexpression of USP10 promoted the protein expression and phosphorylation of FAK, as well as integrins α5, β3, and β4. Knockdown of USP10 had an inhibitory effect on the expression of these proteins (Fig. 3B). However, changes in USP10 did not affect FAK mRNA expression (Fig. 3C). Subsequently, we established stable cell lines with knockdown of FAK and knockdown of FAK followed by overexpression of USP10 and conducted phenotype experiments. The results showed that when overexpressing USP10 on the basis of FAK knockdown, it did not exhibit a significant enhancing effect on cell proliferation and invasion compared to the sole overexpression of USP10 (Fig. 3D-E). Therefore, the activation of the FAK signaling pathway plays a crucial role in USP10's promotion of the malignant phenotype of PDAC cells.
USP10 interacts with FAK and RIOK3
USP10 is a deubiquitinating enzyme, and its basic function is to stabilize the expression of substrate proteins through deubiquitination. The results have shown that overexpression of USP10 upregulates FAK protein expression but not mRNA. Following the finding that RIOK3 binds to and stabilizes FAK protein expression, we speculated that USP10, FAK, and RIOK3 interact. Endogenous co-IP tests with USP10 antibodies in PANC-1, BxPC3, and AsPC-1 cell lines showed that USP10 interacts with FAK and RIOK3 (Fig. 4A). Subsequently, co-IP experiments using RIOK3 antibodies also confirmed interactions between RIOK3 and USP10 as well as FAK (Fig. 4B). To further validate whether the interaction between USP10 and FAK depends on RIOK3, we stably knocked down the expression of RIOK3 in PANC-1 and BxPC-3 cells. The results demonstrated that RIOK3 knockdown eliminated USP10-FAK binding (Fig. 4C). This suggests that the interaction between USP10 and FAK may depend on RIOK3. Subsequently, we performed exogenous immunoprecipitation experiments in 293T cells using tagged antibodies, which verified USP10-RIOK3 interaction (Fig. 4D). Further confocal microscopy tests showed that USP10 and RIOK3 co-localize in the cell membrane and cytoplasm (Fig. 4E). Protein interaction prediction analysis of the CDS sequences of USP10 and RIOK3 also suggested the possibility of their binding (Fig. 4F). In summary, our research suggests that USP10 forms a protein complex with RIOK3 and FAK, with RIOK3 acting as a bridge between USP10 and FAK.
USP10 deubiquitinates and stabilizes RIOK3 expression
To assess whether USP10 deubiquitinated RIOK3, we first examined the protein and mRNA expression of RIOK3 in stable cell lines. The results showed that USP10 knockdown decreased RIOK3 protein expression while overexpression increased it (Fig. 5A). However, RIOK3 mRNA levels did not change (Fig. 5B). This suggests that USP10 may deubiquitinate and stabilize RIOK3 protein. Subsequent deubiquitination assay showed that USP10 knockdown increased RIOK3 ubiquitination, while overexpression decreased it. Notably, overexpressing USP10 with enzymatic active site mutation did not affect the RIOK3 ubiquitination level (Fig. 5C). To determine whether USP10 regulates RIOK3 protein stability, we conducted the cycloheximide (CHX) experiment. Results demonstrated that knocking down USP10 decreased RIOK3's protein half-life, while overexpression extended it. Additionally, overexpression of USP10 with mutated enzymatic activity sites did not affect the protein stability of RIOK3 (Fig. 5D). Furthermore, the IHC study of USP10 and RIOK3 in PDAC patients from Groups 1 and 2 indicated a significant positive correlation in protein expression between USP10 and RIOK3 (Fig. 5E-F). Overall, in PDAC, USP10 deubiquitinates and stabilizes RIOK3.
RIOK3 plays a key role in USP10-driven PDAC tumorigenesis
To further elucidate whether USP10 activates the FAK signaling pathway and promotes pancreatic cancer cell invasion, migration, and tumorigenesis via RIOK3, we overexpressed RIOK3 in stable USP10 knockdown BxPC-3 and AsPC-1 cells. Firstly, Western blot results showed that RIOK3 could re-activate the FAK pathway inhibition caused by the USP10 knockdown (Fig. 6A). Subsequent colony formation and Transwell assays confirmed that the overexpression of RIOK3 could reverse USP10 knockdown-induced cell colony formation and invasion inhibition (Fig. 6B-C). The results of subcutaneous tumor formation and lung metastasis models further clarified the critical role of RIOK3 in reversing the inhibition of tumorigenesis and invasion by USP10 knockdown in PDAC cells (Fig. 6D-E). Additionally, IHC staining results of subcutaneous tumor tissues showed that the protein expression of both RIOK3 and FAK decreased after USP10 knockdown, but the expression of FAK was upregulated in the RIOK3 overexpression group compared to the USP10 knockdown group (Fig. 6F). The above results suggest that RIOK3 plays an important role in the activation of the FAK pathway and tumorigenesis mediated by USP10 in PDAC cells.
USP10 and RIOK3 co-expression is significantly associated with poor prognosis in PDAC
In order to further reveal the correlation between USP10-RIOK3 protein expression and the development of pancreatic cancer, we further detected RIOK3 protein expression in the TMAs. The results showed that high expression of RIOK3 in PDAC patients in both Group 1 and 2 was associated with poorer OS (MST: 13.3m vs. 20.0m, P = 0.0002; 14.0m vs. 20.0m, P = 0.019) (Fig. 7A-B). Importantly, patients with high expression of both USP10 and RIOK3 had significantly worse OS (MST: 12.0m vs. 21.0m, P = 0.027; 14.0m vs. 20.0m, P = 0.005) (Fig. 7C-D). Decision curve analysis (DCA) for other factors that may affect the OS of PDAC patients showed that the net benefit of the combined expression of USP10 and RIOK3 was the highest (Fig. 7E). Therefore, we constructed a Nomogram prognostic model based on the Cox multivariate regression results of these factors and conducted Calibration analysis. The model's predicted survival rates were highly consistent with the observed survival rates, indicating that this prognostic model has a high level of accuracy (Fig. 7F-G). In conclusion, high expression of the USP10-RIOK3 protein can serve as an adverse prognostic factor for PDAC.
USP10 is the biomaker of FAK inhibitor efficacy in PDAC
Our study found that USP10 activated the FAK pathway in PDAC, so we speculated that USP10 could predict the efficacy of FAKi. We performed subcutaneous tumors in nude mice with PANC-1 cell line. USP10 overexpression and Control group received FAKi treatment on 15th d (Fig. 8A). The results showed that USP10 high expression mice changed more significantly compared with the control group receiving FAKi treatment (Fig. 8B-C). Moreover, the tumors size of the high USP10 expression group mice continued to shrink, and the tumors size of the Control group begined slowly to increase at 35th dt (Fig. 8B-C). In conclusion,USP 10 high expression may be an effective biomarker for FAKi treatment of PDAC.