- Research article
- Open Access
A pseudogene-signature in glioma predicts survival
Journal of Experimental & Clinical Cancer Researchvolume 34, Article number: 23 (2015)
Pseudogene was recognized as a potential tumor suppressor or oncogene in varies of diseases, however its roles in glioma have not been investigated. Our study was to identify the pseudogene-signature that predicted glioma survival. Using a pseudogene-mining approach, we performed pseudogene expression profiling in 183 glioma samples from the Chinese Glioma Genome Atlas (CGGA) and set it as the training set. We found a six-pseudogene signature correlated with patients’ clinical outcome via bioinformatics analyses (P ≤ 0.01), and validated it in the Repository of Molecular Brain Neoplasia Data (REMBRANDT) containing 350 cases. A formula calculating the risk score based on the six-pseudogene signature was introduced and the patients of CGGA set were classified into high-risk group and low-risk group with remarkably different survival (P < 0.001) based on their scores. The prognostic value of the signature was confirmed in the REMBRANDT set. Though the function of these pseudogenes is not clear, the identification of the prognostic pseudogenes indicated the potential roles of pseudogenes in glioma pathogenesis and they may have clinical implications in treating glioma.
Gliomas are the most frequent primary tumors of the CNS (central nervous system) , half of which are represented by glioblastoma multiforme (GBM, WHO grade IV), associated with a poor prognosis (median survival less than one year). Identification of markers predicting the survials of gliomas is required for appropriate follow-up and treatment. Over the past few decades, varies of molecular markers were introduced in predicting survival including microRNAs , lncRNAs  and mutations of unique genes . Accordingly, oligodendrogliomas that show frequent 1p19q co-deletions and mutations of the IDH1 gene are associated with a longer survival than astrocytomas . However, the prognostic significance of pseudogenes in glioma has not been investigated.
In the present study, we identified a prognostic six-pseudogene signature from the CGGA patients set, and validated it in the RAMBRANDT set.
Materials and methods
Glioma data sets and corresponding clinical data were downloaded from the publicly available databases, including 183 cases from the Chinese Glioma Genome Atlas (CGGA, http://www.cgcg.org.cn/) and 350 cases from the Repository of Molecular Brain Neoplasia Data (REMBRANDT; http://caintegrator.nci.nih.gov/rembrandt/). Besides, 21 cases without tumors and 102 cases with insufficient data from the REMBRANDT were excluded. Psuedogene gene database was downloaded from the HGNC (HUGO Gene Nomenclature Committee, www.genenames.org), including 12440 pseudogenes, which were cross-matched with the CGGA and the REMBRANDT data set using Microsoft Excel (Microsoft Inc. Redmond, USA). Figure 1 depicts the flow chart of the study.
Pseudogenes represented in both the CGGA set and the REMBRANDT set were identified. Gene expression values represented by multiple probes were collapsed by taking the mean value of the probes [6,7]. The association between the pseudogene expression and patient’s overall survival was assessed by univariable Cox regression analysis along with a permutation test using Biometric Research Branch-Array (BRB-Array) Tools . Genes were considered statistically significant if their permutation P values were less than or equal to 0.01. A risk score formula was established by including each selected pseudogene weighting by its coefficient obtained from multivariable Cox regression analysis [9,10]. Kaplan-Meier survival analysis was used to estimate the survival distributions. The log-rank test was used to assess the statistical significance between stratified survival groups using the median value as the cutoff. A 2-sided P value <0.05 was regarded as significant.
Identification of prognostic pseudogenes from the CGGA data set
A total of 88 pseudogenes were identified, which were included in both the CGGA and REMBRANDT gene set. The 183 CGGA patients was set as the training set, which was used for the detection of the prognostic pseudogenes. By using the BRB-Array tools, the univariable Cox proportional hazards regression analysis was done on the pseudogenes expression data, and a six-pseudogene signature relating to overall survival was identified (p ≤ 0.01). Of these, a hazard ratio above 1 indicated that a high level of expression of a gene was associated with shorter survival (SP3P, ANXA2P3, PTTG3P, LPAL2, CLCA3P). A hazard ratio below 1 indicated that a high level of expression of a gene was associated with longer survival (TDH). The results were shown in Table 1.
The six-pseudogene signature correlates with patients’ survival
A risk-score formula was established according to the expressions of these selected pseudogenes and their coefficients, as follows: Risk Score = (0.097*expressing values of SP3P + 0.722*expressing values of PTTG3P + 0.149*expressing values of LPAL2 + (-0.185*expressing values of CLCA2P) + (-1.052*expressing values of TDH) + 0.656*expressing values of ANXA2P3). We then calculated the six-pseudogene signature risk score for each patient in the CGGA set, and ranked them due to their scores. The 183 patients were divided into a low-risk group and a high-risk group using the median risk score as the cutoff point. Patients in the low-risk group had remarkable longer overall survival time than those in the high-risk group (P < 0.001). Our result was shown in Figure 2.
Validation of the six-pseudogene signature in predicting survival in the REMBRANDT dataset
To confirm our findings, we validated our six-pseudogene signature in the REMBRANDT set. After calculating risk score of all the patients using the formula in the REMBRANDT set, the patients were assigned into a low-risk group and a high-risk group using the same cutoff point as for the CGGA set. Patients in the low-risk group had significant longer over survival time than those in the high-risk group (P < 0.001), which also maintained the consistence with the CGGA set.
In both of the CGGA set and the RAMBRANDT set, we found that the patients of the high-risk group tended to express high levels of risky pseudogenes (SP3P, ANXA2P3, PTTG3P, LPAL2, CLCA3P), whereas the patients of the low-risk group tended to express high level of protective pseudogenes (TDH). Results were shown in Figure 3.
The six-pseudogene signature was independent of age and gender
We stratified the entire CGGA patients (n = 183) into a younger group (age ≤ 50) or an elder one (age > 50). This analysis showed that within each age group, the six-pseudogene risk score could further subdivide the patients into those likely to have longer survival and those likely to have shorter survival. Then, the entire CGGA group was subdivided into 2 groups based on gender. Both in the male group and female group, the patients with lower risk scores had longer overall survival. These results were shown in Figure 4(A-D).
Besides, we made an analysis on relationship between the six-pseudogene signature and WHO Grade, which was showed in Figure 5. We can conclude that mean risk score ascended with the tumor malignance in both the CGGA set and the REMBRANDT set.
We also performed receriver operating characteristic (ROC) analysis to compare the sensitivity and specificity of survival prediction between the six-pseudogene model, MGMT expression status, age, gender, and IDH1 status. The area under receiver operating characteristic (AUROC) was determined and compared between these five prognostic factors. As shown in Figure 6, the AUROC of the six-pseudogene risk score was 0.807, which was significantly larger than that of age (AUROC = 0.619, P = 0.020) and gender (AUROC = 0.508, P = 0.881); when compared with the MGMT status and IDH1 status, the AUROC of the six-pseudogene risk score was larger than theirs’. (For IDH1, 0.807 versus 0.377, P = 0.049; for MGMT, 0.807 versus 0.533, P = 0.051). These results indicated that six-pseudogene signature might have a better predictive ability in predicting worse prognosis (Figure 6).
Pseudogenes are believed to be dysfunctional genes that have lost their protein-coding ability or are otherwise no longer expressed in the cell . However, some of them do play important roles in regulating their parent genes, and some even can be transcribed into RNA. It is also recognized that pseudogenes may regulate tumor suppressors and oncogenes. Pseudogenes were reported to correlate with varies of diseases, such as CYP4Z2P related to breast cancer , while it has not been investigated in glioma.
To verify our hypothesis, we performed pseudogene profiling in large cohorts of glioma patients from CGGA and REMBRANDT. By analyzing the association between gene expression profiling and clinical outcome of glioma patients, we identified a six-pseudogene signature significantly related to the overall survival of glioma patients.
Functional characteristics of the six pseudogenes
As for the functional characteristics of the six pseudogenes, one pseudogene called TDH, or L-threonine dehydrogenase, was believed to be protective, whereas the other five were risky ones.
TDH is a pseudogene of L-threonine 3-dehydrogenase (TDH), whose expression was found in human tissues including heart, brain, placenta, lung, liver, skeletal muscle, kidney, pancreas, spleen, thymus, prostate, testis, ovary and small intestine. TDH mRNA was present in most cell types examined, but was below the level of detection in endothelial cells, glioma cell lines and some leukaemia cell lines according to the study of Alasdair J Edgar . Though TDH is an expressed pseudogene, its function is not clear. So far, we can’t find a single research on TDH. However, in our study, the high level expression of TDH related with longer overall survival, which indicated TDH might be a protect pseudogene. As the expression of TDH mRNA was low in glioma cell line , we conjectured that TDH might serve as an endo-generate siRNA, or work as a molecular sponge combining microRNAs. The mechanism is subject to be investigated.
The present study demonstrated the associations between the high expressions of these five pseudogenes and shortened overall survival. The roles of these pseudogenes in glioma pathogenesis are presently unclear, and our findings suggest that they deserve further studies. We might conjecture some information from their parent genes that were studied before. The parent genes of these pseudogenes (SP3P, ANXA2P3, PTTG3P, LPAL2, CLCA3P) are SP3 transcription factor (Sp3), Annexin A2 (ANXA2), pituitary tumor-transforming 1(PTTG1), Lp(a)(LPA lipoprotein) and chloride channel accessory 3 (CLCA3) respectively.
ANXA2 was found expressed significant higher in NSCLC (Non-small cell lung cancer) tissue compared to that in adjacent non-cancerous tissue according to the study of Jia et al. , and Zhang et al. found annexin A2 silencing inhibits invasion, migration, and tumorigenic potential of hepatoma cells , which indicated that annexin A2 (ANXA2) might serve as an important mediator of malignant transformation and development of hepatocellular carcinoma. The pseudogene of ANXA2 may paly the same role in glioma.
PTTG is an oncogene that plays diverse roles in the occurrence, proliferation, and invasion of a variety of tumors, and PTTG has already served as one of the markers of proliferative activity progress in many tumors . Study were carried out on the association between PTTG expression level and human pituitary macroadenomas , which suggested PTTG may promote invasive tumor growth by stimulating pituitary adenomas proliferation. As the pseudogene of PTTG, PTTG3P may play as an oncogene in glioma, which played as a risky gene.
CLCA may serve as a tumor suppressor, which cound be inferred from the result that hCLCA2 is expressed in normal breast epithelium but not in 29 breast cancer lesions of different stages [17,18]. Our research showed CLCA3P was a risky pseudogene. We conjectured that CLCA3P might be a molecular sponge same as PTENP1. However, the mechanism is subject to be investigated.
The encoded protein of Lp(a) constitutes a substantial portion of lipoprotein(a) and is proteolytically cleaved, resulting in fragments that attach to atherosclerotic lesions and promote thrombogenesis [19,20] which seemed to have no relationship with tumor.
Limitations of the study
Our study does have some small limitations. In our study, the data sets were not designed for detecting pseudogenes, so, only a small minority of human pseudogenes (88 genes) were included in the analysis. Because of the different platform using in the CGGA and REMBRANDT, we only retained the pseudogenes that appeared on both platforms for survival analysis. So, these prognostic pseudogenes identified here might not represent all the pseudogenes correlated with glioma overall survival. Besides, the mechanisms of the pseudogenes are subject to be investigated on which we intended to study in our later research. However, the significant and consistent correlation of our six-pseudogene signature with overall survival in several independent data sets indicates that it is a potentially powerful prognostic marker for glioma.
In conclusion, by employing two independent patients cohorts, our study revealed the prognostic values of pseudogenes in glioma for the first time. Our findings strongly prompt that pseudogene signatures may be of use in predicting the treating outcome and may be novel biomarkers in glioma prognoses. As the updating of common databases, pseudogenes will show more value in future glioma studies, while this is what we continue to explore in the future. We will focus on the function of the pseudogenes we found and try to validate our findings in our clinical outcome.
Adamson C, Kanu OO, Mehta AI, Di C, Lin N, Mattox AK, et al. Glioblastoma multiforme: a review of where we have been and where we are going. Expert Opin Investig Drugs. 2009;18(8):1061–83.
Li R, Gao K, Luo H, Luo H, Wang X, Shi Y, et al. Identification of intrinsic subtype-specific prognostic microRNAs in primary glioblastoma. J Exp Clin Cancer Res. 2014;33:9.
Zhang XQ, Sun S, Lam KF, Kiang KMY, Pu JKS, Ho ASW, et al. A long non-coding RNA signature in glioblastoma multiforme predicts survival. Neurobiol Dis. 2013;58:123–31.
Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson M, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010;17(1):98–110.
Ricard D, Idbaih A, Ducray F, Lahutte M, Hoang-Xuan K, Delattre JY. Primary brain tumours in adults. Lancet. 2012;379(9830):1984–96.
Liu C, Sage JC, Miller MR, Verhaak RG, Hippenmeyer S, Vogel H, et al. Mosaic analysis with double markers reveals tumor cell of origin in glioma. Cell. 2011;146(2):209–21.
Murat A, Migliavacca E, Gorlia T, Lambiv WL, Shay T, Hamou MF, et al. Stem cell-related “self-renewal” signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma. J Clin Oncol. 2008;26(18):3015–24.
Simon R, Lam A, Li MC, Ngan M, Menenzes S, Zhao Y. Analysis of gene expression data using BRB-ArrayTools. Cancer Informat. 2007;3:11–7.
Lossos IS, Czerwinski DK, Alizadeh AA, Wechser MA, Tibshirani R, Botstein D, et al. Prediction of survival in diffuse large B-cell lymphoma based on the expression of 2 genes reflecting tumor and microenvironment. Blood. 2011;118(5):1350–8.
Lossos IS, Czerwinski DK, Alizadeh AA, Wechser MA, Tibshirani R, Botstein D, et al. Prediction of survival in diffuse large-B-cell lymphoma based on the expression of six genes. N Engl J Med. 2004;350(18):1828–37.
Lossos IS, Czerwinski DK, Alizadeh AA, Wechser MA, Tibshirani R, Botstein D, et al. Identification of a novel mammary-restricted cytochrome P450, CYP4Z1, with overexpression in breast carcinoma. Cancer Res. 2004;64(7):2357–64.
Edgar AJ. The human L-threonine 3-dehydrogenase gene is an expressed pseudogene. BMC Genet. 2002;3:18.
Jia JW, Li KL, Wu JX, Guo SL. Clinical significance of annexin II expression in human non-small cell lung cancer. Tumour Biol. 2013;34(3):1767–71.
Zhang HJ, Yao DF, Yao M, Huang H, Wang L, Yan MJ, et al. Annexin A2 silencing inhibits invasion, migration, and tumorigenic potential of hepatoma cells. World J Gastroenterol. 2013;19(24):3792–801.
Demeure MJ, Coan KE, Grant CS, Komorowski RA, Stephan E, Sinari S, et al. PTTG1 overexpression in adrenocortical cancer is associated with poor survival and represents a potential therapeutic target. Surgery. 2013;154(6):1405–16. discussion 16.
Jia W, Lu R, Jia G, Ni M, Xu Z. Expression of pituitary tumor transforming gene (PTTG) in human pituitary macroadenomas. Tumour Biol. 2013;34(3):1559–67.
Sasaki Y, Koyama R, Maruyama R, Hirano T, Tamura M, Sugisaka J, et al. CLCA2, a target of the p53 family, negatively regulates cancer cell migration and invasion. Cancer Biol Ther. 2012;13(14):1512–21.
Walia V, Yu Y, Cao D, Sun M, McLean JR, Hollier BG, et al. Loss of breast epithelial marker hCLCA2 promotes epithelial-to-mesenchymal transition and indicates higher risk of metastasis. Oncogene. 2012;31(17):2237–46.
Nestel PJ, Barnes EH, Tonkin AM, Simes J, Fournier M, White HD, et al. Plasma lipoprotein(a) concentration predicts future coronary and cardiovascular events in patients with stable coronary heart disease. Arterioscler, Thromb, Vasc Biol. 2013;33(12):2902–8.
Anderson JL, Knight S, May HT, Horne BD, Bair TL, Huntinghouse JA, et al. Validation and quantification of genetic determinants of lipoprotein-a levels and predictive value for angiographic coronary artery disease. Am J Cardiol. 2013;112(6):799–804.
Analyses were performed using BRB-ArrayTools developed by Dr. Richard Simon and BRB-ArrayTools Development Team.
This work was supported by grants from the National High Technology Research and Development Program of China (863) (2012AA02A508), the Research Special Fund For Public Welfare Industry of Health (201402008), National Natural Science Foundation of China (91229121江涛, 81272792尤, 81472362尤, 81172389刘, 81372709刘, 81302185张, 81302184王协锋), Jiangsu Province’s Natural Science Foundation (20131019), Jiangsu Province’s Key Provincial Talents Program (RC2011051), Jiangsu Province’s Key Discipline of Medicine (XK201117), Jiangsu Provincial Special Program of Medical Science (BL2012028), and Program for Development of Innovative Research Team in the First Affiliated Hospital of NJMU, and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).
The authors declare that they have no competing interests.
GKM, CXC performed the experiments. GKM, ZJX, WYY, YW analyzed the data. GKM wrote the manuscript. YYP, GKM conceived and designed the experiments. All authors read and approved the final manuscript.
Kai-Ming Gao, Xin-cheng Chen and Jun-xia Zhang contributed equally to this work.