Characteristic of HPV Integration in the Genome and Transcriptome of Cervical Cancer Tissues

High-risk HPV is clearly associated with cervical cancer. HPV integration has been confirmed to promote carcinogenesis in the previous studies. In our study, a total of 285 DNA breakpoints and 287 RNA breakpoints were collected. We analyzed the characteristic of HPV integration in the DNA and RNA samples. The results revealed that the patterns of HPV integration in RNA and DNA samples differ significantly. FHIT, KLF5, and LINC00392 were the hotspot genes integrated by HPV in the DNA samples. RAD51B, CASC8, CASC21, ERBB2, TP63, TEX41, RAP2B, and MYC were the hotspot genes integrated by HPV in RNA samples. Breakpoints of DNA samples were significantly prone to the region of INTRON (P < 0.01, Chi-squared test), whereas in the RNA samples, the breakpoints were prone to EXON. Pathway analysis had revealed that the breakpoints of RNA samples were enriched in the pathways of transcriptional misregulation in cancer, cancer pathway, and pathway of adherens junction. Breakpoints of DNA samples were enriched in the pathway of cholinergic synapse. In summary, our data helped to gain insights into the HPV integration sites in DNA and RNA samples of cervical cancer. It had provided theoretical basis for understanding the mechanism of tumorigenesis from the perspective of HPV integration in the HPV-associated cervical cancers.

[1]  B. Oostra,et al.  A long-range Shh enhancer regulates expression in the developing limb and fin and is associated with preaxial polydactyly. , 2003, Human molecular genetics.

[2]  S. Vinokurova,et al.  Systematic Review of Genomic Integration Sites of Human Papillomavirus Genomes in Epithelial Dysplasia and Invasive Cancer of the Female Lower Genital Tract , 2004, Cancer Research.

[3]  P. Nowell,et al.  Sequence analysis of the MYC oncogene involved in the t(8;14)(q24;q11) chromosome translocation in a human leukemia T-cell line indicates that putative regulatory regions are not altered. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[4]  Li Ding,et al.  Activating HER2 mutations in HER2 gene amplification negative breast cancer. , 2013, Cancer discovery.

[5]  S. Nomoto,et al.  DeltaNp63alpha and TAp63alpha regulate transcription of genes with distinct biological functions in cancer and development. , 2003, Cancer research.

[6]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[7]  M. Machiela,et al.  Genomic characterization of viral integration sites in HPV‐related cancers , 2016, International journal of cancer.

[8]  T. Ried,et al.  Genome-wide analysis of HPV integration in human cancers reveals recurrent, focal genomic instability , 2014, Genome research.

[9]  T. Enomoto,et al.  Aberrant FHIT transcripts in squamous cell carcinoma of the uterine cervix , 1998, International journal of cancer.

[10]  S. Wolf,et al.  Multiplex Identification of Human Papillomavirus 16 DNA Integration Sites in Cervical Carcinomas , 2013, PloS one.

[11]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[12]  S. Nomoto,et al.  ΔNp63α and TAp63α regulate transcription of genes with distinct biological functions in cancer and development , 2003 .

[13]  S. Wacholder,et al.  Rapid clearance of human papillomavirus and implications for clinical focus on persistent infections. , 2008, Journal of the National Cancer Institute.

[14]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[15]  Keegan D. Korthauer,et al.  Chromosomal copy number alterations and HPV integration in cervical precancer and invasive cancer. , 2016, Carcinogenesis.

[16]  X. Fang,et al.  Genome-wide profiling of HPV integration in cervical cancer identifies clustered genomic hot spots and a potential microhomology-mediated integration mechanism , 2015, Nature Genetics.

[17]  Wendy Winckler,et al.  Functional analysis of receptor tyrosine kinase mutations in lung cancer identifies oncogenic extracellular domain mutations of ERBB2 , 2012, Proceedings of the National Academy of Sciences.

[18]  H. Kueh,et al.  A far downstream enhancer for murine Bcl11b controls its T-cell specific expression. , 2013, Blood.

[19]  M. Tarsounas,et al.  RAD51 paralogs: roles in DNA damage signalling, recombinational repair and tumorigenesis. , 2011, Seminars in cell & developmental biology.

[20]  Jay Shendure,et al.  The haplotype-resolved genome and epigenome of the aneuploid HeLa cancer cell line , 2013, Nature.

[21]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[22]  M. Polz-Dacewicz,et al.  Human Papillomavirus (HPV) - structure, epidemiology and pathogenesis. , 2014, Otolaryngologia polska = The Polish otolaryngology.

[23]  Trevor J Pugh,et al.  Landscape of genomic alterations in cervical carcinomas , 2013, Nature.

[24]  C. Crum,et al.  Natural history of cervical neoplasia: defining progression and its consequence. , 2000, Clinical obstetrics and gynecology.

[25]  A. Kimura,et al.  Positive and Negative Regulation of the Cardiovascular Transcription Factor KLF5 by p300 and the Oncogenic Regulator SET through Interaction and Acetylation on the DNA-Binding Domain , 2003, Molecular and Cellular Biology.