- Open Access
Prognostic and clinicopathological insights of phosphodiesterase 9A gene as novel biomarker in human colorectal cancer
BMC Cancer volume 21, Article number: 577 (2021)
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.
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 . Therefore, CRC is described as a malignant tumor that arises from the inner surface of the large intestine . Over 50% of patients die from colon cancer despite having intensive investigations and therapeutic improvements . However, the main reason behind the death can be pointed as a result of spreading the uncontrollable cell division tendency to other organs spontaneously . 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 (https://www.worldlifeexpectancy.com/). 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 . 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 . 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 . So, early identification of CRC is extremely important for the prevention of CRC . Perrectum examination is the simplest way to CRC recognition . However, the most efficient and commonly applied method of CRC recognition is endoscopy . 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) . 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 . PDE9A is a member of eleven PDE isoforms located at q22.3 position of the 21st chromosome in the human cell . 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 . Many hematological malignancies and carcinomas have been connected with decreased levels of cAMP and/or cGMP by the PDE gene family . PDE2A, PDE8B, and PDE11A expression are involved with adrenocortical tumors , high expression of PDE4 has been detected in brain tumors , in glioblastoma multiforme PDE5 is strongly expressed , and increased PDE7B expression found in Chronic lymphocytic leukemia . The cDNA of PDE9A was coined in 1998 and noted as the ninth member of the PDEs family . Only from the PDE9A gene, twenty-eight splice variants can arise (https://asia.ensembl.org/). These can create more proteins through alternative splicing of mRNA or multiple promoters and transcription start sites . 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 . The protein expression of PDE9A is highly conserved and also associated with prostate cancer . DNA hypermethylation of the PDE9A gene results in reduced mRNA expression in tumorous tissues . 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 . 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 . Using these onco-informatics processes it was found that DNA topoisomerases expression in lung cancer was higher than the normal tissues . 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 .
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 (http://www.oncomine.org/) 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 (http://ualcan.path.uab.edu/index.html) 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 (http://www.oncolnc.org/) online tools evaluate survival correlations with expression data and retrieving clinical data for mRNAs, miRNAs, and lncRNAs (long non-coding) . 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 . 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 (http://dna00.bio.kyutech.ac.jp/PrognoScan/) is a cancer microarray dataset that provides information for analyzing potential therapeutic value that could accelerate in silico cancer research by clinical annotation .
Clinicopathological insights of PDE9A
Gene Expression database of Normal and Tumor tissues 2 (GENT2) (http://gent2.appex.kr/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 (https://xena.ucsc.edu/), 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  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 . 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 (https://www.cbioportal.org/) 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) (http://gepia.cancer-pku.cn/index.html) was used to explore the differential expression analysis and also to determine the prognostic effects on COAD cancer patients . 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 .
Co-expression and interacting network analysis
GeneCards, the human gene database (https://www.genecards.org/) 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 (http://r2platform.com) 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 (http://www.genemania.org/) 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 . 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 (https://string-db.org/) 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.
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.
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.
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.
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 . 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 . 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.
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.
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) . cGMP performs crucial roles in cell proliferation, differentiation, and apoptosis via intracellular cGMP signaling [58, 59] and regulates diverse physiological functions covering platelet aggregation , neurotransmission , and vascular smooth-muscle modulation . 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 . 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 . 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 . 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  and STRING arranges the physical and functional interactions of genes (Fig. 11) . 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 . 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 . Besides, GUCY1A2 had been identified as biologically significant to rushing castration-resistant prostate cancer . 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 , Human Cell Line/Animal Model experiment , Gene knockout/knockdown and gene knock-in experiment .
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 https://asia.ensembl.org/; STRING database (https://string-db.org/); geneMANIA https://genemania.org/; R2 platform (http://r2platform.com);OncoLnc (http://www.oncolnc.org/); Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/index.html); Gene Expression database of Normal and Tumor tissues 2 (GENT2) (http://gent2.appex.kr/gent2/); Oncomine database (http://www.oncomine.org/).
American Joint Committee on Cancer
cyclic adenosine 3′,5′-monophosphate and
cyclic guanosine 3′,5′-monophosphate
Colon Mucinous Adenocarcinoma
copy number alterations
Cyclic nucleotide phosphodiesterase
Rectal Mucinous Adenocarcinoma
Serial Analysis of Gene Expression
The Cancer Genome Atlas Program
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. https://doi.org/10.1016/j.canep.2017.05.005.
The, I., T.P.-C.A. of Whole, and G. Consortium. Pan-cancer analysis of whole genomes. Nature. 2020;578(7793):82.
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. https://doi.org/10.3329/jcamr.v3i1.29388.
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.
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. https://doi.org/10.1158/0008-5472.CAN-16-1883.
Tauriello DV, et al. TGFβ drives immune evasion in genetically reconstituted colon cancer metastasis. Nature. 2018;554(7693):538–43. https://doi.org/10.1038/nature25492.
O'keefe SJ. Diet, microorganisms and their metabolites, and colon cancer. Nat Rev Gastroenterol Hepatol. 2016;13(12):691–706. https://doi.org/10.1038/nrgastro.2016.165.
Hussain SA, Sullivan R. Cancer control in Bangladesh. Jpn J Clin Oncol. 2013;43(12):1159–69. https://doi.org/10.1093/jjco/hyt140.
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. https://doi.org/10.3390/medicina55010020.
Haggar FA, Boushey RP. Colorectal cancer epidemiology: incidence, mortality, survival, and risk factors. Clin Colon Rectal Surg. 2009;22(4):191–7. https://doi.org/10.1055/s-0029-1242458.
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. https://doi.org/10.1136/gutjnl-2016-312712.
Wong JJL, Hawkins NJ, Ward RL. Colorectal cancer: a model for epigenetic tumorigenesis. Gut. 2007;56(1):140–8. https://doi.org/10.1136/gut.2005.088799.
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. https://doi.org/10.3748/wjg.v10.i23.3441.
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. https://doi.org/10.1053/j.gastro.2008.02.002.
Świderska M, et al. The diagnostics of colorectal cancer. Contemp Oncol. 2014;18(1):1–6. https://doi.org/10.5114/wo.2013.39995.
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. https://doi.org/10.3390/ijms16022472.
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. https://doi.org/10.1016/j.biopha.2016.12.064.
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. https://doi.org/10.1210/er.2013-1053.
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. https://doi.org/10.1007/s004390050838.
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. https://doi.org/10.1038/nature14332.
Singh N, Patra S. Phosphodiesterase 9: insights from protein structure and role in therapeutics. Life Sci. 2014;106(1–2):1–11. https://doi.org/10.1016/j.lfs.2014.04.007.
Dorner-Ciossek C, Kroker K, Rosenbrock H. Role of PDE9 in Cognition. In: Phosphodiesterases: CNS Functions and Diseases: Springer; 2017. p. 231–54.
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. https://doi.org/10.1517/13543780903485642.
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. https://doi.org/10.3390/cancers12020506.
Szarek E, Stratakis C. Phosphodiesterases and adrenal Cushing in mice and humans. Hormone Metab Res. 2014;46(12):863.
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. https://doi.org/10.1158/1078-0432.CCR-08-0827.
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. https://doi.org/10.18632/oncotarget.14656.
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. https://doi.org/10.1073/pnas.0806152105.
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.
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. https://doi.org/10.1016/j.molcel.2007.05.041.
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. https://doi.org/10.1371/journal.pone.0059444.
Krishnan MS, Racsa M, Yu H-HM. Handbook of Supportive and Palliative Radiation Oncology: Academic Press; 2016.
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.
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. https://doi.org/10.1371/journal.pone.0174515.
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. https://doi.org/10.7717/peerj.6092.
Rhodes DR, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia (New York, NY). 2004;6(1):1.
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.
Wu Y, Xu Y. Integrated bioinformatics analysis of expression and gene regulation network of COL12A1 in colorectal cancer. Cancer Med. 2020.
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. https://doi.org/10.1016/j.neo.2017.05.002.
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. https://doi.org/10.3389/fonc.2020.00068.
Anaya J. OncoLnc: linking TCGA survival data to mRNAs, miRNAs, and lncRNAs. PeerJ Comput Sci. 2016;2:e67. https://doi.org/10.7717/peerj-cs.67.
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. https://doi.org/10.1186/1755-8794-2-18.
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. https://doi.org/10.1186/s12920-019-0514-7.
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. https://doi.org/10.3892/ol.2020.11597.
Goldman Mary J, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nature biotechnology. 2020;38(6):675–8.
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. https://doi.org/10.1159/000493956.
Cerami Ethan, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. 2012:401–404.
Gao J, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.
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. https://doi.org/10.1093/nar/gkx247.
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. https://doi.org/10.3892/ol.2019.10665.
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). https://doi.org/10.1093/database/baq020.
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.
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.
Szklarczyk Damian, et al. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic acids research. 2016: gkw937.
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. https://doi.org/10.1093/nar/gky1131.
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.
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. https://doi.org/10.4161/auto.20268.
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. https://doi.org/10.1016/S0014-5793(96)01392-0.
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. https://doi.org/10.1006/excr.1998.4310.
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.
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. https://doi.org/10.1007/PL00005234.
Vaandrager AB, de Jonge HR. Signalling by cGMP-dependent protein kinases. Mol Cell Biochem. 1996;157(1–2):23–30. https://doi.org/10.1007/BF00227877.
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. https://doi.org/10.1111/j.1365-2184.2012.00819.x.
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. https://doi.org/10.3390/cancers6010436.
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. https://doi.org/10.1523/JNEUROSCI.21-22-09068.2001.
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. https://doi.org/10.1186/s40164-019-0147-8.
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. https://doi.org/10.1200/JCO.2009.23.4781.
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. https://doi.org/10.1186/s12920-020-0657-6.
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. https://doi.org/10.3389/fonc.2019.00394.
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.
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. https://doi.org/10.1186/s12967-018-1415-9.
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.
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. https://doi.org/10.3892/ol.2019.10866.
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. https://doi.org/10.1186/s12916-019-1379-5.
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.
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. https://doi.org/10.1093/nar/gky311.
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. https://doi.org/10.1158/1078-0432.CCR-10-2509.
Rotimi S, et al. Gene expression profiling analysis reveals putative Phytochemotherapeutic target for castration-resistant prostate Cancer. Front Oncol. 2019;9:714. https://doi.org/10.3389/fonc.2019.00714.
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.
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.
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.
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. https://doi.org/10.1016/j.surg.2009.10.056.
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.
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.
Sun HD, et al. Down-regulation of circPVRL3 promotes the proliferation and migration of gastric cancer cells. Sci. Rep. 2018;8(1):1–13.
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.
Ethics approval and consent to participate
Consent for publication
The authors declare that they don’t have competing interests on this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: 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.
About this article
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). https://doi.org/10.1186/s12885-021-08332-3
- PDE9A gene
- Colorectal cancer