Investigation of differentially expressed genes in nasopharyngeal carcinoma by integrated bioinformatics analysis

Zou,Zhenning; Gan,Siyuan; Liu,Shuguang; Li,Rujia; Huang,Jian

doi:10.3892/ol.2019.10382

Investigation of differentially expressed genes in nasopharyngeal carcinoma by integrated bioinformatics analysis

Authors:
- Zhenning Zou
- Siyuan Gan
- Shuguang Liu
- Rujia Li
- Jian Huang
View Affiliations

Affiliations: Department of Pathology, Guangdong Medical University, Zhanjiang, Guangdong 524023, P.R. China, Department of Pathology, The Eighth Affiliated Hospital of Sun Yat‑sen University, Shenzhen, Guangdong 518033, P.R. China
Published online on: May 21, 2019 https://doi.org/10.3892/ol.2019.10382
Pages: 916-926
Copyright: © Zou et al. This is an open access article distributed under the terms of Creative Commons Attribution License.

Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )

Metrics: Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )

Cited By (CrossRef): 0 citations Loading Articles...

Abstract

Nasopharyngeal carcinoma (NPC) is a common malignancy of the head and neck. The aim of the present study was to conduct an integrated bioinformatics analysis of differentially expressed genes (DEGs) and to explore the molecular mechanisms of NPC. Two profiling datasets, GSE12452 and GSE34573, were downloaded from the Gene Expression Omnibus database and included 44 NPC specimens and 13 normal nasopharyngeal tissues. R software was used to identify the DEGs between NPC and normal nasopharyngeal tissues. Distributions of DEGs in chromosomes were explored based on the annotation file and the CYTOBAND database of DAVID. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were applied. Additionally, a protein‑protein interaction (PPI) network, constructed using the STRING database and visualized by Cytoscape, was used to identify hub genes, key modules and important transcription factors (TFs). A total of 906 DEGs were identified; 434 (47.90%) DEGs were upregulated and 472 (52.10%) were downregulated. The DEGs were demonstrated to be enriched in chromosome 7p15‑p14, 2q31, 1q21‑q22, 1q21, 4q21 and 1p31‑p22. DEGs were mainly enriched for the following GO terms: ‘Cilium movement’, ‘microtubule bundle formation’ and ‘axoneme assembly’. KEGG pathway enrichment analysis revealed that pathways for ‘cell cycle’, ‘DNA replication’, ‘interleukin‑17 signaling’, ‘amoebiasis’ and ‘glutathione metabolism’ were enriched. In addition, a PPI network comprising 867 nodes and 1,241 edges was constructed. Finally, five hub genes (aurora kinase A, cell division cycle 6, mitotic arrest deficient 2‑like 1, DNA topoisomerase 2α and TPX2 microtubule nucleation factor), 8 modules, and 14 TFs were identified. Modules analysis revealed that cyclin‑dependent kinase 1 and exportin 1 were involved in the pathway of Epstein‑Barr virus infection. In summary, the hub genes, key modules and TFs identified in this study may promote our understanding of the pathogenesis of NPC and require further in‑depth investigation.

View Figures

View References

1	Chua MLK, Wee JTS, Hui EP and Chan ATC: Nasopharyngeal carcinoma. Lancet. 387:1012–1024. 2016. View Article : Google Scholar : PubMed/NCBI
2	Fu ZT, Guo XL, Zhang SW, Zeng HM, Sun KX, Chen WQ and He J: Incidence and mortality of nasopharyngeal carcinoma in China, 2014. Zhonghua Zhong Liu Za Zhi. 40:566–571. 2018.(In Chinese; Abstract available in Chinese from the Publisher). PubMed/NCBI
3	Wei WI and Sham JS: Nasopharyngeal carcinoma. Lancet. 365:2041–2054. 2005. View Article : Google Scholar : PubMed/NCBI
4	Tsao SW, Yip YL, Tsang CM, Pang PS, Lau VM, Zhang G and Lo KW: Etiological factors of nasopharyngeal carcinoma. Oral Oncol. 50:330–338. 2014. View Article : Google Scholar : PubMed/NCBI
5	Bose S, Yap LF, Fung M, Starzcynski J, Saleh A, Morgan S, Dawson C, Chukwuma MB, Maina E, Buettner M, et al: The ATM tumour suppressor gene is down-regulated in EBV-associated nasopharyngeal carcinoma. J Pathol. 217:345–352. 2009. View Article : Google Scholar : PubMed/NCBI
6	Cai L, Ye Y, Jiang Q, Chen Y, Lyu X, Li J, Wang S, Liu T, Cai H, Yao K, et al: Epstein-Barr virus-encoded microRNA BART1 induces tumour metastasis by regulating PTEN-dependent pathways in nasopharyngeal carcinoma. Nat Commun. 6:73532015. View Article : Google Scholar : PubMed/NCBI
7	Wen X, Tang X, Li Y, Ren X, He Q, Yang X, Zhang J, Wang Y, Ma J and Liu N: Microarray expression profiling of long non-coding RNAs involved in nasopharyngeal carcinoma metastasis. Int J Mol Sci. 17:E19562016. View Article : Google Scholar : PubMed/NCBI
8	Sengupta S, den Boon JA, Chen IH, Newton MA, Dahl DB, Chen M, Cheng YJ, Westra WH, Chen CJ, Hildesheim A, et al: Genome-wide expression profiling reveals EBV-associated inhibition of MHC class I expression in nasopharyngeal carcinoma. Cancer Res. 66:7999–8006. 2006. View Article : Google Scholar : PubMed/NCBI
9	Hu C, Wei W, Chen X, Woodman CB, Yao Y, Nicholls JM, Joab I, Sihota SK, Shao JY, Derkaoui KD, et al: A global view of the oncogenic landscape in nasopharyngeal carcinoma: An integrated analysis at the genetic and expression levels. PLoS One. 7:e410552012. View Article : Google Scholar : PubMed/NCBI
10	Neuwirth E: RColorBrewer: ColorBrewer Palettes. R package version 1.1-2. https://cran.r-project.org/web/packages/RColorBrewer/index.html2014
11	Bolstad BM, Collin F, Brettschneider J, Simpson K, Cope L, Irizarry RA and Speed TA: Quality assessment of Affymetrix GeneChip data. Gentleman R, Carey VJ, Huber W, Irizarry RA and Dudoit S: Bioinformatics and computational biology solutions using R and bioconductor. Statistics for biology and health, Springer. (New York, NY). 2005. View Article : Google Scholar
12	Leek JT, Johnson WE, Parker HS, Jaffe AE and Storey JD: The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 28:882–883. 2012. View Article : Google Scholar : PubMed/NCBI
13	Gautier L, Cope L, Bolstad BM and Irizarry RA: Affy-analysis of affymetrix genechip data at the probe level. Bioinformatics. 20:307–315. 2004. View Article : Google Scholar : PubMed/NCBI
14	Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D and Altman RB: Missing value estimation methods for DNA microarrays. Bioinformatics. 17:520–525. 2001. View Article : Google Scholar : PubMed/NCBI
15	Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W and Smyth GK: Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e472015. View Article : Google Scholar : PubMed/NCBI
16	Huang da W, Sherman BT and Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 4:44–57. 2009. View Article : Google Scholar : PubMed/NCBI
17	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. The Gene Ontology Consortium. Nat Genet. 25:25–29. 2000. View Article : Google Scholar : PubMed/NCBI
18	Kanehisa M and Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28:27–30. 2000. View Article : Google Scholar : PubMed/NCBI
19	Yu G, Wang LG, Han Y and He QY: clusterProfiler: An R package for comparing biological themes among gene clusters. OMICS. 16:284–287. 2012. View Article : Google Scholar : PubMed/NCBI
20	Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, et al: The STRING database in 2017: Quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 45:D362–D368. 2017. View Article : Google Scholar : PubMed/NCBI
21	Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B and Ideker T: Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13:2498–2504. 2003. View Article : Google Scholar : PubMed/NCBI
22	Bader GD and Hogue CW: An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 4:22003. View Article : Google Scholar : PubMed/NCBI
23	Chin CH, Chen SH, Wu HH, Ho CW, Ko MT and Lin CY: CytoHubba: Identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 4 (Suppl 8):S112014. View Article : Google Scholar
24	Janky R, Verfaillie A, Imrichova H, Van de Sande B, Standaert L, Christiaens V, Hulselmans G, Herten K, Naval Sanchez M, Potier D, et al: iRegulon: From a gene list to a gene regulatory network using large motif and track collections. PLoS Computat Biol. 10:e10037312014. View Article : Google Scholar
25	Tulalamba W, Larbcharoensub N, Sirachainan E, Tantiwetrueangdet A and Janvilisri T: Transcriptome meta- analysis reveals dysregulated pathways in nasopharyngeal carcinoma. Tumor Biol. 36:5931–5942. 2015. View Article : Google Scholar
26	Liao KM, Chao TB, Tian YF, Lin CY, Lee SW, Chuang HY, Chan TC, Chen TJ, Hsing CH, Sheu MJ and Li CF: Overexpression of the PSAT1 gene in nasopharyngeal carcinoma is an indicator of poor prognosis. J Cancer. 7:1088–1094. 2016. View Article : Google Scholar : PubMed/NCBI
27	Chen F, Shen C, Wang X, Wang H, Liu Y, Yu C, Lv J, He J and Wen Z: Identification of genes and pathways in nasopharyngeal carcinoma by bioinformatics analysis. Oncotarget. 8:63738–63749. 2017.PubMed/NCBI
28	Dodd LE, Sengupta S, Chen IH, den Boon JA, Cheng YJ, Westra W, Newton MA, Mittl BF, McShane L, Chen CJ, et al: Genes involved in DNA repair and nitrosamine metabolism and those located on chromosome 14q32 are dysregulated in nasopharyngeal carcinoma. Cancer Epidemiol Biomarkers Prev. 15:2216–2225. 2006. View Article : Google Scholar : PubMed/NCBI
29	Lowe JS and Anderson PG: Chapter 10-Respiratory System. Stevens & Lowe's Human Histology (Fourth Edition). Lowe JS and Anderson PG: Mosby, Philadelphia: pp. 166–185. 2015, View Article : Google Scholar
30	Wang L, Ma R, Kang Z, Zhang Y, Ding H, Guo W, Gao Q and Xu M: Effect of IL-17A on the migration and invasion of NPC cells and related mechanisms. PLoS One. 9:e1080602014. View Article : Google Scholar : PubMed/NCBI
31	Gavet O and Pines J: Progressive activation of CyclinB1-Cdk1 coordinates entry to mitosis. Dev Cell. 18:533–543. 2010. View Article : Google Scholar : PubMed/NCBI
32	Zhu J, Liao G, Shan L, Zhang J, Chen MR, Hayward GS, Hayward SD, Desai P and Zhu H: Protein array identification of substrates of the Epstein-Barr virus protein kinase BGLF4. J Virol. 83:5219–5231. 2009. View Article : Google Scholar : PubMed/NCBI
33	Kato K, Yokoyama A, Tohya Y, Akashi H, Nishiyama Y and Kawaguchi Y: Identification of protein kinases responsible for phosphorylation of Epstein-Barr virus nuclear antigen leader protein at serine-35, which regulates its coactivator function. J Gen Virol. 84:3381–3392. 2003. View Article : Google Scholar : PubMed/NCBI
34	Cook ID, Shanahan F and Farrell PJ: Epstein-Barr virus SM protein. Virology. 205:217–227. 1994. View Article : Google Scholar : PubMed/NCBI
35	Boyle SM, Ruvolo V, Gupta AK and Swaminathan S: Association with the cellular export receptor CRM 1 mediates function and intracellular localization of Epstein-Barr virus SM protein, a regulator of gene expression. J Virol. 73:6872–6881. 1999.PubMed/NCBI
36	Vardabasso C, Hasson D, Ratnakumar K, Chung CY, Duarte LF and Bernstein E: Histone variants: Emerging players in cancer biology. Cell Mol Life Sci. 71:379–404. 2014. View Article : Google Scholar : PubMed/NCBI
37	Yang HD, Kim PJ, Eun JW, Shen Q, Kim HS, Shin WC, Ahn YM, Park WS, Lee JY and Nam SW: Oncogenic potential of histone-variant H2A.Z.1 and its regulatory role in cell cycle and epithelial-mesenchymal transition in liver cancer. Oncotarget. 7:11412–11423. 2016.PubMed/NCBI
38	Bui KH, von Appen A, DiGuilio AL, Ori A, Sparks L, Mackmull MT, Bock T, Hagen W, Andrés-Pons A, Glavy JS and Beck M: Integrated structural analysis of the human nuclear pore complex scaffold. Cell. 155:1233–1243. 2013. View Article : Google Scholar : PubMed/NCBI
39	Alanee S, Delfino K, Wilber A, Robinson K, Brard L and Semaan A: Single nucleotide variant in Nucleoporin 107 may be predictive of sensitivity to chemotherapy in patients with ovarian cancer. Pharmacogenet Genomics. 27:264–269. 2017. View Article : Google Scholar : PubMed/NCBI
40	Imaoka H, Toiyama Y, Saigusa S, Kawamura M, Kawamoto A, Okugawa Y, Hiro J, Tanaka K, Inoue Y, Mohri Y and Kusunoki M: RacGAP1 expression, increasing tumor malignant potential, as a predictive biomarker for lymph node metastasis and poor prognosis in colorectal cancer. Carcinogenesis. 36:346–354. 2015. View Article : Google Scholar : PubMed/NCBI
41	Han P, Lin ZR, Xu LH, Zhong Q, Zhu XF, Liang FY, Cai Q, Huang XM and Zeng MS: Ribonucleotide reductase M2 subunit expression and prognostic value in nasopharyngeal carcinoma. Mol Med Rep. 12:401–409. 2015. View Article : Google Scholar : PubMed/NCBI
42	Liu Q, Yang P, Tu K, Zhang H, Zheng X, Yao Y and Liu Q: TPX2 knockdown suppressed hepatocellular carcinoma cell invasion via inactivating AKT signaling and inhibiting MMP2 and MMP9 expression. Chin J Cancer Res. 26:410–417. 2014.PubMed/NCBI
43	Jiang L, Wang P and Chen H: Overexpression of FOXM1 is associated with metastases of nasopharyngeal carcinoma. Ups J Med Sci. 119:324–332. 2014. View Article : Google Scholar : PubMed/NCBI
44	Jiang L, Wang P, Chen L and Chen H: Down-regulation of FoxM1 by thiostrepton or small interfering RNA inhibits proliferation, transformation ability and angiogenesis, and induces apoptosis of nasopharyngeal carcinoma cells. Int J Clin Exp Pathol. 7:5450–5460. 2014.PubMed/NCBI
45	Yu C, Chen L, Yie L, Wei L, Wen T, Liu Y and Chen H: Targeting FoxM1 inhibits proliferation, invasion and migration of nasopharyngeal carcinoma through the epithelialto-mesenchymal transition pathway. Oncol Rep. 33:2402–2410. 2015. View Article : Google Scholar : PubMed/NCBI
46	Liu Z, Long X, Chao C, Yan C, Wu Q, Hua S, Zhang Y, Wu A and Fang W: Knocking down CDK4 mediates the elevation of let-7c suppressing cell growth in nasopharyngeal carcinoma. BMC Cancer. 14:2742014. View Article : Google Scholar : PubMed/NCBI
47	Ding L, Li L, Yang J, Zhou S, Li W, Tang M, Shi Y, Yi W and Cao Y: Latent membrane protein 1 encoded by Epstein-Barr virus induces telomerase activity via p16INK4A/Rb/E2F1 and JNK signaling pathways. J Med Virol. 79:1153–1163. 2007. View Article : Google Scholar : PubMed/NCBI
48	Hsu HP, Li CF, Lee SW, Wu WR, Chen TJ, Chang KY, Liang SS, Tsai CJ and Shiue YL: Overexpression of stathmin 1 confers an independent prognostic indicator in nasopharyngeal carcinoma. Tumour Biol. 35:2619–2629. 2014. View Article : Google Scholar : PubMed/NCBI
49	Chung GT, Lou WP, Chow C, To KF, Choy KW, Leung AW, Tong CY, Yuen JW, Ko CW, Yip TT, et al: Constitutive activation of distinct NF-κB signals in EBV-associated nasopharyngeal carcinoma. J Pathol. 231:311–322. 2013. View Article : Google Scholar : PubMed/NCBI
50	Yi M, Cai J, Li J, Chen S, Zeng Z, Peng Q, Ban Y, Zhou Y, Li X, Xiong W, et al: Rediscovery of NF-κB signaling in nasopharyngeal carcinoma: How genetic defects of NF-κB pathway interplay with EBV in driving oncogenesis? J Cell physiol. 233:5537–5549. 2018. View Article : Google Scholar : PubMed/NCBI