Abstract
Infection with human papillomavirus (HPV) is a major risk factor for head and neck squamous cell carcinoma (HNSCC). The objective of this study is to investigate the gene expression profiles and signaling pathways that are specific to HPV-positive HNSCC (HPV+ HNSCC). Moreover, a competing endogenous RNA (ceRNA) network analysis was utilized to identify the core gene of HPV+ HNSCC and potential targeted therapeutic drugs. Transcriptome sequencing analysis identified 3,253 coding RNAs and 3,903 non-coding RNAs (ncRNAs) that exhibited preferentially expressed in HPV+ HNSCC. Four key signaling pathways were selected through pathway enrichment analysis. By combining ceRNA network and protein-protein interaction (PPI) network topology analysis, RNA Polymerase II Associated Protein 2 (RPAP2), which also exhibited high expression in HPV+ HNSCC based on the TCGA database, was identified as the hub gene. Gene set enrichment analysis (GSEA) results revealed RPAP2's involvement in various signaling pathways, encompassing basal transcription factors, ubiquitin-mediated proteolysis, adherens junction, other glycan degradation, ATP-binding cassette (ABC) transporters, and oglycan biosynthesis. Five potential small molecule targeted drugs (enzastaurin, brequinar, talinolol, phenylbutazone, and afuresertib) were identified using the cMAP database, with enzastaurin showing the highest affinity for RPAP2. Cellular functional experiments confirmed the inhibitory effect of enzastaurin on cell viability of HPV+ HNSCC and RPAP2 expression levels. Additionally, enzastaurin treatment suppressed the expression levels of the top-ranked long non-coding RNA (lncRNA), circular RNA (circRNA), and microRNA (miRNA) in the ceRNA network. This study based on the ceRNA network provides valuable insights into the molecular mechanisms and potential therapeutic strategies for HPV+ HNSCC, and provide theoretical basis for the exploration of HPV+ HNSCC biomarkers and the development of targeted drugs.