IMUP is upregulated in PDAC and associated with poor prognosis
We detected the transcriptomics of tumor and adjacent normal tissues of ten paired primary PDAC patients by microarray scanning using GeneChip. DEGs between tumor and adjacent tissues were analyzed (Fig. S1a, b). In addition, we obtained DEGs and survival-related genes (SRGs) from the TCGA database (http://www.cbioportal.org/index.do), the intersecting DEGs between GeneChip (fold change > 4), TCGA, and SRGs identified IMUP as the most notable gene (Fig. 1a).
The data from the MEXPRESS online tool indicated that tumor stage and OS of PDAC patients were associated with IMUP expression (Fig. 1b) (Koch et al. 2019). The general characteristics of patients, including age, sex, years smoked, and alcohol history, were not significantly different in terms of IMUP expression (P > 0.05; Fig. 1b). Consistent with these results, bioinformatics analysis of data from TCGA and GTEx databases by GEPIA demonstrated that the mRNA level of IMUP was upregulated in tumors (Fig. 1c) (Tang et al. 2019).
Protein levels were investigated in PDAC tumor tissues and adjacent normal tissues using WB and IHC. Results showed that the IMUP abundance in tumor tissues was higher than that in adjacent tissues (Fig. 1d, Fig. S2a). Moreover, data from GEPIA demonstrated that patients with high IMUP abundance had a poor OS, DFS, and an advanced clinical stage (Fig. 1e-g) (Tang et al. 2019), which was consistent with the results of Kaplan–Meier analysis of our clinical tissues (P < 0.001) (Fig. 1h). Univariate and multivariate Cox regression analysis confirmed that high IMUP abundance was a significant predictor of OS (hazard ratio = 3.695, P < 0.001; Table S4). Finally, we verified that IMUP levels in PDAC cell lines, including BxPC-3, SW1990, and PANC-1, was significantly elevated in comparison with that in HPDE cells (Fig. 1i). Collectively, these data suggest that IMUP is upregulated in PDAC and negatively affects patient prognosis. Hence, IMUP may play a critical role in PDAC progression.
IMUP Knockdown inhibits cell growth in vitro and in vivo
Three shRNAs were used to interfere with IMUP expression in BxPC-3 and SW1990 cells. The mRNA and protein levels of IMUP were markedly depleted by shRNA-1 and shRNA-2 (Fig. 2a) and IMUP knockdown significantly inhibited the proliferative and colony formation capacity of BxPC-3 and SW1990 cells (Fig. 2b-d). Flow cytometry results revealed that the cells were blocked at the S phase by IMUP depletion (Fig. 2e, f). Moreover, IMUP depletion inhibited tumor xenograft formation of BxPC-3 cells in BALB/c nude mice (Fig. 2g). The volume and weight of IMUP-sh1 and IMUP-sh2 xenograft tumors also decreased significantly in comparison with normal control group (shNC) xenograft tumors (Fig. 2h, i). IHC was used to assess IMUP downregulation in the xenograft tumors (Fig. S2b). Results indicate that IMUP knockdown significantly suppressed tumor growth in vitro and in vivo, consequently confirming the potent role of IMUP in PDAC tumorigenicity.
IMUP regulates cyclin A2/cyclin E1/CDK2 proteins via FHL1
To investigate the mechanisms of IMUP in cell-cycle arrest, we performed RNA-seq for BxPC-3 cells with or without IMUP knockdown. A total of 621 genes were differentially upregulated, while 153 genes were downregulated following IMUP depletion (Fig. 3a). In addition, we analyzed IMUP-correlated genes from the TCGA database using LinkedOmics (Vasikar et al. 2018). Gene Ontology (GO) analysis results showed that IMUP-associated genes were enriched in regulation of cyclin-dependent protein kinase activity (Fig. 3b). The positively- and negatively-correlated genes (PCGs and NCGs, respectively) associated with IMUP were analyzed using the results from GeneChip, RNA-seq (fold-change > 2) and the TCGA database. The intersection of PCGs from them provided wingless-type MMTV integration site family member 10A (WNT10A) (Fig. S3a), whereas that of NCGs from them provided G protein subunit gamma 2 and FHL1 (Fig. 3c), which might be regulated by IMUP. Among these DEGs, FHL1 is reportedly associated with negative regulation of cell growth. The data from GeneChip and TCGA verified that the FHL1 expression was significantly negatively correlated with IMUP (Fig. 3d, e). Moreover, high expression of FHL1 correlated with longer OS (P = 0.023) in patients with PDAC (Fig. 3f) (Tang et al. 2019). Thus, we hypothesized that FHL1 may mediate the inhibition of tumor growth induced by IMUP depletion.
The WB results showed that knockdown of IMUP decreased the abundance of cyclin A2, cyclin E1, and CDK2 in BxPC-3 and SW1990 cells, which are the key protein kinases of the S phase, and increased the mRNA and protein levels of FHL1 (Fig. 4a; Fig. S3b). The silencing of FHL1 via siRNA restored the protein inhibition induced by IMUP-shRNA (Fig. 4b). In addition, the results of PDAC tumor tissue IHC revealed that high expression of IMUP was associated with a low level of FHL1, while Pearson’s correlation analysis suggested that their IHC scores were negatively correlated (Fig. 4c). These results confirm that IMUP regulates cell-cycle protein activity of PDAC cells by FHL1.
Knockdown of FHL1 rescues the phenotype inhibited by IMUP depletion
To determine whether FHL1 is a potential target of IMUP that promotes tumor growth. FHL1-shRNA was co-transfected with IMUP-shRNA into BxPC-3 and SW1990 cells. Results show that FHL1-shRNA restored cell proliferation, colony formation, S phase arrest, and tumor xenograft growth in mice, which were inhibited by IMUP knockdown (Fig. 5a-c, g). The data showed significant statistical differences between IMUP-shRNA cells, control shRNA cells, and cells co-transfected with IMUP-shRNA and FHL1-shRNA (Fig. 5d-f, h). Thus, we confirmed that IMUP promotes tumor progression of PDAC by regulating FHL1 expression.
IMUP inhibits the transcription of FHL1 by NPM1-induced promoter methylation
To determine the mechanism by which IMUP regulates FHL1 expression, GFP-tagged IMUP overexpression and endogenous IMUP were used for IP (Fig. S4a). The intersection of GFP-IMUP binding proteins and endogenous IMUP binding proteins strongly supports that NPM1 is a probable target of IMUP according to LC/MS analysis (Table S5). Co-IP assay further assessed the direct interaction between IMUP and NPM1 (Fig. 6a). Confocal imaging revealed a partial colocalization of IMUP and NPM1 in PDAC cell nucleus (Fig. 6c).
In addition, we found that knockdown of IMUP reduced the protein level of NPM1 (Fig. 6b). However, the protein level of IMUP was not affected by NPM1 depletion induced by NPM1-siRNAs (Fig. S5a). The mRNA of NPM1 was not affected by IMUP knockdown (Fig. S5b). We believed that IMUP sustained NPM1 stability. CHX was used to block protein synthesis in BxPC-3 and SW1990 cells transfected with control and IMUP siRNAs. The results showed that NPM1 protein degraded slowly after CHX treatment within 48 h in control group, whereas the IMUP knockdown obviously facilitated NPM1 protein degradation (Figs. 6d, e, and S5c). These data indicate that IMUP stabilizes nuclear protein NPM1.
Many studies have reported that FHL1 undergoes epigenetic regulation by promoter methylation in various cancers (Koike et al. 2013), and NPM1 plays an important role in epigenetic regulation (Karimi et al. 2021). Thus, to determine whether IMUP regulated methylation of the FHL1 promoter through NPM1, we assessed the promoter methylation status of two promoter exon fragments. Pyrosequencing results confirmed that the methylation of 16 CpG sites decreased significantly in IMUP-depleted cells (Fig. 6f, g). However, NPM1 overexpression restored promoter methylation. To investigate whether FHL1 promoter methylation directly affected its transcription, we identified two putative specific protein 1 (SP1)-binding sites in FHL1 methylated promoter region: site1 (AGGGCAGGGT) and site2 (GGGGAGGGGT) (Fig. 6h). ChIP qPCR results revealed that SP1 was enriched at two sites of FHL1 promoter. Moreover, IMUP knockdown increased SP1 enrichment in FHL1 promoter. NPM1 overexpression inhibited SP1 binding to the FHL1 promoter (Fig. 6i). Meanwhile, qPCR demonstrated that FHL1 expression was inhibited by the NPM1 vector following IMUP knockdown (Fig. S6). However, FHL1 expression was restored to previous levels when the cells infected with shIMUP and GFP-NPM1 vectors were treated with DNA methyltransferase inhibitor 5′-aza (Fig. S6).
FHL1 regulates cell-cycle protein kinases by interacting with CHK1/CDC25A/14-3-3ξ
IP was performed to explore the mechanisms by which FHL1 regulates cell-cycle arrest. After flag-tagged FHL1 vectors were transfected into BxPC-3 cells, FHL1-binding proteins were analyzed by LC/MS (Fig. S4b, Table S6). Multiple subtypes of 14-3-3 proteins have been identified as candidate binding proteins, which are related to cell-cycle regulation (Gardino et al. 2011). Co-IP was used to validate the interaction between FHL1 and 14-3-3 proteins. We also verified the binding of FHL1 with CDC25A and CHK1 proteins, which reportedly interact with FHL1 in other cancer cells (Xu et al. 2017). The results showed that Flag-FHL1 interacted with 14-3-3ξ, CDC25A, CDC25C, and CHK1 (Fig. 7a). The interactions between FHL1 and 14-3-3ξ, CDC25A, CDC25C, and CHK1 were determined to direct as endogenous 14-3-3ξ, CDC25A, CDC25C, and CHK1 bound to FHL1 (Fig. 7b).
CDC25 subtypes, including CDC25A, CDC25B, and CDC25C, play critical roles in the regulation of cell-cycle checkpoints (Donzelli et al. 2003). Among them, phosphorylation of CDC25A by CHK1 generally contributes to cell-cycle arrest in the S phase (Zhao et al. 2002). FHL1 may regulate CDC25A activity by forming a protein complex with CHK1. Thus, we further investigated whether FHL1 regulated the phosphorylation of CDC25A and found that FHL1 overexpression promoted CDC25A degradation by phosphorylation, thereby inhibiting the protein levels of cyclin A2, cyclin E1, and CDK2 (Fig. 7c). However, FHL1 did not phosphorylate CDC25A directly, as the phosphorylation events were observed to significantly decrease following treatment of cells with a CHK1 inhibitor, which inhibited CHK1 phosphorylation (Fig. 7c). Moreover, overexpression of FHL1 increased the interaction between CDC25A and CHK1 (Fig. 7d). These results indicated that FHL1 increased CHK1-mediated phosphorylation of CDC25A.
It has been reported that 14-3-3ξ binds to CDC25A and prevents CDC25A from activating cell-cycle kinase in HeLa cells (Kahama et al. 2019). Binding between 14-3-3ξ and CDC25A in the cytoplasm inhibited CDC25A activity in cell-cycle progression (Al-Matouq et al. 2017). Further, we determined whether FHL1 affected CDC25A location. WB analysis of nuclear-cytoplasmic isolated protein showed that FHL1 overexpression increased the cytoplasmic distribution of CDC25A in BxPC-3 and SW1990 cells. However, 14-3-3ξ knockdown via siRNA neutralized the modulation of FHL1 on CDC25A localization (Fig. 8a). We also verified that CDC25A was blocked in cytoplasm by IMUP knockdown and that 14-3-3ξ siRNA could restore this effect (Fig. 8b). The mechanism underlying CDC25A distribution regulation by FHL1 was that FHL1 overexpression facilitated the interaction between 14-3-3ξ and CDC25A (Fig. 8c). In addition, immunofluorescence (IF) of xenograft tumors showed that CDC25A was transferred from the nucleus to the cytoplasm after IMUP knockdown (Fig. 8d). This should be attributed to the upregulation of FHL1 induced by IMUP depletion. These results suggest that FHL1 causes localization of CDC25A in the cytoplasm by forming a complex with 14-3-3ξ.