Identification of genes and pathways potentially related to PHF20 by gene expression profile analysis of glioblastoma U87 cell line

Cancer Cell International, Oct 2017

Glioblastoma is the most common and aggressive brain tumor associated with a poor prognosis. Plant homeodomain finger protein 20 (PHF20) is highly expressed in primary human gliomas and its expression is associated with tumor grade. However, the molecular mechanism by which PHF20 regulates glioblastoma remains poorly understood. Genome wide gene expression analysis was performed to identify differentially expressed genes (DEGs) in U87 cells with PHF20 gene knockdown. Gene ontology (GO) and pathway enrichment analyses were performed to investigate the functions and pathways of DEGs. Pathway-net and signal-net analyses were conducted to identify the key genes and pathways related to PHF20. Expression of 540 genes, including FEN1 and CCL3, were significantly altered upon PHF20 gene silencing. GO analysis results showed that DEGs were significantly enriched in small molecule metabolic and apoptotic processes. Pathway analysis indicated that DEGs were mainly involved in cancer and metabolic pathways. The MAPK, apoptosis and p53 signaling pathways were identified as the hub pathways in the pathway network, while PLCB1, NRAS and PIK3 s were hub genes in the signaling network. Our findings indicated that PHF20 is a pivotal upstream regulator. It affects the occurrence and development of glioma by regulating a series of tumor-related genes, such as FEN1, CCL3, PLCB1, NRAS and PIK3s, and activation of apoptosis signaling pathways. Therefore, PHF20 might be a novel biomarker for early diagnosis, and a potential target for glioblastoma therapies.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://cancerci.biomedcentral.com/track/pdf/10.1186/s12935-017-0459-x?site=cancerci.biomedcentral.com

Identification of genes and pathways potentially related to PHF20 by gene expression profile analysis of glioblastoma U87 cell line

Liu et al. Cancer Cell Int Identification of genes and pathways potentially related to PHF20 by gene expression profile analysis of glioblastoma U87 cell line Tianlong Liu 0 2 Tiejun Zhang Feng Zhou Jitao Wang 1 Xiaohu Zhai 0 2 Nan Mu Jongsun Park Minna Liu Wenxing Liu 0 2 Peijin Shang 0 2 Yi Ding 0 2 Aidong Wen 0 2 Yuwen Li 0 1 2 0 Department of Pharmacy, Xijing Hospital, Fourth Military Medical University , Xi'an , China 1 Department of Pharmacy, The First Affiliated Hospital 2 Department of Pharmacy, Xijing Hospital, Fourth Military Medical University , Xi'an , China Background: Glioblastoma is the most common and aggressive brain tumor associated with a poor prognosis. Plant homeodomain finger protein 20 (PHF20) is highly expressed in primary human gliomas and its expression is associated with tumor grade. However, the molecular mechanism by which PHF20 regulates glioblastoma remains poorly understood. Methods: Genome wide gene expression analysis was performed to identify differentially expressed genes (DEGs) in U87 cells with PHF20 gene knockdown. Gene ontology (GO) and pathway enrichment analyses were performed to investigate the functions and pathways of DEGs. Pathway-net and signal-net analyses were conducted to identify the key genes and pathways related to PHF20. Results: Expression of 540 genes, including FEN1 and CCL3, were significantly altered upon PHF20 gene silencing. GO analysis results showed that DEGs were significantly enriched in small molecule metabolic and apoptotic processes. Pathway analysis indicated that DEGs were mainly involved in cancer and metabolic pathways. The MAPK, apoptosis and p53 signaling pathways were identified as the hub pathways in the pathway network, while PLCB1, NRAS and PIK3 s were hub genes in the signaling network. Conclusions: Our findings indicated that PHF20 is a pivotal upstream regulator. It affects the occurrence and development of glioma by regulating a series of tumor-related genes, such as FEN1, CCL3, PLCB1, NRAS and PIK3s, and activation of apoptosis signaling pathways. Therefore, PHF20 might be a novel biomarker for early diagnosis, and a potential target for glioblastoma therapies. PHF20; Glioblastoma; U87 cell; Gene expression profile; Bioinformatics Background Glioblastoma is the most common and lethal tumor of the central nervous system [ 1 ], owing to poor prognosis and repercussions on cognitive function [ 2 ]. Despite advances in knowledge and therapies over several decades, survival has not significantly improved, only 5.1% of patients with glioblastoma have a 5-year survival rate [ 3 ]. Thus, understanding the mechanisms that regulate glioblastoma progression is critical to developing novel therapies to improve patient outcome. One particular protein of interest in glioblastoma regulation is plant homeodomain finger protein 20 (PHF20). PHF20 is a potent transcriptional activator, which binds to methylated lysine residues on the histone tail [ 4 ]. PHF20 is overexpressed in various cancer tissues compared to adjunct normal tissues, including advanced small-cell lung cancers and advanced adenocarcinomas [ 5 ]. Besides, PHF20 is highly expressed in primary human glioma specimens [ 6 ], and functions as an immunogenic antigen in glioblastoma [ 7, 8 ]. Auto-antibodies against PHF-20 were also detected in hepatocellular carcinoma [9] and meduloblastoma [ 10 ]. PHF20 expression levels have also been associated with the pathological tumor grade of gliomas [ 6 ]. To elucidate the mechanisms regulated by PHF20 in glioma as well as identify potential prognostic biomarkers and targets for drug discovery and immunotherapy, a microarray analysis was conducted to harness the systematic gene expression profile related to genomic and phenotypic information on glioblastoma in U87 cells. Methods Cell culture Human glioblastoma cell lines U87, U251 and A172 originated from the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China). Cell lines LN229, HS683 and HEB were kindly provided by the department of  neurosurgery at  The First Affiliated Hospital of SooChow University. The cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) (Corning, NY, USA) containing 10% fetal bovine serum (FBS), 50 U/mL penicillin and 50 μg/mL streptomycin at 37 °C with 5% CO2 incubator. The cell lines tested negative for any mycoplasma contamination. Western blotting 1 × 106 cultured cells were lysed with lysis buffer as previously described [ 5 ]. Protein concentration was measured using the BCA protein assay kit (Beyotime, Shanghai, China). The same amount of protein was separated by 10% sodium dodecyl sulfate–polyacrylamide (SDS-PAGE). A polyvinylidene difluoride (PVDF) membrane (Millipore, Bedford, MA, USA) was then used for electro-transfer. The membrane was blocked with 5% nonfat milk at room temperature for 1 h and incubated in primary antibodies against PHF20 (1:500, Cell Signaling Technology, USA), overnight at 4 °C. Subsequently, the membrane was incubated in the appropriate secondary antibody at room temperature for 1 h. In addition, β-actin was used as the loading control. Protein bands were visualized through enhanced chemiluminescence (ECL) reagent and detected using BioImaging Systems (UVP, Upland, CA, USA). The relative protein levels were calculated with Image J software (National Institutes of Health, USA). All experiments were performed in triplicate. Lentivirus‑based shRNA infection GFP-Lentiviral particles with PHF20-specific shRNA (shPHF20) were purchased from Genechem Co., Ltd. (Shanghai, China). The target sequence was TGACT TGGTTGTATCAGAT. Random sequence, TTCTCCG AACGTGTCACGT, was used as a negative control (shCON). U87 cells in 6-well plates were infected with lentiviral particles containing either shCON or shPHF20 to generate negative control (NC) or PHF20 knockdown (KD) U87 cells, respectively. 12  h after infection, the virus containing culture medium was replaced with fresh DMEM supplemented with 10% FBS for 72  h. The lentiviral infection efficiency was demonstrated by observing the presence of green fluorescent protein within the U87 cells using Olympus-IX71 fluorescence microscope (Tokyo, Japan) and RT-PCR assay. RNA extraction and quantitation Total RNA was isolated using Trizol Reagent (Pufei, Shanghai, China) according to the manufacturer’s protocol. The RNA content was examined by identifying A260 and A280 values by using the Nanodrop 2000 (Thremo Scientific, Waltham, MA, USA). RNA integrity was assessed using a 2100 Bioanalyzer (Agilent Technologies) and an RNA 6000 Nano Kit (Agilent Technologies). RNA with A260/A280  nm values over than 1.9, concentrations over 300 ng/μL and 28S/18S ratios over than 1.4 were used. Quantitative real time PCR analysis Total RNA isolated was processed for cDNA synthesis using M-MLV reverse transcriptase (Promega Corporation, Madison, WI, USA). cDNA was amplified by PCR in StratageneMX3000p (Agilent Technologies, Santa Clara, CA, USA) using SYBR Master Mixture (TaKaRa, Tokyo, Japan). The expression levels of target genes were standardized against the GAPDH, an internal control, and calculated using the 2−△△Ct method. The sequences of the primers used in PCRs are listed in Additional file  1. All the assays were performed in triplicate. mRNA microarray Total RNA was processed for double-strand cDNA synthesis, IVT and amplified RNA fragmentation using the GeneChip 3′IVT Express Kit (Affymetrix, Santa Clara, CA, USA) according to the manufacturer’s instruction. RNA was then processed for hybridization at 45  °C for 17  h using The PrimeView™ Human Gene Expression Array (Affymetrix), which contains 49,395 probes covering more than 36,000 transcripts and variants. The arrays were washed in the GeneChip Fluidic Station 450 (Affymetrix), and scanned by the GeneChip Scanner 3000 (Affymetrix). These microarray data have been deposited in NCBI Gene Expression Omnibus (GEO) under accession number GSE93680. Data processing The raw data, expressed as CEL files, were normalized by the log scale robust multi-array analysis (RMA) method with the Expression Console software version 1.1 (Affymetrix). The screening standard for a distinctly significant gene was an absolute fold change (|FC|)  >  2 and a corrected p < 0.05. Gene ontology and pathway analysis Gene ontology (GO) analysis was applied to analyze the main function of differentially expression genes (DEGs) according to the gene ontology, the key functional classification of National Center of Biotechnology Information (NCBI) [ 11 ]. Two-side Fisher’s exact test and χ2 tests were used to classify the GO category. The false discovery rate (FDR) [ 12 ] was calculated to correct the p value. The standard of difference screening was FDR < 0.05. Pathway analysis was used to find out the significant pathway of the DEGs according to Kyoto Encyclopedia of Genes and Genomes (KEGG) [ 13 ]. The data analysis method and filter criteria were similar to the GO analysis. Pathway‑net analysis Pathway-net analysis was built according to the interaction among pathways of the KEGG database to directly and systemically determine the interaction among the significant pathways [ 14 ]. Signaling processes analysis Based on the KEGG pathway map (http://www.genome. jp/kegg/pathway.html) [ 15 ] ,DEGs involved in key pathways were labeled to clearly visualize the position of specific genes in the signaling processes and determine the regulatory role of DEGs involved in key pathways. Signal‑net analysis A gene–gene interaction network was constructed using the source of the interaction database from KEGG. For instance, if there is confirmative evidence that two genes interact with each other, an interaction edge is assigned between the two genes. The networks are stored and presented as graphs, where nodes represent main genes (protein, compound, etc.) and edges represent the relationship between the nodes, such as activation or phosphorylation. The algorithms and construction of the network were achieved using published methods [ 16 ]. Results PHF20 is highly expressed in glioma cell lines We first examined the expression of PHF20 in various glioma cell lines. Expression of PHF20 protein was significantly higher in glioma cell lines (A172, LN229, U251, HS683 and U87) than that in human astrocyte cell line HEB (p < 0.05, Fig. 1a). The relative expression level of PHF20 in U87 cells was higher than in other glioma cell lines (p < 0.05, Fig. 1b). PHF20 mRNA level was verified by qPCR (Fig. 1c). Thus, the U87 cell line was used to establish PHF20 knockdown cells in following studies. The U87 cells were successfully infected with shCON and shPHF20 72  h after transfection, resulting in a 72% decrease in PHF20 expression (Fig. 1d). Identification of DEGs that regulated by PHF20 A genome wide gene expression analysis was carried out to identify DEGs between NC and KD U87 cells. A total of 540 DEGs were identified, including 175 up-regulated genes and 365 down-regulated genes (Fig. 2a, Additional file 2). A subset of DEGs was verified by qPCR. Expression of BBOF1, FBXO36 and SPARC increased, while the expression of TPM4, FEN1, AGPS, BCAT1 and CCL3 decreased in PHF20 knockdown cells (Fig. 2b). GO analysis of PHF20 associated DEGs To determine the primary functions regulated by PHF20 in glioma cells, a comprehensive gene ontology (GO) analysis was performed. A total of 540 DEGs were assigned to 236 GO terms (111 were up-regulated and 125 were down-regulated, Additional file  3). The analysis revealed that the up-regulated genes were significantly involved in homophilic cell adhesion, protein transport, and ER to Golgi vesicle-mediated transport. Down-regulated genes were involved in small molecule metabolic process, transcription (DNA-dependent), signal and transduction (Fig.  3). Additionally, both up and downregulated genes were enriched in small molecule metabolic, transcription and apoptotic processes. KEGG pathway analysis of PHF20 associated DEGs Pathway enrichment analysis of DEGs was conducted on the basis of the KEGG pathway database. This analysis yielded 147 significant pathways, including 41 up-regulated pathways and 106 down-regulated pathways (see Additional file 4). The most significant up-regulated and down-regulated pathways were shown in Fig. 4. In addition, pathways in cancer were identified as up-regulated pathways as well as down-regulated pathways. Pathway‑net analysis of significant PHF20‑regulated pathways In order to define functional relationships among pathways, an interaction net of the significant pathways associated with PHF20 was built (Fig.  5). 46 key pathways Fig. 2 Differentially expressed genes (DEGs) between the U87 cells infected with the negative control (NC) and those infected with shPHF20 (KD). a Hierarchical clustering for DEGs. Green represents down-regulated genes, red represent up-regulated genes (p<0.05). b A subset of genes differentially expressed between NC and KD cells were validated by qPCR. White bars represent the fold change in expression level between KD and NC as indicated by microarray analysis. Black bars represent the mean fold change of gene expression calculated by qPCR method. GAPDH was used for normalization. Each grey bar is the mean ± SEM of three independent biological replicates and the 147 connections between them were represented by nodes and edges, respectively. Primary interactions occurred between the MAPK, apoptosis, cancer, p53, ErbB, cytokine–cytokine receptor interaction, focal adhesion, and JAK-STAT signaling pathways (Table 1). Signaling processes analysis of key pathways Signaling processes analysis, with the KEGG pathway map, was performed using DEG expression data to determine the regulatory role of DEGs involving in key pathways. As shown in Fig.  6, the MAPK signaling pathway included 7 up-regulated and 8 down-regulated genes while pathways in cancer contained 13 up-regulated and 20 down-regulated PHF20-regulated DEGs. Signal‑net analysis of DEGs Finally, we identified the key gene interactions between PHF20-related DEGs to construct a regulatory network map.  78 genes were included in the signaling network, and 106 potential direct interactions were identified (Fig. 7). PLCB1, NRAS, PIK3CD, PIK3CA, PIK3R1, HDAC4, HDAC8, CRKL, RAB7A and ITGB3 were the most significantly expressed genes according to the degree size (Table 2). Discussion PHF20 was originally identified in glioma patients [ 8 ] and is significantly associated with glioma pathological tumor grade [ 6 ]. In recent years, a growing number of studies have shown that PHF20 is closely related to the development of various tumors [ 17, 18 ] and plays important roles in tumor suppression and progression. However, the underlying molecular pathways regulated by PHF20 in glioma remain largely undetermined. Therefore, further in-depth investigations are essential for better understanding of the biological roles of PHF20 in cancer. In the present study, gene expression profile analysis was performed to identify differentially expressed genes (DEGs) between PHF20 knockdown U87 cells and negative control cells. A total of 540 genes (175 up-regulated genes and 365 down-regulated genes) were differentially expressed following knockdown of PHF20, which suggests that PHF20 may be a key regulator in glioblastoma. Multiple DEGs, including FEN1, BCAT1, AGPS and CCL3, have been implicated in the progression of various cancers. For example, FEN1 is overexpressed in glioblastoma [ 19 ]. FEN1 polymorphisms and variant genotypes are associated with glioma susceptibility [ 20, 21 ]. CCL3 is also highly expressed in glioma, and may promote glioblastoma cell proliferation and migration [22]. Gene ontology enrichment analysis revealed that highly enriched biological functions were related to PHF20, such as hemophilic cell adhesion, protein transport, metabolic process, transcription and apoptotic process. Thus, PHF20 may influence glioma progression by altering these biological processes. Pathway name Several DEG enriched pathways associated with tumorigensis were identified including protein processing in endoplasmic reticulum, metabolic pathways, ubiquitin mediated proteolysis, pathways in cancer, and thyroid hormone signaling pathways. Furthermore, pathwaynet analysis revealed that multiple pathways participate in the occurrence and development of cancer including the p53 signaling pathway, apoptosis, pathways in cancer, and the TLR signaling pathway. The p53 signaling pathway was also enriched as a significant pathway by an array comparative genomic hybridization analysis in pilocytic astrocytoma [ 23 ]. Furthermore, our findings are in line with previous studies that found that PHF20 could stabilize and activate p53 by promoting p53 methylation [ 24 ], and that PHF20 inhibits p53 transcriptional activity via PKB mediated PHF20 phosphorylation [ 5 ]. A recent study also showed that PHF20 inhibits tumorigenicity by inducing apoptosis mediated by p53 and Bax [ 17 ]. Moreover, accumulating evidences suggested that PHF20 was expressed in a number of tumors, including glioma [ 6 ], lung cancer [ 25 ] and myeloid malignancies [ 26 ]. In addition, elevated expression of PHF20 could cause constitutive NF-B activation [ 6 ], which is a key downstream gene of TLR signaling pathway [ 27 ]. Finally, signal-net analysis revealed the interactions between 78 PHF20-regulated genes. Core genes PLCB1, PIK3CD/CA/R1, CRKL, RAB7A and ITGB3 were downregulated while NRAS and HDAC4/8 were up-regulated. PLCB1 plays critical roles in intracellular transduction and regulating signal activation [ 28 ], which are important to tumorigenesis. As one of the RAS oncogene family, NRAS have been reported to be involved in development of leukemia [ 29 ], melanoma [ 30 ] and glioma [ 31, 32 ]. Members of the PIK3 family are frequently detected in a wide range of cancers and have been proposed as biomarkers for patient survival and drug response [ 33, 34 ]. PHF20 has been suggested as a substrate of PKB [ 5 ]. HDACs regulate various nuclear and cytoplasmic processes [ 35 ], which are common in various human neoplasms [ 36, 37 ]. In addition, synergistic anti-tumor actions between HDAC and PIK3 inhibitors have been validated [38]. Conclusions Overall, this study indicated that PHF20 is a pivotal upstream gene that influences the occurrence and development of glioma by regulating a series of tumorrelated genes, like FEN1, CCL3, PLCB1, NRAS and PIK3s, and involved in apoptosis signaling pathways. Thus, PHF20 might be a novel biomarker for early diagnosis and therapeutic target for treatment of glioma. Nevertheless, further studies in molecular pathogenesis and large scale clinical tumor specimen validation are still needed. Additional file 2. The differentially expressed genes (DEGs) between common U87 cells and PHF20 knockdown U87 cells. Additional file 3. Significant gene ontology (GO) analysis of differentially expressed genes (DEGs) related to PHF20. Additional file 4. Significant pathway analysis of differentially expressed genes (DEGs) related to PHF20. Abbreviations DEGs: differentially expressed genes; FDR: false discovery rate; GLEA2: gliomaexpressed antigen 2; GO: gene ontology; HCA58: hepatocellular carcinoma associated antigen 58; KEGG: Kyoto Encyclopedia of Genes and Genomes; PHF20: plant homeodomain finger protein 20. Authors’ contributions LYW and WAD designed and directed the study. DY performed the statistical analysis. LTL collected background information, drafted the manuscript and performed the microarray assay. ZTJ and ZF conducted the networks. ZXH performed the cell culture. MN performed the western blotting. JP participated in study design and helped to draft the manuscript. LMN finished the experiments of RNA extraction. LWX performed the lentivirus infection. SPJ and WJT finished the real-time PCR array. All authors read and approved the final manuscript. of SooChow University, Suzhou, China. 3 State Key Laboratory of Cancer Biology, Department of Biopharmaceutics, School of Pharmacy, Fourth Military Medical University, Xi’an, China. 4 Department of Pharmacology, Chungnam National University, Daejon, South Korea. 5 Department of Nephrology, Xijing Hospital, Fourth Military Medical University, Xi’an, China. 6 Department of Neurosurgery, The First Affiliated Hospital of SooChow University, Suzhou, China. Acknowledgements Not applicable. Competing interests The authors declare that they have no competing interests. Availability of data and materials The datasets generated and analyzed during the current study are available in the Gene Expression Omnibus repository (https://www.ncbi.nlm.nih.gov/geo/ query/acc.cgi?acc=GSE93680). Consent for publication Not applicable. Ethics approval and consent to participate Not applicable. Funding The work has been supported by National Natural Science Foundation of China (Nos. 81501003, 81673631, 81201985) and China Postdoctoral Science Foundation (2015M580465). Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 1. Patrick YW , Santosh K. Malignant gliomas in adults . N Engl J Med . 2008 ; 359 : 492 - 507 . 2. Omuro A , DeAngelis LM . Glioblastoma and other malignant gliomas: a clinical review . JAMA . 2013 ; 310 ( 17 ): 1842 - 50 . 3. Ostrom QT , Gittleman H , Fulop J , Liu M , Blanda R , Kromer C , Wolinsky Y , Kruchko C , Barnholtz-Sloan JS . CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2008-2012 . Neuro-Oncology. 2015 ; 17 ( suppl 4 ): v1 - 62 . 4. Badeaux AI , Yang Y , Cardenas K , Vemulapalli V , Chen K , Kusewitt D , Richie E , Li W , Bedford MT . Loss of the methyl lysine effector protein PHF20 impacts the expression of genes regulated by the lysine acetyltransferase MOF . J Biol Chem . 2011 ; 287 ( 1 ): 429 - 37 . 5. Li Y , Park J , Piao L , Kong G , Kim Y , Park KA , Zhang T , Hong J , Hur GM , Seok JH , et al. PKB-mediated PHF20 phosphorylation on Ser291 is required for p53 function in DNA damage . Cell Signal . 2013 ; 25 ( 1 ): 74 - 84 . 6. Zhang T , Park KA , Li Y , Byun HS , Jeon J , Lee Y , Hong JH , Kim JM , Huang SM , Choi SW , et al. PHF20 regulates NF-kappaB signalling by disrupting recruitment of PP2A to p65 . Nat Commun . 2013 ;4: 2062 . 7. Pallasch CP , Struss A , Munnia A , König J , Steudel W , Fischer U , Meese E. Autoantibodies against GLEA2 and PHF3 in glioblastoma: tumor-associated autoantibodies correlated with prolonged survival . Int J Cancer . 2005 ; 117 ( 3 ): 456 - 9 . 8. Fischer U , Struss AK , Hemmer D , Pallasch CP , Steudel WI , Meese E. Gliomaexpressed antigen 2 (GLEA2): a novel protein that can elicit immune responses in glioblastoma patients and some controls . Clin Exp Immunol . 2001 ; 126 ( 2 ): 206 - 13 . 9. Wang Y , Han KJ , Pang XW , Vaughan HA , Qu W , Dong XY , Peng JR , Zhao HT , Rui JA , Leng XS , et al. Large scale identification of human hepatocellular carcinoma-associated antigens by autoantibodies . J Immunol . 2002 ; 169 ( 2 ): 1102 - 9 . 10. Behrends U , Schneider I , Rossler S , Frauenknecht H , Golbeck A , Lechner B , Eigenstetter G , Zobywalski C , Muller-Weihrich S , Graubner U , et al. Novel tumor antigens identified by autologous antibody screening of childhood medulloblastoma cDNA libraries . Int J Cancer . 2003 ; 106 ( 2 ): 244 - 51 . 11. Ashburner M , Ball CA , Blake JA , Botstein D , Butler H , Cherry JM , Davis AP , Dolinski K , Dwight SS , Eppig JT , et al. Gene ontology: tool for the unification of biology . Nat Genet . 2000 ; 25 ( 1 ): 25 - 9 . 12. Dupuy D , Bertin N , Hidalgo CA , Venkatesan K , Tu D , Lee D , Rosenberg J , Svrzikapa N , Blanc A , Carnec A , et al. Genome-scale analysis of in vivo spatiotemporal promoter activity in Caenorhabditis elegans . Nat Biotechnol . 2007 ; 25 ( 6 ): 663 - 8 . 13. Draghici S , Khatri P , Tarca AL , Amin K , Done A , Voichita C , Georgescu C , Romero R . A systems biology approach for pathway level analysis . Genome Res . 2007 ; 17 ( 10 ): 1537 - 45 . 14. Yi M , Horton JD , Cohen JC , Hobbs HH , Stephens RM . Whole pathway scope: a comprehensive pathway-based analysis tool for high-throughput data . BMC Bioinform . 2006 ; 7 : 30 . 15. Kanehisa M , Goto S , Sato Y , Kawashima M , Furumichi M , Tanabe M . Data, information, knowledge and principle: back to metabolism in KEGG . Nucleic Acids Res . 2014 ; 42 ( D1 ): D199 - 205 . 16. Wei Z , Li H . A Markov random field model for network-based analysis of genomic data . Bioinformatics . 2007 ; 23 ( 12 ): 1537 - 44 . 17. Tang N , Ma L , Lin XY , Zhang Y , Yang DL , Wang EH , Qiu XS . Expression of PHF20 protein contributes to good prognosis of NSCLC and is associated with Bax expression . Int J Clin Exp Pathol . 2015 ; 8 ( 10 ): 12198 - 206 . 18. Taniwaki M , Daigo Y , Ishikawa N , Takano A , Tsunoda T , Yasui W , Inai K , Kohno N , Nakamura Y . Gene expression profiles of small-cell lung cancers: molecular signatures of lung cancer . Int J Oncol . 2006 ; 29 ( 3 ): 567 - 75 . 19. Nikolova T , Christmann M , Kaina B. FEN1 is overexpressed in testis, lung and brain tumors . Anticancer Res . 2009 ; 29 ( 7 ): 2453 - 9 . 20. Ying N , Wang S , Xu H , Wang Y . Association between FEN1 polymorphisms -69G> A and 4150G> T with susceptibility in human disease: a metaanalysis . Iran J Public Health . 2015 ; 44 ( 12 ): 1574 - 9 . 21. Chen Y , Zhang X , Qiu X , Li J , Yuan Q , Jiang T , Yang M. Functional FEN1 genetic variants and haplotypes are associated with glioma risk . J NeuroOncol . 2013 ; 111 ( 2 ): 145 - 51 . 22. Xu BJ , An QA , Srinivasa GS , Yan W , Pierce LA , Abel TW , Rush SZ , Cooper MK , Ye F , Shyr Y , et al. Identification of blood protein biomarkers that aid in the clinical assessment of patients with malignant glioma . Int J Oncol . 2012 ; 40 ( 6 ): 1995 - 2003 . 23. Pećina-Šlaus N , Gotovac K , Kafka A , Tomas D , Borovečki F. Genetic changes observed in a case of adult pilocytic astrocytoma revealed by array CGH analysis . Mol Cytogenet . 2014 ; 7 : 95 . 24. Cui G , Park S , Badeaux AI , Kim D , Lee J , Thompson JR , Yan F , Kaneko S , Yuan Z , Botuyan MV , et al. PHF20 is an effector protein of p53 double lysine methylation that stabilizes and activates p53 . Nat Struct Mol Biol . 2012 ; 19 ( 9 ): 916 - 24 . 25. Bankovic J , Stojsic J , Jovanovic D , Andjelkovic T , Milinkovic V , Ruzdijic S , Tanic N. Identification of genes associated with non-small-cell lung cancer promotion and progression . Lung Cancer . 2010 ; 67 ( 2 ): 151 - 9 . 26. Lasho TL , Pardanani A , Tefferi A , Smith DI . Altered expression of tumor suppressor PHF20 in myeloproliferative neoplasms . Leukemia . 2014 ; 28 ( 8 ): 1762 - 4 . 27. Verstrepen L , Bekaert T , Chau TL , Tavernier J , Chariot A , Beyaert R . TLR-4, IL-1R and TNF-R signaling to NF-κB: variations on a common theme . Cell Mol Life Sci . 2008 ; 65 ( 19 ): 2964 - 78 . 28. Spyridakis S , Leondaritis G , Nakos G , Lekka ME , Galanopoulou D. A specific phospholipase C activity regulates phosphatidylinositol levels in lung surfactant of patients with acute respiratory distress syndrome . Am J Resp Cell Mol . 2010 ; 42 ( 3 ): 357 - 62 . 29. Kong G , Chang YI , You X , Ranheim EA , Zhou Y , Burd CE , Zhang J . The ability of endogenous Nras oncogenes to initiate leukemia is codondependent . Leukemia . 2016 ; 30 ( 9 ): 1935 - 8 . 30. Jakob JA , Bassett RL , Ng CS , Curry JL , Joseph RW , Alvarado GC , Rohlfs ML , Richard J , Gershenwald JE , Kim KB , et al. NRAS mutation status is an independent prognostic factor in metastatic melanoma . Cancer Am Cancer Soc . 2012 ; 118 ( 16 ): 4014 - 23 . 31. Bleeker FE , Lamba S , Zanon C , Molenaar RJ , Hulsebos TJ . Mutational profiling of kinases in glioblastoma . BMC Cancer . 2014 ; 14 : 718 . 32. Tsurushima H , Tsuboi K , Yoshii Y , Ohno T , Meguro K , Nose T. Expression of N-ras gene in gliomas . Neurol Med Chir (Tokyo). 1996 ; 36 ( 10 ): 704 - 8 . 33. Herrero-Gonzalez S , Di Cristofano A . New routes to old places: PIK3R1 and PIK3R2 Join PIK3CA and PTEN as endometrial cancer genes . Cancer Discov . 2011 ; 1 ( 2 ): 106 - 7 . 34. Kita D , Yonekawa Y , Weller M , Ohgaki H. PIK3CA alterations in primary (de novo) and secondary glioblastomas . Acta Neuropathol . 2007 ; 113 ( 3 ): 295 - 302 . 35. Li L , Yang XJ . Molecular and functional characterization of histone deacetylase 4 (HDAC4) . Methods Mol Biol . 2016 ; 1436 : 31 - 45 . 36. Oehme I , Deubzer HE , Wegener D , Pickert D , Linke JP , Hero B , KoppSchneider A , Westermann F , Ulrich SM , von Deimling A , et al. Histone deacetylase 8 in neuroblastoma tumorigenesis . Clin Cancer Res . 2009 ; 15 ( 1 ): 91 - 9 . 37. Zeng LS , Yang XZ , Wen YF , Mail SJ , Wang MH , Zhang MY , Zheng XF , Wang HY . Overexpressed HDAC4 is associated with poor survival and promotes tumor progression in esophageal carcinoma . Aging (Albany NY) . 2016 ; 8 ( 6 ): 1236 - 49 . 38. Pei Y , Liu K , Wang J , Garancher A , Tao R , Esparza LA , Maier DL , Udaka YT , Murad N , Morrissy S , et al. HDAC and PI3K antagonists cooperate to inhibit growth of MYC-driven medulloblastoma . Cancer Cell . 2016 ; 29 ( 3 ): 311 - 23 .


This is a preview of a remote PDF: https://cancerci.biomedcentral.com/track/pdf/10.1186/s12935-017-0459-x?site=cancerci.biomedcentral.com

Tianlong Liu, Tiejun Zhang, Feng Zhou, Jitao Wang, Xiaohu Zhai, Nan Mu, Jongsun Park, Minna Liu, Wenxing Liu, Peijin Shang, Yi Ding, Aidong Wen, Yuwen Li. Identification of genes and pathways potentially related to PHF20 by gene expression profile analysis of glioblastoma U87 cell line, Cancer Cell International, 2017, 87, DOI: 10.1186/s12935-017-0459-x