Integration of human papilloma virus (HPV) DNA into the human genome may progressively contribute to cervical carcinogenesis. To explore how HPV integration affects gene expression by altering DNA methylation during carcinogenesis, we analyzed a multiomics dataset for cervical cancer. We obtained multiomics data by HPV-capture sequencing, RNA sequencing, and Whole Genome Bisulfite Sequencing from 50 patients with cervical cancer. We detected 985 and 485 HPVintegration sites in matched tumor and adjacent paratumor tissues. Of these, LINC00486 (n = 19), LINC02425 (n = 11), LLPH (n = 11), PROS1 (n = 5), KLF5 (n = 4), LINC00392 (n = 3), MIR205HG (n = 3) and NRG1 (n = 3) were identified as highfrequency HPV-integrated genes, including five novel recurrent genes. Patients at clinical stage II had the highest number of HPV integrations. E6 and E7 genes of HPV16 but not HPV18 showed significantly fewer breakpoints than random distribution. HPV integrations occurring in exons were associated with altered gene expression in tumor tissues but not in paratumor tissues. A list of HPV-integrated genes regulated at transcriptomic or epigenetic level was reported. We also carefully checked the candidate genes with regulation pattern correlated in both levels. HPV fragments integrated at MIR205HG mainly came from the L1 gene of HPV16. RNA expression of PROS1 was downregulated when HPV integrated in its upstream region. RNA expression of MIR205HG was elevated when HPV integrated into its enhancer. The promoter methylation levels of PROS1 and MIR205HG were all negatively correlated with their gene expressions. Further experimental validations proved that upregulation of MIR205HG could promote the proliferative and migrative abilities of cervical cancer cells. Our data provides a new atlas for epigenetic and transcriptomic regulations regarding HPV integrations in cervical cancer genome. We demonstrate that HPV integration may affect gene expression

No abstract

No abstract