Abstract
Background/Aim: Amplification of chromosome 7p (Ch.7p) is common in colorectal cancer (CRC). The aim of this study was to identify potential driver genes on Ch.7p that are overexpressed due to DNA copy number amplification and determine their clinical significance in CRC. Materials and Methods: We identified phosphoserine phosphatase (PSPH) as a potential driver gene using a CRC dataset from The Cancer Genome Atlas (TCGA) using a bioinformatics approach. The expression of PSPH in 124 primary CRCs was examined by quantitative reverse transcription polymerase chain reaction (PCR) and immunohistochemistry. The biological effect of PSPH expression was explored by Gene Set Enrichment Analysis (GSEA) using the TCGA dataset. Results: PSPH was overexpressed in tumor tissues and PSPH positively correlated with depth of invasion and distant metastasis. On multivariate analysis, high PSPH expression was an independent poor prognostic factor. These results were supported by GSEA. Conclusion: PSPH could be a novel prognostic biomarker with malignant potential on Ch.7p in CRC.
Colorectal cancer (CRC) is one of the most common types of cancer worldwide and the second most common cause of cancer mortality (1). Relapse of CRC frequently occurs despite multidisciplinary therapy, such as surgical resection, chemotherapy and molecular-targeted therapy (2).
Intra-tumor heterogeneity (ITH) is the existent molecular and cellular heterogeneity within a single tumor. ITH likely contributes to the occurrence of drug-resistant cells and increases the chance of therapeutic failure, since the presence of distinct subpopulations of cells with different sensitivities to anti-cancer drugs increases the chance of resistance and recurrence (3). Therefore, it is important to identify oncogenes that exist in tumor cells, as such oncogenes could be promising therapeutic targets despite ITH.
Amplification of chromosome 7 occurs frequently in CRC and colorectal adenoma tissues (4-6). By performing multiregional copy number analysis of CRC tumor tissues, we showed that amplification of the short arm of chromosome 7 (Ch.7p) occurs in all regions of an individual tumor (7). These data suggest that amplification of Ch.7p is a fundamental and predominant event in tumorigenesis of CRC. We hypothesized that Ch.7p harbors driver genes that promote tumorigenesis or tumor progression through gain-of-function in CRC.
In this study, we aimed to identify a novel potential driver gene on Ch.7p by a bioinformatics approach and determine the clinical significance of its expression in CRC.
Materials and Methods
Selection of candidate genes. We obtained the RNA sequencing data and the DNA copy number data of 615 CRC patients in The Cancer Genome Atlas (TCGA) from the Broad Institute's Firehose (http://gdac.broadinstitute.org/runs/stddata__2015_08_21/data/COADREAD/20150821/). The RNA sequencing data also included expression profiles of 51 paired normal colon samples. Using these data, we extracted the candidate genes from 426 genes on Ch.7p that satisfied the following two criteria. First, the DNA copy number and mRNA expression levels had to be positively correlated with each other. Cut-off of the correlation coefficient was determined as 0.4. Second, the gene in question was overexpressed in tumor tissues compared to normal tissues (>2-fold change).
Patients and clinical samples. All protocols used in this study were approved by the local ethics review board of Kyushu University. A total of 124 patients with CRC who underwent surgical resection of a primary tumor at Kyushu University Beppu Hospital and the affiliated hospitals between 1992 and 2007 were enrolled in this study. Clinicopathological factors and clinical stage were classified using the tumor-node-metastasis (TNM) system of classification. All patients were treated in accordance with the Japanese Society of Cancer of the Colon and Rectum Guidelines for the Treatment of Colorectal Cancer. Written informed consent was obtained from all patients. Resected tumor tissues and paired normal colon tissues were immediately stored in RNAlater (Ambion, Palo Alto, CA, USA), frozen in liquid nitrogen and kept at −80°C until RNA extraction.
RNA extraction and reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Total RNA was isolated from the frozen tissue specimens through use of ISOGEN (Nippon Gene, Tokyo, Japan), according to the manufacture's protocol. Complementary DNA (cDNA) was synthesized from 8 μg of total RNA using M-MLV reverse transcriptase (Invitrogen, Carlsbad, CA, USA). qPCR was performed using a LightCycler 480 and a SYBR Green I Master kit (Roche Applied Science, Basel, Switzerland) following the manufacturer's protocol. Gene expression was quantified using the following specific oligonucleotide primers: phosphoserine phosphatase (PSPH): 5’-GCAGTGCCTTTCAAAGCTGC-3’ (sense) and 5’-TAGGCGACTTACCAGCTCCC-3’ (antisense); glyceraldehyde-3-phosphate dehydrogenase (GAPDH): 5’-AGCCACATCGCTCAGACAC-3’ (sense) and 5’-GCCCAATACGACCAAATCC-3’ (antisense). The expression levels were normalized to the expression of GAPDH as an internal control in each sample. The expression levels were calculated as values relative to the expression level of Human Universal Reference Total RNA (Clontech Laboratories, Palo Alto, CA, USA).
Immunohistochemical analysis. Immunohistochemistry of PSPH in CRC cases was performed on formalin-fixed, paraffin-embedded tissues as previously described (8). A polyclonal rabbit anti-PSPH antibody (HPA020376; Sigma-Aldrich, St. Louis, MO, USA) was used as the primary antibody diluted at 1:200. Tumor histology was independently reviewed by an experienced pathologist.
Gene set enrichment analysis (GSEA). Using the RNA sequencing data of CRC tissues in TCGA, we performed GSEA (9) to elucidate the correlation between the candidate genes and gene sets that are associated with tumor progression.
Statistical analysis. Associations between the variables were tested by the Mann-Whitney U-test or Fisher's exact test. Overall survival (OS) curves were plotted according to the Kaplan-Meier method and compared using the log-rank test. Univariate and multivariate analyses were performed using the Cox proportional hazards model to identify independent variables predictive of OS. Statistical analyses were performed using JMP Pro 11 software (SAS Institute, Cary, NC, USA) and R version 3.2.0 (The R Foundation for Statistical Computing, Vienna, Austria). The differences were considered significant when the p-value was lower than 0.05.
Results
PSPH was identified as a potential driver gene in CRC. We identified PSPH, the gene for the rate-limiting enzyme in the serine biosynthesis pathway (SSP), as a potential driver gene on Ch.7p in CRC (Figure 1a). In this pathway, PSPH converts 3-phosphoserine to L-serine, which fuels glycine biosynthesis (10). Glycine is a major source of glutathione (GSH) and purine nucleotides (10). In the TCGA dataset, PSPH mRNA expression in CRC tumor tissues was 2.64-fold higher than in normal tissues (p<0.001) (Figure 1b). PSPH mRNA expression and PSPH copy number were positively correlated (R=0.51, p<0.001). Furthermore, copy number amplification of PSPH was observed in most of the CRC tissues in the TCGA dataset (Figure 1c).
PSPH was overexpressed in CRC. Next, we conducted qRT-PCR to measure the expression levels of PSPH mRNA in tumor tissues and paired normal colon tissues from 124 CRC patients from our Hospital. Consistent with the result from the TCGA dataset, PSPH mRNA expression was significantly higher in CRC tissues than in normal colon tissues (p<0.001) (Figure 2a). We also performed immunohistochemical staining to observe PSPH protein expression in CRC tissues and normal colon tissues. PSPH was strongly stained in tumor cells, whereas it was weak to moderate in normal colon cells (Figure 2b). Higher magnification of the image revealed cytoplasmic staining of PSPH in tumor cells.
Prognostic significance of PSPH mRNA expression in patients with CRC. We divided the 124 patients with CRC into a high PSPH mRNA expression group (n=57) and a low PSPH mRNA expression group (n=67) using the minimum p-value approach (11). The high PSPH mRNA expression group had poorer OS than the low expression group in our dataset (Figure 2c, left). Furthermore, we found that the high PSPH mRNA expression group had poorer prognosis than the low expression group in the TCGA dataset (Figure 2c, right). Multivariate analysis demonstrated that high PSPH mRNA expression was an independent prognostic factor in our dataset (5-year survival rate; low/high: 89.6/66.7 %, p<0.01) (Table I).
Clinicopathological characteristics of PSPH mRNA expression in patients with CRC. Clinicopathological analysis of our dataset revealed that high expression of PSPH mRNA was correlated with depth of invasion (p<0.05) and distant metastasis (p<0.05) (Table II).
Correlation between PSPH mRNA expression and gene sets associated with cell-cycle progression and apoptosis in CRC. GSEA revealed a positive correlation between PSPH mRNA expression and the expression of a gene set involved in cell-cycle progression (p<0.05) (Figure 3a, left). Negative correlation between PSPH mRNA expression and the expression of a gene set up-regulated in apoptosis progression was also detected (p<0.05) (Figure 3a, right).
Identification of candidate driver genes on chromosome 7p in CRC. (a) Schematic diagram of the strategy for candidate genes selection as described in Materials and Methods. (b) Boxplot of PSPH mRNA expression in 615 CRC tissues and 51 normal colon tissues in the TCGA dataset. (c) Correlation between PSPH copy number and PSPH mRNA expression in the TCGA dataset. R is the Pearson correlation coefficient. ***p<0.001.
Univariate and multivariate analyses of clinicopathological factors affecting 5-year overall survival.
Discussion
Here, we identified PSPH as a candidate of driver gene on Ch.7p in CRC. We demonstrated that overexpression of PSPH was associated with depth of invasion and distant metastasis and, also, that it was an independent poor prognostic factor in CRC patients. Moreover, overexpression of PSPH was likely caused by amplification of Ch.7p. These findings indicate that PSPH expression is a novel biomarker for poor prognosis and could play an important role in tumor progression of CRC.
Clinical significance of PSPH expression in CRC. (a) PSPH mRNA expression in 124 CRC tissues and paired normal colon tissues by qRT-PCR. (b) Immunohistochemical staining for PSPH in normal colon tissues (left) and CRC tissues (right). (c) Kaplan-Meier overall survival curves for patients with CRC according to PSPH mRNA expression in two independent datasets. (left; Kyushu University Beppu Hospital, n=124) (right; TCGA dataset, n=615). Original magnification; ×200. ***p<0.001.
As shown in Figure 3b, PSPH plays a crucial role in the production of L-serine in SSP (10). L-serine is converted to glycine, a major source of GSH and purine nucleotides (10). Importantly, glycine consumption and activation of mitochondrial glycine biosynthesis correlate with rapid proliferation across several types of cancer cell lines, while starvation of glycine inhibits proliferation by prolonging the G1 phase of the cell cycle (12). It is reported that the accumulation of GSH increases the ability to scavenge reactive oxygen species (ROS), which results in promoting proliferation of cancer cells and resistance against drug-induced apoptosis (13, 14). Consistent with these data, our GSEA with the TCGA dataset indicated that high PSPH expression is significantly associated with cell cycle progression and acquisition of an anti-apoptotic characteristic of CRC. Moreover, SSP enzymes, such as PSPH, are overexpressed in several types of cancer (15-18) and overexpression of SSP enzymes has been associated with poor prognosis (18-22). These data support our hypothesis that overexpression of PSPH due to amplification of Ch.7p enhances its function in serine synthesis and accelerates the production of glycine and GSH, which leads to an aggressive phenotype of tumor cells and poor prognosis in CRC patients.
GSEA for PSPH mRNA expression in the TCGA dataset. (a) GSEA using the TCGA dataset showed that PSPH mRNA expression was positively correlated with the expression of a gene set involved in cell cycle progression (left), whereas negatively correlated with a gene set up-regulated in apoptosis progression (right). (b) Schematic model indicating that overexpression of PSPH promotes proliferation and cancer cell survival. ES, Enrichment score; 3-PG, 3-phosphoglycerate; 3-PHP, 3-phosphohydroxy pyruvate; PHGDH, phosphoglycerate; PSAT1, phosphoserine aminotransferase-1; PSPH, phosphoserine phosphatase; SHMT, serine hydroxymethyltransferase; GSH, glutathione; ROS, reactive oxygen species.
Furthermore, knockdown of PSPH reduces GSH and impairs anti-oxidant activity of CRC cell lines, thereby enhancing the anticancer efficacy of 5-fluorouracil (5-FU) in vitro and in vivo (14). Of note, preclinical studies are ongoing for small molecules targeting SSP enzymes in the treatment of cancer (23). These reports strongly support our clinical findings that overexpression of PSPH is associated with tumor aggressiveness and poor prognosis in CRC.
In conclusion, we found that PSPH copy number amplification accompanies alterations in Ch.7p in CRC. Overexpression of PSPH is likely a ubiquitous driving event in all CRC cells. Thus, PSPH expression constitutes a novel prognostic biomarker for the malignant potential of CRC. PSPH may offer a therapeutic target to overcome ITH in CRC patients.
PSPH mRNA expression and clinicopathological factors in patients with CRC (n=124).
Acknowledgements
This research used the super-computing resource provided by the Human Genome Center, The Institute of Medical Science, The University of Tokyo (http://sc.hgc.jp/shirokane.html). We thank K. Oda, M. Kasagi, M. Sakuma, M. Oshiumi, N. Mishima and T. Kawano for their technical assistance. This work was supported in part by the following grants and foundations: Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research (grant number 16K10543, 16H01576, 16K07177, 16K10397, 16K19197 and 16K19107) and the OITA Cancer Research Foundation.
Footnotes
Disclosure
None of the Authors has a conflict of interest to disclose.
- Received March 15, 2017.
- Revision received March 28, 2017.
- Accepted March 31, 2017.
- Copyright© 2017, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved
References
In this issue
Jump to section
Related Articles
Cited By...
- Programmed Cell Death 10 (PDCD10) Is a Candidate Tumor-associated Gene in Esophageal Squamous Cell Carcinoma
- Exploring tumor-normal cross-talk with TranNet: role of the environment in tumor progression
- PRKRIP1, A Splicing Complex Factor, Is a Marker of Poor Prognosis in Colorectal Cancer
- CRMP5-associated GTPase (CRAG) Is a Candidate Driver Gene for Colorectal Cancer Carcinogenesis
- Cytolytic Activity (CYT) Score Is a Prognostic Biomarker Reflecting Host Immune Status in Hepatocellular Carcinoma (HCC)
- A Scoring System to Predict the Development of Bone Metastasis After Radical Resection of Colorectal Cancer