High-risk human papillomavirus (HPV) is a causal factor for cervical cancer, of which HPV16 is the predominant genotype, but the detailed mechanism remains to be elucidated. In this study, we performed transcriptome sequencing in cervical cancer tissues with HPV16-positive and normal tissues with HPV16-negative, and SiHa cells with or without HPV16 E6/E7 knockdown, and identified 140 differential expressed genes (DEGs) in two data sets. We carried out a series of bioinformatic analyses to learn more about the 140 DEGs, and found that 140 DEGs were mostly enriched in cell cycle and DNA repair through Kyoto Encyclopedia of Genes and Genomes pathway enrichment, Gene Ontology annotation, and gene set enrichment analysis. A total of 20 genes including RMI1, MKI67, FANCB, KIF14, CENPI, RACGAP1, EXO1, KIF4A, FOXM1, C19orf57, PSRC1, NUSAP1, CIT, NDC80, MCM7, GINS2, MCM6, ORC1, TLX2, and UHRF1 were screened by co-expression analysis; of those, the expressions of 6 (CENPI, FANCB, KIF14, ORC1, RACGAP1, and RMI1) were verified by qRT-PCR. Further, we found that E2F family, NF-Y, AhR:Arnt, and KROX family may be involved in modulating DEGs by TransFind prediction. TF2DNA database and co-expression analysis suggested that 12 TFs (ZNF367, TLX2, DEPDC1B, E2F8, ZNF541, EGR2, ZMAT3, HES6, CEBPA, MYBL2, FOXM1, and RAD51) were upstream modulators of DEGs. Our findings may provide a new understanding for effects of HPV oncogenes in the maintenance of cancerous state at the transcriptional level.