Viral E6 is overexpressed via high viral load in invasive cervical cancer with episomal HPV16
© The Author(s). 2017
Received: 1 August 2016
Accepted: 8 February 2017
Published: 15 February 2017
The integration of HR-HPV genome into host DNA is regarded as a key step for the development of cervical cancer. However, HR-HPV genome indeed exists as episome except for integrant. It may be alternative mechanisms in episome-associated carcinogenesis, although, by which HPV 16 episome induces cervical carcinogenesis is unclear now.
Ninety-three invasive cervical cancer tissues with HPV16 positive were collected. Viral physical status was calculated from comparing E2 to E6-copies and detection of viral load was made with realtime-PCR using copy numbers of E6. HPV16 E6 mRNA transcript levels were measured by realtime-PCR. The methylation frequency of HPV16 promoter was detected by PCR and pyrosequencing.
In 93 samples, 21.5% (20/93) presented purely integrated viral genome, 53.8% (50/93) mixed viral genome, and 24.7% (23/93) purely episomal viral genome. Mean E6 expression in samples with purely episomal viral genomes was 7.13-fold higher than that with purely integrated viral genomes. Meanwhile, viral load in samples with purely episomal viral genomes was 4.53-fold higher than that with purely integrated viral genomes. E6 mRNA expression increased with the viral load in purely episomal cases. There were no differences of mean methylation frequency between purely episomal and integrated virus and among five CpG positions of HPV16 promoter for all samples. And there also was no correlation between E6 mRNA expression and methylation of HPV16 promoter among all samples with purely HPV16 episomal virus.
HPV16 with the purely episomal viral genomes exists in a definite proportion of invasive cervical cancer, and episomal HPV16 also overexpresses E6 mRNA, probably through a high level of viral load.
KeywordsCervical cancer HPV16 Episome Integration Oncogene expression Viral load
Cervical cancer is one of the most common cancers among women worldwide. Infection with high-risk human papillomaviruses (HPV) is a causal factor for cervical intraepithelial neoplasia and cervical cancer. Among all the HR-HPV genotypes, HPV 16 is the most prevalent, reaching 65.2% of all genotypes in cervical cancer . Human papillomavirus (HPV) genome integration into the host chromosome is considered as a crucial event during the life cycle of the virus and a major step towards carcinogenesis . The integration of HPV16 DNA promotes a constitutive high expression level of E6 and E7 oncoproteins, resulting in the extensive proliferation of the infected epithelial cells. Typical integrants have complete or partial disruption of the open reading frame (ORF) in E2 ORF . An important consequence of viral integration is an abolishment of E2 gene inhibiting the expression of oncogene E6 and E7, resulting in the extensive proliferation and malignant transformation of infected epithelial cells. Higher steady-state levels of viral oncogene transcripts in precancers and invasive cancers were expected as a consequence of virus integration. Thus, the integration of HR-HPV genome into host DNA is regarded as a key, even prerequisite, step for the development of cervical cancer.
However, the actual physical status of HPV 16 genome in cervical cancer cells appears far more complicated. For instance, Li et al.  reported the existence of HPV 16 integration in all of 15 cases of cervical cancer tissues, but Dutta S et al.  detected 82% cases contained HPV16 integrant in cervical cancer samples, and Mazumder D and his collogues found that 70.3% samples harbored integration HPV16 in cervical cancer tissue . Further, Vinokurova and collaborators found that as high as 45% samples contained purely episomal virus in HPV16-positive cervical cancer tissues , and Cheung et al.  found pure episomal HPV16 genomes in 14 of 29 (48.3%) cervical cancer tissues. Therefore, HR-HPV genome indeed exists as episome except for integrant. This phenomenon points toward the biological plausibility of cervical carcinogenesis under the impact of HPV16 episome, in addition to E2 disruption due to viral genome integration into the host genome. It may be alternative mechanisms in episome-associated carcinogenesis, although, by which HPV 16 episome induces cervical carcinogenesis is unclear now.
Here, we detected HPV16 physical status, E6 mRNA transcript level, viral load, and E6 promoter methylation in HPV16 positive invasive cervical cancer tissues, and analyzed the association of HPV16 E6 expression with viral load and promoter methylation. The aim of the study was to understand the potential mechanism in cervical cancer pathogenesis with HPV16 episome infection.
Tissue sample collection
A total of 93 biopsied cervical cancer tissues were obtained from female patients (range 23–71 years, median 39 years) treated at the Department of Gynecology Oncology, Women’s Hospital, School of Medicine, Zhejiang University, China. All samples were histology-confirmed as invasive cervical cancer and had been confirmed to HPV16 infection by type specific PCR. Eighty-four cases of squamous cancer, three of adenosquamous cancer and six of adenocarcinoma were included. Two cases were stage IA1, 66 stage IB1, 13 stage IB2, ten stage IIA and two stage IIB. The study was done in accordance with the guidelines of the local ethical committee.
DNA extraction, RNA extraction and reverse transcription
Biopsy tissues were physically disrupted by magnetic beads. DNA was extracted using UniversalGen DNA Kit (Cowin Biotech, Beijing, China) and stored at −80 °C until use. Total RNAs were extracted with TRIZOL reagent (Invitrogen Carlsbad, CA, USA) following manufacturer's instructions. cDNA was synthesized with PrimeScriptTM RT Master Mix (TaKaRa Otsu, Shiga, Japan) and stored at −80 °C until use.
Detection of HPV16 physical status
To determine the physical status of the virus, E2 and E6 gene of HPV 16 DNA were quantified and the E2/E6 ratio was calculated. This approach is based on the fact that episomes present an identical amount of E2 and E6, whereas integration induces loss of E2.
The primer sequences for PCRs
Detection of HPV16 viral load
The strategy involved the measurement of the total viral load of HPV16 DNA by quantification of the HPV16 E6 gene and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) with real-time PCR. The sequence information of GAPDH DNA primers (GAPDH-DF and GAPDH-DR)  were seen in Table 1. Standard curves for GPADH was obtained by amplification of a 2-fold dilution series of female human DNA between 112.5 and 1.75 ng/μL (Promega).
The amount of genomic DNA (ng) presented in each sample was divided by the weight of 1 genome equivalent (6.6 pg/cell) to obtain the number of cells in the sample . Viral load was expressed as the copies of E6 per cell.
Detection of E6 mRNA transcripts
E6 mRNA transcript expression was quantified by using the same primer. Human Endogenous Control GAPDH (GAPDH-cDF and GAPDH-cDR) was used as normalizer, listed in Table 1. Real-time PCR was detected under the same conditions with an initial denaturation at 95 °C for 15 s, reaction mixtures underwent 40 cycles at 95 °C for 5 s and followed by 60 °C for 30 s. Samples were also analyzed in triplicates in each trial.
Detection the methylation frequency of HPV16 promoter
The Methylation frequency of HPV16 promoter was detected with PCR and pyrosequencing. The promoter P97 of HPV 16 contains potential methylation sites with five CpG dinucleotides located at 31, 37, 43, 52 and 58 nt (reference sequence NC_001526). Samples were treated with bisulphite using EZ DNA Methylation-Gold Kit according to manufacturer’s instructions (ZYMO Research, Irvine, USA). HPV16MP-F and HPV16MP-Rbio primers (Seen in Table 1) were used in PCR . PCR conditions were as follows: preheating at 94 °C for 5 min, 40 cycles at 94 °C for 45 s, 55 °C for 45 s, 72 °C for 45 s, and a final extension at 72 °C for 10 min. One of the primers was labeled with 5-biotin. Then PCR amplifiers were detected by pyrosequencing with the primer HPV16MP-S. Samples were prepared for pyrosequencing using the Vacuum prep Workstation (Qiagen). Single-strand sequencing template was transferred to a 96-well sequencing plate containing sequencing primer HPV16MP-S (Seen in Table 1). The plate was incubated in 85 °C for 2 min. Pyrosequencing was performed in a PSQ 96 MA system using PyroMark Q96 Reagents (Qiagen). The results were analyzed in the Pyro Q-CpG Software to determine the proportion of C/T at the targeted position(s). Methylated samples were further divided into different methylation degree groups of highly methylated (>50%), medium methylated (11–50%) and low methylated (<10%).
Data were analyzed with independent-samples t-test, Mann–Whitney test, Kruskal–Wallis test and Chi-Square test using IBM SPSS statistics 20. All P-values less than 0.05 were considered statistically significant. All reported P-values were bilateral.
Physical status of HPV16 in invasive cervical cancer tissues
The expression of E6 mRNA and viral loads in different physical status from invasive cervical cancer tissues
We quantified E6 expression (normalized by GAPDH) by qRT-PCR of cDNA products generated directly from mRNA. Relative quantification based on comparative CT method revealed the significant difference (P = 0.003; Kruskal-Wallis H test) among ∆CT (E6 CT - GPADH CT) values of cases with integrated viral genome (mean ∆CT = 3.28), mixed viral genome (mean ∆CT = 2.27) and episomal viral genome (mean ∆CT = 0.45). The fold-change analysis (using 2-∆∆CT, where ∆∆CT = mean ∆CT of purely episomal - mean ∆CT of purely integrated), depicted that E6 expression in cases with purely episomal viral genomes was 7.13 folds higher than that in cases with purely integrated viral genomes.
Samples were further divided into two groups of high viral load group (≧50copies/cell) and low viral load group (<50copies/cell) on the basis of viral load. E6 mRNA expression was higher in high viral load group than that in low viral load group in purely episomal cases (P = 0.007; Mann–Whitney Test), but E6 mRNA expression was not significantly different between high and low viral load group in purely integrated cases.
E6 mRNA expression and the methylation of HPV16 promoter from invasive cervical cancer tissues
The correlation between the expression of E6 mRNA and methylation of HPV16 promoter was investigated. There were no differences in E6 mRNA expression among different methylation degrees of promoter in all 90 cases of invasive cervical cancer tissues. And there also was no correlation between E6 mRNA expression and methylation of HPV16 promoter among all samples with purely HPV16 episomal virus.
When high-risk human papillomavirus infects host cells, some of viral genomes integrated into host genome. As a consequence, deregulated expression of the HPV E6 and E7 genes in epithelial stem cells leads to malignant transformation in the respective cells at early stages of dysplasia. Integration of viral genome into the host genome, chiefly at fragile sites [2, 12], is considered as a critical event in the pathogenesis of cervical neoplasia due to the loss of negative feedback control of oncogene expression from viral itself through E2 disruption [13, 14]. In addition to viral oncogene over-expression, integrated HPV also may be regarded as a selectable form of the virus because this form not only is resistant to clearance from host, but also enables infected cells to maintain viral oncogene expression and avoids cell death .
However, HPV 16 also can exist as episome in clinical invasive cervical cancer tissues actually. Our results showed that a substantial proportion of invasive cervical cancer cases (24.7%) contained purely episomal viral genome, only 21.5% of invasive cervical cancer samples contained fully integrated genome, and 53.8% of samples had a mixed genome (integrated and episomal viral genome). Our study identified that a substantial proportion of individuals of invasive cervical cancer carried episomal virus, in line with previous studies [5, 8, 16]. Thus, our study suggest, together with other studies , that purely episomal viral genome can exist in invasive cervical cancer and integration is not a prerequisite step for invasive cervical cancer development.
It is known that high levels of the viral oncoproteins E6 and E7 in proliferating cells of epithelia are a necessity for oncogenic transformation. Up-regulation of E6/E7 expression also is constitutively required to maintain the transformed phenotype of invasive cervical cancer regardless of integration or episome status . An in vitro experiment has revealed similarities between the episome and integrant associated routes of neoplastic progression . Some reports have also showed that integration of the HPV16 genome does not invariably result in high levels of viral oncogene transcripts , and HR-HPV integration per se dose not necessarily lead to increased oncogene expression . Our study found that purely episomal viral genome of HPV 16 in invasive cancer tissues also had E6 mRNA overexpression, which even higher than that of integrated virus, suggesting that HPV oncogenes are also overexpressed in invasive cervical cancer with viral episomal infection, but the mechanism may be involved in a way other than viral integration and E2 disruption.
It has been shown that viral load is one of the risk factors in invasive cervical cancer development . Therefore, we detected viral load in different viral physical status. In our study, purely HPV 16 episome contained 238.31 copies per cell in invasive cervical cancer tissues, but purely integrated HPV 16 contained only 52.56 copies per cell. Mean viral load in cases with purely episomal viral genomes was 4.53-fold compared with purely integrated viral genomes. We further analyzed viral load and E6 mRNA transcript levels of HPV16 in different viral physical status, and found that mean E6 expression in cases harboring purely episomal viral genomes was 7.13-fold higher than those harboring purely integrated forms. Furthermore, our result showed that E6 mRNA expression increased simultaneously with the viral loads in HPV 16 episomal genomes, but not in HPV 16 integrated genomes. Similarly, Marongiu L and his collogues previously found that samples harboring solely episomal HPV16 DNA had a higher viral load than samples with solely integrated forms . Thus, our and previously findings suggest that the overexpression of HPV oncogene may be mediated through a high viral load under HPV episomal status without E2 disruption . In fact, intact E2 protein can enhance viral DNA replication by interacting with the viral replication factor E1 and recruiting it to the origin of replication [21, 22], and consequently facilitates viral genome segregation by tethering the viral genomes to host mitotic chromosomes .
Gene expression is influenced by promoter DNA methylation. HPV 16 E6 and E7 genes are transcribed from the promoter P97. The promoter P97 region contains several binding sites for cellular transcription factors, such as Sp1, TFIID, and E2BS. DNA methylation within the binding sites of transcription factors, such as Sp1, might block binding indirectly, either by changing the conformation of chromatin or by interacting with methyl-CpG-specific repressor proteins . The capacity of the full-length E2 ORF gene product E2 to bind E2BSs in vitro is inhibited by methylation of cytosines within its binding site . It is known that the interaction of E2 with various E2BSs mediates repression or activation of transcription and, through its interaction with E1, viral DNA replication . However to date, few studies have examined the E6 mRNA transcript levels with respect to the promoter methylation of HPV16 in clinical invasive cervical tissues. Cheung JL and his collogues found that higher proportion of samples harboring the pure episomal form had high-level methylation at the 4 CpGs of two E2BSs in the promoter, compared to those of integrated forms (including mixed forms) . However, our study did not find the differences in E6 mRNA expression between purely episomal and integrated virus and among different methylation degrees of promoter in the cases that containing purely HPV16 episomal virus. Thus, the impact of HPV promoter methylation on viral ocogene expression appears to be limited, but a further investigation with larger samples is needed.
In summary, we found that HPV16 with the purely episomal viral genomes exists in a definite proportion of invasive cervical cancer, and episomal HPV16 also overexpresses E6 mRNA, probably through a high level of viral load. The impact of HPV16 E6/E7 promoter methylation on E6 expression may be limited in invasive cervical cancer with episomal HPV16 infection. Our results may provide another mechanism by which HR-HPV with episomal viral genomes induces the development of cervical cancer.
High-risk human papillomaviruses
Messenger Ribonucleic Acid
Open reading frame
We are grateful to members of the Department of Gynecologic Oncology, Women’s Hospital, School of Medicine, Zhejiang University, for access to and assistance with the sample collection.
This work was supported by the National Natural Science Foundation of China (81202065, 81202066, 81472421 and 81301687) and Natural Science Foundation of Zhejiang Province (Y12H160091).
Availability of data and materials
The datasets generated and/or analyzed during the current study are not publicly available due complicated calculations but are available from the corresponding author on reasonable request.
Study conceptualization and design: DH; data acquisition: DH, JL and YH; quality control of data and algorithms: DH, XNL and BHL; data analysis and interpretation: DH, WGL, XX and XDC; statistical analysis: DH, YL and DXH; manuscript preparation: DH, JL, YH and XDC; manuscript editing: DH, XX and XDC; and manuscript review: DH, XX and XDC. All authors read and approved the final manuscript.
The authors declare no conflicts of interest.
Consent for publication
Ethics approval and consent to participate
This study was approved by Ethics Committee of Woman’ Hospital, School of Medicine, Zhejiang University. All patients had written informed consent for tissue sampling used to scientific research prior to surgical treatment.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hong D, Ye F, Chen H, Lu W, Cheng Q, Hu Y, Xie X. Distribution of human papillomavirus genotypes in the patients with cervical carcinoma and its precursors in Zhejiang province, china. Int J Gynecol Cancer. 2008;18(1):104–9.View ArticlePubMedGoogle Scholar
- Wentzensen N, Vinokurova S, von Knebel DM. Systematic review of genomic integration sites of human papillomavirus genomes in epithelial dysplasia and invasive cancer of the female lower genital tract. Cancer Res. 2004;64(11):3878–84.View ArticlePubMedGoogle Scholar
- Luft F, Klaes R, Nees M, Durst M, Heilmann V, Melsheimer P, Von Knebel DM. Detection of integrated papillomavirus sequences by ligation-mediated PCR (DIPS-PCR) and molecular characterization in cervical cancer cells. Int J Cancer. 2001;92(1):9–17.View ArticlePubMedGoogle Scholar
- Li H, Yang Y, Zhang R, Cai Y, Yang X, Wang Z, Li Y, Cheng X, Ye X, Xiang Y, et al. Preferential sites for the integration and disruption of human papillomavirus 16 in cervical lesions. J Clin Virol. 2013;56(4):342–7.View ArticlePubMedGoogle Scholar
- Dutta S, Chakraborty C, Dutta AK, Mandal RK, Roychoudhury S, Basu P, Panda CK. Physical and methylation status of human papillomavirus 16 in asymptomatic cervical infections changes with malignant transformation. J Clin Pathol. 2015;68(3):206–11.View ArticlePubMedGoogle Scholar
- Mazumder Indra D, Singh RK, Mitra S, Dutta S, Chakraborty C, Basu PS, Mondal RK, Roychoudhury S, Panda CK. Genetic and epigenetic changes of HPV16 in cervical cancer differentially regulate E6/E7 expression and associate with disease progression. Gynecol Oncol. 2011;123(3):597–604.View ArticlePubMedGoogle Scholar
- Vinokurova S, Wentzensen N, Kraus I, Klaes R, Driesch C, Melsheimer P, Kisseljov F, Durst M, Schneider A, Von Knebel DM. Type-dependent integration frequency of human papillomavirus genomes in cervical lesions. Cancer Res. 2008;68(1):307–13.View ArticlePubMedGoogle Scholar
- Cheung JL, Cheung TH, Yu MY, Chan PK. Virological characteristics of cervical cancers carrying pure episomal form of HPV16 genome. Gynecol Oncol. 2013;131(2):374–9.View ArticlePubMedGoogle Scholar
- Boulet GA, Benoy IH, Depuydt CE, Horvath CA, Aerts M, Hens N, Vereecken AJ, Bogers JJ. Human papillomavirus 16 load and E2/E6 ratio in HPV16-positive women: biomarkers for cervical intraepithelial neoplasia > or = 2 in a liquid-based cytology setting? Cancer Epidemiol Biomarkers Prev. 2009;18(11):2992–9.View ArticlePubMedGoogle Scholar
- Hong D, Lu W, Ye F, Hu Y, Xie X. Gene silencing of HPV16 E6/E7 induced by promoter-targeting siRNA in SiHa cells. Br J Cancer. 2009;101(10):1798–804.View ArticlePubMedPubMed CentralGoogle Scholar
- Hong D, Ye F, Lu W, Hu Y, Wan X, Chen Y, Xie X. Methylation status of the long control region of HPV 16 in clinical cervical specimens. Mol Med Rep. 2008;1(4):555–60.PubMedGoogle Scholar
- Kalantari M, Blennow E, Hagmar B, Johansson B. Physical state of HPV16 and chromosomal mapping of the integrated form in cervical carcinomas. Diagn Mol Pathol. 2001;10(1):46–54.View ArticlePubMedGoogle Scholar
- Tan SH, Leong LE, Walker PA, Bernard HU. The human papillomavirus type 16 E2 transcription factor binds with low cooperativity to two flanking sites and represses the E6 promoter through displacement of Sp1 and TFIID. J Virol. 1994;68(10):6411–20.PubMedPubMed CentralGoogle Scholar
- Jeon S, Lambert PF. Integration of human papillomavirus type 16 DNA into the human genome leads to increased stability of E6 and E7 mRNAs: implications for cervical carcinogenesis. Proc Natl Acad Sci U S A. 1995;92(5):1654–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Pett M, Coleman N. Integration of high-risk human papillomavirus: a key event in cervical carcinogenesis? J Pathol. 2007;212(4):356–67.View ArticlePubMedGoogle Scholar
- Li W, Wang W, Si M, Han L, Gao Q, Luo A, Li Y, Lu Y, Wang S, Ma D. The physical state of HPV16 infection and its clinical significance in cancer precursor lesion and cervical carcinoma. J Cancer Res Clin Oncol. 2008;134(12):1355–61.View ArticlePubMedGoogle Scholar
- Gray E, Pett MR, Ward D, Winder DM, Stanley MA, Roberts I, Scarpini CG, Coleman N. In vitro progression of human papillomavirus 16 episome-associated cervical neoplasia displays fundamental similarities to integrant-associated carcinogenesis. Cancer Res. 2010;70(10):4081–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Hafner N, Driesch C, Gajda M, Jansen L, Kirchmayr R, Runnebaum IB, Durst M. Integration of the HPV16 genome does not invariably result in high levels of viral oncogene transcripts. Oncogene. 2008;27(11):1610–7.View ArticlePubMedGoogle Scholar
- Das D, Bhattacharjee B, Sen S, Mukhopadhyay I, Sengupta S. Association of viral load with HPV16 positive cervical cancer pathogenesis: causal relevance in isolates harboring intact viral E2 gene. Virology. 2010;402(1):197–202.View ArticlePubMedGoogle Scholar
- Marongiu L, Godi A, Parry JV, Beddows S. Human Papillomavirus 16, 18, 31 and 45 viral load, integration and methylation status stratified by cervical disease stage. BMC Cancer. 2014;14:384.View ArticlePubMedPubMed CentralGoogle Scholar
- Gillitzer E, Chen G, Stenlund A. Separate domains in E1 and E2 proteins serve architectural and productive roles for cooperative DNA binding. EMBO J. 2000;19(12):3069–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Yasugi T, Benson JD, Sakai H, Vidal M, Howley PM. Mapping and characterization of the interaction domains of human papillomavirus type 16 E1 and E2 proteins. J Virol. 1997;71(2):891–9.PubMedPubMed CentralGoogle Scholar
- Oliveira JG, Colf LA, McBride AA. Variations in the association of papillomavirus E2 proteins with mitotic chromosomes. Proc Natl Acad Sci U S A. 2006;103(4):1047–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu WG, Srinivasan K, Dai Z, Duan W, Druhan LJ, Ding H, Yee L, Villalona-Calero MA, Plass C, Otterson GA. Methylation of adjacent CpG sites affects Sp1/Sp3 binding and activity in the p21(Cip1) promoter. Mol Cell Biol. 2003;23(12):4056–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Thain A, Jenkins O, Clarke AR, Gaston K. CpG methylation directly inhibits binding of the human papillomavirus type 16 E2 protein to specific DNA sequences. J Virol. 1996;70(10):7233–5.PubMedPubMed CentralGoogle Scholar
- Kim K, Garner-Hamrick PA, Fisher C, Lee D, Lambert PF. Methylation patterns of papillomavirus DNA, its influence on E2 function, and implications in viral infection. J Virol. 2003;77(23):12450–9.View ArticlePubMedPubMed CentralGoogle Scholar