Abstract
Circular RNAs (circRNAs) are regulatory molecules that participate in the occurrence, development and progression of tumors. To obtain a complete blueprint of cervical carcinogenesis, we analyzed the temporal transcriptomic landscapes of mRNAs and circRNAs. Microarrays were performed to identify the circRNA and mRNA expression profiles of cervical squamous cell carcinoma (CSCC) and high-grade squamous intraepithelial lesion (HSIL) patients compared with normal controls (NC). Short time-series expression miner (STEM) was utilized to characterize the time-course expression patterns of circRNAs and mRNAs from NC to HSIL and CSCC. A total of 3 circRNA profiles and 3 mRNA profiles with continuous upregulated patterns were identified and selected for further analysis. Furthermore, functional annotation showed that the mRNAs were associated with DNA repair and cell division. The protein-protein interaction (PPI) network analysis revealed that the ten highest-degree genes were considered to be hub genes. Subsequently, a competing endogenous RNA (ceRNA) network analysis and real-time PCR validation indicated that hsa_circ_0001955/hsa-miR-6719-3p/CDK1, hsa_circ_0001955/hsa-miR-1277-5p/NEDD4L and hsa_circ_0003954/hsa-miR-15a-3p/SYCP2 were highly correlated with cervical carcinogenesis. Silencing of hsa_circ_0003954 inhibited SiHa cell proliferation and perturb the cell cycle in vitro. This study provides insight into the molecular events regulating cervical carcinogenesis, identifies functional circRNAs in CSCC, and improves the understanding of the pathogenesis and molecular biomarkers of CSCC and HSIL.