Next Article in Journal
Chemokine Profile in Psoriasis Patients in Correlation with Disease Severity and Pruritus
Next Article in Special Issue
Complexities of Prostate Cancer
Previous Article in Journal
Extracellular Vesicles in Diffuse Large B Cell Lymphoma: Characterization and Diagnostic Potential
Previous Article in Special Issue
ALDH3A2, ODF2, QSOX2, and MicroRNA-503-5p Expression to Forecast Recurrence in TMPRSS2-ERG-Positive Prostate Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 21
10.3390/ijms232113326
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Association between Cyclin Dependent Kinase 2 Associated Protein 1 (CDK2AP1) and Molecular Subtypes of Lethal Prostate Cancer

by
Yaser Gamallat
1,2,3,
Andrea Bakker
1,
Ealia Khosh Kish
1,2,
Muhammad Choudhry
1,2,
Simon Walker
1,
Saood Aldakheel
1,
Sima Seyedi
2,
Kuo-Cheng Huang
1,
Sunita Ghosh
4,
Geoffrey Gotto
5 and
Tarek A. Bismar
1,2,3,*
1
Department of Pathology and Laboratory Medicine, Cumming School of Medicine, University of Calgary, Calgary, AB T2N 1N4, Canada
2
Departments of Oncology, Biochemistry and Molecular Biology, Cumming School of Medicine, Calgary, AB T2N 4N1, Canada
3
Arnie Charbonneau Cancer Institute and Tom Baker Cancer Center, Calgary, AB T2N 4N1, Canada
4
Departments of Mathematical and Statistical Sciences and Medical Oncology, Faculty of Medicine and Dentistry, University of Alberta, Edmonton, AB T6G 2R7, Canada
5
Prostate Cancer Center, Calgary, AB T2V 1P9, Canada
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(21), 13326; https://doi.org/10.3390/ijms232113326
Submission received: 22 September 2022 / Revised: 25 October 2022 / Accepted: 28 October 2022 / Published: 1 November 2022
(This article belongs to the Special Issue Advances in Agressive Prostate Cancer)
Browse Figures
Review Reports Versions Notes

Abstract

:
Prostate cancer (PCa) is one of the most commonly diagnosed types of malignancy and is the second leading cause of cancer-related death in men in developed countries. Cyclin dependent kinase 2 associate protein 1(CDK2AP1) is an epigenetic and cell cycle regulator gene which has been downregulated in several malignancies, but its involvement in PCa has not yet been investigated in a clinical setting. We assessed the prognostic value of CDK2AP1 expression in a cohort of men diagnosed with PCa (n = 275) treated non-surgically by transurethral resection of the prostate (TURP) and studied the relationship between CDK2AP1 expression to various PCa molecular subtypes (ERG, PTEN, p53 and AR) and evaluated the association with clinical outcome. Further, we used bioinformatic tools to analyze the available TCGA PRAD transcriptomic data to explore the underlying mechanism. Our data confirmed increased expression of CDK2AP1 with higher Gleason Grade Group (GG) and metastatic PCa (p <0.0001). High CDK2AP1 expression was associated with worse overall survival (OS) (HR: 1.62, CI: 1.19–2.21, p = 0.002) and cause-specific survival (CSS) (HR: 2.012, CI 1.29–3.13, p = 0.002) using univariate analysis. When compared to each sub-molecular type. High CDK2AP1/PTEN-loss, abnormal AR or p53 expression showed even worse association to poorer OS and CCS and remained significant when adjusted for GG. Our data indicates that CDK2AP1 directly binds to p53 using the Co-Immunoprecipitation (Co-IP) technique, which was validated using molecular docking tools. This suggests that these two proteins have a significant association through several binding features and correlates with our observed clinical data. In conclusion, our results indicated that the CDK2AP1 overexpression is associate with worse OS and CSS when combined with certain PCa molecular subtypes; interaction between p53 stands out as the most prominent candidate which directly interacts with CDK2AP1.

1. Introduction

Prostate cancer (PCa) is a heterogeneous disease and remains as one of the most common cancers worldwide [1]. As of yet, there are few clinical biomarkers used to estimate the likelihood of lethal PCa and predict disease outcomes. Evaluation of molecular profiles allows a more precise risk classification of PCa [2], and within this classification alteration in ETS-related gene (ERG) and phosphatase tensin homologue (PTEN) are among the most common genomic alterations to occur in PCa [3,4]. ERG rearrangements occur in approximately half of PCa cases and may be associated with poor prognosis in some studies [5]. PTEN deletions, have been associated with a worse prognosis [6]. Androgen receptor (AR) expression plays a pivotal role in castration-resistant prostate cancer (CRPC) [7]. TP53 mutation affecting 50% of metastatic PCa cases are well known to serve as a prognostic biomarker of PCa, especially CRPC [8]. However, there remains conflicting results regarding the implementation of biomarker signatures, especially given variance between different ethnicities and the search for an effective biomarker continues [9,10].
Cyclin-dependent Kinase 2-associated Protein 1 (CDK2AP1) was originally called p12DOC−1 or doc-1 (deleted in oral cancer-1) is a highly conserved, ubiquitously expressed gene on chromosome 12q24, which encodes a 115-amino acid protein [11,12]. CDK2AP1 has been reported to function as an epigenetic and cell cycle regulator [13]. Additionally, CDK2AP1 has been suggested to have a role in TGF-Beta mediated growth suppression and may also be involved in epigenetic gene regulation [14]. Varying levels of CDK2AP1 have been observed in tumors and adjacent non-cancerous normal tissues. Decreased expression has also been found in cancers such as breast [15], oral squamous cell carcinoma [16] and lung cancer [17]. On the contrary, a few studies reported high CDK2AP1 mRNA levels associate with tumorigenesis. Xu et al. found that CDK2AP1 knockdown by RNA interference (RNAi) resulted in reduction in cell proliferation, glioma growth and tumorigenesis. [18]. Furthermore, Gera et al. [15] found 39-fold increased expression of CDK2AP1 in adjacent non-cancerous normal tissues when compared to breast tumors.
In the present study, we are the first to investigate the role and prognostic value of CDK2AP1 in a clinical setting of large cohort of PCa patients treated non-surgically by transurethral resection of the prostate (TURP) including incidental, advanced and castrate-resistant prostate cancer (CRPC). We evaluated the possible relevance to known molecular subtypes in PCa including ERG, PTEN, AR and p53 and performed early in vitro studies as well as expanded bioinformatic analysis into molecular associations between CDK2AP1 and lethal PCa.

2. Materials and Methods

2.1. Study Population and Tissue Microarray Construction

The study cohort consisted of a group of men diagnosed with PCa by TURP (n = 275). Patients within this cohort were either not treated actively or treated with Androgen deprivation therapy (ADT) pre or post TURP. Patients with treatments post TURP were categorized as advanced group and those with treatments prior to TURP sample who had advanced local disease with obstructive symptoms while on ADT, were categorized as CRPC group. The incidental group was categorized as patients with no prior ADT therapy and who had Gleason grade group 1–3. The cohort samples were collected between 1999 and 2014, with an overall median follow-up of 39.28 (1.48–101.22) months. Each sample diagnosis and GG was confirmed by two of the study pathologists. Gleason scoring was assessed according to the 2018 WHO and ISUP grade groups. Table 1 demonstrates patient demographics of the study cohort relative to the proposed subgroups, GG and biomarker groups expression. Clinical follow-up was approved by the University of Calgary, Cumming School of Medicine Ethics Review Board from the Alberta Tumor Registry and included dates of therapy, overall survival (OS), and prostate cancer specific survival (CSS). Tissue samples of the cohort were assembled on two tissue microarrays (TMAs) with an average of two cores per patient using a manual tissue arrayer (Beecher Instruments, Silver Spring, MD, USA).

2.2. Immunohistochemistry

CDK2AP1 protein expression was assessed using immunohistochemistry (IHC) on a Dako Omnis autostainer. Briefly, about 4 µm formalin-fixed paraffin-embedded sections were pretreated with citrate pH 6.0 epitope retrieval buffer. Mouse monoclonal CDK2AP1 antibody (Cat# SC-390283, Santa Cruz, Dallas, TX, USA) (1:200) and Androgen receptor clone AR27 (1:25) (mouse monoclonal AR (Cat# NCL-L-AR-318, Leica Biosystems Inc. Buffalo Grove, IL 60089, USA). The FLEX DAB+ Substrate Chromogen system was used as the post-incubation detection reagent. PTEN and ERG were evaluated using Fluorescence in situ hybridization (FISH), as previously described [19,20]. We assessed p53 expression using a previously validated method [21].

2.3. Pathological Analysis

Histological diagnoses of individual TMA cores were confirmed by the two pathologists. on the initial slides. Gleason scoring was assessed according to the 2018 World Health Organization/International Society of Urological Pathology GGs. In each patient, the two predominant patterns of PCa were sampled and included on the TMAs for analysis. CDK2AP1 IHC expression was assessed using a four-tiered system (0, negative; 1, weak; 2, moderate; and 3, high expression). PTEN and ERG IHC were evaluated as binary values (negative vs. positive) reflective of ERG gene rearrangements and homozygous PTEN deletions as detected by FISH, as previously described [19,20]. We assessed p53 expression using a previously validated method reflective of TP53 sequencing mutations [22]. In brief, score 1: wild type; nuclear staining (strong or weak) with internal control; score 0: absent nuclear staining with positive control; score 2: Overexpressed nuclear staining and score 3: Cytoplasmic staining. In short, score 1 = normal and score 0, 2 and 3 = abnormal (but each represents a different type of alteration to p53). AR IHC was assessed using three-tiered system; 0 negative, 1; weak; 2 moderate and 3 high intensities. The AR groups were then combined into low (0–1) and high (2–3) for simplifying analysis.

2.4. Cell line and Western Blotting

Cell lines Hek293, Hela, RWPE1, PC3, DU-145, LNCAP and VCAP were obtained from ATCC. The cells were grown and maintained according to ATCC recommended protocols. PC3-ERG stable cell lines received from Felix Feng, University of Michigan and were maintained in our lab in DMEM/F12 with 10% FBS [23].

2.5. Protein–Protein Interaction Using Co-Immunoprecipitation and Molecular Docking

For investigating the protein–protein interaction relevance between the PCa molecular subtypes and CDK2AP1, we used Pierce Direct Magnetic IP/CO-IP kit (Cat # 88828, Thermo Scientific, Waltham, MA, USA). Briefly, HEK293 cells were grown in DMEM media supplemented with 10% FBS and 100 U/mL penicillin/streptomycin. Cells were harvested and lysed with IP lyses buffer (Ref#1861603, Thermo Scientific, USA) containing Protease/Phosphatase Inhibitor Cocktail (Cell signaling Cat# 5872S). Pierce NHS-Activated magnetic beads incubated with 5 μg of CDK2AP1 monoclonal mouse antibody (Santa Cruz Cat# sc-390283) or IgG (as a negative control) for 2 h at 4 °C on rotator, followed by incubation with 1500 μg of total protein lysate for 2 h or overnight at 4 °C. Finally, the mixture was washed and eluted. The whole cell lysate (Input) and the immune-precipitant or IgG (negative control) was subjected to Western blotting to confirm the interaction.
To construct protein–protein interaction model, active residues involved in protein structure of p53 (PDB: 2OCJ) [24] and CDK2AP1 (PDB: 2KW6) [25] were utilized. HADDOCK 2.4 developed by the Bonvin lab utilized for protein docking (https://www.bonvinlab.org/software/haddock2.4) accessed on 10 July 2022 [26]. Chain A of p53 and chain B for CDK2AP1 were used to model this interaction. For both molecules, a value of 15% was used as the minimum percentage of relative solvent accessibility (RSA) to consider a residue as accessible. Additionally, 40% was used as the Minimum percentage of RSA to automatically define surface neighbors of active residues as passive. The clustering method used was Fraction of Common Contacts (FCC), with a RMSD cut-off of 0.6 and a minimum cluster size of 4. PyMol software used to model the protein structure and interactions. For the physical contacts, a distance of 4 angstroms were used as the cut-off for classifying an interaction as physical interaction. Furthermore, for modeling the interaction we used the following p53 active residues, amino acids specified in the two transactivation’s domains, TAD (residues 6–30) and TAD2 (residues 35–59), p53 domain (residues 99 to 289) and tetramerization domain (residues 319–358).

2.6. CDK2AP1 mRNA Expression in TCGA PRAD and Associated Gene Set Enrichment Analysis

The gene signature analysis was obtained using TCGA PRAD transcriptomics data [27]. Pan-cancer analysis used to analyze the expression range of CDK2AP1 gene in tumor vs. normal across TCGA, GTex and TARGET databases and includes 22 tissue types simultaneously. This tool based on RNA-seq- rapid analysis servers to furnish comparative information for a selected gene of interest. We further compared the tumor gene array data paired with either normal noncancerous or adjacent normal tissue. Then, the results were blotted and presented with the proportions of tumor samples that show higher expression of the selected gene compared to normal samples at each of the quantile cut-off values (minimum, 1st quartile, median, 3rd quartile, maximum) [27]. Then, we used TCGA PRAD mutations profiling of PTEN, TP53 and AR to assess CDK2AP1 expression and compared CDK2AP1 between mutated and un-mutated samples.
To further elaborate the CDK2AP1 gene potential mechanism in PCa, we used linkedOmics tools (http://www.linkedomics.org (accessed on 10 July 2022)) to analyze the cohort available in TCGA PRAD [28]. Both CDK2AP1 overrepresented (ORA) and gene set enrichment (GSEA) were analyzed based on FDR and categorized as molecular function, cellular components and biological process using WEB-based GEne SeT AnaLysis Toolkit and Explorer [29].

2.7. Statistical Analysis

SPSS used to conduct clinical data statistical analysis (IBM SPSS Statistics for Windows, version 25.0, released 2017; IBM Corp., Armonk, NY, USA). Frequency and proportions were reported for categorical data. Mean and standard deviations were reported for normally distributed continuous data, and median and range were reported for non-normally continuous data. Chi-square test was used to compare two categorical variables and Fisher’s exact test was used where the cell frequencies were <5. Overall survival (OS) was calculated from the time of diagnosis to death; patients who were alive at the end of the study period were censored. Prostate cancer–specific survival (CSS) was defined as death due to PCa; patients who died due to any other reason or who were alive at the end of the study period were censored. OS and CSS were analyzed using the Kaplan–Meier method. Median time and the corresponding 95% confidence intervals (CI) were reported. Log-rank tests were used to compare two or more survival curves. Univariate Cox’s proportional hazard models were used to determine the factors associated with OS and CSS; hazard ratios (HR) and the corresponding 95% CIs were reported. Adjusted Cox’s models for Gleason score variable were reported as well, in addition to the univariate association. A p value of <0.05 was used for statistical significance, and two-sided tests were utilized.

3. Results

3.1. CDK2AP1 Expression in PCa

In our cohort, high CDK2AP1 expression was noted in 141/275 (51.3%) of cases and 79.4% of these cases were either advanced or castrate resistant. CDK2AP1 IHC intensity had mean expression of 0.71 ± 0.93 in incidental, 1.49 ± 1.17 in advanced and, 1.63 ± 1.06 in castrate resistant group (p < 0.0001) (Figure 1A,B).

3.2. High CDK2AP1 Expression Is Associated with Higher Gleason Group (GG) and Metastatic Disease

Separately, high CDK2AP1 expression, PTEN-loss, ERG negativity, and wild type p53 expression were associated with higher Gleason Groups in this cohort (p < 0.0001, p < 0.0001, p < 0.0001, p < 0.001, respectively). Furthermore, a significant association was observed between higher Gleason Groups when combining CDK2AP1 expression with other PCa biomarkers, (high CDK2AP1 expression/PTEN loss, high CDK2AP1/AR expression, high CDK2AP1/ERG positive, and high CDK2AP1/abnormal p53 expression (p < 0.0001, p < 0.0001, p < 0.0001, p < 0.0001, respectively) (Table 1).

3.3. High CDK2AP1 Expression Is Associated with Metastatic Prostate Cancer

Separately in this cohort, high CDK2AP1 and high AR expression were significantly associated with metastatic disease (p value 0.034 and 0.019, respectively) and when combined, there is a higher rate of metastatic disease observed in high CDK2AP1 and high AR expression cases (p value 0.013). However, this relationship was not significant for CDK2AP1 when combined with other PCa subtypes. Notably, the mean intensity of CDK2AP1 was significantly higher for metastatic (1.56 ± 1.07) versus non-metastatic (1.19 ± 1.14) disease (p value 0.033). Supplementary Table S1.

3.4. High CDK2AP1 Expression Is Associated with Poor Overall and Cause Specific Survival in PCa

High CDK2AP1 expression was significantly associated with poor OS (HR: 1.62, CI: 1.19–2.21, p = 0.002) and poor CSS (HR: 2.01, CI: 1.29–3.13, p = 0.002) (Figure 1C,D). Moreover, high CDK2AP1 expression was noted in 24% of PTEN negative cases, 18.2% of ERG positive cases, 3.6% AR low expression cases and in 14.2% of abnormal p53 cases (Table 1).

3.5. High CDK2AP1 Expression in Combination with ERG, PTEN, p53 or AR Signify Poorer Overall and Cause Specific Survival

When high CDK2AP1 expression was combined with other PCa molecular subtypes (Figure 2), PTEN-loss and high CDK2AP1 expression were even more significantly associated with poorer OS (HR: 2.83, CI: 1.89–4.24, p < 0.0001) and CSS (HR: 4.43, CI: 2.28–7.93, p < 0.0001) compared to each marker individually or any other combination of the two specifically in relation to CCS in univariate analysis (Table 2). However, PTEN-loss and high CDK2AP1 expression had a higher CSS HR compared to PTEN-loss and low CDK2AP1 expression cases (HR: 4.18, CI: 2.07–8.42, p < 0.0001). In addition, our data further indicate that combining CDK2AP1 expression with other biomarkers such as AR, ERG or p53 status is associated with added prognostic value compared to individual markers.
Interestingly, abnormal p53 combined with high CDK2AP1 expression was the most significant combination, significantly associated with poorer OS (HR: 3.81, CI: 2.46–5.90, p < 0.0001) and CSS (HR: 6.01, CI: 3.36–10.76, p < 0.0001). More importantly, this relationship was still significant in multivariate analysis, after adjusting for Gleason grade grouping; OS (HR: 2.29, CI: 1.42–3.70, p = 0.001) and CSS (HR: 2.36, CI: 1.28–4.34, p = 0.006). However, other combined biomarkers signatures showed comparable prognostic values, depending on outcome measured (Table 3).

3.6. CDK2AP1 Interacts with p53

Co-Immunoprecipitation (Co-IP) was used to explore possible direct interactions between CDK2AP1 and different PCa molecular subtypes. Initially, we screened several PCa cell line (Hela and Hek293) to check whether they express CDK2AP1 and the other molecular subtypes. We found that Hek293 and Hela cells expressed both CDK2AP1 and p53 (Figure 3A). Co-IP results revealed that p53 potentially immunoprecipitated with CDK2AP1 but PTEN, ERG or AR did not (Figure 3B). Further, we used protein–protein molecular docking tools to validate the potential binding sites between CDK2AP1 and p53.
When analyzing the docking results obtained from the HADDOCK 2.4 webs server, the lowest score was selected as the best model (55.8 ± 18.0, this cluster was also shown to have a z-score of −1.5) which interpreted the highest possible interaction.
CDK2AP1 protein structure is comprised of 115 amino acid residues, which come together to form a structure with two alpha helices. Each monomer of CDK2AP1 undergoes dimerization to form an active protein. This protein only contains one active domain which lies between residue 3 and 114, this is hypothesized as being the active site domain that interacts with a variety of other proteins. The p53 monomer structure utilized is modelled in the absence of DNA consists of 2 alpha helices and 8 beta-strands, the monomers undergo tetramerization to form an active p53 protein. The active residues of this p53 protein include two transactivation’s domains (TADs), the first TAD is between residues 6–30, while the second TAD is between residues 35–59. Moreover, the p53 active site domain is situated between residues 99 and 28. The protein also contains a tetramerization domain which is found between residues 319–358.
Visualizing the modelled interaction (Figure 3C), we were able to detect the presence of hydrogen bonding and physical contact interactions between CDK2AP1 and p53. From the in silico derivation, it was observed that CDK2AP1 exhibited 15 hydrogen bonding contacts with p53 residues. Additionally, the modelled interactions also demonstrated 29 physical contacts within 4 angstroms of interacting p53 residues (Supplementary Tables S3 and S4). In the visualization of the docking, we found hydrogen bonding to be present between the N-terminal histidine residues (H3, H8, H10) with side chains of glutamate, glutamine, serine, and histidine residues found on p53 (Supplementary Tables S3 and S4). The remainder of the hydrogen bonding interactions were between the residues found on the alpha-helixes of CDK2AP1 (between S11–G44) with various residue sidechains spread throughout the p53 molecule. Next, a significant number of physical contacts were found to be present when analyzing this protein interaction model, the exact interactions have been listed in supplementary Table S3. The M1 and K25 residues of CDK2AP1 were shown to have physical contacts with 6 of the p53 residues. Furthermore, two arginine residues (R40 and R43) on the second alpha-helix motif were shown to have significant physical interaction with p53 residues. Additionally, moderate amounts of physical contact were also observed between histidine residues 3, 8, and 10 with each residue interacting with 4–5 other residues of p53.
Moreover, this modelled protein–protein interaction had a Van der Waals energy value that was calculated to be −72.1 ± 14.0, along with electrostatics energy value of −521.8 ± 6.8. In addition to this, the modelled interaction also had a desolvation and restraints violation energy value of −15.4 ± 0.8 and 2476.1 ± 34.4, respectively.

3.7. CDK2AP1 mRNA Expression in Pan-Cancer Data

Pan-Cancer data analysis for CDK2AP1 mRNA expression in 22 tumors revealed CDK2AP1 was significantly upregulated in all 22 tumors including PCa (Figure 4A). Interestingly, the significant difference was more apparent when the tumor was compared to noncancerous normal samples vs. adjacent normal tissue in prostate samples (Figure 4B,C). PRAD TCGA mutation profiling for the PTEN, TP53 and AR, we found CDK2AP1 is highly expressed in PTEN and TP53 but not AR mutations (Figure 4D).

3.8. CDK2AP1 Gene Set Enrichment Analysis in PCa TCGA PRAD

CDK2AP1 Gene Set Enrichment analysis data revealed a significant and distinctive gene set profiling as presented by volcano blots (Figure 5A). Heat map (Figure 5B) showing the top 50 positively and top 50 negatively correlated genes (Figure 5C). CDK2AP1 overexpression GSEA the genes set enrichment analysis (Figure 5D) indicated upregulating biological processes including cell proliferation and growth, upregulation in metabolic process, biological regulations, cellular components organization, localization, development process, cell communication. The Cellular component upregulated genes appear mostly in the nucleus (1686 genes) and then followed by the cytoplasmic membrane (1593 genes). The positive correlated genes for molecular functions show a wide range of potential functions from protein and nucleic acid binding to chromatin binding. When we analyzed the overrepresented (ORA) genes (Figure 5E), many cellular functions such as mitochondria respiratory chain complex assembly, mitochondria gene expression, mRNA processing, ribonucleoprotein complex biogenesis, protein localization to endoplasmic reticulum, DNA-template transcription, and termination, tRNA metabolic process, protein-DNA complex subunit organization, RNA localization and protein localization to mitochondria were overrepresented. Interestingly, the negative correlated gene sets were mostly involved declined immune and inflammatory responses which were downregulated. All together, these data indicated a potential function of CDK2AP1 in tumorigenesis and immune-evading mechanism.

4. Discussion

CDK2AP1 has diverse functions ranging from epigenetics moderator to cell cycle regulator [13]. Earlier studies suggest that CKD2AP1 may function as a tumor suppressor, however, we saw a unique CDK2AP1 overexpression pattern in prostatic adenocarcinoma along with most other malignancies. Our data indicates poor survival associated with increasing expression of CDK2AP1 in patients with PCa. Interestingly, this association was variable when implicated with other molecular subtypes of PCa (i.e.,: when combined with PTEN, ERG, AR and p53). Previous studies have shown the association between PTEN-loss and high Gleason scores [30], however the role of CDK2AP1 with PTEN and other various PCa subtypes remains a mystery. In this cohort, high CDK2AP1 expression with PTEN-loss, high AR, positive ERG and abnormal p53 expression were significantly associated with higher Gleason Groups and poor overall survival. Furthermore, there was a significant increase in the mean intensity of CDK2AP1 expression in metastatic disease when compared to non-metastatic disease. Overall, the higher expression of CDK2AP1 in high Gleason Group and metastatic tumors demonstrate its role as a potential biomarker for lethal PCa. Our data suggested the role of CDK2AP1 as a possible oncogene in this malignancy.
Even though PCa is still trailing other cancers in biomarker guided treatment, genes such as PTEN, ERG, AR and TP53 have been noted as clinically relevant biomarkers to predict patient’s outcomes and suggest therapeutic strategies [4,31]. p53 inactivation or abnormal expression has been demonstrated to be associated with adverse patient clinical outcomes and metastatic disease while high expression of ERG and PTEN-loss have been demonstrated to be clinically relevant to shorter patient survival, higher Gleason scores and unfavorable pathological and clinical outcomes [32]. We found an interesting association between CDK2AP1 and p53 expression in our clinical data which were confirmed by Co-IP and protein–protein interaction docking tools. Alsayegh et al. demonstrated CDK2AP1 knockdown downregulated p53 in primary human fibroblasts [33]. Although p53 is classically thought of as a tumor suppressor, recent studies have demonstrated that some types of mutant TP53 display gain-of-function with increased expression of malignant properties [34].
We saw a significant overexpression of CDK2AP1 in PCa tumors through mRNA-seq data when compared to normal tissue, however, this significance was not persistent when compared to adjacent normal tissues (Figure 4C). More interestingly, this significance was observed with PTEN and TP53 mutated samples. In breast cancer, Gera et al. [15] reports that they found 39-fold increase of normalized mRNA expression of CDK2AP1 in noncancerous adjacent tissues [15]. Moreover, Xu et al. [18] studied CDK2AP1 overexpression association with tumorigenesis in glioma in vitro and in vivo and found that, CDK2AP1 knockdown inhibits cells proliferation, and colony formation. Furthermore, CDK2AP1 knockdown induced cell cycle arrest and apoptosis. Furthermore, an increased CDK2AP1 expression was reported in precursor lesions of invasive tumors in the oral epithelium [35]. This data supports with our findings in PCa.
Many cellular functions were found to be significantly over and under-represented. Mostly, pathways directly related to DNA-templated transcription, termination, mRNA processing, RNA localization, and mitochondrial respiratory chain complex assembly and mitochondrial gene expression were over-represented. Such pathways correspond to the oncogenic role of CDK2AP1 in cellular proliferation as mitochondrial dysfunction and dysregulation is very common in PCa behavior [36]. Furthermore, pathways directly related to the immune response against tumor cells such as leukocyte activation involved in inflammatory response, regulation of leukocyte activation, interleukin-6 production were under-represented in CDK2AP1 gene set enrichment analysis. Previous studies have shown that leukocyte-to-monocyte ratio is higher in patients with prostate cancer which correspond to the under-represented regulation of leukocyte activation [37]. Furthermore, it has been described that interleukin-6 is down regulated in human prostatic carcinoma cells which provides further evidence for this enrichment gene analysis and role of CDK2AP1 as an oncogene. It appears that CDK2AP1 oncogenic overexpression may enhance the ability of tumor cells to evade the immune system [38].
Some limitations of this study, relate to the cohort investigated being a heterogenous population ranging from incidental PCa to advanced and CRPC in addition to samples being TURP in nature. Nonetheless, this represents a real sample of patients at a tertiary institution in Western country. Further in vitro studies are needed to elucidate more of CDK2AP1 function in tumorigenesis.

5. Conclusions

In conclusion, our study confirms a significant role of CDK2AP1 in PCa. CDK2AP1 upregulation was shown to be associated with poor OS and CSS in PCa. Furthermore, when combined with other PCa molecular subtypes, p53 and PTEN were the most prevalent mutations associated with CDK2AP1.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232113326/s1, Table S1: CDK2AP1 expression in metastatic and non-metastatic PCa. Table S2: Number of patients at risk. Table S3: Hydrogen Bonding Contacts, Table S4: Physical contacts within 4 angstroms.

Author Contributions

Conceptualization, T.A.B.; Data curation, S.G., A.B., G.G. and T.A.B.; Formal analysis, S.G., Y.G., S.W., K.-C.H. and E.K.K.; Funding acquisition, T.A.B.; Methodology, Y.G., M.C., E.K.K., S.S., and S.A.; Supervision, T.A.B.; Writing—original draft, A.B., Y.G. and E.K.K.; Writing—review and editing, E.K.K., M.C., G.G., S.W., K.-C.H. and T.A.B. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported in part by the Prostate Cancer Foundation Young Investigator Award (T.A.B). This work is also supported by Prostate Cancer Canada and The Canadian Cancer Society and is proudly funded by the Canadian Cancer Society #400391.

Institutional Review Board Statement

The study was approved by the University of Calgary Cumming School of Medicine Ethics Review Board and in accordance with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Sandra Zhuang and Jennifer Pollonais for technical support in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CDK2AP1Cyclin Dependent Kinase 2 Associated Protein 1
PCaProstate cancer
TURPTransurethral resection of the prostate
CSSCause-specific survival
OSOverall survival
GGGleason Grade Group
Co-IPCo-Immunoprecipitation
ERGETS-related gene
PTENPhosphatase tensin homologue
ARAndrogen receptor
CRPCCastration-resistant prostate cancer
ADTAndrogen deprivation therapy
TMAsTissue microarrays
IHCImmunohistochemistry
FISHFluorescence in situ hybridization
FCCFraction of Common Contacts
RSARelative solvent accessibility
KMKaplan–Meier
CIConfidence intervals
HRHazard ratios
TADTransactivation’s domains

References

  1. Siegel, R.L.; Miller, K.D.; Fuchs, H.E.; Jemal, A. Cancer Statistics, 2021. CA Cancer J. Clin. 2021, 71, 7–33. [Google Scholar] [CrossRef]
  2. Arora, K.; Barbieri, C.E. Molecular Subtypes of Prostate Cancer. Curr. Oncol. Rep. 2018, 20, 58. [Google Scholar] [CrossRef]
  3. Imada, E.L.; Sanchez, D.F.; Dinalankara, W.; Vidotto, T.; Ebot, E.M.; Tyekucheva, S.; Franco, G.R.; Mucci, L.A.; Loda, M.; Schaeffer, E.M.; et al. Transcriptional Landscape of PTEN Loss in Primary Prostate Cancer. BMC Cancer 2021, 21, 856. [Google Scholar] [CrossRef]
  4. Bismar, T.A.; Hegazy, S.; Feng, Z.; Yu, D.; Donnelly, B.; Palanisamy, N.; Trock, B.J. Clinical Utility of Assessing PTEN and ERG Protein Expression in Prostate Cancer Patients: A Proposed Method for Risk Stratification. J. Cancer Res. Clin. Oncol. 2018, 144, 2117–2125. [Google Scholar] [CrossRef]
  5. Khosh Kish, E.; Choudhry, M.; Gamallat, Y.; Buharideen, S.M.; Dhananjaya, D.; Bismar, T.A. The Expression of Proto-Oncogene ETS-Related Gene (ERG) Plays a Central Role in the Oncogenic Mechanism Involved in the Development and Progression of Prostate Cancer. Int. J. Mol. Sci. 2022, 23, 4772. [Google Scholar] [CrossRef]
  6. Abou-Ouf, H.; Assem, H.; Ghosh, S.; Karnes, R.J.; Stoletov, K.; Palanisamy, N.; Lewis, J.D.; Bismar, T.A. High Serine-arginine Protein Kinase 1 Expression with PTEN Loss Defines Aggressive Phenotype of Prostate Cancer Associated with Lethal Outcome and Decreased Overall Survival. Eur. Urol. Open Sci. 2021, 23, 1–8. [Google Scholar] [CrossRef]
  7. Wasmuth, E.V.; Hoover, E.A.; Antar, A.; Klinge, S.; Chen, Y.; Sawyers, C.L. Modulation of Androgen Receptor Dna Binding Activity through Direct Interaction with the ETS Transcription Factor ERG. Proc. Natl. Acad. Sci. USA 2020, 117, 8584–8592. [Google Scholar] [CrossRef] [Green Version]
  8. Liu, Z.; Guo, H.; Zhu, Y.; Xia, Y.; Cui, J.; Shi, K.; Fan, Y.; Shi, B.; Chen, S. TP53 Alterations of Hormone-Naïve Prostate Cancer in the Chinese Population. Prostate Cancer Prostatic Dis. 2021, 24, 482–491. [Google Scholar] [CrossRef]
  9. Stern, N.; Ly, T.L.; Welk, B.; Chin, J.; Ballucci, D.; Haan, M.; Power, N. Association of Race and Ethnicity with Prostate Cancer–Specific Mortality in Canada. JAMA Netw. Open 2021, 4, e2136364. [Google Scholar] [CrossRef]
  10. Lotan, T.L.; Tomlins, S.A.; Bismar, T.A.; Van Der Kwast, T.H.; Grignon, D.; Egevad, L.; Kristiansen, G.; Pritchard, C.C.; Rubin, M.; Bubendorf, L. Report from the International Society of Urological Pathology (ISUP) Consultation Conference on Molecular Pathology of Urogenital Cancers. I. Molecular Biomarkers in Prostate Cancer. Am. J. Surg. Pathol. 2020, 44, e15–e29. [Google Scholar] [CrossRef]
  11. Todd, R.; McBride, J.; Tsujl, T.; Donoff, R.B.; Nagai, M.; Chou, M.Y.; Chiang, T.; Wong, D.T.W. Deleted in Oral Cancer-1 (doc-l), a Novel Oral Tumor Suppressor Gene. FASEB J. 1995, 9, 1362–1370. [Google Scholar] [CrossRef] [PubMed]
  12. Daigo, Y.; Suzuki, K.; Maruyama, O.; Miyoshi, Y.; Yasuda, T.; Kabuto, T.; Imaoka, S.; Fujiwara, T.; Takahashi, E.; Fujino, M.A.; et al. Isolation, Mapping and Mutation Analysis of a Human cDNA Homologous to the doc-1 Gene of the Chinese Hamster, a Candidate Tumor Suppressor for Oral Cancer. Genes Chromosom. Cancer 1997, 20, 204–207. [Google Scholar] [CrossRef]
  13. Wong, D.T.; Kim, J.J.; Khalid, O.; Sun, H.-H.; Kim, Y. Double Edge: CDK2AP1 in Cell-Cycle Regulation and Epigenetic Regulation. J. Dent. Res. 2012, 91, 235–241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Spruijt, C.G.; Bartels, S.J.J.; Brinkman, A.B.; Tjeertes, J.V.; Poser, I.; Stunnenberg, H.G.; Vermeulen, M. CDK2AP1/DOC-1 is a Bona Fide Subunit of the Mi-2/NuRD Complex. Mol. Biosyst. 2010, 6, 1700–1706. [Google Scholar] [CrossRef] [PubMed]
  15. Gera, R.; Mokbel, L.; Jiang, W.G.; Mokbel, K. mRNA Expression of CDK2AP1 in Human Breast Cancer: Correlation with Clinical and Pathological Parameters. Cancer Genom. Proteom. 2018, 15, 447–452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Shintani, S.; Ohyama, H.; Zhang, X.; McBride, J.; Matsuo, K.; Tsuji, T.; Hu, M.G.; Hu, G.; Kohno, Y.; Lerman, M.; et al. p12DOC-1 Is a Novel Cyclin-Dependent Kinase 2-Associated Protein. Mol. Cell. Biol. 2000, 20, 6300–6307. [Google Scholar] [CrossRef]
  17. Sun, M.; Jiang, R.; Wang, G.; Zhang, C.; Li, J.; Jin, C.; Zhang, X. Cyclin-Dependent Kinase 2-Associated Protein 1 Suppresses Growth and Tumorigenesis of Lung Cancer. Int. J. Oncol. 2013, 42, 1376–1382. [Google Scholar] [CrossRef] [Green Version]
  18. Xu, Y.; Wang, J.; Fu, S.; Wang, Z. Knockdown of CDK2AP1 by RNA Interference Inhibits Cell Growth and Tumorigenesis of Human Glioma. Neurol. Res. 2014, 36, 659–665. [Google Scholar] [CrossRef]
  19. Huang, K.-C.; Bégin, L.R.; Palanisamy, N.; Donnelly, B.; Bismar, T.A. SPINK1 Expression in Relation to PTEN and ERG in Matched Primary and Lymph Node Metastatic Prostate Cancer: Implications for Biomarker Development. Urol. Oncol. Semin. Orig. Investig. 2016, 34, 235.e1–235.e10. [Google Scholar] [CrossRef]
  20. Bismar, T.A.; Yoshimoto, M.; Duan, Q.; Liu, S.; Sircar, K.; Squire, J.A. Interactions and Relationships of PTEN, ERG, SPINK1 and AR in Castration-Resistant Prostate Cancer. Histopathology 2012, 60, 645–652. [Google Scholar] [CrossRef]
  21. Yemelyanova, A.; Vang, R.; Kshirsagar, M.P.; Lu, D.; Marks, M.A.; Shih, I.M.; Kurman, R.J. Immunohistochemical Staining Patterns of p53 Can Serve as a Surrogate Marker for TP53 Mutations in Ovarian Carcinoma: An Immunohistochemical and Nucleotide Sequencing Analysis. Mod. Pathol. 2011, 24, 1248–1253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Guedes, L.B.; Almutairi, F.; Haffner, M.C.; Rajoria, G.; Liu, Z.; Klimek, S.; Zoino, R.; Yousefi, K.; Sharma, R.; De Marzo, A.M.; et al. Analytic, Preanalytic, and Clinical Validation of p53 IHC for Detection of TP53 Missense Mutation in Prostate Cancer. Clin. Cancer Res. 2017, 23, 4693–4703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Han, S.; Brenner, J.C.; Sabolch, A.; Jackson, W.; Speers, C.; Wilder-Romans, K.; Knudsen, K.E.; Lawrence, T.S.; Chinnaiyan, A.M.; Feng, F.Y. Targeted Radiosensitization of ETS Fusion-Positive Prostate Cancer through PARP1 Inhibition. Neoplasia 2013, 15, 1207–1217. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, Y.; Rosengarth, A.; Luecke, H. Structure of the Human p53 Core Domain in the Absence of DNA. Acta Crystallogr. D Biol. Crystallogr. 2007, 63 Pt. 3, 276–281. [Google Scholar] [CrossRef]
  25. Ertekin, A.; Aramini, J.M.; Rossi, P.; Leonard, P.G.; Janjua, H.; Xiao, R.; Maglaqui, M.; Lee, H.-W.; Prestegard, J.H.; Montelione, G.T. Human Cyclin-dependent Kinase 2-associated Protein 1 (CDK2AP1) Is Dimeric in Its Disulfide-reduced State, with Natively Disordered N-terminal Region. J. Biol. Chem. 2012, 287, 16541–16549. [Google Scholar] [CrossRef] [Green Version]
  26. Honorato, R.V.; Koukos, P.I.; Jiménez-García, B.; Tsaregorodtsev, A.; Verlato, M.; Giachetti, A.; Rosato, A.; Bonvin, A.M.J.J. Structural Biology in the Clouds: The WeNMR-EOSC Ecosystem. Front. Mol. Biosci. 2021, 8, 729513. [Google Scholar] [CrossRef]
  27. Bartha, Á.; Győrffy, B. TNMplot.com: A Web Tool for the Comparison of Gene Expression in Normal, Tumor and Metastatic Tissues. Int. J. Mol. Sci. 2021, 22, 2622. [Google Scholar] [CrossRef]
  28. Vasaikar, S.V.; Straub, P.; Wang, J.; Zhang, B. LinkedOmics: Analyzing Multi-Omics Data within and across 32 Cancer Types. Nucleic Acids Res. 2018, 46, D956–D963. [Google Scholar] [CrossRef] [Green Version]
  29. Liao, Y.; Wang, J.; Jaehnig, E.J.; Shi, Z.; Zhang, B. WebGestalt 2019: Gene Set Analysis Toolkit with Revamped UIs and APIs. Nucleic Acids Res. 2019, 47, W199–W205. [Google Scholar] [CrossRef] [Green Version]
  30. Jamaspishvili, T.; Berman, D.M.; Ross, A.E.; Scher, H.I.; De Marzo, A.M.; Squire, J.A.; Lotan, T.L. Clinical Implications of PTEN Loss in Prostate Cancer. Nat. Rev. Urol. 2018, 15, 222–234. [Google Scholar] [CrossRef]
  31. Garcia-Marques, F.; Liu, S.; Totten, S.M.; Bermudez, A.; Bs, C.T.; Hsu, E.; Ba, R.N.; Hembree, A.; Stoyanova, T.; Brooks, J.D.; et al. Protein Signatures to Distinguish Aggressive from Indolent Prostate Cancer. Prostate 2022, 82, 605–616. [Google Scholar] [CrossRef] [PubMed]
  32. Gesztes, W.; Schafer, C.; Young, D.; Fox, J.; Jiang, J.; Chen, Y.; Kuo, H.-C.; Mwamukonda, K.B.; Dobi, A.; Burke, A.P.; et al. Focal p53 Protein Expression and Lymphovascular Invasion in Primary Prostate Tumors Predict Metastatic Progression. Sci. Rep. 2022, 12, 5404. [Google Scholar] [CrossRef] [PubMed]
  33. Alsayegh, K.N.; Gadepalli, V.S.; Iyer, S.; Rao, R.R. Knockdown of CDK2AP1 in Primary Human Fibroblasts Induces p53 Dependent Senescence. PLoS ONE 2015, 10, e0120782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Hafner, A.; Bulyk, M.L.; Jambhekar, A.; Lahav, G. The Multiple Mechanisms that Regulate p53 Activity and Cell Fate. Nat. Rev. Mol. Cell Biol. 2019, 20, 199–210. [Google Scholar] [CrossRef]
  35. Katarkar, A.; Patel, L.; Mukherjee, S.; Ray, J.G.; Haldar, P.K.; Chaudhuri, K. Association of Oral Tumor Suppressor Gene Deleted in Oral Cancer-1 (DOC-1) in Progression of Oral Precancer to Cancer. Oral Sci. Int. 2015, 12, 15–21. [Google Scholar] [CrossRef] [Green Version]
  36. Xiao, J.; Cohen, P.; Stern, M.C.; Odedina, F.; Carpten, J.; Reams, R. Mitochondrial Biology and Prostate Cancer Ethnic Disparity. Carcinogenesis 2018, 39, 1311–1319. [Google Scholar] [CrossRef]
  37. Cihan, Y.B.; Arslan, A.; Ergul, M.A. Subtypes of White Blood Cells in Patients with Prostate Cancer or Benign Prostatic Hyperplasia and Healthy Individuals. Asian Pac. J. Cancer Prev. 2013, 14, 4779–4783. [Google Scholar] [CrossRef] [Green Version]
  38. Tsui, K.-H.; Feng, T.-H.; Hsieh, W.-C.; Chang, P.-L.; Juang, H.-H. Expression of interleukin-6 is Downregulated by 17-(allylamino)-17-demethoxygeldanamycin in Human Prostatic Carcinoma Cells. Acta Pharmacol. Sin. 2008, 29, 1334–1341. [Google Scholar] [CrossRef]
Ijms 23 13326 g001 550
Figure 1. The expression of CDK2AP1 in Prostate tissues and associated survival analysis (A) IHC staining of CDK2AP1 in benign and PCa tissues (20× magnification). (B) Box plots represents the mean expression of CDK2AP1 in Incidental, Advanced, and CRPC in our study cohort samples. The error bars indicate the standard deviation of the mean (SD). Student t-test was performed, p value < 0.05 was considered significant between Advanced or CRPC compared with Incidental (Benign). (C) Kaplan–Meier (KM) curves outlining the expression of CDK2AP1 in relation to overall survival (OS) and (D) Cause-specific survival (CSS).
Figure 1. The expression of CDK2AP1 in Prostate tissues and associated survival analysis (A) IHC staining of CDK2AP1 in benign and PCa tissues (20× magnification). (B) Box plots represents the mean expression of CDK2AP1 in Incidental, Advanced, and CRPC in our study cohort samples. The error bars indicate the standard deviation of the mean (SD). Student t-test was performed, p value < 0.05 was considered significant between Advanced or CRPC compared with Incidental (Benign). (C) Kaplan–Meier (KM) curves outlining the expression of CDK2AP1 in relation to overall survival (OS) and (D) Cause-specific survival (CSS).
Ijms 23 13326 g001
Ijms 23 13326 g002 550
Figure 2. KM curves demonstrating the relationship between CDK2AP1 expression and OS and CSS (A,B) CDK2AP1 and PTEN expression. (C,D) CDK2AP1 and p53 expression. (E,F) CDK2AP1 and ERG expression. Samples were scores semi-quantitatively using a four-tiered system for CDK2AP1 (negative—0; weak—1; moderate—2, strong—3; or in short low risk = 0, 1 and high risk = 2, 3) and p53 (wildtype—1; absent nuclear staining—0; overexpressed nuclear staining—2; cytoplasmic staining—3; or in short, low risk = 1 and normal, high risk = 0, 2, 3 and abnormal). PTEN low risk = PTEN positive staining; PTEN high risk = PTEN loss and negative staining. ERG low risk = ERG negative staining; ERG high risk = ERG positive staining. Number of patients at risk given in Supplementary Table S2.
Figure 2. KM curves demonstrating the relationship between CDK2AP1 expression and OS and CSS (A,B) CDK2AP1 and PTEN expression. (C,D) CDK2AP1 and p53 expression. (E,F) CDK2AP1 and ERG expression. Samples were scores semi-quantitatively using a four-tiered system for CDK2AP1 (negative—0; weak—1; moderate—2, strong—3; or in short low risk = 0, 1 and high risk = 2, 3) and p53 (wildtype—1; absent nuclear staining—0; overexpressed nuclear staining—2; cytoplasmic staining—3; or in short, low risk = 1 and normal, high risk = 0, 2, 3 and abnormal). PTEN low risk = PTEN positive staining; PTEN high risk = PTEN loss and negative staining. ERG low risk = ERG negative staining; ERG high risk = ERG positive staining. Number of patients at risk given in Supplementary Table S2.
Ijms 23 13326 g002
Ijms 23 13326 g003 550
Figure 3. CDK2AP1 potential binding target for p53 (A) Western blot screening of CDK2AP1 and p53 protein expression. GAPDH used as loading control. (B) The CO-Immunoprecipitation analysis of CDK2AP1 with PCa molecular subtypes (p53, AR, PTEN and ERG). (C) Molecular docking module shows CDK2AP1(RED) and p53 (Blue) Protein–protein interaction.
Figure 3. CDK2AP1 potential binding target for p53 (A) Western blot screening of CDK2AP1 and p53 protein expression. GAPDH used as loading control. (B) The CO-Immunoprecipitation analysis of CDK2AP1 with PCa molecular subtypes (p53, AR, PTEN and ERG). (C) Molecular docking module shows CDK2AP1(RED) and p53 (Blue) Protein–protein interaction.
Ijms 23 13326 g003
Ijms 23 13326 g004 550
Figure 4. Differential expression of CDK2AP1 in tumors and normal tissues (A) Boxplots shows the CDK2AP1 genes expression in most of common tumor types Normal (left) Tumor(right). Significant differences by a Mann–Whitney U test are marked with red color (* p < 0.01). (B) Boxplot shows CDK2AP1 RNA-seq expression in PCa tumor vs. noncancerous normal tissue. (C) and PCa tumor vs. Paired adjacent-normal tissue and (D) CDK2AP1 expression in PTEN/TP53 or AR mutations (green is wildtype and red is mutated).
Figure 4. Differential expression of CDK2AP1 in tumors and normal tissues (A) Boxplots shows the CDK2AP1 genes expression in most of common tumor types Normal (left) Tumor(right). Significant differences by a Mann–Whitney U test are marked with red color (* p < 0.01). (B) Boxplot shows CDK2AP1 RNA-seq expression in PCa tumor vs. noncancerous normal tissue. (C) and PCa tumor vs. Paired adjacent-normal tissue and (D) CDK2AP1 expression in PTEN/TP53 or AR mutations (green is wildtype and red is mutated).
Ijms 23 13326 g004
Ijms 23 13326 g005 550
Figure 5. CDK2AP1 gene set enrichment association in PCa tumorigenesis (A) Volcano blots showed the CDK2AP1 associated genes in TCGA PRAD. Significant up-regulation genes (dark red dots, FDR < 0.01) whereas genes showed significant down-regulation (dark green dots, FDR < 0.01). (B,C) Heatmap represents the top 50 positively or negatively correlated genes associated with CDK2AP1 overexpression. (D) Bar blots represents GSEA categories; Biological Process, Cellular Component and Molecular Function (red, blue and green bars, respectively). The height of the bar represents the number of IDs in the user list and in the category. The FDR is calculated using Benjamini-Hochberg method. (E) Bars blot showed the positively correlated GO (Molecular Function) for Overrepresentation Enrichment Analysis (ORA). Data were analyzed using Pearson correlation coefficient and ranked according to FDR value.
Figure 5. CDK2AP1 gene set enrichment association in PCa tumorigenesis (A) Volcano blots showed the CDK2AP1 associated genes in TCGA PRAD. Significant up-regulation genes (dark red dots, FDR < 0.01) whereas genes showed significant down-regulation (dark green dots, FDR < 0.01). (B,C) Heatmap represents the top 50 positively or negatively correlated genes associated with CDK2AP1 overexpression. (D) Bar blots represents GSEA categories; Biological Process, Cellular Component and Molecular Function (red, blue and green bars, respectively). The height of the bar represents the number of IDs in the user list and in the category. The FDR is calculated using Benjamini-Hochberg method. (E) Bars blot showed the positively correlated GO (Molecular Function) for Overrepresentation Enrichment Analysis (ORA). Data were analyzed using Pearson correlation coefficient and ranked according to FDR value.
Ijms 23 13326 g005
Table
Table 1. Demographics of TURB cohort: clinical and pathological biomarkers.
Table 1. Demographics of TURB cohort: clinical and pathological biomarkers.
VariablesTotal (n = 275) (%)
Gleason Score
Grade group 1 82 (29.8)
Grade group 2 (3 + 4)29 (10.5)
Grade group 3 (4 + 3)17 (6.2)
Grade group 4 (8)20 (7.3)
Grade group 5 (9/10)121 (44.0)
Missing6 (2.2)
Deceased
Yes165 (60.0)
No108 (39.3)
Missing2 (0.7)
Prostate cancer specific mortality
Yes85 (30.9)
No188 (68.4)
Missing2 (0.7)
Cancer Subgroup
Incidental94 (34.2)
Advanced109 (39.6)
Castrate Resistant72 (26.2)
CDK2AP1 Score (Score by Cancer subgroup)
Score 2 or 3141 (51.3%)
Incidental Score29 (10.5)
Advanced Score66 (24.0)
Castrate Resistant Score46 (16.7)
Score 0 or 1134 (48.7)
Incidental Score65 (23.6)
Advanced Score43 (15.6)
Castrate Resistant Score26 (9.5)
PTEN and CDK2AP1 combined
PTEN negative and CDK2AP1 Score 2, 366 (24.0)
PTEN negative and CDK2AP1 Score 0, 128 (10.2)
PTEN positive and CDK2AP1 score 2, 373 (26.5)
PTEN positive and CDK2AP1 Score 0, 1102 (37.1)
Missing6 (2.2)
ERG and CDK2AP1 Combined
ERG positive and CDK2AP1 Score 2, 350 (18.2)
ERG positive and CDK2AP1 Score 0, 120 (7.3)
ERG negative and CDK2AP1 Score 2, 389 (32.4)
ERG negative and CDK2AP1 Score 0, 1110 (40.0)
Missing6 (2.2)
AR and CDK2AP1 Combined
AR Score 0, 1 and CDK2AP1 Score 2, 310 (3.6)
AR Score 0, 1 and CDK2AP1 Score 0, 112 (4.4)
AR Score 2, 3 and CDK2AP1 Score 2, 3126 (45.8)
AR Score 2, 3 and CDK2AP1 Score 0, 1115 (41.8)
Missing12 (4.4)
p53 and CDK2AP1 Combined
p53 Score 0, 2, 3 and CDK2AP1 Score 2, 339 (14.2)
p53 Score 0, 2, 3 and CDK2AP1 Score 0, 110 (3.6)
p53 Score 1 and CDK2AP1 Score 2, 398 (35.6)
p53 Score 1 and CDK2AP1 Score 0, 199 (36.0)
Missing29 (10.5)
CDK2AP1, AR scores as 0; negative, 1; weak, 2; moderate and 3, high intensity. p53 scores 0; absent nuclear staining, 1 wild type, 2; overexpressed nuclear staining, 3; cytoplasmic staining; PTEN scores 0; negative, 1; positive.
Table
Table 2. Univariate analysis of CDK2AP1 expression in relation to OS and CSS.
Table 2. Univariate analysis of CDK2AP1 expression in relation to OS and CSS.
VariablesOverall SurvivalProstate Cancer Specific Survival
HR (95% CI)p-ValueHR (95% CI)p-Value
CDK2AP1
(Score 0, 1) (n = 129)
Score 2, 3 (n = 135)1.62 (1.19–2.21)0.0022.01 (1.29–3.13)0.002
PTEN
(Positive) (n = 170)
PTEN negative (n = 89)2.44 (1.79–3.33)<0.00013.26 (2.12–5.01)<0.0001
CDK2AP1 and PTEN
(PTEN positive and CDK Score 0, 1) (n = 98)
PTEN negative and CDK score 2, 3 (n = 62)2.83 (1.89–4.24)<0.00014.43 (2.28–7.93)<0.0001
PTEN negative and CDK score 0, 1 (n = 28)3.07 (1.87–5.01)<0.00014.18 (2.07–8.42)<0.0001
PTEN positive and CDK score 2, 3 (n = 72)1.64 (1.08–2.51)0.0211.97 (1.05–3.71)<0.0001
CDK2AP1 and ERG
(ERG negative and CDK score 0, 1) (n = 106)
ERG positive and CDK score 2, 3 (n = 48)2.00 (1.31–3.06)0.0012.42 (1.33–4.40)0.004
ERG positive and CDK score 0, 1 (n = 20)2.62 (1.51–4.53)0.0012.63 (1.18–5.88)0.018
ERG negative and CDK score 2, 3 (n = 86)1.85 (1.27–2.69)0.0012.28 (1.34–3.87)0.002
CDK2AP1 and AR
(AR score 2, 3 and CDK score 0, 1) (n = 110)
AR score 0, 1 and CDK score 2, 3 (n = 10)3.53 (1.73–7.17)0.0015.41 (2.19–13.37)<0.0001
AR score 0, 1 and CDK score 0, 1 (n = 12)2.09 (1.06–4.12)0.0332.81 (1.14–6.92)0.025
AR score 2, 3 and CDK score 2, 3 (n = 120)1.74 (1.24–2.44)0.0012.20 (1.34–3.62)0.002
CDK2AP1 and p53
(p53 score 1 and CDK score 0, 1) (n = 96)
p53 score 0, 2, 3 and CDK score 2, 3 (n = 37)3.81 (2.46–5.90)<0.00016.01 (3.36–10.76)<0.0001
p53 score 0, 2, 3 and CDK score 0, 1 (n = 9)2.37 (1.08–5.24)0.0331.49 (0.35–6.35)0.589
p53 score 1 and CDK score 2, 3 (n = 94)1.34 (0.92–1.97)0.1301.44 (0.82–2.54)0.203
CDK2AP1, AR scores as 0; negative, 1; weak, 2; moderate and 3; high intensity. p53 scores 0; absent nuclear staining, 1; wild type, 2; overexpressed nuclear staining, 3; cytoplasmic staining.
Table
Table 3. Multivariate analysis of CDK2AP1 expression in relation to OS and CSS adjusted for Gleason Score.
Table 3. Multivariate analysis of CDK2AP1 expression in relation to OS and CSS adjusted for Gleason Score.
VariablesOverall SurvivalProstate Cancer Specific Mortality
HR (95% CI)p-ValueHR (95% CI)p-Value
CDK2AP1 and PTEN
(PTEN positive and CDK2AP1 Score 0, 1) (n = 91)
PTEN negative and CDK2AP1 score 2, 3 (n = 60)1.74 (1.12–2.71)0.0141.84 (0.99–3.42)0.054
PTEN negative and CDK2AP1 score 0, 1 (n = 26)2.36 (1.41–3.93)0.0012.39 (1.14–5.00)0.021
PTEN positive and CDK2AP1 score 2, 3 (n = 68)1.44 (0.93 –2.22)0.1021.54 (0.80–2.94)0.193
ERG and CDK2AP1
(ERG negative and CDK2AP1 score 0, 1) (n = 98)
ERG positive and CDK2AP1 score 2, 3 (n = 46)1.22 (0.77–1.93)0.3901.06 (0.56–2.01)0.861
ERG positive and CDK2AP1 score 0, 1 (n = 19)2.30 (1.32–4.01)0.0032.08 (0.92–4.71)0.078
ERG negative and CDK2AP1 score 2, 3 (n = 82)1.58 (1.07–2.33)0.0211.79 (1.04–3.10)0.037
AR and CDK2AP1
(AR score 2, 3 and CDK2AP1 score 0, 1) (n = 104)
AR score 0, 1 and CDK2AP1 score 2, 3 (n = 10)2.07 (1.00–4.29)0.0502.17 (0.87–5.42)0.098
AR score 0, 1 and CDK2AP1 score 0, 1 (n = 11)1.40 (0.68–2.87)0.3611.27 (0.48–3.40)0.629
AR score 2, 3 and CDK2AP1 score 2, 3 (n = 115)1.32 (0.92–1.89)0.1281.32 (0.79–2.22)0.294
P53 and CDK2AP1
(p53 score 1 and CDK2AP1 score 0, 1) (n = 91)
p53 score 0, 2, 3 and CDK2AP1 score 2, 3 (n = 37)2.29 (1.42–3.70)0.0012.36 (1.28–4.34)0.006
p53 score 0, 2, 3 and CDK2AP1 score 0, 1 (n = 9)1.51 (0.67–3.40)0.3250.63 (0.15–2.72)0.631
p53 score 1 and CDK2AP1 score 2, 3 (n = 88)1.07 (0.71–1.60)0.7570.89 (0.49–1.61)0.697
CDK2AP1, AR scores as 0; negative, 1; weak, 2; moderate and 3; high intensity. p53 scores 0; absent nuclear staining, 1; wild type, 2; overexpressed nuclear staining, 3; cytoplasmic staining.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gamallat, Y.; Bakker, A.; Khosh Kish, E.; Choudhry, M.; Walker, S.; Aldakheel, S.; Seyedi, S.; Huang, K.-C.; Ghosh, S.; Gotto, G.; et al. The Association between Cyclin Dependent Kinase 2 Associated Protein 1 (CDK2AP1) and Molecular Subtypes of Lethal Prostate Cancer. Int. J. Mol. Sci. 2022, 23, 13326. https://doi.org/10.3390/ijms232113326

AMA Style

Gamallat Y, Bakker A, Khosh Kish E, Choudhry M, Walker S, Aldakheel S, Seyedi S, Huang K-C, Ghosh S, Gotto G, et al. The Association between Cyclin Dependent Kinase 2 Associated Protein 1 (CDK2AP1) and Molecular Subtypes of Lethal Prostate Cancer. International Journal of Molecular Sciences. 2022; 23(21):13326. https://doi.org/10.3390/ijms232113326

Chicago/Turabian Style

Gamallat, Yaser, Andrea Bakker, Ealia Khosh Kish, Muhammad Choudhry, Simon Walker, Saood Aldakheel, Sima Seyedi, Kuo-Cheng Huang, Sunita Ghosh, Geoffrey Gotto, and et al. 2022. "The Association between Cyclin Dependent Kinase 2 Associated Protein 1 (CDK2AP1) and Molecular Subtypes of Lethal Prostate Cancer" International Journal of Molecular Sciences 23, no. 21: 13326. https://doi.org/10.3390/ijms232113326

APA Style

Gamallat, Y., Bakker, A., Khosh Kish, E., Choudhry, M., Walker, S., Aldakheel, S., Seyedi, S., Huang, K. -C., Ghosh, S., Gotto, G., & Bismar, T. A. (2022). The Association between Cyclin Dependent Kinase 2 Associated Protein 1 (CDK2AP1) and Molecular Subtypes of Lethal Prostate Cancer. International Journal of Molecular Sciences, 23(21), 13326. https://doi.org/10.3390/ijms232113326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop