Integrated clinicopathological features and gene microarray analysis of pancreatic neuroendocrine tumors

Pancreatic neuroendocrine tumors are relatively rare pancreatic neoplasms over the world. Investigations about molecular biology of PNETs are insufficient for nowadays. We aimed to explore the expression of messenger RNA and regulatory processes underlying pancreatic neuroendocrine tumors from different views. The expression profile of GSE73338 were downloaded, including samples with pancreatic neuroendocrine tumors. First, the Limma package was utilized to distinguish the differentially expressed messenger RNA. Gene Ontology classification and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were performed to explore the functions and pathways of target genes. In addition, we constructed a protein-protein interaction network. NEK2, UBE2C, TOP2A and PPP1R1A were revealed with continuous genomic alterations in higher tumor stage. 91 up-regulated and 36 down-regulated genes were identified to be differentially expressed in malignant PNETs. Locomotory behavior was significantly enriched for biological processes of metastasis PNETs. GCGR and GNAS were identified as the hub of proteins in the protein-protein interaction sub-network of malignant PNETs. We showed the gene expression differences in PNETs according to different clinicopathological aspects. NEK2, UBE2C, TOP2A are positively associated with high tumor grade, and PPP1R1A negatively. GCGR and GNAS are regarded as the hub of the PPI sub-network. CXCR4 may affect the progression of PNETs via the CXCR4-CXCL12-CXCR7 chemokine receptor axis. However, more studies are required.