Background:Multiple genes were previously identified to be associated with cervical cancer; however, the genetic architecture of cervical cancer remains unknown and many causal genes have yet been discovered.Methods: To explore causal genes related to cervical cancer, a two-stage causal inference approach was proposed within the framework of Mendelian randomization, where the gene expression was treated as exposure, with methylations located within that gene serving as instrumental variables. Five prediction models were first utilized to characterize the relationship between the expression and methylations for each gene; then the methylation-regulated gene expression (MReX) was obtained and the association was evaluated via Cox mixed-effects model based on MReX. We further implemented the harmonic mean p-value (HMP) combination to take advantage of respective strengths of these prediction models while accounting for dependency among the p-values.Results: A total of 14 causal genes were discovered to be associated with the survival risk of cervical cancer in TCGA when the five prediction models were separately employed. The total number of causal genes was brought to 23 when conducting HMP. Some of the newly discovered genes may be novel (e.g. YJEFN3, SPATA5L1, IMMP1L, C5orf55, PPIP5K2, ZNF330, CRYZL1, PPM1A, ESCO2, ZNF605, ZNF225, ZNF266, FICD and OSTC). Functional analyses showed these genes were enriched in tumor-associated pathways. Additionally, four genes (i.e. COL6A1, SYDE1, ESCO2 and GIPC1) were differentially expressed.Conclusion: Overall, our study discovered promising candidate genes that are causally associated with the survival risk of cervical cancer and thus provided new insights into the genetic etiology of cervical cancer.