Skip to main content

Prognostic and clinicopathological insights of phosphodiesterase 9A gene as novel biomarker in human colorectal cancer



PDE9A (Phosphodiesterase 9A) plays an important role in proliferation of cells, their differentiation and apoptosis via intracellular cGMP (cyclic guanosine monophosphate) signaling. The expression pattern of PDE9A is associated with diverse tumors and carcinomas. Therefore, PDE9A could be a prospective candidate as a therapeutic target in different types of carcinoma. The study presented here was designed to carry out the prognostic value as a biomarker of PDE9A in Colorectal cancer (CRC). The present study integrated several cancer databases with in-silico techniques to evaluate the cancer prognosis of CRC.


The analyses suggested that the expression of PDE9A was significantly down-regulated in CRC tissues than in normal tissues. Moreover, methylation in the DNA promoter region might also manipulate PDE9A gene expression. The Kaplan–Meier curves indicated that high level of expression of PDE9A gene was associated to higher survival in OS, RFS, and DSS in CRC patients. PDE9A demonstrated the highest positive correlation for rectal cancer recurrence with a marker gene CEACAM7. Furtheremore, PDE9A shared consolidated pathways with MAPK14 to induce survival autophagy in CRC cells and showed interaction with GUCY1A2 to drive CRPC.


Overall, the prognostic value of PDE9A gene could be used as a potential tumor biomarker for CRC.

Peer Review reports


Cancer is the second most common cause of death worldwide which occurs due to the uncontrolled growth of cells that can migrate to other parts of the body [1, 2]. According to WHO-Cancer Report-2020-Global Profile, around 18 million new cancer cases and 9.6 million deaths are estimated worldwide from cancer in 2018. Therein, Colon cancer (CRC) is imputed as the third most diagnosed cancer in the world [3]. Therefore, CRC is described as a malignant tumor that arises from the inner surface of the large intestine [4]. Over 50% of patients die from colon cancer despite having intensive investigations and therapeutic improvements [5]. However, the main reason behind the death can be pointed as a result of spreading the uncontrollable cell division tendency to other organs spontaneously [6]. CRC is known as westernized diseases and it shows comparatively lower lethality in the South Asian population than the western countries [3, 7]. WHO reported in 2018 that, Colon cancer deaths in Bangladesh are 3598, or 0.46% of total deaths ( According to Syed Akram Hussain and Richard Sullivan, colon cancer is one of the 10 common causes of deaths from cancer in both males and females [8]. But the affected case of CRC is expected to increase by 60% to more than 2.2 million new cases and 1.1 million deaths by 2030 [9]. Besides, epidemiological studies reveal that CRC is related to familial and hereditary factors, besides some other factors like age, environment, and lifestyle also affect [10,11,12]. Experiments on CRC patients show that DNA methylation and covalent modification of histones results in carcinogenesis [13]. So, early identification of CRC is extremely important for the prevention of CRC [14]. Perrectum examination is the simplest way to CRC recognition [15]. However, the most efficient and commonly applied method of CRC recognition is endoscopy [15]. Molecular diagnosis methods for CRC, based on genetic and epigenetic tests, also have limited applicability [16, 17]. Stool DNA (sDNA) is very stable which can indicate the changes in DNA in colorectal adenocarcinomas (COAD) [15]. All of these instigate to find the unknown mechanisms contributing to CRC malignancy and the necessity for the development of potential biomarkers for the CRC prognosis.

PDEs (Cyclic nucleotide phosphodiesterases) consisted of a large superfamily of 11 PDE gene families (PDE1 to PDE11) and involves in regulating the intracellular levels of the second messengers cAMP and cGMP [18]. PDE9A is a member of eleven PDE isoforms located at q22.3 position of the 21st chromosome in the human cell [19]. PDE9A encodes high-affinity cGMP-Specific 3′,5′-Cyclic Phosphodiesterase 9A protein which hydrolyzes the second messenger cGMP, that regulates many important physiological processes [20,21,22]. The regulation and expression of PDEs (cyclic nucleotide phosphodiesterases) play a significant role in tumor progression and inhibition [23]. Many hematological malignancies and carcinomas have been connected with decreased levels of cAMP and/or cGMP by the PDE gene family [24]. PDE2A, PDE8B, and PDE11A expression are involved with adrenocortical tumors [25], high expression of PDE4 has been detected in brain tumors [26], in glioblastoma multiforme PDE5 is strongly expressed [27], and increased PDE7B expression found in Chronic lymphocytic leukemia [28]. The cDNA of PDE9A was coined in 1998 and noted as the ninth member of the PDEs family [29]. Only from the PDE9A gene, twenty-eight splice variants can arise ( These can create more proteins through alternative splicing of mRNA or multiple promoters and transcription start sites [18]. Although PDE9 mRNA has been detected in many organs; but, the highest expression has been noted in prostate, colon, small intestine, brain, kidney, spleen, thymus, and hematopoietic cells [22]. The protein expression of PDE9A is highly conserved and also associated with prostate cancer [30]. DNA hypermethylation of the PDE9A gene results in reduced mRNA expression in tumorous tissues [31]. Thus, our investigation was designed to determine the prognostic and clinicopathological importance of the PDE9A gene in CRC since there is no obvious report regarding its surveillance.

The most difficult job for an oncologist is to evaluate the prognostic significance of a gene which is also an essential skill for treatment information [32]. Cancer bioinformatics is believed to be a way of early evaluation of a cancer which is used to identify the role of a gene for any cancer before going laboratories. Various database provides gene expression and related information for analysis in various cancer types like UALCAN provides information on pan-cancer molecular subtypes [33]. Using these onco-informatics processes it was found that DNA topoisomerases expression in lung cancer was higher than the normal tissues [34]. Ping Yan et al., recommended that COL1A1, MMP2, FN1, TIMP1, SPARC, COL4A1, and ITGA5 may be potential biomarkers and therapeutic targets for GC by adopting in silico process [35].

In this study, our target is to assess and identify the role of phosphodiesterase 9A (PDE9A) gene which can portray a complete concept to understand the expression and clinical significance in CRC. Therefore, we investigate the patterns of expression, promoter methylation patterns of our target gene. We also explore the combined prognostic relevance, the association of genes co-expressed with PDE9A, and the interaction network of PDE9A in CRC. Moreover, this multi-omics data mining approach can provide useful intimation to instigate the researcher for finding new approaches in anti-cancer therapies.

Materials and methods

Expression analysis of PDE9A

To retrieve the mRNA expression of the PDE9A gene in various cancer types, we have employed Oncomine database ( that investigates differential expression patterns of a gene in diverse clinical cancer specimens and analogous normal controls, incorporates 715 independent datasets and 86,733 samples [36, 37]. The threshold parameters were settled as reflects: keywords, PED9A gene, primary filter, cancer vs. normal analysis; cancer type, Colorectal Cancer, data type: mRNA, fold change: 2, P-value: 1E− 4, gene ranking: 10%. Co-express genes with PDE9A were also investigated for further analysis. Thereafter, we used UALCAN database ( to analyze expression, survival, methylation, pan-cancer view, and correlation data [38, 39]. PDE9A gene expression and methylation were analyzed based on sample types, individual cancer stages, patient’s race, patient’s weight, patient’s gender, patient’s age, histological subtype, nodal metastasis status, and TP53 mutation status. The positive correlation between PDE9A and CEACAM7 genes for Colon adenocarcinoma (COAD) was assessed and Pearson’s correlation analysis value was 0.54. The p-value was considered statistically significant when it was less than 0.05 (p < 0.05) for all the results. The screening parameters were set: “Gene: PDE9A”; “Analysis Type: TCGA Gene analysis”; “Cancer Type: Colon adenocarcinoma”; “Data Type: TCGA dataset”. OncoLnc ( online tools evaluate survival correlations with expression data and retrieving clinical data for mRNAs, miRNAs, and lncRNAs (long non-coding) [40]. OncoLnc database gives the results of the Cox analysis based on 21 cancer studies for 8647 patients from that link to TCGA survival data or MiTranscriptome beta lncRNA data [41]. This webtool analysis Cox regression results for PDE9A gene and provides information on cox coefficient, Rank, Median expression, P-value, Mean expression, FDR corrected, and Kaplan-Meier plot against Colon adenocarcinoma (COAD). The value 25 was inputted for Lower percentile and Upper percentile to generate Kaplan Meier plot. PrognoScan Database ( is a cancer microarray dataset that provides information for analyzing potential therapeutic value that could accelerate in silico cancer research by clinical annotation [42].

Clinicopathological insights of PDE9A

Gene Expression database of Normal and Tumor tissues 2 (GENT2) ( was utilized to retrieve the information of PDE9A in tissue wide expression pattern, meta-survival analysis, and prognostic significance of a gene of interest-based on tumor subtypes [43, 44]. Accordingly, we put the PDE9A gene for tissue wide gene expression patterns across 72 paired tissue. Subtype profiling along with durability analysis in colon tissue for the PDE9A gene was also explored where results were shown in the box and dot-plot of individual subtype. Moreover, progression-free survival in colon tissue for the PDE9A gene was evaluated by meta-analysis of analytical tests (t-test, log2 fold changes, etc.). UCSC Xena (, a high-performance data mining platform to visualize and analyze functional genomic data sets for both enormous public archives (TCGA, GDC, etc) and private datasets between genomic and phenotypic variables [45] was employed to set up a correlation heat map between PDE9A and CEACAM7 gene expression in the same patient cohort through data processing in TCGA (The Cancer Genome Atlas Program) colon cancer (COAD). Subsequently, PDE9A DNA methylation status is also examined by adopting this UCSC Xena portal [46]. The matched 551 TCGA Colon Cancer specimens were wielded for all analyses. The clinical features of PDE9A was explored by the utilization of cBioPortal for Cancer Genomics ( which is an open-source, open-access web-based data repositories that explore, visualize, and analyze multidimensional cancer genomics data from TCGA [47, 48]. Two TCGA datasets of COAD, particularly “TCGA Nature 2012 (276 samples)” and “TCGA PanCancer Atlas (594 samples) were culled for further investigation of PDE9A gene alterations or copy number alterations (CNA) using the cBioPortal database. The OncoPrint, cancer types summary, plots, mutations, survival, copy number segment, pathways, expression tabs were guided following the default settings of the cBioPortal.

Determination of prognostic effects

Gene Expression Profiling Interactive Analysis (GEPIA) ( was used to explore the differential expression analysis and also to determine the prognostic effects on COAD cancer patients [49]. GEPIA facilitates user’s data mining in diverse areas including differential expression analysis in cancer vs. normal samples, profiling plotting according to cancer categories or various pathological phases, patient survival analysis, correlation analysis, identical gene disclosure, and dimensionality reduction analysis [50].

Co-expression and interacting network analysis

GeneCards, the human gene database ( is a searchable, speedy, and sophisticated user-friendly integrative database that furnishes information of all annotated and predicted human genes [51, 52]. For the, GeneCards was scrutinized to retrieve the information of PDE9A gene regarding the gene functions, pathways, and interactions, mRNA expression in normal human tissues, SNPs, disorders associated with this gene, etc. Later, the R2 platform ( was utilized the R2 web interface to explore the PDE9A gene expression, correlation analysis between PDE9A and CEACAM7 genes as well as to generate a Kaplan-Meier survival plot. R2 platform generates Kaplan-Meier plot (overall survival) for a specific dataset and the Kaplan-Meier plot was generated for PDE9A gene against Tumor Colon Adenocarcinoma - TCGA - 286 - rsem - tcgars dataset with the optimum cut-off values. The co-expression between PDE9A and CEACAM7 genes was identified in Tumor Colon Adenocarcinoma - TCGA - 286 - rsem - tcgars dataset. This dataset provides data analyzing 286 samples and source: TCGA ID: COAD with date: 2000-01-01.

Functional insights of PDE9A

GeneMANIA ( web address predicts the function of a gene or gene lists and identifies the physical interaction, genetic interactions, co-expression, pathway, co-localization, and shared protein domain [53]. We used GeneMANIA database to obtain an interaction network of related connected genes for PDE9A. Genes can be connected by the interacted network based on different attributes. Here nodes represent genes whereas links show networks. The screening status was set: Organism: Homo sapiens; Query Gene: PDE9A; Maximum Resultant Gene: 25; Maximum Resultant Attributes: 10; Query-Dependent Weighting: Automatically and rest are used as default. Furtheremore, STRING database ( was also employed to predict the protein-protein interactions including direct (physical) and indirect (functional) associations [54, 55]. The amino acid sequence used in STRING database was retrieved from KEGG database (593 aa) and Homo sapiens were selected as initial parameters.


The workflow of the study is shown in Fig. 1.

Fig. 1

Graphical demonstration of PDE9A gene expression and their relations with patients survivability in colorectal cancer tissues

PDE9A mRNA expression level evaluation in various cancer

For analyzing the expression of PDE9A in various cancers, we compared its transcriptional levels in various cancer tissues with their corresponding normal tissues by adopting various databases like Oncomine, GEPIA, Gent2, UALCAN, GeneCard. Using the differential analysis tool of the Oncomine database, we analyzed cDNA microarray data for PDE9A in 20 different human cancer types. Significant mRNA upregulation (red) or downregulation (blue) was depicted by the Oncomine graphical results. According to Oncomine, PDE9A was discovered to be downregulated in various colorectal cancer including Colorectal Adenoma, Rectal Mucinous Adenocarcinoma, Cecum Adenocarcinoma, Colon Adenoma, Rectal Adenocarcinoma, Colorectal Carcinoma, Colon Mucinous Adenocarcinoma, Colon Adenocarcinoma, Colon Carcinoma (Fig. 2, Table 1). The expression of PDE9A was also estimated by the GENT2 database where results are depicted in boxplot across 72 paired cancer vs. normal tissues. GENT2 database reveals that PDE9A expression is downregulated in the bladder, blood, breast, cervix, colon, skin, esophagus, head and neck, oral, stomach, and tongue cancers, and is upregulated in various cancer including adrenal gland, bone, brain, ovary, prostate, teeth, thyroid, uterus, and vulva cancer through the GLP570 platform (Fig. 3). UALCAN website also provides expression information of PDE9A across TCGA cancers in which bladder, breast, colon, kidney, rectal, and stomach cancers expression are lower and cholangiocarcinoma, glioblastoma, lung, prostate, thyroid cancers expression are higher following their corresponding normal samples (Supplementary Fig. 1A and 1B). To confirm the result attained from Oncomine, GENT2, UALCAN database, we performed PDE9A expression analysis using another platform, GEPIA. It provides results in the dot plot and bar plot. GEPIA reveals the expression profile across all tumor tissues and corresponding normal tissues (Supplementary Fig. 1C). mRNA expression of PDE9A in COAD was significantly lower in COAD compared to normal tissues. GeneCards database shows the location of the PDE9A gene in 21q22.3 cytogenetic band along with mRNA expression in normal human tissues based on RNAseq, Microarray, and SAGE (Serial Analysis of Gene Expression) (Supplementary Fig. 1D). The immunohistochemical staining images of PDE9A protein expression in CRC are given in (Fig. 4). Systematic mRNA expression analysis carried out here from multiple databases was enough to prove that PDE9A expression across a wide range of cancer types in which it is significantly downregulated in CRC.

Fig. 2

a Expression of PDE9A in different cancers from Oncomine database with significant mRNA depending on different types of cancer and corresponding normal tissue. Here, upregulation is marked in red and downregulation is blue. The threshold parameters were set as follows: P-value of 1E− 4, fold change 2, gene ranking 10%. b-m box plot comparing PDE9A expression in different colorectal cancer types (normal tissue in left and cancerous tissue in right) depicted significant downregulation of PDE9A in b Colorectal Adenoma (COA), c Rectal Mucinous Adenocarcinoma (REMAD), d Cecum Adenocarcinoma (CEAD), e Rectal Adenocarcinoma (READ), f Colon Mucinous Adenocarcinoma (COMAD), g Colon Adenocarcinoma (COAD), h Colon Adenoma (COA), i Colon Carcinoma (COCA), j Colon Adenoma (COA), k Colorectal Carcinoma, l Rectal Adenocarcinoma (READ), m Colorectal Carcinoma

Table 1 Expression of PDE9A among different subtypes of colorectal cancer and normal individuals using the Oncomine database
Fig. 3

PDE9A mRNA expression level across various cancer types was obtained from GENT2 database in normal and tumor tissues. Red color designates a boxplot of cancer samples whereas blue specifies a boxplot of normal samples

Fig. 4

The figure represents antibody staining immunohistochemical images of PDE9A protein expression in colorectal cancer tissues by the Human Protein Atlas website (a) sample of a male patient aged 63, (b) sample of a female patient aged 57

mRNA expression and methylation status of the PDE9A gene with various clinicopathological parameters in CRC patients

Based on different clinicopathological parameters, the UALCAN database has been employed to analyze the expression pattern and methylation status of the PDE9A gene in COAD. Statistical analysis from this database revealed that the expression of PDE9A is significantly attenuated in colon cancer samples compared with the normal counterparts. The connection between PDE9A expression and different clinicopathological characteristics in CRC samples exhibits that PDE9A is down-regulated in all the variables (sample types, individual cancer stages, patient’s race, patient’s weight, patient’s gender, patient’s age, histological subtype, nodal metastasis status, and TP53 mutation status) used here compared to the normal (Fig. 5 and Supplementary Table 1). Next, promoter methylation of PDE9A had been analyzed as it is thought that epigenetic modulation is the key factor to reduce mRNA expression and methylation levels in the promoter regions. Thus methylation is closely associated with the development of tumors. The PDE9A promoter methylation had been analyzed in colon cancer tissues. It was found that the promoter methylation level of PDE9A in sample types, individual cancer stages, patient race, patient gender, patient age, patient weight, tumor histology, and T53 mutation status was reduced than normal tissues (Fig. 6 and Supplementary Table 2). Furthermore, the correlation analysis between PDE9A expression and DNA methylation heat map in the TCGA colon cancer sample types depicted that PDE9A expression positively related to PDE9A promoter methylation at some CpG sites (blank frame) (Supplementary Fig. 2). Hence, the results implied that the status of methylation in the PDE9A DNA promoter region can be related to PDE9A expression in cancer tissues. Besides, the GEPIA database was also employed to identify mRNA expression of PDE9A between COAD samples and normal tissue samples. The PDE9A expression is significantly lower or downregulated (significance; P > 0.05) in COAD than compared with normal tissues (Fig. 7a). The relatedness between PDE9A mRNA expression levels and different tumor stages of COAD was also evaluated, where PDE9A is significantly downregulated in stage II (Fig. 7b). To establish our expression concept strongly in COAD, we further used the Gent2 database to add COAD subtype profiling for PDE9A. The database provides five Subtypes (Molecular subtype, AJCC Stage, Duke Stage, Grade, and Histology) for colon cancer analyzing 1146 samples. The PDE9A gene expression in colon cancer tissue based on different subtypes are shown in Supplementary Fig. 3 via box and dot-plot.

Fig. 5

mRNA expression analysis of PDE9A based on different variables for Colon Adenocarcinoma (COAD) patients. Expression of PDE9A in COAD based on (a) sample types, (b) individual cancer stages, (c) patient race, (d) patient gender, (e) patient weight, (f) patient weight, (g) tumor histological subtype, (h) nodal metastasis status, and (i)T53 mutation status

Fig. 6

PDE9A promoter methylation profile based on different variables of TCGA colon adenocarcinoma (COAD) from UALCAN database. The variables are (a) sample types, (b) individual cancer stages, (c) patient’s race, (d) patient’s gender, (e) patient’s age, (f) patient’s weight, (g) nodal metastasis status, (h) histological subtype, and (i) TP53 mutation status. According to UALCAN, the Beta value indicates the level of DNA methylation ranging from 0 (unmethylated) to 1 (fully methylated). Different beta value cut-off has been considered to indicate hypermethylation [Beta value: 0.7–0.5] or hypo-methylation [Beta-value: 0.3–0.25]

Fig. 7

a Box plot of PDE9A mRNA expression using GEPIA database in Colon Adenocarcinoma (COAD) patients, b The connection between PDE9A expression with different tumor stages in patients with COAD from GEPIA database

Survival plot analysis and prognostic value

The survival plot provides basic statistical concepts to analyze the patient’s survival period on-time event for cancer patients. Survival plot can be Kaplan-Meier plots to visualize survival curves, Log-rank test to compare the survival curves of two or more groups, Cox proportional hazards regression to describe the effect of variables on survival. The most important measures for survival plot in cancer studies include the overall survival (OS) curve, Relapse-free survival (RFS) or Disease-free survival (DFS) curve, and Disease-specific survival (DSS) curve. The relapse-free survival time represents the time between the response of patients to treatment and recurrence of the disease and Disease-free survival is corresponds with the percentage of people in a study who have not died from a specific disease in a defined period. We used different databases namely Gent2 (Fig. 8), PrognoScan (Fig. 9, supplementary Table 3), GEPIA (Fig. 10), UALCAN (Supplementary Fig. 4), OncoLnc (Supplementary Fig. 4), and R2 (Supplementary Fig. 4) to analyze the survival plot or Kaplan-Meier plot for the PDE9A gene against colon cancer. OncoLnc and GENT2 provide the Log-rank test curve and other databases provide Kaplan-Meier plots. The Kaplan–Meier curves disclosed that high expression of the PDE9A gene was related to encouraging conditions in OS (overall survival), RFS (relapse-free survival), and DSS (disease-specific survival) in colon cancer patients. Besides, the Log-rank test curve also depicted that a higher survival rate was detected with high levels of PDE9A expression in colon cancer patients. Conversely, low PDE9A expression is correlated with poor survival, but high PDE9A expression is associated with a higher survival rate. Excluding little exceptions, a low prognostic value was observed by analyzing the survival plot in the case of OS, RFS, and DSS in colon cancer patients. Thus, this study suggests that high expression of the PDE9A gene is positively correlated with a good prognosis in colon cancer patients.

Fig. 8

The prognostic outcome of PDE9A in Colon Adenocarcinoma (COAD) patients using GENT2 database. Relation of PDE9A expression level and survival rate based on different variables were analyzed for (a) AJCC stage plot divided by median cut off Disease-specific survival; (b) AJCC stage plot divided by median cutt off Overall survival; (c) AJCC stage plot divided by subtype Disease-specific survival; (d) AJCC stage plot divided by subtypes Overall Survival; (e) Grade plot divided by median cut off Disease-specific survival; (f) Grade plot divided by median cut off Overall survival; (g) Grade plot divided by subtype Disease-specific survival; (h) Grade plot divided by subtype Overall survival; (i) Histology plot divided by median cut off Disease-specific survival; (j) Histology plot divided by median cut off Overall survival; (k) Histology plot divided by subtypes Disease-specific survival; (l) Histology plot divided by subtypes Overall survival; (m) Duke stage plot divided by median cutt off Disease-specific survival; (n) Duke stage plot divided by median cutt off Overall survival; (o) Subtype plot divided by median cut off Disease-specific survival; (p) Subtype plot divided by median cut-off Overall Survival

Fig. 9

Kaplan-Meier plot of PDE9A from PrognoScan Web server. a Disease-free survival; b Disease-specific survival; c Overall survival

Fig. 10

Survival assay of PDE9A in Colon adenocarcinoma (COAD) patient using survival plot from GEPIA database. The survival plots are shown for (a) overall survival (Days) and (b) overall survival (Months); (c) Dease-free survival (Days) and (d) Dease-free survival (Months)

Interacting network analysis

To investigate the interaction network, we implemented GeneMANIA and STRING, two different web-based network analysis tools. As we all know that Protein-protein interactions (PPIs) perform a crucial role in cellular functions and biological signaling in all organisms, this may facilitate extensive information regarding interactions and pathways [56]. GeneMANIA provides an interaction network based on the parameters of physical interaction, genetic interactions, co-expression, pathway, co-localization, and shared protein domain. While STRING arranges the physical and functional interactions of genes. To identify the PPIs involving with PDE9A, GeneMANIA provided predicted protein partners are POU2F1, PDE5A, PDE3A, PDE3B, PDE2A, PDE11A, PDE1A, PDE10A, PDE6B, PDE1B, PDE6A, PDE6G, GUCY1A2, KCNMB2, MAPK14, MRVI1, KCNMB4, KCNMB3, PDE7A, KCNMB1, GUCY1A3, PRKG2, NOS1, NOS2, and KCNMA1. KCNMA1, KCNMB1, and KCNMB2 share physical interaction; PDE5A, PDE3A, PDE3B, PDE2A, PDE1A, PDE10A, PDE6B, PDE1B, PDE6A, PDE6G, GUCY1A2, MRVI1, KCNMB4, KCNMB3, PDE7A, GUCY1A3, PRKG2, NOS1, and NOS2 shares co-expression interaction (Fig. 11a). PDE9A contributes consolidated pathways (platelet homeostasis mediated) with MAPK14 which can carry out an important role in colon cancer prognosis. Recent research demonstrates that overexpression of active MAPK14 induces survival autophagy in CRC cells depleted of TP53 and that inhibition of autophagy in such cells enhances the cytotoxic effect [57]. This interaction may lead to PDE9A high expression to a higher survival rate for colon cancer patients. PPI network analysis obtained from the STRING database revealed that PDE9A showed a high confidence interaction (score 0.813) with GUCY1A2 (Guanylate cyclase soluble subunit alpha-2; Has guanylyl cyclase on binding to the beta-1 subunit). It also shares interaction with GUCY1B3, GUCY1A3, ADCY10, NPR1, NPR2, ADCY5, GUCY2D, GUCY2F, PRKG2, GUCA1A, GUCA1B, NPPC, GUCA1C, NPPA, MRVI1, GNAL, GNAS, GIPR, and MTNR1B (Fig. 11b). The network stats are the number of nodes: 21; the number of edges: 89; average node degree:8.48; avg. local clustering coefficient: 0.81 and PPI enrichment p-value: < 1.0e-16. These predicted interacting genes of PDE9A may be involved in the regulation of cancer progression and prognosis.

Fig. 11

a Interaction network of PDE9A gene derived from GeneMANIA. GeneMANIA shows the interactions for genes linked with each other according to physical interaction, genetic interactions, co-expression, pathway, co-localization, and shared protein domain. b Protein-protein interaction network for PDE9A gene from STRING database. The interaction network shows that these 20 proteins significantly interact with PDE9A. According to the website, colored nodes represent query proteins and first shell of interactors; white nodes represent second shell of interactors; empty nodes represent proteins of unknown 3D structure; filled nodes represent some 3D structure is known or predicted

Co-expression analysis

To identify a co-ordinated expression mediated cancer prognosis of PDE9A, we employed the co-expression profile via UALCAN website. PDE9A shows the highest positive correlation with CEACAM7 in COAD (Pearson CC = 0.54) (Fig. 12a and b). This positive interaction was confirmed by adopting the Oncomine database which ensured that CEACAM7 was also under-expressed in colon cancer tissues. To ensure this positive co-expression pattern between PDE9A and CEACAM7, we also look into TCGA colon cancer patient data via UCSC Xena (Fig. 12c). R2 database was also used to explore the relations of the PDE9A gene with the CEACAM7 gene using correlation statistics. The co-expression between these two genes was identified in Tumor Colon Adenocarcinoma - TCGA - 286 - rsem - tcgars dataset. The analysis reveals r-value = − 0.306; p-value = 1.31e-07; T-value = 5.415; degrees of freedom = 284 which declares the significance of the correlation in colon cancer sample. Figure 12d shows their positive correlation and the expression of CEACAM7 follows the expression of PDE9A. The figure portrays the expression of both genes in Tumor Colon Adenocarcinoma dataset which also indicates that the increasing PDE9A expression was followed by increasing expression of CEACAM7 in a good way.

Fig. 12

The co-expression profile between the two correlated gene PDE9A and CEACAM7. Gene expression between PDE9A and CEACAM7 using UALCAN database a based on tumor; b based on patient’s race- Caucasian, African- American, Asian; c Correlation heat map of PDE9A and CEACAM7 across sample types in the TCGA colon cancer database using UCSC Xena. d The correlated expression of the PDE9A gene with the CEACAM7 gene in Tumor Colon Adenocarcinoma - TCGA - 286 - rsem - tcgars dataset collected from R2 web tool


Phosphodiesterase 9A (PDE9A) is a protein-coding gene, which catalyzes the hydrolysis of the intracellular second messenger levels cyclic adenosine 3′,5′-monophosphate (cAMP) and cyclic guanosine 3′,5′-monophosphate (cGMP) [18]. cGMP performs crucial roles in cell proliferation, differentiation, and apoptosis via intracellular cGMP signaling [58, 59] and regulates diverse physiological functions covering platelet aggregation [60], neurotransmission [61], and vascular smooth-muscle modulation [62]. Several studies have been reported that dysregulation of PDE activities is associated with diverse tumors and carcinomas like chronic lymphocytic leukemia, glioblastoma, prostate cancer, breast cancer, adrenocortical tumor, etc. [24, 63, 64]. However, the distinct role of PDE9A in CRC is yet to be explored. In the present investigation, the mRNA expression levels, methylation status, prognostic values, PPI interactions, co-expressions in COAD patients were investigated through different databases. Here, we have utilized several bioinformatics databases including Oncomine, UALCAN, GEPIA, GENT2, R2, OncoLnc, PrognoScan, cBioportal, UCSC Xena, GeneCards, GeneMania, and STRING to detect the role of PDE9A in cancer prognosis and progression.

Here, the mRNA expression patterns of the PDE9A in CRC among specific parameters were evaluated adopting Oncomine, UALCAN, GENT2, GEPIA databases. mRNA expression analysis from all these databases demonstrates that PDE9A was conspicuously downregulated in CRC tissues with their corresponding normal tissues. PDE9A expression was downregulated in various colorectal cancer types including Colorectal Adenoma, Colon Mucinous Adenocarcinoma, Colon Adenocarcinoma, Colon Adenoma, Colon Carcinoma, Colorectal Carcinoma, etc. (Fig. 2 and Table 1). Moreover, mRNA expression among various clinicopathological criteria from the UALCAN database also ensured the authenticity that PDE9A was downregulated in all the variables used here (Fig. 5). mRNA expression of cGMP-Specific Phosphodiesterase 9A in various tissues was also reported (Figs. 3 and 7 and Supplementary Fig. 1, 3) [65, 66].

Next, we had scrutinized promoter DNA methylation status through the UALCAN database, where promoter methylation of PDE9A in COAD was attenuated than normal tissues based on sample types, individual cancer stages, patient race, patient gender, patient age, patient weight, tumor histology, and T53 mutation status (Fig. 6). Furthermore, we correlated PDE9A gene expression and DNA methylation heat map in TCGA colon cancer sample types by adopting the UCSC Xena database (Supplementary Fig. 2). Promoter DNA methylation status analysis from both these databases revealed that variations in the DNA methylation status are considered to be associated in the PDE9A gene expression in different stages of CRC. Diversified genetic and epigenetic modifications are responsible to arise CRC among which DNA methylation is thought one of the most influential epigenetic events occurs during the early stages of such oncogenic transformation [67]. Here, we noticed that DNA hypomethylation is correlated with down-regulated gene expression in CRC. Alternatively, in prostate cancer DNA hypermethylation is demonstrated to be associated with upregulated gene expression resembles the diversity of epigenetic regulation [68]. Thus, the analysis of promoter methylation standards and its clinicopathological attributes in CRC may act as a potential biomarker for tracing the tumor [69,70,71].

Further, we investigated the prognostic relevance of PDE9A in CRC by analyzing the survival plot or Kaplan-Meier plot from different databases like cBioportal, R2, UALCAN, PrognoScan, OncoLnc, GEPIA, and GENT2, etc. (Figs. 8, 9 and 10 and Supplementary Fig. 4). The Kaplan–Meier curves analysis depicted that elevated PDE9A gene expression was associated with favorable conditions in OS, RFS, and DSS. On the contrary, lower PDE9A expression was correlated with a poor survival rate. Several studies have conferred to investigate the prognostic value of various genes and their survival outcomes in CRC to recognize genes that could potentially perform as novel prognostic predictors [72,73,74]. In Malaysia, survival study and prognostic circumstances for CRC patients were analyzed for effective initial detection and advancements in cancer medication [75]. Therefore, from the survival plot analysis, we can confer that PDE9A expression might be associated with CRC progression and prognosis.

Protein-protein interaction networks were investigated from GeneMANIA and STRING databases, where GeneMANIA facilitates an enormous set of functional association data [76] and STRING arranges the physical and functional interactions of genes (Fig. 11) [55]. GeneMANIA contributes to predicting the PDE9A interaction patterns with other physically interacting proteins. These physical interactions depicted that PDE9A shares consolidated pathways (platelet homeostasis mediated) with MAPK14. MAPK14 had been reported in the progress of irinotecan resistance in HCT116 cells in which TP53 was wiped out that hinders cell proliferation and brings about survival-autophagy in CRC [57]. Such interactions may promote higher PDE9A expression to a higher survival rate in CRC patients. STRING network analysis for PDE9A demonstrated a higher confidence interaction with GUCY1A2. GUCY1A2 displayed a higher level (> = 50%) DNA hypermethylation that resembles a prognostic signature in colon cancer [77]. Besides, GUCY1A2 had been identified as biologically significant to rushing castration-resistant prostate cancer [78]. Thus, PDE9A expression might have significant relevance to compelling CRC that can perform as inhibitors for the identified targets.

For understanding, the molecular mechanism of CRC prognosis and progression with PDE9A expression must be scrutinized. We explored the co-expression analysis of the PDE9A gene adopting UALCAN, R2, and UCSC Xena databases (Fig. 12). The co-expression profile analysis adopting TCGA data through the UALCAN database demonstrated the highest positive correlation with CEACAM7 in COAD. R2 database and TCGA colon cancer patients data from UCSC Xena also ensured the positive correlation between PDE9A and CEACAM7. TCGA data analysis and Oncomine database analysis exhibited that it is under-expressed in CRC. CEACAM7 is a member of the carcinoembryonic antigen (CEA) protein family and encodes a cell surface glycoprotein. Down-regulated gene expression of this gene may be inspected in the colon and rectal cancer [79, 80]. CEACAM7 plays important role in cancer pathology and may provide valuable information to predict promising prognostic markers of CRC and to unravel the molecular mechanism of CRC [73, 81, 82].

In summary, through a standardized data mining process, PDE9A gene expression patterns and methylation status were examined handling openly accessible expression and clinicopathological data. These analyses suggest that PDE9A expression was downregulated in various CRC tissues with their corresponding normal tissues and promoter DNA methylation might have significant relevance with PDE9A expression. Besides, KM plotter analysis revealed prognostic significance in which elevated PDE9A gene expression was combined with positive conditions in OS, RFS, and DSS, alternatively lower PDE9A expression corresponded with a poor survival rate. In addition, PPIs and co-expression analysis demonstrated that PDE9A shares interaction networks with other genes that had a promising role in cancer pathology and may act as a potential biomarker or prognostic marker. Some functional proteins were co-expressed with PDE9A like CEACAM7, a prognostication marker for rectal cancer reappearance. Consequently, it is speculated that PDE9A is a promising prognosis predictor and curative target of CRC.

As the current research centered solely on in silico analysis, a large scale clinical experiment is needed to scrutinize the molecular mechanism of PDE9A in CRC both in vitro and in vivo. Therefore the present study could be validated through multiple wet lab experiment including Gene expression analysis [83], Human Cell Line/Animal Model experiment [84], Gene knockout/knockdown and gene knock-in experiment [85].


The present research was conducted to explore the role of the PDE9A gene in CRC using bioinformatics analysis methods. We executed a significant onco-informatics analysis to explore the expression profile and various clinicopathological parameters of the PDE9A gene in CRC. Extensive data mining from various publicly available databases revealed that PDE9A was downregulated in various CRC corresponding with their normal tissues. Besides, promoter DNA methylation status, survival plot analysis, PPIs, and co-expression analysis demonstrated that PDE9A has significant prognostic and clinicopathological value in CRC. Thus, PDE9A might have a potential function as a prognostic biomarker or tumor marker.

Availability of data and materials

All data generated and analyzed during this study are included in this article including: - Ensembl database; STRING database (; geneMANIA; R2 platform (;OncoLnc (; Gene Expression Profiling Interactive Analysis (GEPIA) (; Gene Expression database of Normal and Tumor tissues 2 (GENT2) (; Oncomine database (

UALCAN database (; cBioPortal for Cancer Genomics (;



American Joint Committee on Cancer


cyclic adenosine 3′,5′-monophosphate and


cyclic guanosine 3′,5′-monophosphate


carcinoembryonic antigen


Colorectal Adenoma


Colon Adenoma


colorectal adenocarcinoma


Colon Carcinoma


Colon Mucinous Adenocarcinoma


copy number alterations


Cecum Adenocarcinoma


Colorectal cancer


Colorectal Carcinoma


Disease-free survival


Disease-specific survival




Long non-coding


overall survival


Cyclic nucleotide phosphodiesterase


phosphodiesterase 9A


Rectal Adenocarcinoma


Rectal Mucinous Adenocarcinoma


Relapse-free survival


Serial Analysis of Gene Expression


Stool DNA


The Cancer Genome Atlas Program


  1. 1.

    Rivu SF, Apu MNH, Shabnaz S, Nahid NA, Islam MR, al-Mamun MMA, et al. Association of TP53 codon 72 and CDH1 genetic polymorphisms with colorectal cancer risk in Bangladeshi population. Cancer Epidemiol. 2017;49:46–52.

    Article  PubMed  Google Scholar 

  2. 2.

    The, I., T.P.-C.A. of Whole, and G. Consortium. Pan-cancer analysis of whole genomes. Nature. 2020;578(7793):82.

    Article  Google Scholar 

  3. 3.

    Raza AM, Kamal M, Begum F, Yusuf MA, Mohammad D, Begum M, et al. Clinico-demographic characteristics of colorectal carcinoma in Bangladeshi patients. J Curr Adv Med Res. 2016;3(1):22–5.

    Article  Google Scholar 

  4. 4.

    Kobaek-Larsen M, Thorup I, Diederichsen A, Fenger C, Hoitinga MR. Review of colorectal cancer and its metastases in rodent models: comparative aspects with those in humans. Comparative Med. 2000;50(1):16–26.

    CAS  Google Scholar 

  5. 5.

    Hsiao K-Y, Lin YC, Gupta SK, Chang N, Yen L, Sun HS, et al. Noncoding effects of circular RNA CCDC66 promote colon cancer growth and metastasis. Cancer Res. 2017;77(9):2339–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Tauriello DV, et al. TGFβ drives immune evasion in genetically reconstituted colon cancer metastasis. Nature. 2018;554(7693):538–43.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    O'keefe SJ. Diet, microorganisms and their metabolites, and colon cancer. Nat Rev Gastroenterol Hepatol. 2016;13(12):691–706.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Hussain SA, Sullivan R. Cancer control in Bangladesh. Jpn J Clin Oncol. 2013;43(12):1159–69.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Rahman M, et al. Identification of prognostic biomarker signatures and candidate drugs in colorectal cancer: insights from systems biology analysis. Medicina. 2019;55(1):20.

    Article  PubMed Central  Google Scholar 

  10. 10.

    Haggar FA, Boushey RP. Colorectal cancer epidemiology: incidence, mortality, survival, and risk factors. Clin Colon Rectal Surg. 2009;22(4):191–7.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Doubeni CA, Corley DA, Quinn VP, Jensen CD, Zauber AG, Goodman M, et al. Effectiveness of screening colonoscopy in reducing the risk of death from right and left colon cancer: a large community-based study. Gut. 2018;67(2):291–8.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Wong JJL, Hawkins NJ, Ward RL. Colorectal cancer: a model for epigenetic tumorigenesis. Gut. 2007;56(1):140–8.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Xu X-L, Yu J, Zhang HY, Sun MH, Gu J, du X, et al. Methylation profile of the promoter CpG islands of 31 genes that may contribute to colorectal carcinogenesis. World J Gastroenterol: WJG. 2004;10(23):3441–54.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Levin B, Lieberman DA, McFarland B, Andrews KS, Brooks D, Bond J, et al. Screening and surveillance for the early detection of colorectal cancer and adenomatous polyps, 2008: a joint guideline from the American Cancer Society, the US multi-society task force on colorectal Cancer, and the American College of Radiology. Gastroenterology. 2008;134(5):1570–95.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Świderska M, et al. The diagnostics of colorectal cancer. Contemp Oncol. 2014;18(1):1–6.

    Article  Google Scholar 

  16. 16.

    Ng JM-K, Yu J. Promoter hypermethylation of tumour suppressor genes as potential biomarkers in colorectal cancer. Int J Mol Sci. 2015;16(2):2472–96.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Das V, Kalita J, Pal M. Predictive and prognostic biomarkers in colorectal cancer: a systematic review of recent advances and challenges. Biomed Pharmacother. 2017;87:8–19.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Azevedo MF, Faucz FR, Bimpaki E, Horvath A, Levy I, de Alexandre RB, et al. Clinical and molecular genetics of the phosphodiesterases (PDEs). Endocr Rev. 2014;35(2):195–233.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Guipponi M, Scott HS, Kudoh J, Kawasaki K, Shibuya K, Shintani A, et al. Identification and characterization of a novel cyclic nucleotide phosphodiesterase gene (PDE9A) that maps to 21q22. 3: alternative splicing of mRNA transcripts, genomic structure and sequence. Hum Genet. 1998;103(4):386–92.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Lee DI, Zhu G, Sasaki T, Cho GS, Hamdani N, Holewinski R, et al. Phosphodiesterase 9A controls nitric-oxide-independent cGMP and hypertrophic heart disease. Nature. 2015;519(7544):472–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Singh N, Patra S. Phosphodiesterase 9: insights from protein structure and role in therapeutics. Life Sci. 2014;106(1–2):1–11.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Dorner-Ciossek C, Kroker K, Rosenbrock H. Role of PDE9 in Cognition. In: Phosphodiesterases: CNS Functions and Diseases: Springer; 2017. p. 231–54.

    Chapter  Google Scholar 

  23. 23.

    Savai R, Pullamsetti SS, Banat GA, Weissmann N, Ghofrani HA, Grimminger F, et al. Targeting cancer with phosphodiesterase inhibitors. Expert Opin Investig Drugs. 2010;19(1):117–31.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Pinto EM, Faucz FR, Paza LZ, Wu G, Fernandes ES, Bertherat J, et al. Germline variants in Phosphodiesterase genes and genetic predisposition to pediatric adrenocortical tumors. Cancers. 2020;12(2):506.

    CAS  Article  PubMed Central  Google Scholar 

  25. 25.

    Szarek E, Stratakis C. Phosphodiesterases and adrenal Cushing in mice and humans. Hormone Metab Res. 2014;46(12):863.

    CAS  Article  Google Scholar 

  26. 26.

    Goldhoff P, Warrington NM, Limbrick DD Jr, Hope A, Woerner BM, Jackson E, et al. Targeted inhibition of cyclic AMP phosphodiesterase-4 promotes brain tumor regression. Clin Cancer Res. 2008;14(23):7717–25.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Cesarini V, Martini M, Vitiani LR, Gravina GL, di Agostino S, Graziani G, et al. Type 5 phosphodiesterase regulates glioblastoma multiforme aggressiveness and clinical outcome. Oncotarget. 2017;8(8):13223–39.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Zhang L, Murray F, Zahno A, Kanter JR, Chou D, Suda R, et al. Cyclic nucleotide phosphodiesterase profiling reveals increased expression of phosphodiesterase 7B in chronic lymphocytic leukemia. Proc Natl Acad Sci. 2008;105(49):19532–7.

    Article  PubMed  Google Scholar 

  29. 29.

    Soderling SH, Bayuga SJ, Beavo JA. Identification and characterization of a novel family of cyclic nucleotide phosphodiesterases. J Biol Chem. 1998;273(25):15553–8.

    CAS  Article  Google Scholar 

  30. 30.

    Wang Q, Li W, Liu XS, Carroll JS, Jänne OA, Keeton EK, et al. A hierarchical network of transcription factors governs androgen receptor-dependent prostate cancer growth. Mol Cell. 2007;27(3):380–92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Sato T, Arai E, Kohno T, Tsuta K, Watanabe SI, Soejima K, et al. DNA methylation profiles at precancerous stages associated with recurrence of lung adenocarcinoma. PLoS One. 2013;8(3):e59444.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Krishnan MS, Racsa M, Yu H-HM. Handbook of Supportive and Palliative Radiation Oncology: Academic Press; 2016.

    Google Scholar 

  33. 33.

    Chen F, et al. Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers. Nat Commun. 2019;10(1):1–15.

    Article  Google Scholar 

  34. 34.

    Hou G-X, Liu P, Yang J, Wen S. Mining expression and prognosis of topoisomerase isoforms in non-small-cell lung cancer by using Oncomine and Kaplan–Meier plotter. PLoS One. 2017;12(3):e0174515.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Yan P, He Y, Xie K, Kong S, Zhao W. In silico analyses for potential key genes associated with gastric cancer. PeerJ. 2018;6:e6092.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Rhodes DR, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia (New York, NY). 2004;6(1):1.

    CAS  Article  Google Scholar 

  37. 37.

    Rhodes DR, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia (New York, NY). 2007;9(2):166.

    CAS  Article  Google Scholar 

  38. 38.

    Wu Y, Xu Y. Integrated bioinformatics analysis of expression and gene regulation network of COL12A1 in colorectal cancer. Cancer Med. 2020.

  39. 39.

    Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Zheng H, Zhang G, Zhang L, Wang Q, Li H, Han Y, et al. Comprehensive review of web servers and bioinformatics tools for cancer prognosis analysis. Front Oncol. 2020;10:68.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Anaya J. OncoLnc: linking TCGA survival data to mRNAs, miRNAs, and lncRNAs. PeerJ Comput Sci. 2016;2:e67.

    Article  Google Scholar 

  42. 42.

    Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genet. 2009;2(1):18.

    CAS  Article  Google Scholar 

  43. 43.

    Park S-J, Yoon BH, Kim SK, Kim SY. GENT2: an updated gene expression database for normal and tumor tissues. BMC Med Genet. 2019;12(5):101.

    Article  Google Scholar 

  44. 44.

    Ouyang S, Liu J‑H, Ni Z, Ding G‑F, Wang Q‑Z. Downregulation of ST3GAL5 is associated with muscle invasion, high grade and a poor prognosis in patients with bladder cancer. Oncol Lett. 2020;20(1):828–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Goldman Mary J, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nature biotechnology.  2020;38(6):675–8.

  46. 46.

    Cui X, Yi Q, Jing X, Huang Y, Tian J, Long C, et al. Mining prognostic significance of MEG3 in human breast cancer using bioinformatics analysis. Cell Physiol Biochem. 2018;50(1):41–51.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Cerami Ethan, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. 2012:401–404.

  48. 48.

    Gao J, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.

    Article  Google Scholar 

  49. 49.

    Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Xue C, Zhang J, Zhang G, Xue Y, Zhang G, Wu X. Elevated SPINK2 gene expression is a predictor of poor prognosis in acute myeloid leukemia. Oncol Lett. 2019;18(3):2877–84.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Safran M, Dalah I, Alexander J, Rosen N, Iny Stein T, Shmoish M, et al. GeneCards version 3: the human gene integrator. Database. 2010;2010(0).

  52. 52.

    Stelzer G, et al. The GeneCards suite: from gene data mining to disease genome sequence analyses. Curr Protocols Bioinform. 2016;54(1):1.30 1–1.30. 33.

    Article  Google Scholar 

  53. 53.

    Warde-Farley D, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(suppl_2):W214–20.

    CAS  Article  Google Scholar 

  54. 54.

    Szklarczyk Damian, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic acids research. 2016: gkw937.

  55. 55.

    Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

    CAS  Article  Google Scholar 

  56. 56.

    De Las Rivas J, Fontanillo C. Protein–protein interactions essentials: key concepts to building and analyzing interactome networks. PLoS Comput Biol. 2010;6(6):e1000807.

    Article  Google Scholar 

  57. 57.

    Paillas S, Causse A, Marzi L, de Medina P, Poirot M, Denis V, et al. MAPK14/p38α confers irinotecan resistance to TP53-defective cells by inducing survival autophagy. Autophagy. 2012;8(7):1098–112.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Loweth AC, Williams GT, Scarpello JHB, Morgan NG. Evidence for the involvement of cGMP and protein kinase G in nitric oxide-induced apoptosis in the pancreatic B-cell line, HIT-T15. FEBS Lett. 1997;400(3):285–8.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Shimojo T, Hiroe M, Ishiyama S, Ito H, Nishikawa T, Marumo F. Nitric oxide induces apoptotic death of cardiomyocytes via a cyclic-GMP-dependent pathway. Exp Cell Res. 1999;247(1):38–47.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Eigenthaler M, et al. Signal transduction by cGMP-dependent protein kinases and their emerging roles in the regulation of cell adhesion and gene expression. Rev Physiol Biochem Pharmacol. 1999;135:173–209.

    CAS  Article  Google Scholar 

  61. 61.

    Smolenski A, Burkhardt AM, Eigenthaler M, Butt E, Gambaryan S, Lohmann SM, et al. Functional analysis of cGMP-dependent protein kinases I and II as mediators of NO/cGMP effects. Naunyn Schmiedeberg's Arch Pharmacol. 1998;358(1):134–9.

    CAS  Article  Google Scholar 

  62. 62.

    Vaandrager AB, de Jonge HR. Signalling by cGMP-dependent protein kinases. Mol Cell Biochem. 1996;157(1–2):23–30.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Saravani R, Karami-Tehrani F, Hashemi M, Aghaei M, Edalat R. Inhibition of phosphodiestrase 9 induces c GMP accumulation and apoptosis in human breast cancer cell lines, MCF-7 and MDA-MB-468. Cell Prolif. 2012;45(3):199–206.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Fajardo AM, Piazza GA, Tinsley HN. The role of cyclic nucleotide signaling pathways in cancer: targets for prevention and treatment. Cancers. 2014;6(1):436–58.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Andreeva SG, Dikkes P, Epstein PM, Rosenberg PA. Expression of cGMP-specific phosphodiesterase 9A mRNA in the rat brain. J Neurosci. 2001;21(22):9068–76.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Razmkhah F, Ghasemi S, Soleimani M, Amini Kafi-abad S. LY86, LRG1 and PDE9A genes overexpression in umbilical cord blood hematopoietic stem progenitor cells by acute myeloid leukemia (M3) microvesicles. Exp Hematol Oncol. 2019;8(1):23.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Shen L, Kantarjian H, Guo Y, Lin E, Shan J, Huang X, et al. DNA methylation predicts survival and response to therapy in patients with myelodysplastic syndromes. J Clin Oncol. 2010;28(4):605–13.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Rauluseviciute I, Drabløs F, Rye MB. DNA hypermethylation associated with upregulated gene expression in prostate cancer demonstrates the diversity of epigenetic regulation. BMC Med Genet. 2020;13(1):6.

    CAS  Article  Google Scholar 

  69. 69.

    Li D, Bai Y, Feng Z, Li W, Yang C, Guo Y, et al. Study of promoter methylation patterns of HOXA2, HOXA5, and HOXA6 and its clinicopathological characteristics in colorectal cancer. Front Oncol. 2019;9:394.

    Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Kerachian MA, et al. Crosstalk between DNA methylation and gene expression in colorectal cancer, a potential plasma biomarker for tracing this tumor. Sci Rep. 2020;10(1):1–13.

    Article  Google Scholar 

  71. 71.

    Freitas M, Ferreira F, Carvalho S, Silva F, Lopes P, Antunes L, et al. A novel DNA methylation panel accurately detects colorectal cancer independently of molecular pathway. J Transl Med. 2018;16(1):45.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Idos GE, et al. The prognostic implications of tumor infiltrating lymphocytes in colorectal Cancer: a systematic review and meta-analysis. Sci Rep. 2020;10(1):1–14.

    Article  Google Scholar 

  73. 73.

    Bian Q, Chen J, Qiu W, Peng C, Song M, Sun X, et al. Four targeted genes for predicting the prognosis of colorectal cancer: a bioinformatics analysis case. Oncol Lett. 2019;18(5):5043–54.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  74. 74.

    Yu Y, Carey M, Pollett W, Green J, Dicks E, Parfrey P, et al. The long-term survival characteristics of a cohort of colorectal cancer patients and baseline variables associated with survival outcomes with or without time-varying effects. BMC Med. 2019;17(1):150.

    Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Hassan MRA, et al. Survival analysis and prognostic factors for colorectal cancer patients in Malaysia. Asian Pac J Cancer Prev. 2016;17(7):3575–81.

    PubMed  Google Scholar 

  76. 76.

    Franz M, Rodriguez H, Lopes C, Zuberi K, Montojo J, Bader GD, et al. GeneMANIA update 2018. Nucleic Acids Res. 2018;46(W1):W60–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  77. 77.

    Yi JM, Dhir M, van Neste L, Downing SR, Jeschke J, Glöckner SC, et al. Genomic and epigenomic integration identifies a prognostic signature in colon cancer. Clin Cancer Res. 2011;17(6):1535–45.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  78. 78.

    Rotimi S, et al. Gene expression profiling analysis reveals putative Phytochemotherapeutic target for castration-resistant prostate Cancer. Front Oncol. 2019;9:714.

    Article  PubMed  PubMed Central  Google Scholar 

  79. 79.

    Thompson J, Seitz M, Chastre E, Ditter M, Aldrian C, Gespach C, et al. Down-regulation of carcinoembryonic antigen family member 2 expression is an early event in colorectal tumorigenesis. Cancer Res. 1997;57(9):1776–84.

    CAS  PubMed  Google Scholar 

  80. 80.

    Thompson J, et al. CGM2, a member of the carcinoembryonic antigen gene family is down-regulated in colorectal carcinomas. J Biol Chem. 1994;269(52):32924–31.

    CAS  Article  Google Scholar 

  81. 81.

    Gemei M, et al. Carcinoembryonic antigen family cell adhesion molecules (CEACAM) as colorectal cancer biomarkers. In: Biomarkers in Disease: Methods, Discoveries and Applications: Biomarkers in Cancer: Springer Netherlands; 2015. p. 685–705.

    Google Scholar 

  82. 82.

    Messick CA, Sanchez J, DeJulius KL, Hammel J, Ishwaran H, Kalady MF. CEACAM-7: a predictive marker for rectal cancer recurrence. Surgery. 2010;147(5):713–9.

    Article  PubMed  Google Scholar 

  83. 83.

    Song W, et al. Rhomboid domain containing 1 promotes colorectal cancer growth through activation of the EGFR signalling pathway Nat. Commun. 2015;6(1):1–13.

    Google Scholar 

  84. 84.

    Feng Hailiang, et al. ALDH1A3 affects colon cancer in vitro proliferation and invasion depending on CXCR4 status. British journal of cancer. 2018;118(2):224–32.

  85. 85.

    Sun HD, et al. Down-regulation of circPVRL3 promotes the proliferation and migration of gastric cancer cells. Sci. Rep. 2018;8(1):1–13.

Download references


All the authors are grateful to Mahmudul Hasan (Assistant Professor, Department of Pharmaceuticals and Industrial Biotechnology, Faculty of Biotechnology and Genetic Engineering, Sylhet Agricultural University, Sylhet-3100, Bangladesh) for assisting in plagiarism checking and literature review.


This study is supported by King Saud University (KSU), Researchers Supporting Project number (RSP-2020/257), King Saud University, Riyadh, Saudi Arabia.

Author information




Conceptualization and Methodology was performed by TFS, AR, MMRK, FY, MSI, MUH; Formal analysis and Data curation was performed by TFS, AR, MMRK and FY; Writing-original draft prepared by TFS, AR, MMRK, FY and MSI; Writing-review and editing performed by MUH, GB, ON, SAA, and MSI; Supervised by MUH. All authors have read and agreed to submit the final version of the manuscript.

Corresponding author

Correspondence to Md. Shariful Islam.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they don’t have competing interests on this study.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table 01

: Differential expression analysis of PDE9A in COAD based on different variables using UALCAN database. Supplementary Table 02: PDE9A promoter methylation level based on different variables in COAD from UALCAN database. Supplementary Table 03: The biological relationship between gene expression and survival from PrognoScan database. Supplementary Fig. 01: (A + B) Expression of PDE9A across TCGA cancer with tumor and normal samples from UALCAN database, (C) PDE9A expression profile across all tumor samples and paired normal tissues from GEPIA database (Dot plot, each dots represent expression of samples), (D) PDE9A mRNA expression in normal human tissues based on RNAseq, Microarray, and SAGE using GeneCards database. Supplementary Fig. 02: (A) Heat map of PDE9A expression and DNA methylation status across TCGA Colon cancer sample types from UCSC Xena genome browser, (B) PDE9A expression in different colon cancer DNA methylation clusters from UCSC Xena. Supplementary Fig. 03: Colon cancer tissue subtype profiling analysis for PDE9A gene based on molecular subtype, AJCC stage, Duke stage, Grade, and Histology from Gent2 database. Supplementary Fig. 04: The prognostic value of PDE9A expression on colon cancer. Survival plot from UALCAN database (A) Effect of PDE9A expression level on COAD patient survival; (B) Effect of PDE9A expression level and body weight on COAD patient survival; (C) Effect of PDE9A expression level & Race on COAD patient survival; (D) Effect of PDE9A expression level & gender on COAD patient survival; (E) Survival curve from OncoLnc database, and (F) Kaplan–Meier plot from R2 database in Tumor Colon Adenocarcinoma - TCGA - 286 - rsem - tcgars dataset.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Susmi, T.F., Rahman, A., Khan, M.M.R. et al. Prognostic and clinicopathological insights of phosphodiesterase 9A gene as novel biomarker in human colorectal cancer. BMC Cancer 21, 577 (2021).

Download citation


  • PDE9A gene
  • Colorectal cancer
  • Methylation
  • Pathways
  • Biomarker