Potential Biomarkers and Signaling Pathways Associated with the Pathogenesis of Primary Ameloblastoma: A Systems Biology Approach

Abstract

Objective Ameloblastoma is a benign odontogenic tumor that may lead to ameloblastic carcinoma. This study aimed to determine potential signaling pathways and biological processes, critical genes and their regulating transcription factors (TFs), and miRNAs, as well as protein kinases involved in the etiology of primary ameloblastoma. Methods The dataset GSE132472 was obtained from the GEO database, and multivariate statistical analyses were applied to identify differentially expressed genes (DEGs) in primary ameloblastoma tissues compared to the corresponding normal gingiva samples. A protein-protein interaction (PPI) map was built using the STRING database. The Cytoscape software identified significant modules and the hub genes within the PPI network. Gene Ontology annotation and signaling pathway analyses were executed by employing the DAVID and Reactome databases, respectively. Significant TFs and miRNAs acting on the hub genes were identified using the iRegulon plugin and MiRWalk 2.0 database, respectively. A protein kinase enrichment analysis was conducted using the online Kinase Enrichment Analysis 2 (KEA2) web server. The approved drugs acting on the hub genes were also found. Results A total of 1,629 genes were differentially expressed in primary ameloblastoma (P value <0.01 and |Log2FC| > 1). HRAS, CDK1, MAPK3, ERBB2, COL1A1, CYCS, and BRCA1 demonstrated high degree and betweenness centralities in the PPI network. E2F4 was the most significant TF acting on the hub genes. BTK was the protein kinase significantly enriched by the TFs. Cholesterol biosynthesis was considerably involved in primary ameloblastoma. Conclusions This study provides an intuition into the potential mechanisms involved in the etiology of ameloblastoma.