MAPK8IP2 is a potential prognostic biomarker and promote tumor progression in prostate cancer
BMC Cancer volume 22, Article number: 1162 (2022)
MAPK8IP2 is one of the JNK-interacting proteins (JIPs) family members, and is involved in the regulation of the JNK and P38 MAPK signaling pathways. MAPK8IP2 has been reported to be closely associated with several cancers. However, the biological function of MAPK8IP2 in prostate cancer (PCa) remains unclear.
MAPK8IP2 expression in PCa and subgroups of PCa was analyzed by public databases. The prognostic role of MAPK8IP2 in prostate cancer was analyzed using the Cox regression method. The potential mechanism by which MAPK8IP2 affects PCa progression was investigated by utilizing public data, including genetic alteration, DNA methylation, m6A methylation, and immune infiltration data. We further performed in vitro assays to validate the effect of MAPK8IP2 on PCa cell proliferation, migration and invasion.
MAPK8IP2 is highly expressed in PCa tissues. Overexpression of MAPK8IP2 is associated with adverse clinicopathological factors and a poor prognosis in PCa. Receiver operating curve analysis showed that MAPK8IP2 can distinguish PCa tissues from non-PCa tissues with a certain accuracy (AUC = 0.814). The MAPK8IP2 genetic alteration rate was 2.6% and MAPK8IP2 alterations correlated with a poor prognosis. We also found that CDK12 and TP53 mutations were associated with MAPK8IP2 expression. The DNA methylation level of MAPK8IP2 was higher in primary tumors than in normal tissues, and the high MAPK8IP2 DNA methylation group of PCa patients had poor survival. Enrichment analysis indicated that MAPK8IP2 was involved in the MAPK signaling pathway. In vitro, knockdown of MAPK8IP2 inhibited PCa cell proliferation, migration and invasion.
MAPK8IP2 is a potential target for PCa treatment and can serve as a novel biomarker for PCa diagnosis and prognosis evaluation.
Prostate cancer (PCa) is the second most common malignant disease and the fifth leading cause of death in male cancer patients, accounting for nearly 3.8% of deaths and 7.3% of new cancer cases worldwide in 2020 . PCa can be effectively treated by surgery, radiotherapy and hormone therapy . However, tumor metastasis or castration resistance leads to a poor prognosis. Although the prognosis of prostate cancer patients has improved with the development of prostate cancer diagnosis and treatment strategies, it remains unsatisfactory . Thus, there is a need for further research to explore new biomarkers to predict the occurrence and progression of prostate cancer.
The mitogen-activated protein kinase (MAPK) signaling pathway is one of the most widely studied signal transduction pathways. As one of the core members of the MAPK family, c-Jun NH2-terminal kinase (JNK) is implicated in cell migration, differentiation, proliferation, and apoptosis, and is also closely related to tumor generation and progression [4, 5]. JNK-interacting protein (JIP) plays a critical role in the MAPK signaling pathway .
Mitogen-activated protein kinase 8 interacting protein 2 (MAPK8IP2) (also known as JIP2) is located on chromosome 22 (22q13), and is involved in the regulation of the JNK and P38 MAPK signaling pathways . MAPK8IP2 has a similar structure as JIP1, including an SH3 domain, a PTB domain, and a JNK-binding domain [6, 7]. Northern blot analysis demonstrated that MAPK8IP2 is expressed not only in the brain but also in other human tissues, such as the pancreas, ovary, and prostate . Although studies have shown that the JNK and P38 MAPK pathways are associated with tumors , only a few studies have reported on the role of MAPK8IP2 in cancer. Recently, it was founded that MAPK8IP2 may be related to the prognosis of the glioblastoma and pancreatic cancer [9, 10]. Furthermore, MAPK8IP2 was identified to be associated with the progression of cervical cancer . However, the function of MAPK8IP2 in PCa is unclear.
In this study, we analyzed the expression of MAPK8IP2 in PCa samples from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Then, we evaluated the association between the expression level of MAPK8IP2 and clinicopathological parameters in PCa and assessed the diagnostic and prognostic value of MAPK8IP2 for PCa. In addition, to determine the potential molecular mechanism of MAPK8IP2 in prostate cancer, we used various databases to analyze the relationship of MAPK8IP2 expression with immune infiltration, gene mutation, m6A modification, and DNA methylation and explored the gene and functional network correlated with MAPK8IP2 expression. Finally, we confirmed that knockdown of MAPK8IP2 inhibited the proliferation, migration and invasion of PCa cell lines in vitro. Our study uncovered the important role of MAPK8IP2 in PCa and offers a novel potential biomarker for improving the diagnosis and prognosis of PCa.
Materials and methods
Expression and clinical correlation analysis of MAPK8IP2
Then, the expression of MAPK8IP2 and the relationship between MAPK8IP2 expression level and clinicopathological factors of prostate cancer patients were analyzed by using PRAD data from TCGA database (https://portal.gdc.cancer.gov/). The GSE70768 dataset from the GEO database (https://www.ncbi.nlm.nih.gov/gds) was also used to analyze the expression of MAPK8IP2 in PCa. In addition, we compared the expression of MAPK8IP2 in tumor tissues with normal tissues using the Human Protein Atlas (HPA) database(https://www.proteinatlas.org/) [13, 14].
Genetic alteration and methylation analysis of MAPK8IP2
The genetic alterations of MAPK8IP2 in PCa were assessed using three datasets (MCTP, Nature 2012; SU2C/PCF Dream Team, PNAS 2019; TCGA, PanCancer Atlas) containing survival data in cBioPortal (www.cbioportal.org) [15, 16]. We used the muTarget (https://www.mutarget.com/)  platform to identify gene mutations altering MAPK8IP2 expression level and to identify the expression changes of genes associated with MAPK8IP2 mutation.
We analyzed the correlation between MAPK8IP2 and DNA methylation using the TCGA(PRAD) dataset from the UCSC Xena database. The MAPK8IP2 promoter methylation profile based on sample types, patient age, patient race, and nodal metastasis status was assessed using the TCGA(PRAD) dataset in the UALCAN database (http://ualcan.path.uab.edu/) . Then, the value of MAPK8IP2 DNA methylation level as a survival prediction and diagnosis tool for prostate cancer was assessed based on the DNMIVD database (http://www.unimd.org/dnmivd/) .
Additionally, the relationship between MAPK8IP2 and m6A RNA methylation -related genes was also analyzed in PRAD samples from the TCGA database, and the differential expression of m6A RNA methylation-related genes between the low and high MAPK8IP2 expression groups was assessed. Twenty-one m6A RNA methylation-related genes  were included in the study.
Immune infiltration analysis of MAPK8IP2
The single sample gene set enrichment algorithm (ssGSEA) from the R package “GSVA”  was used for immune infiltration analysis for PCa in the TCGA (PRAD) dataset. In addition, the association of MAPK8IP2 expression with the infiltration levels of twenty-four immune cells was evaluated. These 24 immune cell types are based on the literature reports .
Protein–protein interaction (PPI) analysis and enrichment analysis of MAPK8IP2-related genes
We employed the STRING database (https://cn.string-db.org/)  to construct the PPI network. An interaction score > 0.9 was set to obtain the protein that interacts with MAPK8IP2. Nine genes interacting with MAPK8IP2 were obtained for Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis [24, 25]. The R package “ClusterProfiler” was used to GO and KEGG analysis . Then, we used the TCGA (PRAD) dataset to analyze the correlation between MAPK8IP2 and nine partner genes and performed Kaplan–Meier analysis for patients grouped according to the expression of the nine partner genes. The differentially expressed genes correlated with MAPK8IP2 and MAPK8IP3 were analyzed via the LinkedOmics database (http://www.linkedomics.org/login.php) .
Cell culture and transfection
The human normal prostate epithelial cell lines (RWPE-1 and HPrEC) and the human PCa cell lines (PC3, DU145, 22RV1, LNCaP, and VCaP) were obtained from ATCC. PC3 cells were grown with F12K medium (Gibco, USA). HPrECs were maintained in prostate epithelial cell basal medium (ATCC, USA). RWPE-1 and VCaP cells were cultured in DMEM (Gibco, USA). DU145, 22RV1, and LNCaP cells were maintained in RPMI-1640 medium (Gibco, USA). All the cells were grown in media with 10% FBS (Gibco, Australia) at 37 °C in a 5% CO2 incubator. MAPK8IP2 siRNA was synthesized by RiboBio (Guangzhou, China). The siRNA sequence was as follows: siMAPK8IP2: 5′-GCCATTTCTTCCAGATGAA-3′. Transfection was performed with Lipofectamine 2000 (Invitrogen, Thermo Fisher Scientific) according to the manufacturer’s instructions.
Quantitative real-time PCR (qRT–PCR)
Total RNA was extracted by TRIzol reagent (CWBIO, China) and reverse transcribed into cDNA using the First-Strand cDNA Synthesis kit (Transgen, AT341-01) according to the manufacturer’s protocols. TransStart Green qPCR SuperMix (Transgen, AQ101-01) was used for qRT–PCR. Based on the 2-∆∆CT method and using β-actin as a reference gene, relative expression was calculated. The primer sequences were as follows: β-actin: forward, 5′-TCTCCCAAGTCCACACAGG-3′ and reverse, 5′-GGCACGAAGGCTCATCA-3′; MAPK8IP2: forward, 5′-CGCTGCAGCCATTT.
CTTCC-3′ and reverse, 5′-ACTCCTGGGAGACAAAGACG-3′.
Cell proliferation assays
The proliferation of cells was tested with 5-ethynyl-2′-deoxyuridine (EdU) and Cell Counting Kit 8 (CCK-8). Two days after siRNA transfection, PCa cells were seeded in 96-well plates for EdU and CCK-8 assays. PC3 (1 × 104 cells/well), DU145 (1× 104 cells/well), and 22RV1 (3 × 104 cells/well) cells were seeded in 96-well plates, After the cells were cultured for 1 day, EdU staining was performed following the EdU kit (RiboBio, China) manufacturer’s guidelines. Cell proliferation was observed using fluorescence microscopy, and the percentage of EdU-positive cells was calculated. PC3 (8 × 103 cells/well), DU145 (6× 103 cells/well), and 22RV1 (5 × 104 cells/well) cells were seeded in 96-well plates, 100 μl of medium containing 10% CCK-8 (APExBIO, USA) was added to each well at 24 h, 48 h, 72 h, and 96 h. After incubation for 2 h at 37 °C, the OD values at 450 nm were measured by a microplate reader.
Transwell chambers (8 μm pores, Corning, USA) were used for migration and invasion experiments. Then, 600 μl of complete medium containing 20% FBS was added to the bottom chamber. The transfected PC3 (6 × 104 cells/well), DU145 (5 × 104 cells/well), and 22RV1 (1.2 × 105 cells/well) cells were seeded in the upper chamber in 200 μl FBS-free medium for cell migration (without Matrigel) and invasion (with Matrigel). After being cultured for 48 h (for PC3 and 22RV1 cells) or 24 h (for DU145 cells), cells that passed through the transwell membrane were fixed, stained and observed under a microscope.
Western blot assay
Total protein was extracted with RIPA lysis buffer (APExBIO, USA) containing protease and phosphatase inhibitors and quantified with a BCA kit (CWBIO, China). Then, 10% SDS-PAGE was used to separate the proteins, and the transfer of proteins to PVDF membranes was performed. The membranes were incubated overnight at 4 °C with primary antibodies. The primary antibodies included GAPDH (Servicebio, China, 1:1000) and MAPK8IP2 (Santa Cruz, USA, 1:200).
Statistical analysis of data from the TCGA and GEO databases was performed using R (v 3.6.3). MAPK8IP2 functional experimental data were statistically analyzed using GraphPad Prism (v 8.0.1). Continuous variables are summarized as the mean (SD), and categorical variables are shown as the number (%). The chi-squared test, Fisher’s test, and Mann-Whitney U test were used to analyze the association between MAPK8IP2 expression levels and clinicopathological factors. To identify significant differences between groups of MAPK8IP2 functional experiments, Student’s t-test, and ANOVA were used. Spearman’s rank test was used to assess the correlation between the expression of two genes. Univariate and multivariate Cox regression analyses revealed significant predictors of poor overall survival (OS) and short progression-free interval (PFI). Estimating OS and PFI were performed using Kaplan–Meier analysis, and the statistical significance was evaluated by log-rank test. The predictive ability of MAPK8IP2 to predict PCa outcomes was assessed by receiver operating characteristic (ROC) curve analysis. p < 0.05 was considered to indicate statistical significance.
The expression level of MAPK8IP2 in PCa
We used TCGA (tumor and normal data) and GTEx (normal data) cohorts to analyze the expression of MAPK8IP2 in 33 cancer types (tumor vs. normal) (Fig. 1A). The results showed significantly higher expression of MAPK8IP2 in 15 types of human cancer, including prostate adenocarcinoma (PRAD), breast invasive carcinoma (BRCA), and bladder urothelial carcinoma (BLCA) (cancer tissues vs. normal tissues). Then, we further evaluated MAPK8IP2 expression in prostate cancer by analyzing TCGA (PRAD) and GEO (GSE70768) datasets, [the TCGA included unpaired samples (52 normal tissues vs. 499 tumor tissues) (Fig. 1B) and paired samples (52 pairs) (Fig. 1C); GSE70786 included unpaired samples (73 normal tissues vs. 125 tumor tissues) (Fig. 1D) and paired samples (73 pairs) (Fig. 1E)]. In the analysis of paired and unpaired samples, MAPK8IP2 expression was significantly lower tin normal tissues than in tumor tissues. Moreover, MAPK8IP2 protein expression in tumor and normal tissues was analyzed via the HPA database. We found positive staining for MAPK8IP2 in most PCa tissues. As shown Fig. 1F (Patient ID: 2932) and Fig. 1G (Patient ID: 525), MAPK8IP2 staining was significantly higher tin PRAD tissue than in normal tissue.
High expression of MAPK8IP2 was associated with adverse clinicopathological factors and a poor prognosis in PCa
According to the expression levels of MAPK8IP2, we categorized the 499 PCa patients in the TCGA (PRAD) dataset into low (n = 249) and high (n = 250) expression groups according to the median value. The expression of MAPK8IP2 was significantly correlated with Gleason score (P < 0.001), T stage (P < 0.001), N stage (P = 0.002), primary therapy outcome (P = 0.048), and residual tumor (P = 0.009), but not with age, prostate specific antigen (PSA), M stage, race, and zone of origin, as shown in Table 1. Furthermore, patients were grouped according to clinicopathological factors, and we found that patients with PSA (≥ 4 ng/ml), Gleason score (9-10), primary treatment outcome (progressive disease (PD)/stable disease (SD)/partial response (PR)), T stage (T3-T4), N stage (N1) and residual tumor (R1-R2) showed higher levels of MAPK8IP2 expression than patients with other features, as shown in Fig. 2A-F.
We further performed Kaplan–Meier analysis to estimate the relationship between MAPK8IP2 expression and survival. The results indicated that high expression of MAPK8IP2 was significantly associated with short PFI (Fig. 2G) and OS (Fig. 2H) in PCa patients.
In addition, Cox regression analysis was used to determine whether MAPK8IP2 was an independent prognostic factor for predicting PFI and OS. Univariate Cox analysis associated with PFI showed that T stage, N stage, Gleason score, PSA, primary therapy outcome, residual tumor, and MAPK8IP2 affected PFI (Table 2). Multivariate Cox analysis demonstrated that Gleason score 9-10 (HR, 2.375; 95% CI 1.382-4.082; P = 0.002), primary treatment outcome PD/SD/PR (HR, 3.586; 95% CI 2.026-6.348; P < 0.001), and high MAPK8IP2 expression group (HR, 1.690; 95% CI 1.045-2.733; P = 0.032) were significant risk factors for PFI (Table 2). As shown in Table 3, the univariate Cox analysis showed that M stage, Gleason score, PSA, primary therapy outcome, and MAPK8IP2 were associated with OS.Multivariate Cox analysis also showed that high MAPK8IP2 expression group (HR, 1.690; 95% CI 1.045-2.733; P = 0.032) was a significant independent risk factors for OS, besides M1 stage.
Next, we used ROC curves to evaluate the ability of MAPK8IP2 to discriminate between PCa sample and non-PCa samples. Figure 2I shows that MAPK8IP2 can serve as a diagnostic biomarker to distinguish PCa from non-PCa tissues with a certain accuracy (AUC = 0.814, CI: 0.779-0.850).
MAPK8IP2 genetic alteration in PCa patients
A total of 999 PCa patients from three datasets containing survival data were evaluated for MAPK8IP2 genetic alterations. The proportion of MAPK8IP2 genetic alteration in PCa was 2.6% (Fig. 3A), which consisted of missense mutation (Fig. 3B), amplification, and deep deletion were 0.2% (2/999), 0.8% (8/999), and 1.6% (16/999), respectively. In three datasets (MCTP, Nature 2012; SU2C/PCF Dream Team, PNAS 2019; and TCGA, PanCancer Atlas), the MAPK8IP2 genetic alteration rates were 0.81% (4/494), 4.28% (19/444) 4.92% (3/61), respectively (Fig. 3C). Kaplan–Meier plotters showed that MAPK8IP2 altered group was correlated with poor progression-free survival (PFS) (log-rank test: P = 0.0358) (Fig. 3F), but was not significantly correlated with disease specific survival (DSS) (Fig. 3D) and OS (Fig. 3E). In addition, through the muTarget platform, MAPK8IP2 mutations did not affect the expression of other genes, but the CDK12 and TP53 mutation groups were found to have higher MAPK8IP2 expression levels than the wild-type group (Fig. 3G, H).
DNA methylation of MAPK8IP2 in PCa patients
We used the UCSC Xena and UALCAN databases to compar the DNA methylation levels of MAPK8IP2 between normal and tumor tissues in PCa. The heat map (Fig. 4A) suggested that the methylation levels of most CpG sites in primary tumors were higher than those in normal tissues, and the box plot (Fig. 4B) analysis also yielded the same result that primary tumor tissues had higher MAPK8IP2 DNA methylation levels than normal tissues. Subgroup analysis indicated that MAPK8IP2 DNA methylation levels were also upregulated in various groups based on race, age, and nodal metastasis status (Fig. 4C-E).
We further used the DNMIVD database to investigate the diagnostic value of CpG sites located within MAPK8IP2 to distinguish tumor samples from normal samples and to perform survival analysis of PCa patients in the low and high groups of DNA methylation of MAPK8IP2. The diagnostic model (Fig. 4F-H) results showed that four CpG sites (cg09877744, cg20995573, cg17852224, and cg09764221) were identified as potential markers for the diagnosis of PCa (the AUC of the ROC curve was 0.810, Fig. 4G), and cg09877744 had the highest importance score (bar plot, Fig. 4F) and DNA methylation levels (clustering heatmap, Fig. 4H). Kaplan–Meier survival analysis showed that the high MAPK8IP2 methylation group had a shorter disease free interval (DFI) (P = 4.26e-03, Fig. 4I) and PFI (P = 3.70e-03, Fig. 4J) than the low MAPK8IP2 methylation group.
Relationship between MAPK8IP2 expression and the expression of m6A RNA methylation related genes in PCa
The TCGA PRAD dataset was used to analyze the association between MAPK8IP2 expression and 21 m6A RNA methylation-related genes in PCa. The results showed that 9 m6A-related genes correlated with MAPK8IP2 expression (Fig. 5A). The scatter plot (Fig. 5B) shows 9 significantly positively correlated m6A-related genes including METTL3 (r = 0.218, P < 0.001), RBM15B (r = 0.151, P < 0.001), VIRMA (r = 0.122, P = 0.007), YTHDF1 (r = 0.220, P < 0.001), YTHDF2 (r = 0.119, P = 0.008), YTHDF3 (r = 0.107, P = 0.017), HNRNPC (r = 0.119, P = 0.008), HNRNPA2B1 (r = 0.226, P < 0.001), and IGF2BP3 (r = 0.138, P = 0.002). Then, a total of 499 samples were divided into low (n = 249) and high (n = 250) MAPK8IP2 expression group. We compared the expression of 21 m6A-related genes between the low and high MAPK8IP2 expression groups, and found significantly increased expression of METTL3, RBM15B, VIRMA, YTHDF1, YTHDF2, HNRNPC, HNRNPA2B1, and IGF2BP3 in the high expression group (P < 0.05) (Fig. 5C).
Correlation between MAPK8IP2 expression and immune infiltration in PCa
We analyzed the correlation of MAPK8IP2 expression with immune infiltration by using the ssGSEA method. As shown in Fig. 6A, there was a negatively correlation between MAPKJ8IP2 expression and the infiltration levels of most immune cells. Then, a total of 499 samples patients from the TCGA PRAD cohort were divided into low (n = 249) and high (n = 250) MAPK8IP2 expression groups. We determined the correlation coefficients between the levels of different immune cell types and MAPK8IP2 expression with |r| > 0.1 as the cut off and found that the infiltration levels of B cells, mast cells, Th1 cells and T helper cells were significantly lower in the high expression group (P < 0.05) (Fig. 6B-E). A scatter plot was used to display the correlations between MAPK8IP2 expression and the infiltration levels of 4 types of immune cells (Fig. 6F-I).
MAPK8IP2-interacting protein analysis and GO and KEGG enrichment analysis
We constructed the PPI network of MAPK8IP2 via the STRING database, and the top 9 functional partner genes were identified: MAPK8IP1, MAP2K7, MAP3K11, MAPK8IP3, MAPK10, MAPK8, MAPK9, FGF12, and LRP8 (Fig. 7A). GO analysis showed that MAPK8IP2 and partner genes are mainly involved in kinesin binding, neuronal cell body, growth cone, and stress-activated protein kinesin signaling cascade (Fig. 7B). KEGG analysis revealed that MAPK8IP2 and partner genes were mainly correlated with the MAPK signaling pathway (Fig. 7C). Then, the correlation between MAPK8IP2 and 9 partner genes was analyzed. The results showed that MAPK8IP2 expression was significantly positively correlated with the expression of FGF12, LRP8, MAP2K7, MAP3K11, MAPK8IP1, MAPK8IP3, and MAPK9 (Fig. 7D-J). We further performed Kaplan–Meier survival analysis for MAPK8IP2-related genes, grouped by median value. The results showed that the MAPK8IP2 related genes not significant difference in survival, besides MAP2K7, MAP3K11, MAPK8IP1 and MAPK8IP3. As shown in Fig. 7K-O and Fig. S1A-D, the high MAPK8IP1 expression group had poor OS (Fig. 7K) and PFI (Fig. 7M), and the high MAP2K7, MAP3K11, and MAPK8IP3 expression groups had a short PFI (Fig. S1A, C; Fig. 7O) but no significant difference in OS (Fig. S1B, D; Fig. 7N). In addition, we further analyzed the expression of MAPK8IP2 and its partner genes through public databases. Univariate and multivariate Cox regression analyses were performed on MAPK8IP2 and its partner genes (MAPK8IP1, MAPK8IP3, MAP2K7, and MAP3K11), and the results showed that MAPK8IP2 and MAPK8IP3 were significantly associated with the progression-free interval of PCa patients (Table S4). Then, the differentially expressed genes correlated with MAPK8IP2 and MAPK8IP3 were analyzed via the LinkedOmics database. We identified MAPK8IP2 and MAPK8IP3 related genes by Spearman correlation analysis (Fig. S4A, D). The heatmaps show the top 50 positively correlated genes and the top 50 negatively correlated genes (Fig. S4B, C, E, F). We selected those with positive Spearman correlation coefficients > 0.4 for the analysis and found that only PDIA2 was positively correlated with both MAPK8IP1 and MAPK8IP3 (Fig. S4G).
SiRNA-MAPK8IP2 transfection inhibited the proliferation, migration, and invasion of PCa cells in vitro
The expression levels of MAPK8IP2 were detected by qRT–PCR, and we found that MAPK8IP2 mRNA expression levels were significantly higher in PCa cell lines (PC3, DU145, 22RV1, LNCAP, and VCAP) than in normal prostate epithelial cells (HPrEC and RWPE1) (Fig. 8A). Western blot analysis demonstrated that MAPK8IP2 protein was highly expressed in PC3, DU145, and 22RV1 cell lines compared to HPrEC cell line (Fig. 8B).
Therefore, we further evaluated the biological function of MAPK8IP2 in PCa cells, and MAPK8IP2 specific siRNA was transfected into PC3, DU145, and 22RV1 cells. MAPK8IP2 mRNA and protein expression levels were significantly downregulated after transfection with siRNA-MAPK8IP2 (Fig. 8C, D). EdU and CCK-8 assays showed that downregulation of MAPK8IP2 efficiently inhibited the growth of PC3, DU145, and 22RV1 cells (Fig. 8E, F). Then, transwell assays were performed to assess the effect of MAPK8IP2 on PC3, DU145, and 22RV1 cells migration and invasion. The results revealed that knockdown of MAPK8IP2 significantly suppressed the migration and invasion capabilities of PC3, DU145, and 22RV1 cells (Fig. 8G).
Studies have demonstrated that JNK and P38 MAPK signaling are closely related to human cancers. MAPK8IP2 is member of the JIP group of MAPK scaffolding proteins and can interact with JNK and p38 [6, 28, 29]. We therefore considered that MAPK8IP2 may play an important role in the development of cancer.
In this study, we found that MAPK8IP2 was upregulated in most cancer tissues by analyzing the TCGA and GTEx cohorts of the UCSC database. Based on the TCGA-PRAD and GSE70768 datasets, we confirmed that MAPK8IP2 expression in PCa tissues was significantly higher than that in normal tissues (paired or unpaired). High expression of MAPK8IP2 was associated with adverse clinicopathological factors and with poor survival in PCa patients. Multivariate Cox analysis demonstrated that MAPK8IP2 expression level was a significant independent prognostic factor terms of OS and PFI. The ROC curve analysis showed that MAPK8IP2 had high diagnostic value for PCa.
Genetic alterations are closely correlated with cancers. Several genetic alterations have been shown to be associated with the initiation and development of PCa . We found that the MAPK8IP2 genetic alteration frequency was only 2.6% in PCa, mainly deep deletion. Although MAPK8IP2 genetic alteration was not significantly correlated with DSS and OS, MAPK8IP2 altered group had significantly poorer PFS than the unaltered group. Considering the small sample size of the MAPK8IP2 altered group, more data are needed for validation. In addition, we found that the CDK12 and TP53 mutation groups had higher expression levels of MAPK8IP2 than the wild-type group. A previous study demonstrated that advanced PCa patients with CDK12 mutations had worse clinical characteristics and were more likely to progress . Patients with metastatic castration-resistant prostate cancer with CDK12 mutations had a worse prognosis than those without CDK12 mutations . TP53 is also one of the most commonly mutated genes in prostate cancer. TP53 mutations were correlated with a poor outcome in PCa patients . The detection of TP53 alterations has clinical value for guiding the precise treatment of hormone-native prostate cancer . These results showed that the prognosis of PCa patients was closely associated with MAPK8IP2 genetic mutation and MAPK8IP2-related genetic mutation.
DNA methylation is a common epigenetic mechanism, and aberrant DNA methylation is strongly associated with cancer. In PCa, aberrant methylation patterns are correlated with tumorigenesis and progression [35, 36]. In this study, we discovered that MAPK8IP2 methylation levels in prostate cancer tissues were higher than in normal tissues. In PCa, survival was worse in the high methylation level group than in the low methylation level group of MAPK8IP2. The ROC analysis showed that four CpG sites located within MAPK8IP2 have better reliability for diagnosing PCa. A recent study found that the expression level of MAPK8IP2 was significantly negatively correlated with DNA methylation in pancreatic cancer (R = -0.507) . However, our results suggest that the expression level of MAPK8IP2 may be positively correlated with DNA methylation in PCa. This regulatory mechanism may differ from the classical transcriptional silencing mechanism of DNA methylation, but is associated with DNA methylation-induced transcriptional activation .
m6A methylation is the most abundant epigenetic modification in eukaryotic mRNA. It plays an important role in cancer pathogenesis and development . Li et al.  demonstrated that YTHDF2 upregulates pAKT by inducing mRNA degradation and promotes prostate cancer progression. Chen et al.  found that silencing METTL3 inhibits the invasion and metastasis of prostate cancer cells. Wen et al.  showed that m6A modification of lncRNA NEAT1 promotes bone metastasis in PCa. In this study, we found that MAPK8IP2 expression was significantly positively associated with METTL3, RBM15B, VIRMA, YTHDF1, YTHDF2, YTHDF3, HNRNPC, HNRNPA2B1, and IGF2BP3 expression; we also found significantly increased expression of METTL3, RBM15B, VIRMA, YTHDF1, YTHDF2, HNRNPC, HNRNPA2B1, and IGF2BP3 in the high MAPK8IP2 expression group. We speculate that MAPK8IP2 is likely to be modified by m6A to promote PCa progression.
In the analysis of immune infiltration, we discovered that MAPK8IP2 expression was negatively related to the infiltration of some immune cells, mainly including B cells, mast cells, Th1 cells, and T helper cells. However, the correlation was not strong enough. The results indicated that immune infiltration may be associated with MAPK8IP2 to promote PCa progression, but is not a major regulatory factor.
In addition, we further performed GO and KEGG analyses on MAPK8IP2 and 9 partner genes. GO and KEGG analyses revealed that these genes are mainly involved in the stress-activated protein kinesin signaling cascade and MAPK signaling pathway. JNK and P38 activation has been shown to promote prostate cancer cell migration and invasion [41,42,43,44]. However, they seem to play a dual role in the proliferation of PCa cells. JNK and P38 activation have been reported to promote PCa cells proliferation [45, 46]. Conversely, Xie et al.  found that activation of JNK inhibits the proliferation of PCa cells, and Zhang et al.  showed that inhibition of P38 promotes PCa cell proliferation. These results indicated that MAPK8IP2 may mediate the JNK and P38 signaling pathways to promote prostate cancer cell proliferation, migration and invasion. Furthermore, Zhao et al.  demonstrated that MAPK8IP2 is regulated by E6 to promote cervical cancer progression via the noncanonical WNT pathway. Therefore, we think that MAPK8IP2 may be involved in other noncanonical pathways to promote PCa progression, besides, via the canonical MAPK signaling pathway. We also found that seven partner genes were significantly positively correlated with MAPK8IP2 expression, but only MAP2K7, MAP3K11, MAPK8IP1 and MAPK8IP3 were significantly associated with survival. Overexpression of MAP2K7, MAP3K11, MAPK8IP1 and MAPK8IP3 was also correlated with worse survival in PCa. Further multivariate Cox regression analysis showed that MAPK8IP2 and MAPK8IP3 were significantly associated with the progression-free interval of PCa patients. MAPK8IP3 like MAPK8IP2 gene belong to the JIPs family and are involved in the regulation of the JNK signaling pathway . Previous studies indicated that MAPK8IP3 may be associated with the malignant progression of prostate cancer . Therefore, we consider that MAPK8IP2and MAPK8IP3 have synergistic effects in promoting PCa progression. The synergistic effect is associated with the PDIA2 gene.
Finally, we confirmed that MAPK8IP2 expression was significantly higher in PCa cell lines than in normal prostate epithelial cells. The proliferation, migration and invasion of prostate cancer cells were inhibited after transfection with siRNA-MAPK8IP2 in vitro.
Several limitations of our research should be recognized. First, potential mechanisms associated with MAPK8IP2 were obtained by bioinformatics analysis, and mechanism validation assays were not performed. Second, we conducted functional assays of PCa cells in vitro, but not in vivo.
Overall, our study revealed that MAPK8IP2 is highly expressed in PCa. High MAPK8IP2 expression is correlated with a poor prognosis in PCa patients. Furthermore, DNA methylation of MAPK8IP2 and MAPK8IP2 genetic alteration were associated with the prognosis of PCa. MAPK8IP2 might be involved in the regulation of the JNK and P38 MAPK signaling pathways and be modified by m6A. MAPK8IP2 is a promising biomarker for PCa diagnosis, treatment and prognosis evaluation, but more experiments are needed for validation.
Availability of data and materials
The datasets generated or analyzed during this study are available in public database. Public databases are all open access (links are included in the methods section of this article). All data supporting this study are included within the article (and its supplementary information files). All data are available from the corresponding author upon reasonable request.
mitogen-activated protein kinase
c-Jun NH2-terminal kinase
mitogen-activated protein kinase 8 interacting protein 2
Human Protein Atlas
The Cancer Genome Atlas
Gene Expression Omnibus
Single Sample Gene Set Enrichment algorithm
Kyoto Encyclopedia of Genes and Genomes
quantitative real-time PCR
Cell Counting Kit 8
Sung H, Ferlay J, Siegel RL, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.
Parker C, Castro E, Fizazi K, et al. Prostate cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2020;31(9):1119–34. https://doi.org/10.1016/j.annonc.2020.06.011.
Litwin MS, Tan HJ. The diagnosis and treatment of prostate Cancer: a review. JAMA. 2017;317(24):2532–42. https://doi.org/10.1001/jama.2017.7248.
Davis RJ. Signal transduction by the JNK group of MAP kinases. Cell. 2000;103(2):239–52. https://doi.org/10.1016/s0092-8674(00)00116-1.
Dhanasekaran DN, Kashef K, Lee CM, et al. Scaffold proteins of MAP-kinase modules. Oncogene. 2007;26(22):3185–202. https://doi.org/10.1038/sj.onc.1210411.
Yasuda J, Whitmarsh AJ, Cavanagh J, et al. The JIP group of mitogen-activated protein kinase scaffold proteins. Mol Cell Biol. 1999;19(10):7245–54. https://doi.org/10.1128/MCB.19.10.7245.
Negri S, Oberson A, Steinmann M, et al. cDNA cloning and mapping of a novel islet-brain/JNK-interacting protein. Genomics. 2000;64(3):324–30. https://doi.org/10.1006/geno.2000.6129.
Wagner EF, Nebreda AR. Signal integration by JNK and p38 MAPK pathways in cancer development. Nat Rev Cancer. 2009;9(8):537–49. https://doi.org/10.1038/nrc2694.
Hsu JB, Lee TY, Cheng SJ, et al. Identification of differentially expressed genes in different glioblastoma regions and their association with Cancer stem cell development and Temozolomide response. J Pers Med. 2021;11(11):1047. https://doi.org/10.3390/jpm11111047.
Jia Y, Shen M, Zhou Y, Liu H. Development of a 12-biomarkers-based prognostic model for pancreatic cancer using multi-omics integrated analysis. Acta Biochim Pol. 2020;67(4):501–8. https://doi.org/10.18388/abp.2020_5225.
Zhao L, Wang L, Zhang C, et al. E6-induced selective translation of WNT4 and JIP2 promotes the progression of cervical cancer via a noncanonical WNT signaling pathway. Signal Transduct Target Ther. 2019;4:32. https://doi.org/10.1038/s41392-019-0060-y.
Goldman MJ, Craft B, Hastie M, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38(6):675–8. https://doi.org/10.1038/s41587-020-0546-8.
Uhlén M, Fagerberg L, Hallström BM, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419. https://doi.org/10.1126/science.1260419.
Uhlen M, Zhang C, Lee S, et al. A pathology atlas of the human cancer transcriptome. Science. 2017;357(6352):eaan2507. https://doi.org/10.1126/science.aan2507.
Sander C, Schultz N. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4. https://doi.org/10.1158/2159-8290.CD-12-0095 Erratum in: Cancer Discov. 2012;2(10):960.
Gao J, Aksoy BA, Dogrusoz U, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1. https://doi.org/10.1126/scisignal.2004088.
Nagy Á, Győrffy B. muTarget: a platform linking gene expression changes and mutation status in solid tumors. Int J Cancer. 2021;148(2):502–11. https://doi.org/10.1002/ijc.33283.
Chandrashekar DS, Bashel B, Balasubramanya SAH, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58. https://doi.org/10.1016/j.neo.2017.05.002.
Ding W, Chen J, Feng G, et al. DNMIVD: DNA methylation interactive visualization database. Nucleic Acids Res. 2020;48(D1):D856–62. https://doi.org/10.1093/nar/gkz830.
Sun T, Wu R, Ming L. The role of m6A RNA methylation in cancer. Biomed Pharmacother. 2019;112:108613. https://doi.org/10.1016/j.biopha.2019.108613.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7. https://doi.org/10.1186/1471-2105-14-7.
Bindea G, Mlecnik B, Tosolini M, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95. https://doi.org/10.1016/j.immuni.2013.10.003.
Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13. https://doi.org/10.1093/nar/gky1131.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.
The Gene Ontology Consortium. The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(D1):D330–8. https://doi.org/10.1093/nar/gky1055.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63. https://doi.org/10.1093/nar/gkx1090.
Buchsbaum RJ, Connolly BA, Feig LA. Interaction of Rac exchange factors Tiam1 and Ras-GRF1 with a scaffold for the p38 mitogen-activated protein kinase cascade. Mol Cell Biol. 2002;22(12):4073–85. https://doi.org/10.1128/MCB.22.12.4073-4085.2002.
Schoorlemmer J, Goldfarb M. Fibroblast growth factor homologous factors and the islet brain-2 scaffold protein regulate activation of a stress-activated protein kinase. J Biol Chem. 2002;277(51):49111–9. https://doi.org/10.1074/jbc.M205520200.
Fontana F, Anselmi M, Limonta P. Molecular mechanisms and genetic alterations in prostate cancer: from diagnosis to targeted therapy. Cancer Lett. 2022;534:215619. https://doi.org/10.1016/j.canlet.2022.215619.
Reimers MA, Yip SM, Zhang L, et al. Clinical outcomes in cyclin-dependent kinase 12 mutant advanced prostate Cancer. Eur Urol. 2020;77(3):333–41. https://doi.org/10.1016/j.eururo.2019.09.036.
Rescigno P, Gurel B, Pereira R, et al. Characterizing CDK12-mutated prostate cancers. Clin Cancer Res. 2021;27(2):566–74. https://doi.org/10.1158/1078-0432.CCR-20-2371.
Nientiedt C, Budczies J, Endris V, et al. Mutations in TP53 or DNA damage repair genes define poor prognostic subgroups in primary prostate cancer. Urol Oncol. 2022;40(1):8.e11–8. https://doi.org/10.1016/j.urolonc.2021.06.024.
Liu Z, Guo H, Zhu Y, et al. TP53 alterations of hormone-naïve prostate cancer in the Chinese population. Prostate Cancer Prostatic Dis. 2021;24(2):482–91. https://doi.org/10.1038/s41391-020-00302-3.
Baylin SB, Esteller M, Rountree MR, et al. Aberrant patterns of DNA methylation, chromatin formation and gene expression in cancer. Hum Mol Genet. 2001;10(7):687–92. https://doi.org/10.1093/hmg/10.7.687.
Aryee MJ, Liu W, Engelmann JC, et al. DNA methylation alterations exhibit intraindividual stability and interindividual heterogeneity in prostate cancer metastases. Sci Transl Med. 2013;5(169):169ra10. https://doi.org/10.1126/scitranslmed.3005211.
Smith J, Sen S, Weeks RJ, et al. Promoter DNA Hypermethylation and paradoxical gene activation. Trends Cancer. 2020;6(5):392–406. https://doi.org/10.1016/j.trecan.2020.02.007.
Li J, Xie H, Ying Y, et al. YTHDF2 mediates the mRNA degradation of the tumor suppressors to induce AKT phosphorylation in N6-methyladenosine-dependent way in prostate cancer. Mol Cancer. 2020;19(1):152. https://doi.org/10.1186/s12943-020-01267-6.
Chen Y, Pan C, Wang X, et al. Silencing of METTL3 effectively hinders invasion and metastasis of prostate cancer cells. Theranostics. 2021;11(16):7640–57. https://doi.org/10.7150/thno.61178.
Wen S, Wei Y, Zen C, et al. Long non-coding RNA NEAT1 promotes bone metastasis of prostate cancer through N6-methyladenosine. Mol Cancer. 2020;19(1):171. https://doi.org/10.1186/s12943-020-01293-4.
Si-Tu J, Cai Y, Feng T, et al. Upregulated circular RNA circ-102004 that promotes cell proliferation in prostate cancer. Int J Biol Macromol. 2019;122:1235–43. https://doi.org/10.1016/j.ijbiomac.2018.09.076.
Chen J, Li L, Yang Z, et al. Androgen-deprivation therapy with enzalutamide enhances prostate cancer metastasis via decreasing the EPHB6 suppressor expression. Cancer Lett. 2017;408:155–63. https://doi.org/10.1016/j.canlet.2017.08.014.
Liu C, Jiang S, Xie H, et al. Long non-coding RNA AC245100.4 contributes to prostate cancer migration via regulating PAR2 and activating p38-MAPK pathway. Med Oncol. 2022;39(5):94. https://doi.org/10.1007/s12032-022-01689-w.
He J, Mao Y, Huang W, et al. Methylcrotonoyl-CoA carboxylase 2 promotes proliferation, migration and invasion and inhibits apoptosis of prostate Cancer cells through regulating GLUD1-P38 MAPK signaling pathway. Onco Targets Ther. 2020;13:7317–27. https://doi.org/10.2147/OTT.S249906.
El-Haibi CP, Singh R, Sharma PK, et al. CXCL13 mediates prostate cancer cell proliferation through JNK signalling and invasion through ERK activation. Cell Prolif. 2011;44(4):311–9. https://doi.org/10.1111/j.1365-2184.2011.00757.x.
Tiwari A, Mukherjee B, Hassan MK, et al. Reduced FRG1 expression promotes prostate cancer progression and affects prostate cancer cell migration and invasion. BMC Cancer. 2019;19(1):346. https://doi.org/10.1186/s12885-019-5509-4.
Xie H, Liu T, Chen J, et al. Activation of PSGR with β-ionone suppresses prostate cancer progression by blocking androgen receptor nuclear translocation. Cancer Lett. 2019;453:193–205. https://doi.org/10.1016/j.canlet.2019.03.044.
Zhang L, Lu SY, Guo R, et al. STK10 knockout inhibits cell migration and promotes cell proliferation via modulating the activity of ERM and p38 MAPK in prostate cancer cells. Exp Ther Med. 2021;22(2):851. https://doi.org/10.3892/etm.2021.10283.
Engström W, Ward A, Moorwood K. The role of scaffold proteins in JNK signalling. Cell Prolif. 2010;43(1):56–66. https://doi.org/10.1111/j.1365-2184.2009.00654.x.
Cheng Y, Li L, Qin Z, et al. Identification of castration-resistant prostate cancer-related hub genes using weighted gene co-expression network analysis. J Cell Mol Med. 2020;24(14):8006–17. https://doi.org/10.1111/jcmm.15432.
This study was sponsored by the National Natural Science Foundation of China (No. 82060467) and the Key Research and Development Program of Jiangxi Province (No. 20171ACB20029).
Ethics approval and consent to participate
Consent for publication
All authors declare that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Kaplan-Meier curves comparing OS (A, C), and PFI (B, D) in low and high expression, groups of MAP2K7 and MAP3K11. OS: overall survival; PFS: progression-free survival.
Kaplan-Meier curves comparing OS (A, C), and PFI (B, D) in low and high expression groups of FGF12, LRP8 and MAPK9. OS: overall survival; PFS: progression-free survival.
The original images for Western blot in artical. The bands were visualized using the ImageQuant LAS 500.
Differentially expressed genes associated with MAPK8IP2 or MAPK8IP3 in PCa. (A) The correlation between MAPK8IP2 and differentially expressed genes. (B-C) Heat map showing genes positively or negatively correlated with MAPK8IP2 (the top 50 genes). (D) The correlation between MAPK8IP3 and differentially expressed genes. (E-F) Heat map showing genes positively or negatively correlated with MAPK8IP3 (the top 50 genes). (G) The Venn results displayed that only PDIA2 was positively associated with both MAPK8IP2 and MAPK8IP3.
(A) the expression level of MAPK8IP2 was significantly downregulated in the siMAPK8IP2 group compared with the siNC group and the blank group after transfection. (B-D) proliferation migration and invasion of PC3 cells were significantly inhibitor in the siMAPK8IP2 group compared with the siNC group and the blank group.
Browse interactions in tabular form.
Gene Ontology term enrichment analysis of the MAPK8IP2 gene and its partner genes.
KEGG enrichment analysis of the MAPK8IP2 gene and its partner genes.
Univariate and multivariate Cox regression analyses of MAPK8IP2 and its partners genes predict progression-free interval in PCa.
About this article
Cite this article
Zeng, Z., He, W., Jiang, Y. et al. MAPK8IP2 is a potential prognostic biomarker and promote tumor progression in prostate cancer. BMC Cancer 22, 1162 (2022). https://doi.org/10.1186/s12885-022-10259-2