Yi Liu1
Juan Xiang
Gang Peng
Chenfu Shen- 1Department of Neurosurgery, Xiangya Hospital, Central South University, Changsha, China
- 2Department of Geriatrics, Xiangya Hospital, Central South University, Changsha, China
PDZ-binding kinase (PBK) is known to regulate tumor progression in some cancer types. However, its relationship to immune cell infiltration and prognosis in different cancers is unclear. This was investigated in the present study by analyzing data from TCGA, GEO, GETx, TIMER, CPTAC, GEPIA2, cBioPortal, GSCALite, PROGNOSCAN, PharmacoDB, STRING, and ENCORI databases. PBK was overexpressed in most tumors including adenocortical carcinoma (hazard ratio [HR] = 2.178, p < 0.001), kidney renal clear cell carcinoma (KIRC; HR = 1.907, p < 0.001), kidney renal papillary cell carcinoma (HR = 3.024, p < 0.001), and lung adenocarcinoma (HR = 1.255, p < 0.001), in which it was associated with poor overall survival and advanced pathologic stage. PBK methylation level was a prognostic marker in thyroid carcinoma (THCA). PBK expression was positively correlated with the levels of BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, and TTK in 32 tumor types; and with the levels of the transcription factors E2F1 and MYC, which regulate apoptosis, the cell cycle, cell proliferation and invasion, tumorigenesis, and metastasis. It was also negatively regulated by the microRNAs hsa-miR-101-5p, hsa-miR-145-5p, and hsa-miR-5694. PBK expression in KIRC, liver hepatocellular carcinoma, THCA, and thymoma was positively correlated with the infiltration of immune cells including B cells, CD4+T cells, CD8+ T cells, macrophages, monocytes, and neutrophils. The results of the functional enrichment analysis suggested that PBK and related genes contribute to tumor development via cell cycle regulation. We also identified 20 drugs that potentially inhibit PBK expression. Thus, PBK is associated with survival outcome in a variety of cancers and may promote tumor development and progression by increasing immune cell infiltration into the tumor microenvironment. These findings indicate that PBK is a potential therapeutic target and has prognostic value in cancer treatment.
Introduction
PDZ-binding kinase (PBK), also known as T lymphokine-activated killer cell-derived protein kinase (TOPK), is a serine-threonine kinase of the mitogen-activated protein kinase (MAPKK) family that regulates cell survival, proliferation, growth, apoptosis, and inflammation (Abe et al., 2000; Gaudet et al., 2000; Zhao et al., 2020). As a mitotic kinase, PBK is activated by the cyclin-dependent kinase 1 (CDK1)/cyclin B1 complex during mitosis and promotes cytokinesis by phosphorylating polycomb inhibitory complex 1 (PRC1) or anaphase-promoting complex (APC) ubiquitin ligase (Kraft et al., 2003; Abe et al., 2007). PBK knockdown was shown to result in cytokinesis defects in cancer cells (Park et al., 2010).
PBK is overexpressed in a variety of active proliferating cells including malignant tumor cells as well as normal sperm cells. The upregulation of PBK has diagnostic and prognostic significance in cancer. For example, increased expression of PBK has been linked to poor prognosis in patients with ovarian serous cystadenocarcinoma (OV), esophageal carcinoma (ESCA), and stomach adenocarcinoma (STAD) (Ikeda et al., 2016; Ohashi et al., 2016; Ohashi et al., 2017). PBK is also related to cell transformation induced by H-Ras, activation of c-Jun-NH2 kinase (JNK), and expression of p53 induced by ultraviolet A, ultraviolet B, and DNA damage, respectively (Oh et al., 2007; Hu et al., 2010). PBK directly phosphorylates extracellular signal-regulated kinase (ERK), histone H3 (Ser10), histone H2AX (Ser139), peroxide-reduced protein 1 (PRX1, Ser32), JNK1 (Thr183/Tyr185), and p53-associated protein kinase (PRPK, Ser250) (Zykova et al., 2010; Roh et al., 2018) and activates downstream signaling pathways including MAPK, ribosomal S6 kinase (RSK), and transcription factors (TFs) such as activating protein 1 (AP-1) or nuclear factor kappa B (NF-κB) through phosphorylation, thereby promoting cell proliferation, migration, and invasion (Aksamitiene et al., 2010; Park et al., 2014; Ishikawa et al., 2018). Additionally, PBK promotes cell migration by regulating phosphatidylinositol 3-kinase (PI3K)/phosphatase and tensin homolog (PTEN)/protein kinase B (AKT) or transforming growth factor beta 1 (TGF-β1)/mothers against decapentaplegic (Smad) signaling pathways (Shih et al., 2012; Lee et al., 2020). PBK is differentially expressed between cancer and normal tissues and modulates radiosensitivity in tumors, suggesting that patients can be treated with radiotherapy (Pirovano et al., 2017); and a recent retrospective analysis showed that PBK expression was negatively correlated with the prognosis of patients with oral squamous cell carcinoma treated with radiotherapy (Yu et al., 2021). The PBK inhibitor OTS964/OTS514 was shown to inhibit cytokinesis and cancer cell growth (Matsuo et al., 2014), and OTS964/OTS514 induced the apoptosis of lung cancer cells including small cell lung cancer (SCLC) cells (Park et al., 2017). Thus, PBK likely plays an important role in cancer development.
Previous studies on PBK were limited to a single cancer type and lacked a comparative analysis of multiple tumors. The aim of the present study was to examine the role of PBK in different cancers and elucidate the molecular mechanisms and functions of PBK and its interacting molecules in carcinogenesis. We also analyzed the relationships between PBK expression, mutation and methylation status, clinical outcome, and tumor immune infiltration in cancers, and screened drugs that can potentially inhibit PBK. Our results demonstrate that PBK plays an important role in carcinogenesis and is a potential therapeutic target in cancer treatment.
Materials and Methods
Gene Expression Analysis
RNA expression and clinical data in The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases were downloaded from the University of California at Santa Cruz (UCSC) Xena database (https://xenabrowser.net/datapages/). The R “limma” package was used to analyze the mRNA expression of PBK in the pan-cancer cohort.
The UALCAN portal (http://ualcan.path.uab.edu/analysis-prot.html) was used to access the Clinical Proteome Oncology Analysis Alliance (CPTAC) dataset (Chen F. et al., 2019), which includes breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), OV, clear cell renal cell carcinoma, endometrial carcinoma (UCEC), lung adenocarcinoma (LUAD), and childhood brain cancer. We used these data to compare total protein expression between primary tumors and normal tissue by entering “PBK” into the search field.
We used the Gene Expression Profiling Interactive Analysis v2 (GEPIA2; http://gepia2.cancer-pku.cn/#index) “Pathological Stage Plot” module to obtain maps of PBK expression in different pathologic stages (stage I–IV) of TCGA tumors. Log2 (transcripts per million+1) expression data and an adjusted p value < 0.05 were used as cutoffs to generate the stage diagram.
Alterations in Gene Expression
We investigated the inheritance of PBK gene modifications in different cancers using the cBioPortal portal (https://www.cbioportal.org/) (Gao et al., 2013) through “TCGA Pan-Cancer Atlas Studies” in the “Quick select” section with “PBK” as query. The results of genetic changes in TCGA tumor samples were displayed in the “Cancer Types Summary” module, including the type of mutation and copy number alterations. The comparison module also returned disease specificity and overall survival (OS), progression-free survival (PFS), and disease-free survival (DFS) in TCGA cancer samples with or without PBK gene mutations. A log-rank p value < 0.05 was used as the cutoff.
PBK Regulatory Networks of TFs, microRNAs (miRNAs), and Methylation
In order to identify regulatory factors that affect PBK expression, we used Harmonizome (https://maayanlab.cloud/Harmonizome) (Rouillard et al., 2016) to predict TFs that can bind to the PBK gene promoter in the Chromatin Immunoprecipitation Enrichment Analysis (ChEA), Encyclopedia of DNA Elements (ENCODE), JASPAR, and TRANScription FACtor (TRANSFAC) databases. TFs that appeared at least twice in these databases were re-screened in the ENCORI database (http://starbase.sysu.edu.cn/panCancer.php) (p < 0.05, R = 1 to −1).
miRNAs targeting PBK were identified using the microRNA database (miRDB) (http://mirdb.org/cgi-bin/search.cgi) and TargetScan (http://www.targetscan.org/vert_72/). miRNAs that appeared at least twice were re-screened in ENCORI to determine their relationship to PBK (p < 0.05, R = 1 to −1).
Gene Set Cancer Analysis (GSCALite; http://bioinfo.life.hust.edu.cn/web/GSCALite/) is a bioinformatics platform that allows methylation, drug sensitivity, genetic, gene variation, and survival analyses of cancer datasets (Liu et al., 2018). We used the GSCALite platform to analyze the differential methylation of PBK and TFs between tumors and paired normal tissues in different types of cancer in TCGA. After logging onto the GSCALite server, we used the “TCGA Cancer–methylation” module to select 32 TCGA cancer types. In the final map generated by GSCALite, only 14 cancer types showed methylation profiles of PBK and related TFs. Using the same platform, we analyzed the relationship between methylation and the expression of PBK and related TFs in different cancers.
Survival Analysis
The PrognoScan database (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html) allows meta-analysis of the prognostic value of genes by comparing the relationship between gene expression and variables of interest in a large number of published reports (Mizuno et al., 2009). We used this platform to evaluate the relationship between PBK expression and patient prognosis in different cancer cohorts. Survival curves were generated with PrognoScan and the “survival” package in R from TCGA data. The survival results are shown as hazard ratio (HR) with 95% confidence interval (CI) and log-rank p value.
Correlation Analysis
Pearson correlation analysis was performed to evaluate the correlation between PBK and immune checkpoints, tumor mutation burden (TMB), microsatellite instability (MSI), DNA methylation, and mismatch repair (MMR). The “heatmap” package in R was used to display the results as a heatmap.
Immune Cell Infiltration Into the Tumor Microenvironment
The Tumor IMmune Estimation Resource (TIMER; http://timer.cistrome.org/) (Li et al., 2017) web server is a comprehensive resource for systematic analysis of immune cell infiltration in different cancer types. We used the immune gene module of TIMER to assess the correlation between PBK expression in six immune cell types and their infiltration in 32 TCGA tumors. We also used the Tumor-Immune System Interactions database (TISIDB) (http://cis.hku.hk/TISIDB/) to evaluate whether there was a significant difference in PBK expression between responders and non-responders to immunotherapies (e.g., anti-programmed death 1 [PD1] or anti-programmed death ligand 1 [PD-L1]).
Enrichment Analysis of PBK-Related Genes
We used the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) website (https://string-db.org/) with “PBK” and organism = Homo sapiens to construct a network of the interaction between PBK and related genes. We then used the GSCALite platform to analyze the function of PBK-related genes by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Using the ggplot2 package in R v3.6.2, a point map was generated to visualize the results (p < 0.05).
Drugs Targeting PBK
We downloaded data on PBK-related drugs from PharmacoDB (https://www.pharmacodb.pmgenomics.ca/) (Smirnov et al., 2018), with p < 0.01 and a standardized coefficient of −1 to 1 as thresholds. The “heatmap” package in R was used to display the results as a heatmap.
Results
Expression of PBK in Human Cancers
We analyzed the differential expression of PBK in 33 cancers as compared to normal tissue using TCGA and GTEx databases. PBK expression was higher in adrenocortical carcinoma (ACC), bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), glioblastoma multiforme (GBM), head and neck squamous cell carcinoma (HNSC), kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), brain lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), ovarian serous cystadenocarcinoma (OV), rectum adenocarcinoma (READ), prostate adenocarcinoma (PRAD), skin cutaneous melanoma (SKCM), stomach adenocarcinoma (STAD), thyroid carcinoma (THCA), uterine corpus endometrial carcinoma (UCEC), and uterine carcinosarcoma (UCS) but lower in acute myeloid leukemia (LAML) compared to normal tissues (all p < 0.01; Figure 1A).
FIGURE 1. PBK expression. (A) Expression of PBK in tumor and normal tissues in TCGA + GETx dataset. (B) PBK protein level in tumor tissues compared to normal tissues. (C) PBK expression in different stages of ACC, COAD, KICH, KIRC, KIRP, LIHC, and LUAD.
Analysis of the CPTAC dataset showed that PBK protein expression was significantly higher in KIRC, COAD, OV, and UCEC tissues than in normal tissues (Figure 1B). The results of the pathologic staging map based on GEPIA2 data revealed significant associations between PBK expression and pathologic stage in ACC (p = 0.0076), COAD (p = 6.32e−09), KICH (p = 8.04e−06), KIRC (p = 0.000294), KIRP (p = 2.19e−15), LIHC (p = 1.64e−05), and LUAD (p = 0.000189) (Figure 1C). In ACC, KICH, KIRC, KIRP, and LUAD, the association was positive whereas in COAD and LIHC, there was an inverse relationship between pathologic stage and PBK expression.
PBK Mutation Profile in Different Tissues
We used cBioPortal to determine the mutation frequency of PBK in the TCGA database (10,967 samples in 32 studies). We found that PBK had a relatively high mutation frequency in PRAD, OV, and UCEC (>6%) (Figures 2A,B). A total of 40 mutations (including 31 missense, 7 truncation, 1 deletion, and 1 fusion) were detected between amino acid 0 and 322 (Figure 2C).
FIGURE 2. Genetic alterations in PBK. (A) PBK mutation frequency in multiple TCGA pan-cancer dataset in the cBioPortal database. (B) Overall mutation count for CD96 in various TCGA cancer types in the cBioPortal database. (C) Mutation diagram of CD96 across protein domains in different cancer types. (D) OS, DSS, DFS, and PFS in cancers with genetic alterations in PBK.
PBK gene mutations were unrelated to survival and prognosis of cancers in TCGA. Using cBioPortal, we determined that there was no significant difference in OS (p = 0.163), DFS (p = 0.883), disease-specific survival (DSS; p = 0.173), or PFS (p = 0.106) between samples with vs those without PBK mutation (Figure 2D).
Association Between PBK Expression and Cancer Patient Prognosis
To evaluate the prognostic value of PBK expression, we used the PrognoScan database to analyze the relationship between PBK and survival outcomes of cancer patients. The results based on 7 cohorts including GSE13507 (BLCA, HR = 1.60, p = 0.001998) (Kim et al., 2010; Lee et al., 2010), GSE4271-GP96 (GBM, HR = 1.33, p = 0.012008) (Phillips et al., 2006; Costa et al., 2010; Lin et al., 2011), GSE4412-GP96 (GBM, HR = 1.32, p = 0.0416) (Freije et al., 2004), GSE1456-GP96 (BRCA, HR = 2.21, p = 0.00088) (Pawitan et al., 2005; Hall et al., 2006), GSE 13213 (LUAD, HR = 1.43,p = 0.000337) (Tomida et al., 2009), GSE31210 (LUAD, HR = 1.96, p = 0.000087) (Okayama et al., 2012; Yamauchi et al., 2012), and GSE19234 (SKCM, HR = 3.44, p = 0.000211) (Bogunovic et al., 2009) showed that high PBK expression was significantly associated with a worse prognosis (Figures 3A–G).
FIGURE 3. Survival analysis comparing high and low expression of PBK in different types of cancer in the Gene Expression Omnibus and TCGA datasets. (A–G) Kaplan–Meier survival curves in 8 cancer cohorts (GSE13507, GSE4271-GP96, GSE4412-GP96, GSE1456-GP96, GSE13213, GSE31210, and GSE19234). (H, I) Relationship between PBK expression and patient prognosis (OS and DSS) of different cancers in TCGA database. (J) Survival curves of OS in 11 cancer types (ACC, DLBC, KIRC, KIRP, LGG, LIHC, LUAD, MESO, PCPG, OV, and THYM) in TCGA.
In order to further assess the prognostic value of PBK in cancer, we performed univariate Cox regression and Kaplan–Meier survival analyses using TCGA pan-cancer clinical data from the UCSC database. The results of the univariate Cox regression analysis showed that PBK expression level was associated with clinical outcome in ACC (HR = 2.178, p < 0.001), CHOL (HR = 1.941, p = 0.029), KICH (HR = 2.344, p < 0.001), KIRC (HR = 1.907, p < 0.001), KIRP (HR = 3.024, p < 0.001), LGG (HR = 1.515, p < 0.001), LIHC (HR = 1.492, p < 0.001), LUAD (HR = 1.255, p < 0.001), mesothelioma (MESO; HR = 1.898, p < 0.001), pancreatic adenocarcinoma (PAAD; HR = 1.714, p < 0.001), pheochromocytoma and paraganglioma (PCPG; HR = 3.843, p = 0.001), and sarcoma (SARC; HR = 1.245) and was a risk factor for OS in PCPG (p = 0.02), but a protective factor in OV (HR = 0.859, p = 0.031) and THYM (HR = 0.531, p = 0.015) (Figure 3H). To avoid deviations caused by non-cancer events, we performed a univariate Cox regression analysis for DSS. PBK was a risk factor for DSS in patients with ACC (HR = 2.166, p < 0.001), KICH (HR = 2.336, p < 0.001), KIRC (HR = 2.215, p < 0.001), KIRP (HR = 3.633, p < 0.001), LGG (HR = 1.559, p < 0.001), LIHC (HR = 1.361, p = 0.002), LUAD (HR = 1.225, p = 0.005), MESO (HR = 1.954, p < 0.001), PAAD (HR = 1.781, p < 0.001), PCPG (HR = 4.683, p = 0.002), and SARC (HR = 1.263, p = 0.022) but was a protective factor in COAD (HR = 0.746, p = 0.017) (Figure 3I).
The Kaplan–Meier survival analysis revealed that high PBK expression predicted worse OS in ACC, KIRC, KIRP, LGG, LIHC, LUAD, MESO, and PCPG but better OS in DLBC, OV and THYM (Figure 3J).
TFs and miRNAs Regulating PBK
Given the importance of PBK in cancer, we explored TFs and miRNAs that regulate PBK expression. We identified 153 transcriptional regulators of PBK in 4 databases including ChEA, NCODE, JASPAR, and TRANSFAC. After screening for those that were represented at least twice in these databases, 11 TFs were ultimately obtained (Figure 4A). We analyzed the expression of these TFs in different cancers using the GSCALite platform. The 11 TFs were highly expressed in most cancers, with E2F1 showing the most significant upregulation (Figure 4B). The heat map indicated that the 11 TFs in 32 cancers were associated with altered PBK expression (Figure 4C): E2F1 in THYM (cor = 0.916), tripartite motif-containing 28 (TRIM28) in THYM (cor = 0.69), MYC in uveal melanoma (UVM) (cor = 0.577), and TATA-binding protein (TBP) in LIHC (cor = 0.497) were positively correlated with PBK level whereas signal transducer and activator of transcription 5A (STAT5A) and Spi-1 proto-oncogene (SPI1) showed negative correlations. E2F1 regulates the cell cycle, DNA repair, apoptosis, and cell proliferation (Fouad et al., 2020). TRIM28 is involved in the regulation of target gene transcription, response to DNA damage, downregulation of p53 activity, stimulation of epithelial–mesenchymal transition (EMT), autophagy induction, and regulation of reverse transcription-mediated transposition (Czerwińska et al., 2017). MYC promotes tumorigenesis by regulating cell proliferation, cell cycle, and genomic instability (Casey et al., 2017). Thus, PBK may contribute to tumorigenesis through interaction with these factors.
FIGURE 4. TF and miRNA regulatory networks. (A) TF and miRNA network. (B) Expression of TFs in 14 cancers. Red and blue dots indicate increased and decreased expression, respectively, in cancers vs normal tissues. (C, D) Heatmap of TFs (C) and miRNAs (D) altered in 32 cancers. Red means TFs or miRNAs is positively correlated with PBK expression, green means TFs or miRNAs is negatively correlated with PBK expression.
MiRNAs control gene expression at the post-transcriptional level and thereby regulate cell proliferation, apoptosis, and invasion in cancer (Li et al., 2018; Zeng et al., 2021). miRNAs targeting PBK were predicted with TargetScan and the miRDB database; a total of 45 miRNAs were obtained in the intersection of the 2 datasets. Correlations between the miRNAs and PBK were calculated using the ENCORI database (p < 0.05). Of the 7 miRNAs with significant correlations (hsa-miR-6505-5p, hsa-miR-192-3p, hsa-miR-101-5p, hsa-miR-3651, hsa-miR-302d-3p, hsa-miR-5694, and hsa-miR-145-5p) (Figure 4A), hsa-miR-101-5p in UCS (cor = −0.61), hsa-miR-145-5p in testicular germ cell tumor (TGCT) (cor = −0.357), and hsa-miR-5694 in MESO (cor = −0.368) showed the strongest inhibitory effect on PBK expression (Figure 4D). miR-101-5p is known to inhibit the progression of KIRC, BRCA, non-SCLC (NSCLC), and CESC (Chen Q. et al., 2019; Shen et al., 2019; Toda et al., 2020; Yamada et al., 2020), while miR-145-5p has this effect in GBM and PRAD (Ozen et al., 2015; Kurogi et al., 2018). Our results indicate that these miRNAs negatively regulate PBK expression, thereby affecting its function in tumorigenesis.
Changes in the Methylation State of PBK Influence Patient Prognosis
Alterations in DNA methylation state are related to gene expression changes in cancer (Masser et al., 2018). We used the GSCALite platform to analyze the methylation profiles of PBK and related TFs in various TCGA cancers. We first compared the methylation status of tumor and normal tissues in 14 cancer types and found that PBK methylation was decreased in THCA, LIHC, and LUAD (Figure 5A). We also examined the correlation between PBK methylation and expression in 32 cancer types and found that the expression levels of PBK and its TFs were negatively correlated with PBK methylation, with positive correlations observed only in a few cases (Figure 5B). Finally, we examined the association between PBK methylation and OS in 22 cancer types and found that in LGG, THCA, GBM, and PCPG, a higher methylation level was associated with a lower risk of death, whereas a positive association was observed between PBK methylation and mortality in UCEC (Figure 5C).
FIGURE 5. PBK gene methylation in TCGA cancer tissues. (A) Bubble map of the differential methylation of PBK and TFs between tumors and corresponding normal tissues in 14 cancer types. Red and blue dots indicate increased and decreased methylation, respectively, in cancers vs normal tissues. (B) Bubble map of the association between methylation and expression of PBK and TFs across different cancer types. Red dots represent increased methylation and expression levels, and blue dots represent increased methylation and decreased expression level. (C) Bubble map showing the association between OS and methylation of PBK and TFs in different cancer types. Red and blue dots indicate high and low risk, respectively. The size of the dot reflects statistical significance, with a larger size indicating greater significance.
High PBK Expression Is Correlated With Immune Cell Infiltration in Cancers
To determine whether PBK is involved in immune cell invasion into the tumor microenvironment, we analyzed the correlation between PBK expression and the infiltration of six immune cell types in 32 cancers in the TIMR database. PBK expression level was positively correlated with the degree of immune cell infiltration in KIRC, LIHC, THCA, and THYM (Figures 6A–D). The correlations between PBK expression and immune cell infiltration in other cancers are shown in Supplementary Figure 1. Our results suggest that in some cancers such as KIRC and LIHC, PBK overexpression can affect patient prognosis by increasing immune cell infiltration.
FIGURE 6. Correlations between PBK expression and immune cell infiltration in cancers. (A) KIRC. (B) LIHC. (C) THCA. (D) THYM.
Relationship Between Immunotherapy, Immune Checkpoints, and PBK
Immune checkpoint molecules maintain immune system activation within a normal range; abnormal expression of these molecules is closely related to the occurrence and development of some tumors. Immune checkpoint inhibitors (ICIs) are a treatment option for certain malignancies. Here we analyzed the correlation between PBK expression and immune checkpoint molecules in different cancers. In CHOL, HNSC, KICH, KIRC, LIHC, and THCA, PBK expression was positively correlated with the expression of cluster of differentiation 86 (CD86), T cell immunoreceptor with Ig and ITIM domains (TIGIT), indoleamine 2,3-dioxygenase 1 (IDO1), cytotoxic T-lymphocyte–associated protein 4 (CTLA4), inducible T cell costimulator (ICOS), lymphocyte-activating 3 (LAG3), and tumor necrosis factor superfamily member 4 (TNFSF4). On the contrary, in GBM, LUSC, and THCA, PBK expression was negatively correlated with the expression of immune checkpoint molecules such as CD86, V-set immunoregulatory receptor (VSIR), hepatitis A virus cellular receptor 2 (HAVCR2), and neuropilin 1 (NRP1) (Figure 7A). These results indicate that PBK plays an important role in tumor immune regulation and may thus impact patient survival. Although PBK expression was significant in urothelial cancer, there was no significant difference in expression level between immunotherapy responders and non-responders (Supplementary Table 1).
FIGURE 7. Correlations between PBK and immune checkpoints and other variables of interest. (A–E) Correlations between PBK expression and confirmed immune checkpoints (A), genes involved in MMR (B), TMB (C), MSI (D), and DNA methylation (E) in multiple cancers. *p < 0.05, **p < 0.01, ***p < 0.001.
Defects in MMR and subsequent MSI can lead to the accumulation of mutations that increase TMB, which in turn stimulates the host’s antitumor immune response (Zhao et al., 2019). We examined the correlation between PBK expression and MMR signatures and found that PBK level was positively correlated with MutL homolog 1 (MLH1), MSH2, MSH6, PMS1 homolog 2 (PMS2), and epithelial cell adhesion molecule (EPCAM) levels in BLCA, BRCA, CESC, LGG, LUAD, LUSC, PAAD, PCPG, PRAD, READ, and UCEC. It was also negatively correlated with EPCAM and positively correlated with MLH1, MSH2, and MSH6 in THYM (Figure 7B). PBK score was positively correlated with TMB in BLCA, BRCA, CHOL, COAD, HNSC, KICH, KIRC, UCEC, LAML, LGG, LUAD, LUSC, MESO, PAAD, PRAD, SARC, SKCM, and STAD; but negatively correlated with TMB in THYM (p < 0.05; Figure 7C). Finally, PBK expression was positively correlated with MSI in COAD, HNSC, LIHC, MESO, READ, SARC, STAD, and UCEC (p < 0.05; Figure 7D) and with the expression of DNA methyltransferases (DNMT1, DNMT2, DNMT3A, and DNMT3B) in ACC, BLCA, BRCA, CESC, ESCA, KICH, KIRC, KIRP, LGG, LIHC, MESO, PCPG, SKCM, THYM, UCEC, and UVM (p < 0.05; Figure 7E).
Enrichment of PBK-Binding Protein Networks
We explored the functional enrichment of PBK in cancer by mapping the protein–protein interaction network of PBK with the STRING database. Ten genes including baculoviral IAP repeat-containing 5 (BIRC5), cyclin B1 (CCNB1), cell division cycle protein 20 (CDC20), cyclin-dependent kinase 1 (CDK1), DLG-associated protein 5 (DLGAP5), mitotic arrest deficient 2-like 1 (MAD2L1), maternal embryonic leucine zipper kinase (MELK), polo-like kinase 1 (PLK1), DNA topoisomerase II alpha (TOP2A), TTK protein kinase (TTK), and PDZ-binding kinase (PBK) formed an interaction network with PBK (Figure 8A). We analyzed the association between these 10 genes and PBK in 32 TCGA cancers and found that the expression of all of the genes in all of the cancers was positively correlated with PBK expression (Figure 8B). The results of the GO functional enrichment analysis showed that PBK is involved in the regulation of cell cycle, serine/threonine kinase activity, and cell division (Figure 8C). The KEGG pathway analysis results revealed enrichment in “cell cycle,” “p53 signal pathway,” and “cell senescence” (Figure 8D). Thus, PBK has potential functions in cell cycle regulation.
FIGURE 8. PBK network. (A) PBK correlation network. (B) Correlation heatmap of 10 genes interacting with PBK in 32 cancers. (C, D) GO (biological process [BP], cellular component [CC], and molecular function [MF]) (C) and KEGG pathway enrichment (D) analyses of the PBK correlation network.
Drug and Pathway Activity Analysis
We examined the relationship between the expression of PBK and its interaction partners and the activation or inhibition of cancer-related signaling pathways using the GSCALite platform. Consistent with results of the functional enrichment analysis, BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, TTK, and PBK were related to the activation of apoptosis, cell cycle, and DNA damage response; and BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, TTK, and PBK were related to the inhibition of hormone receptors and RAS/MAPK pathways. Additionally, PBK was related to the induction of EMT as well as activation of TSC/mammalian target of rapamycin (mTOR) and inhibition of PI3K/AKT and receptor tyrosine kinase (RTK) signaling pathways (Figures 9A,B).
FIGURE 9. Role of PBK in known cancer-related pathways (GSCALite) and analysis of drugs potentially targeting PBK based on PharmacoDB data. (A, B) Pie chart (A) and heatmap (B) of pathway activity. (C) Heatmap of drugs correlated with 10 genes (BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, and TTK) and PBK.
We searched the PharmacoDB database for drugs related to PBK with a cutoff of p < 0.01. Drugs without targets molecules or pathways were excluded based on standard coefficients >0.1 and <−0.1. PBK was found to be resistant to 31 drugs and sensitive to 20 including those targeting the cell cycle such as paclitaxel, docetaxel, doxorubicin, and irinotecan; and apoptosis regulators such as brivanib, TW37, and elesclomol (Figure 9C).
The mTOR signaling pathway regulates cell metabolism, growth, proliferation, survival, and autophagy in many cancers; mTOR inhibitors suppress inflammation, proliferation, autophagy, and apoptosis and are used in cancer treatment (Chen and Zhou, 2020). PBK was found to be sensitive to rapamycin and temsirolimus, which target the mTOR signaling pathway.
PBK overexpression was associated with resistance to linifanib, lisitinib, neratinib, erlotinib, BRD-K04800985, VX-11e, sildenafil, PD-0325901, erlotinib, gefitinib, AZD6244, VNLG/124, WZ8040, PD318088, BMS-536924, FTI-277, canertinib, PD153035, sorafenib, ML083, WZ4002, and BRD-K75293299.
Discussion
PBK has been described in some human cancers and its overexpression is associated with tumor invasion, proliferation, and metastasis in OV, lung cancer, COAD, and PRAD (Shih et al., 2012; Chen et al., 2015; Ikeda et al., 2016; Zykova et al., 2017). However, the function of PBK in many other tumor types remains unclear. In this study, we provided evidence for PBK as a pan-cancer biomarker based on gene expression and survival data and analyses of immune cell infiltration, TFs, miRNAs, pharmacogenomics, and related intracellular signaling pathways.
The TCGA and GTEx data showed that PBK was highly expressed in ACC, BLCA, BRCA, CESC, CHOL, COAD, ESCA, GBM, HNSC, KICH, KIRC, KIRP, LGG, LIHC, LUAD, LUSC, OV, RAAD, PRAD, SKCM, STAD, THCA, UCEC, and UCS and downregulated in LAML compared to normal tissue. The Kaplan–Meier survival analysis revealed that high expression of PBK was a risk factor in ACC, DLBC, KIRC, KIRP, LGG, LIHC, LUAD, MESO, and PCPG and a protective factor in OV and THYM. Additionally, we found that PBK expression was positively correlated with cancer stage in ACC, KICH, KIRC, KIRP, and LUAD. PBK expression level was shown to be significantly associated with the progression of ACC and NSCLC (Lei et al., 2013; Kar et al., 2019). Thus, PBK is a prognostic biomarker in many cancers, especially in ACC and LUAD as well as in KIRC and KIRP, which has not been previously reported.
The expression of TFs related to PBK showed the highest positive correlations with UVM and THYM, while PBK-related miRNAs had the strongest inhibitory effects in UCS, TGCT, and MESO. These factors are involved in the regulation of the cell cycle, apoptosis, and cell proliferation and are linked to genome instability. For example, MYC is an oncogene in many cancers and plays an important role in regulating metabolic pathways, especially glutamine and glutamine breakdown (Tambay et al., 2021). TFs regulating PBK expression may thus influence its function in tumorigenesis.
PBK methylation was decreased in THCA, LIHC, and LUAD, which was negatively correlated with the expression level of PBK. Additionally, PBK methylation level was correlated with OS in LGG, THCA, GBM, PCPG, and UCEC cancer. The correlation was especially strong for THCA; therefore, while PBK expression cannot be used as a prognostic marker for THCA, we show here for the first time that a change in its methylation status may predict survival outcome.
There have been few studies on the function of PBK in the tumor immune microenvironment. We observed a negative correlation between PBK expression and the infiltration of immune cells including B cells, CD4+ T cells, CD8+ T cells, and macrophages in most tumors, but this association was especially strong in KIRC, LIHC, THCA, and THYM. Based on this observation combined with the survival data, we speculate that PBK affects the prognosis of cancers—particularly KIRC and LIHC—by enhancing immune cell infiltration.
ICIs are a type of immunotherapy with demonstrated efficacy and safety in multiple cancers (He and Xu, 2020). Patients with high TMB or MSI are more likely to experience long-term survival benefits from immunotherapy (Samstein et al., 2019; Zhu et al., 2021). To date there have been no studies on the relationship between TMB and MSI and PBK expression in cancers. Our results showed that PBK expression in various types of tumor was related to TMB and MSI, although additional studies are needed in order to determine the significance of this observation in terms of response to immunotherapy.
Functional enrichment analysis showed that PBK and related genes were functionally enriched in cell cycle regulation. Previous studies have shown that PBK binds to the CDK1/cyclin B1 complex on the mitotic spindle and induces its phosphorylation, leading to cell proliferation and cell cycle progression (Côté et al., 2002; Park et al., 2010). Dysregulation of the cell cycle results in abnormal cell proliferation, inhibition of apoptosis, and increased invasion and metastasis (Zaretsky et al., 2016). Thus, PBK may play an oncogenic role by interfering with cell cycle regulation. This finding can guide the development of targeted drugs for cancer treatment.
BIRC5, CCNB1, CDC20, CDK1, DLGAP5, MAD2L1, MELK, PLK1, TOP2A, TTK, and PBK formed an interaction network that was associated with activation of cell apoptosis and the DNA damage response and inhibition of RAS/MAPK signaling. PBK was also related to the induction of EMT, activation of TSC/mTOR signaling, and inhibition of the PI3K/AKT and RTK pathways. Thus, PBK and PBK-related genes promote tumorigenesis through a variety of mechanisms.
The low solubility and toxicity of the PBK inhibitors HI-TOPK-032 and OTS964TOPK limit their clinical application (Matsuo et al., 2014; Ikeda et al., 2016). As such, new PBK inhibitors for cancer treatment are needed. In this study, we screened 31 drugs that potentially target PBK and identified several candidates including linifanib, which is used in the treatment of LIHC, COAD, STAD, NSCLC, and HNSC (Hsu et al., 2014; Chen et al., 2016; Horinouchi, 2016; Chan et al., 2017; Han et al., 2021); neratinib, which is used in BRCA (Saura et al., 2021); erlotinib and gefitinib, which are used in NSCLC (Yang et al., 2017); and BMS-536924, which is used in GBM and ESCA (Adachi et al., 2014; Zhou, 2015). The anticancer effect of these drugs may be related to PBK inhibition.
Conclusion
We carried out a comprehensive evaluation of PBK to investigate its potential roles in promoting cancer and as a prognostic indicator through bioinformatics analyses. We identified 4 TFs and 3 miRNAs that were related to PBK upregulation and downregulation, respectively, in cancer. We also found that PBK expression is related to cancer immunity: in some cancers, PBK may affect tumor progression by increasing immune cell infiltration into the tumor microenvironment, making it a potential target for cancer immunotherapy. Functional enrichment analysis showed that PBK and its related genes regulate tumor-related signaling pathways. Finally, we found that PBK is associated with sensitivity and resistance to various anticancer drugs and is related to drug-targeted genes in tumor cells. These findings indicate that evaluating PBK expression levels may have diagnostic utility in cancer and can guide treatment decisions, although experimental studies are needed to confirm this possibility.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
YL designed the experiment and collected and analyzed the data. YL and XJ wrote the manuscript. GP provided technical guidance. CS and GP revised the manuscript. All authors have read and approved the final manuscript.
Funding
This study was supported by the National Natural Science Foundation of China Youth Fund Project (no. 81801908) and Hunan Provincial Natural Science Committee Project (no. 2021JJ70074).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.785370/full#supplementary-material
References
Abe, Y., Matsumoto, S., Kito, K., and Ueda, N. (2000). Cloning and Expression of a Novel MAPKK-like Protein Kinase, Lymphokine-Activated Killer T-Cell-Originated Protein Kinase, Specifically Expressed in the Testis and Activated Lymphoid Cells. J. Biol. Chem. 275 (28), 21525–21531. doi:10.1074/jbc.M909629199
Abe, Y., Takeuchi, T., Kagawa-Miki, L., Ueda, N., Shigemoto, K., Yasukawa, M., et al. (2007). A Mitotic Kinase TOPK Enhances Cdk1/Cyclin B1-Dependent Phosphorylation of PRC1 and Promotes Cytokinesis. J. Mol. Biol. 370 (2), 231–245. doi:10.1016/j.jmb.2007.04.067
Adachi, Y., Ohashi, H., Imsumran, A., Yamamoto, H., Matsunaga, Y., Taniguchi, H., et al. (2014). The Effect of IGF-I Receptor Blockade for Human Esophageal Squamous Cell Carcinoma and Adenocarcinoma. Tumor Biol. 35 (2), 973–985. doi:10.1007/s13277-013-1131-2
Aksamitiene, E., Kholodenko, B. N., Kolch, W., Hoek, J. B., and Kiyatkin, A. (2010). PI3K/Akt-sensitive MEK-independent Compensatory Circuit of ERK Activation in ER-Positive PI3K-Mutant T47D Breast Cancer Cells. Cell Signal. 22 (9), 1369–1378. doi:10.1016/j.cellsig.2010.05.006
Bogunovic, D., O'Neill, D. W., Belitskaya-Levy, I., Vacic, V., Yu, Y.-L., Adams, S., et al. (2009). Immune Profile and Mitotic index of Metastatic Melanoma Lesions Enhance Clinical Staging in Predicting Patient Survival. Proc. Natl. Acad. Sci. 106 (48), 20429–20434. doi:10.1073/pnas.0905139106
Casey, S. C., Baylot, V., and Felsher, D. W. (2017). MYC: Master Regulator of Immune Privilege. Trends Immunol. 38 (4), 298–305. doi:10.1016/j.it.2017.01.002
Chan, E., Goff, L. W., Cardin, D. B., Ancell, K., Smith, S. J., Whisenant, J. G., et al. (2017). Phase II Study of the Multikinase Inhibitor of Angiogenesis, Linifanib, in Patients with Metastatic and Refractory Colorectal Cancer Expressing Mutated KRAS. Invest. New Drugs 35 (4), 491–498. doi:10.1007/s10637-017-0458-8
Chen, Y., and Zhou, X. (2020). Research Progress of mTOR Inhibitors. Eur. J. Med. Chem. 208, 112820. doi:10.1016/j.ejmech.2020.112820
Chen, J.-H., Liang, Y.-X., He, H.-C., Chen, J.-Y., Lu, J.-M., Chen, G., et al. (2015). Overexpression of PDZ-Binding Kinase Confers Malignant Phenotype in Prostate Cancer via the Regulation of E2F1. Int. J. Biol. Macromolecules 81, 615–623. doi:10.1016/j.ijbiomac.2015.08.048
Chen, J., Guo, J., Chen, Z., Wang, J., Liu, M., and Pang, X. (2016). Linifanib (ABT-869) Potentiates the Efficacy of Chemotherapeutic Agents through the Suppression of Receptor Tyrosine Kinase-Mediated AKT/mTOR Signaling Pathways in Gastric Cancer. Sci. Rep. 6, 29382. doi:10.1038/srep29382
Chen, F., Chandrashekar, D. S., Varambally, S., and Creighton, C. J. (2019a). Pan-cancer Molecular Subtypes Revealed by Mass-Spectrometry-Based Proteomic Characterization of More Than 500 Human Cancers. Nat. Commun. 10 (1), 5679. doi:10.1038/s41467-019-13528-0
Chen, Q., Liu, D., Hu, Z., Luo, C., and Zheng, S. L. (2019b). miRNA-101-5p Inhibits the Growth and Aggressiveness of NSCLC Cells through Targeting CXCL6. Onco. Targets. Ther. 12, 835–848. doi:10.2147/ott.S184235
Côté, S., Simard, C., and Lemieux, R. (2002). Regulation of Growth-Related Genes by Interleukin-6 in Murine Myeloma Cells. Cytokine 20 (3), 113–120. doi:10.1006/cyto.2002.1988
Costa, B. M., Smith, J. S., Chen, Y., Chen, J., Phillips, H. S., Aldape, K. D., et al. (2010). ReversingHOXA9Oncogene Activation by PI3K Inhibition: Epigenetic Mechanism and Prognostic Significance in Human Glioblastoma. Cancer Res. 70 (2), 453–462. doi:10.1158/0008-5472.Can-09-2189
Czerwińska, P., Mazurek, S., and Wiznerowicz, M. (2017). The Complexity of TRIM28 Contribution to Cancer. J. Biomed. Sci. 24 (1), 63. doi:10.1186/s12929-017-0374-4
Fouad, S., Hauton, D., and D'Angiolella, V. (2020). E2F1: Cause and Consequence of DNA Replication Stress. Front. Mol. Biosci. 7, 599332. doi:10.3389/fmolb.2020.599332
Freije, W. A., Castro-Vargas, F. E., Fang, Z., Horvath, S., Cloughesy, T., Liau, L. M., et al. (2004). Gene Expression Profiling of Gliomas Strongly Predicts Survival. Cancer Res. 64 (18), 6503–6510. doi:10.1158/0008-5472.Can-04-0452
Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signal. 6 (269), pl1. doi:10.1126/scisignal.2004088
Gaudet, S., Branton, D., and Lue, R. A. (2000). Characterization of PDZ-Binding Kinase, a Mitotic Kinase. Proc. Natl. Acad. Sci. 97 (10), 5167–5172. doi:10.1073/pnas.090102397
Hall, P., Ploner, A., Bjöhle, J., Huang, F., Lin, C.-Y., Liu, E. T., et al. (2006). Hormone-replacement Therapy Influences Gene Expression Profiles and Is Associated with Breast-Cancer Prognosis: a Cohort Study. BMC Med. 4, 16. doi:10.1186/1741-7015-4-16
Han, Y., Zhi, W.-H., Xu, F., Zhang, C.-B., Huang, X.-Q., and Luo, J.-F. (2021). Selection of First-Line Systemic Therapies for Advanced Hepatocellular Carcinoma: A Network Meta-Analysis of Randomized Controlled Trials. World J. Gastreonterol. 27 (19), 2415–2433. doi:10.3748/wjg.v27.i19.2415
He, X., and Xu, C. (2020). Immune Checkpoint Signaling and Cancer Immunotherapy. Cell Res. 30 (8), 660–669. doi:10.1038/s41422-020-0343-4
Horinouchi, H. (2016). Anti-vascular Endothelial Growth Factor Therapies at the Crossroads: Linifanib for Non-small Cell Lung Cancer. Transl. Lung Cancer Res. 5 (1), 78–81. doi:10.3978/j.issn.2218-6751.2015.06.08
Hsu, H.-W., de Necochea-Campion, R., Williams, V., Duerksen-Hughes, P. J., Simental, A. A., Ferris, R. L., et al. (2014). Linifanib (ABT-869), Enhances Cytotoxicity with Poly (ADP-Ribose) Polymerase Inhibitor, Veliparib (ABT-888), in Head and Neck Carcinoma Cells. Oral Oncol. 50 (7), 662–669. doi:10.1016/j.oraloncology.2014.03.006
Hu, F., Gartenhaus, R. B., Eichberg, D., Liu, Z., Fang, H.-B., and Rapoport, A. P. (2010). PBK/TOPK Interacts with the DBD Domain of Tumor Suppressor P53 and Modulates Expression of Transcriptional Targets Including P21. Oncogene 29 (40), 5464–5474. doi:10.1038/onc.2010.275
Ikeda, Y., Park, J.-H., Miyamoto, T., Takamatsu, N., Kato, T., Iwasa, A., et al. (2016). T-LAK Cell-Originated Protein Kinase (TOPK) as a Prognostic Factor and a Potential Therapeutic Target in Ovarian Cancer. Clin. Cancer Res. 22 (24), 6110–6117. doi:10.1158/1078-0432.Ccr-16-0207
Ishikawa, C., Senba, M., and Mori, N. (2018). Mitotic Kinase PBK/TOPK as a Therapeutic Target for Adult T-cell L-eukemia/lymphoma. Int. J. Oncol. 53 (2), 801–814. doi:10.3892/ijo.2018.4427
Kar, A., Zhang, Y., Yacob, B. W., Saeed, J., Tompkins, K. D., Bagby, S. M., et al. (2019). Targeting PDZ-Binding Kinase Is Anti-Tumorigenic in Novel Preclinical Models of ACC. Endocr. Relat. Cancer 26 (10), 765–778. doi:10.1530/erc-19-0262
Kim, W.-J., Kim, E.-J., Kim, S.-K., Kim, Y.-J., Ha, Y.-S., Jeong, P., et al. (2010). Predictive Value of Progression-Related Gene Classifier in Primary Non-muscle Invasive Bladder Cancer. Mol. Cancer 9, 3. doi:10.1186/1476-4598-9-3
Kraft, C., Herzog, F., Gieffers, C., Mechtler, K., Hagting, A., Pines, J., et al. (2003). Mitotic Regulation of the Human Anaphase-Promoting Complex by Phosphorylation. EMBO J. 22 (24), 6598–6609. doi:10.1093/emboj/cdg627
Kurogi, R., Nakamizo, A., Suzuki, S. O., Mizoguchi, M., Yoshimoto, K., Amano, T., et al. (2018). Inhibition of Glioblastoma Cell Invasion by Hsa-miR-145-5p and Hsa-miR-31-5p Co-overexpression in Human Mesenchymal Stem Cells. J. Neurosurg. 130 (1), 44–55. doi:10.3171/2017.8.Jns1788
Lee, J.-S., Leem, S.-H., Lee, S.-Y., Kim, S.-C., Park, E.-S., Kim, S.-B., et al. (2010). Expression Signature ofE2F1and its Associated Genes Predict Superficial to Invasive Progression of Bladder Tumors. J. Clin.Oncol. 28 (16), 2660–2667. doi:10.1200/jco.2009.25.0977
Lee, Y.-J., Park, J.-H., and Oh, S.-M. (2020). TOPK Promotes Epithelial-Mesenchymal Transition and Invasion of Breast Cancer Cells through Upregulation of TBX3 in TGF-β1/Smad Signaling. Biochem. Biophysical Res. Commun. 522 (1), 270–277. doi:10.1016/j.bbrc.2019.11.104
Lei, B., Liu, S., Qi, W., Zhao, Y., Li, Y., Lin, N., et al. (2013). PBK/TOPK Expression in Non-Small-cell Lung Cancer: its Correlation and Prognostic Significance with Ki67 and P53 Expression. Histopathology 63 (5), a–n. doi:10.1111/his.12215
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.Can-17-0307
Li, Z., Lin, C., Zhao, L., Zhou, L., Pan, X., Quan, J., et al. (2018). Oncogene miR-187-5p Is Associated with Cellular Proliferation, Migration, Invasion, Apoptosis and an Increased Risk of Recurrence in Bladder Cancer. Biomed. Pharmacother. 105, 461–469. doi:10.1016/j.biopha.2018.05.122
Lin, K., Taylor, J. R., Wu, T. D., Gutierrez, J., Elliott, J. M., Vernes, J.-M., et al. (2011). TMEFF2 Is a PDGF-AA Binding Protein with Methylation-Associated Gene Silencing in Multiple Cancer Types Including Glioma. PLoS One 6 (4), e18608. doi:10.1371/journal.pone.0018608
Liu, C.-J., Hu, F.-F., Xia, M.-X., Han, L., Zhang, Q., and Guo, A.-Y. (2018). GSCALite: a Web Server for Gene Set Cancer Analysis. Bioinformatics 34 (21), 3771–3772. doi:10.1093/bioinformatics/bty411
Masser, D. R., Hadad, N., Porter, H., Stout, M. B., Unnikrishnan, A., Stanford, D. R., et al. (2018). Analysis of DNA Modifications in Aging Research. Geroscience 40 (1), 11–29. doi:10.1007/s11357-018-0005-3
Matsuo, Y., Park, J.-H., Miyamoto, T., Yamamoto, S., Hisada, S., Alachkar, H., et al. (2014). TOPK Inhibitor Induces Complete Tumor Regression in Xenograft Models of Human Cancer through Inhibition of Cytokinesis. Sci. Transl. Med. 6 (259), 259ra145. doi:10.1126/scitranslmed.3010277
Mizuno, H., Kitada, K., Nakai, K., and Sarai, A. (2009). PrognoScan: a New Database for Meta-Analysis of the Prognostic Value of Genes. BMC Med. Genomics 2, 18. doi:10.1186/1755-8794-2-18
Oh, S.-M., Zhu, F., Cho, Y.-Y., Lee, K. W., Kang, B. S., Kim, H.-G., et al. (2007). T-lymphokine-activated Killer Cell-Originated Protein Kinase Functions as a Positive Regulator of c-Jun-NH2-kinase 1 Signaling and H-Ras-Induced Cell Transformation. Cancer Res. 67 (11), 5186–5194. doi:10.1158/0008-5472.Can-06-4506
Ohashi, T., Komatsu, S., Ichikawa, D., Miyamae, M., Okajima, W., Imamura, T., et al. (2016). Overexpression of PBK/TOPK Contributes to Tumor Development and Poor Outcome of Esophageal Squamous Cell Carcinoma. Ar 36 (12), 6457–6466. doi:10.21873/anticanres.11244
Ohashi, T., Komatsu, S., Ichikawa, D., Miyamae, M., Okajima, W., Imamura, T., et al. (2017). Overexpression of PBK/TOPK Relates to Tumour Malignant Potential and Poor Outcome of Gastric Carcinoma. Br. J. Cancer 116 (2), 218–226. doi:10.1038/bjc.2016.394
Okayama, H., Kohno, T., Ishii, Y., Shimada, Y., Shiraishi, K., Iwakawa, R., et al. (2012). Identification of Genes Upregulated in ALK-Positive and EGFR/KRAS/ALK-negative Lung Adenocarcinomas. Cancer Res. 72 (1), 100–111. doi:10.1158/0008-5472.Can-11-1403
Ozen, M., Karatas, O. F., Gulluoglu, S., Bayrak, O. F., Sevli, S., Guzel, E., et al. (2015). Overexpression of miR-145-5p Inhibits Proliferation of Prostate Cancer Cells and Reduces SOX2 Expression. Cancer Invest. 33 (6), 251–258. doi:10.3109/07357907.2015.1025407
Park, J.-H., Nishidate, T., Nakamura, Y., and Katagiri, T. (2010). Critical Roles of T-LAK Cell-Originated Protein Kinase in Cytokinesis. Cancer Sci. 101 (2), 403–411. doi:10.1111/j.1349-7006.2009.01400.x
Park, J.-H., Jeong, Y.-J., Won, H. K., Choi, S.-Y., Park, J.-H., and Oh, S.-M. (2014). Activation of TOPK by Lipopolysaccharide Promotes Induction of Inducible Nitric Oxide Synthase through NF-Κb Activity in Leukemia Cells. Cell Signal. 26 (5), 849–856. doi:10.1016/j.cellsig.2014.01.004
Park, J.-H., Inoue, H., Kato, T., Zewde, M., Miyamoto, T., Matsuo, Y., et al. (2017). TOPK (T-LAK Cell-Originated Protein Kinase) Inhibitor Exhibits Growth Suppressive Effect on Small Cell Lung Cancer. Cancer Sci. 108 (3), 488–496. doi:10.1111/cas.13160
Pawitan, Y., Bjöhle, J., Amler, L., Borg, A.-L., Egyhazi, S., Hall, P., et al. (2005). Gene Expression Profiling Spares Early Breast Cancer Patients from Adjuvant Therapy: Derived and Validated in Two Population-Based Cohorts. Breast Cancer Res. 7 (6), R953–R964. doi:10.1186/bcr1325
Phillips, H. S., Kharbanda, S., Chen, R., Forrest, W. F., Soriano, R. H., Wu, T. D., et al. (2006). Molecular Subclasses of High-Grade Glioma Predict Prognosis, Delineate a Pattern of Disease Progression, and Resemble Stages in Neurogenesis. Cancer Cell 9 (3), 157–173. doi:10.1016/j.ccr.2006.02.019
Pirovano, G., Ashton, T. M., Herbert, K. J., Bryant, R. J., Verrill, C. L., Cerundolo, L., et al. (2017). TOPK Modulates Tumour-specific Radiosensitivity and Correlates with Recurrence after Prostate Radiotherapy. Br. J. Cancer 117 (4), 503–512. doi:10.1038/bjc.2017.197
Roh, E., Lee, M.-H., Zykova, T. A., Zhu, F., Nadas, J., Kim, H.-G., et al. (2018). Targeting PRPK and TOPK for Skin Cancer Prevention and Therapy. Oncogene 37 (42), 5633–5647. doi:10.1038/s41388-018-0350-9
Rouillard, A. D., Gundersen, G. W., Fernandez, N. F., Wang, Z., Monteiro, C. D., McDermott, M. G., et al. (2016). The Harmonizome: a Collection of Processed Datasets Gathered to Serve and Mine Knowledge about Genes and Proteins. Database (Oxford) 2016, baw100. doi:10.1093/database/baw100
Samstein, R. M., Lee, C.-H., Shoushtari, A. N., Hellmann, M. D., Shen, R., Janjigian, Y. Y., et al. (2019). Tumor Mutational Load Predicts Survival after Immunotherapy across Multiple Cancer Types. Nat. Genet. 51 (2), 202–206. doi:10.1038/s41588-018-0312-8
Saura, C., Matito, J., Oliveira, M., Wildiers, H., Brufksy, A. M., Waters, S. H., et al. (2021). Biomarker Analysis of the Phase III NALA Study of Neratinib + Capecitabine versus Lapatinib + Capecitabine in Patients with Previously Treated Metastatic Breast Cancer. Clin. Cancer Res. 27 (21), 5818–5827. doi:10.1158/1078-0432.Ccr-21-1584
Shen, W., Xie, X. Y., Liu, M. R., and Wang, L. L. (2019). MicroRNA-101-5p Inhibits the Growth and Metastasis of Cervical Cancer Cell by Inhibiting CXCL6. Eur. Rev. Med. Pharmacol. Sci. 23 (5), 1957–1968. doi:10.26355/eurrev_201903_17234
Shih, M.-C., Chen, J.-Y., Wu, Y.-C., Jan, Y.-H., Yang, B.-M., Lu, P.-J., et al. (2012). TOPK/PBK Promotes Cell Migration via Modulation of the PI3K/PTEN/AKT Pathway and Is Associated with Poor Prognosis in Lung Cancer. Oncogene 31 (19), 2389–2400. doi:10.1038/onc.2011.419
Smirnov, P., Kofia, V., Maru, A., Freeman, M., Ho, C., El-Hachem, N., et al. (2018). PharmacoDB: an Integrative Database for Mining In Vitro Anticancer Drug Screening Studies. Nucleic Acids Res. 46 (D1), D994–d1002. doi:10.1093/nar/gkx911
Tambay, V., Raymond, V.-A., and Bilodeau, M. (2021). MYC Rules: Leading Glutamine Metabolism toward a Distinct Cancer Cell Phenotype. Cancers 13 (17), 4484. doi:10.3390/cancers13174484
Toda, H., Seki, N., Kurozumi, S., Shinden, Y., Yamada, Y., Nohata, N., et al. (2020). RNA‐sequence‐based microRNA Expression Signature in Breast Cancer: Tumor‐suppressive miR‐101‐5p Regulates Molecular Pathogenesis. Mol. Oncol. 14 (2), 426–446. doi:10.1002/1878-0261.12602
Tomida, S., Takeuchi, T., Shimada, Y., Arima, C., Matsuo, K., Mitsudomi, T., et al. (2009). Relapse-related Molecular Signature in Lung Adenocarcinomas Identifies Patients with Dismal Prognosis. Jco 27 (17), 2793–2799. doi:10.1200/jco.2008.19.7053
Yamada, Y., Nohata, N., Uchida, A., Kato, M., Arai, T., Moriya, S., et al. (2020). Replisome Genes Regulation by Antitumor miR‐101‐5p in clear Cell Renal Cell Carcinoma. Cancer Sci. 111 (4), 1392–1406. doi:10.1111/cas.14327
Yamauchi, M., Yamaguchi, R., Nakata, A., Kohno, T., Nagasaki, M., Shimamura, T., et al. (2012). Epidermal Growth Factor Receptor Tyrosine Kinase Defines Critical Prognostic Genes of Stage I Lung Adenocarcinoma. PLoS One 7 (9), e43923. doi:10.1371/journal.pone.0043923
Yang, Z., Hackshaw, A., Feng, Q., Fu, X., Zhang, Y., Mao, C., et al. (2017). Comparison of Gefitinib, Erlotinib and Afatinib in Non‐small Cell Lung Cancer: A Meta‐analysis. Int. J. Cancer 140 (12), 2805–2819. doi:10.1002/ijc.30691
Yu, W.-N., Lin, H.-F., Lee, Y., Shia, W.-C., Sung, W.-W., Yeh, C.-M., et al. (2021). PBK Expression Is Associated with Prognosis of Patients with Oral Squamous Cell Carcinoma Treated with Radiotherapy: A Retrospective Study. Anticancer Res. 41 (4), 2177–2182. doi:10.21873/anticanres.14991
Zaretsky, J. M., Garcia-Diaz, A., Shin, D. S., Escuin-Ordinas, H., Hugo, W., Hu-Lieskovan, S., et al. (2016). Mutations Associated with Acquired Resistance to PD-1 Blockade in Melanoma. N. Engl. J. Med. 375 (9), 819–829. doi:10.1056/NEJMoa1604958
Zeng, Z., Zhu, Q., Zhao, Z., Zu, X., and Liu, J. (2021). Magic and Mystery of microRNA‐32. J. Cel Mol Med 25, 8588–8601. doi:10.1111/jcmm.16861
Zhao, P., Li, L., Jiang, X., and Li, Q. (2019). Mismatch Repair Deficiency/microsatellite Instability-High as a Predictor for Anti-PD-1/pd-L1 Immunotherapy Efficacy. J. Hematol. Oncol. 12 (1), 54. doi:10.1186/s13045-019-0738-1
Zhao, R., Choi, B. Y., Wei, L., Fredimoses, M., Yin, F., Fu, X., et al. (2020). Acetylshikonin Suppressed Growth of Colorectal Tumour Tissue and Cells by Inhibiting the Intracellular Kinase, T‐lymphokine‐activated Killer Cell‐originated Protein Kinase. Br. J. Pharmacol. 177 (10), 2303–2319. doi:10.1111/bph.14981
Zhou, Q. (2015). BMS-536924, an ATP-Competitive IGF-1R/IR Inhibitor, Decreases Viability and Migration of Temozolomide-Resistant Glioma Cells In Vitro and Suppresses Tumor Growth In Vivo. Onco. Targets. Ther. 8, 689–697. doi:10.2147/ott.S80047
Zhu, M., Jin, Z., and Hubbard, J. M. (2021). Management of Non-Colorectal Digestive Cancers with Microsatellite Instability. Cancers 13 (4), 651. doi:10.3390/cancers13040651
Zykova, T. A., Zhu, F., Vakorina, T. I., Zhang, J., Higgins, L. A., Urusova, D. V., et al. (2010). T-LAK Cell-Originated Protein Kinase (TOPK) Phosphorylation of Prx1 at Ser-32 Prevents UVB-Induced Apoptosis in RPMI7951 Melanoma Cells through the Regulation of Prx1 Peroxidase Activity. J. Biol. Chem. 285 (38), 29138–29146. doi:10.1074/jbc.M110.135905
Zykova, T. A., Zhu, F., Wang, L., Li, H., Bai, R., Lim, D. Y., et al. (2017). The T-LAK Cell-Originated Protein Kinase Signal Pathway Promotes Colorectal Cancer Metastasis. EBioMedicine 18, 73–82. doi:10.1016/j.ebiom.2017.04.003
Glossary
ACC adrenocortical carcinoma
BLCA bladder urothelial carcinoma
BRCA breast invasive carcinoma;
CESC cervical squamous cell carcinoma and endocervical adenocarcinoma
CHOL cholangiocarcinoma
COAD colon adenocarcinoma
DLBC lymphoid neoplasm diffuse large B-cell lymphoma
ESCA esophageal carcinoma
GBM glioblastoma multiforme
HNSC head and neck squamous cell carcinoma
KICH kidney chromophobe
KIRC kidney renal clear cell carcinoma
KIRP kidney renal papillary cell carcinoma
LAML acute myeloid leukemia
LGG brain lower grade glioma
LIHC liver hepatocellular carcinoma
LUAD lung adenocarcinoma
LUSC lung squamous cell carcinoma
MESO mesothelioma
NSCLC non-small cell lung cancer
OV ovarian serous cystadenocarcinoma
PAAD pancreatic adenocarcinoma
PCPG pheochromocytoma and paraganglioma
PRAD prostate adenocarcinoma
READ rectum adenocarcinoma
SARC sarcoma
SKCM skin cutaneous melanoma
STAD stomach adenocarcinoma
TGCT testicular germ cell tumor
THCA thyroid carcinoma
THYM thymoma
UCEC uterine corpus endometrial carcinoma
UCS uterine carcinosarcoma
UVM uveal melanoma.
Keywords: PBK, immune infiltration, methylation, cell cycle, tumorigenesis, pan-cancer
Citation: Liu Y, Xiang J, Peng G and Shen C (2021) Omics- and Pharmacogenomic Evidence for the Prognostic, Regulatory, and Immune-Related Roles of PBK in a Pan-Cancer Cohort. Front. Mol. Biosci. 8:785370. doi: 10.3389/fmolb.2021.785370
Received: 29 September 2021; Accepted: 27 October 2021;
Published: 11 November 2021.
Edited by:
Florence Le Calvez-Kelm, International Agency for Research on Cancer (IARC), FranceReviewed by:
Carmelina Antonella Iannuzzi, Istituto Nazionale Tumori Fondazione G. Pascale (IRCCS), ItalyFrancesca Pentimalli, Istituto Nazionale Tumori Fondazione G. Pascale (IRCCS), Italy
Copyright © 2021 Liu, Xiang, Peng and Shen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gang Peng, penggang83403236@126.com; Chenfu Shen, scf680524@csu.edu.cnscf680524@csu.edu.cn