Gene expression profile of THZ1-treated nasopharyngeal carcinoma cell lines indicates its involvement in the inhibition of the cell cycle

Background: The aim of this study was to identify downstream target genes and pathways regulated by THZ1 in nasopharyngeal carcinoma (NPC). Methods: The gene expression profile of GSE95750 in two NPC cell lines, untreated group and treated with THZ1 group, was analyzed. Differentially expressed genes (DEGs) were compared using the R-software. Then Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathways (KEGG) was analyzed using Database for Annotation, Visualization, and Integrated Discovery (DAVID). Cytoscape was used for protein-protein interaction (PPI) analysis. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to verified the gene expression. Results: We identified 25 genes with increased expression and 567 genes with decreased expression in THZ1-treated NPC cells. The top 10 significantly DEGs between untreated group and THZ1 treated group were identified by qRT-PCR and the results were in agreement with RNA-seq. The total 592 DEGs were found enriched in 1,148 GO terms and 38 KEGG pathways. The most important enriched pathways identified were cell cycle related, and several related node genes were identified, such as CDC6, CDC34, CDK7, CDK9, CCNA2, CCNB1, CDT1, KIF11, LIN9, PLK1, and POLR family, which consistent with RNA-seq. Conclusions: Our results emphasize the differential genes and pathways occurring in THZ1-treated NPC cells, which increases our understanding of the anti-tumor mechanisms of THZ1.