mRNA sequencing of novel cell lines from human papillomavirus type-16 related vulval intraepithelial neoplasia: consequences of expression of HPV16 E4 and E5

HPV DNA Integration at Actionable Cancer-Related Genes Loci in HPV-Associated Carcinomas

In HPV-associated carcinomas, some examples of cancer-related genes altered by viral insertion and corresponding to potential therapeutic targets have been described, but no quantitative assessment of these events, including poorly recurrent targets, has been reported to date. To document these occurrences, we built and analyzed a database comprised of 1455 cases, including HPV genotypes and tumor localizations. Host DNA sequences targeted by viral integration were classified as "non-recurrent" (one single reported case; 838 loci), "weakly recurrent" (two reported cases; 82 loci), and highly recurrent (≥3 cases; 43 loci). Whereas the overall rate of cancer-related target genes was 3.3% in the Gencode database, this rate increased to 6.5% in "non-recurrent", 11.4% in "weakly recurrent", and 40.1% in "highly recurrent" genes targeted by integration (p = 4.9 × 10-4). This rate was also significantly higher in tumors associated with high-risk HPV16/18/45 than other genotypes. Among the genes targeted by HPV insertion, 30.2% corresponded to direct or indirect druggable targets, a rate rising to 50% in "highly recurrent" targets. Using data from the literature and the DepMap 23Q4 release database, we found that genes targeted by viral insertion could be new candidates potentially involved in HPV-associated oncogenesis. A more systematic characterization of HPV/host fusion DNA sequences in HPV-associated cancers should provide a better knowledge of HPV-driven carcinogenesis and favor the development of personalize patient treatments.

Integrative analysis of gene expression and DNA methylation to identify biomarkers of non-genital warts induced by low-risk human papillomaviruses infection

Background

Human papillomaviruses have been shown to dysregulate the gene expression and DNA methylation profiles of their host cells over the course of infection. However, there is a lack of information on the impact of low-risk HPV infection and wart formation on host cell's expression and methylation patterns. Therefore, the objective of this study is to analyse the genome and methylome of common warts using an integrative approach.

Methods

In the present study, gene expression (GSE136347) and methylation (GSE213888) datasets of common warts were obtained from the GEO database. Identification of the differentially expressed and differentially methylated genes was carried out using the RnBeads R package and the edgeR Bioconductor package. Next, functional annotation of the identified genes was obtained using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). Network construction and analyses of the gene-gene, protein-protein, and signaling interactions of the differentially expressed and differentially methylated genes was performed using the GeneMANIA web interface, the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database, and the Signaling Network Open Resource 2.0 (SIGNOR 2.0), respectively. Lastly, significant hub genes were identified using the Cytoscape application CytoHubba.

Results

A total of 276 genes were identified as differentially expressed and differentially methylated in common warts, with 52% being upregulated and hypermethylated. Functional enrichment analysis identified extracellular components as the most enriched annotations, while network analyses identified ELN, ITGB1, TIMP1, MMP2, LGALS3, COL1A1 and ANPEP as significant hub genes.

Conclusions

To the best knowledge of the authors, this is the first integrative study to be carried out on non-genital warts induced by low-risk HPV types. Future studies are required to re-validate the findings in larger populations using alternative approaches.

The Major Constituent of Green Tea, Epigallocatechin-3-Gallate (EGCG), Inhibits the Growth of HPV18-Infected Keratinocytes by Stimulating Proteasomal Turnover of the E6 and E7 Oncoproteins

Epigallocatechin-3-gallate (EGCG), the primary bioactive polyphenol in green tea, has been shown to inhibit the growth of human papilloma virus (HPV)-transformed keratinocytes. Here, we set out to examine the consequences of EGCG treatment on the growth of HPV18-immortalised foreskin keratinocytes (HFK-HPV18) and an authentic HPV18-positive vulvar intraepithelial neoplasia (VIN) clone, focusing on its ability to influence cell proliferation and differentiation and to impact on viral oncogene expression and virus replication. EGCG treatment was associated with degradation of the E6 and E7 oncoproteins and an upregulation of their associated tumour suppressor genes; consequently, keratinocyte proliferation was inhibited in both monolayer and organotypic raft culture. While EGCG exerted a profound effect on cell proliferation, it had little impact on keratinocyte differentiation. Expression of the late viral protein E4 was suppressed in the presence of EGCG, suggesting that EGCG was able to block productive viral replication in differentiating keratinocytes. Although EGCG did not alter the levels of E6 and E7 mRNA, it enhanced the turnover of the E6 and E7 proteins. The addition of MG132, a proteasome inhibitor, to EGCG-treated keratinocytes led to the accumulation of the E6/E7 proteins, showing that EGCG acts as an anti-viral, targeting the E6 and E7 proteins for proteasome-mediated degradation.

Proportion of transcriptionally active DNA virus integrants: a meta-analysis

Oncoviruses are collectively responsible for over 1,000,000 new cases of cancer per year; some can integrate into the host's chromosomes. The present work was aimed at assessing the proportion of transcriptionally active viral integrants through a systematic review of the scientific publications present on the MedLine database. From the articles screened, 628 viral integrants overall were retrieved, of which 530.84 were transcriptionally active (84.53%); among the clinical samples, 264 of 323 integrants were active (81.73%). The causes for the silencing were not addressed in the articles analyzed. These findings might highlight a possible risk factor for the insurgence of cancer since some oncovirus integrants could be reactivated by stimuli of disparate nature. Further studies should address such possibility.

Human Papillomavirus DNA Methylation Predicts Response to Treatment Using Cidofovir and Imiquimod in Vulval Intraepithelial Neoplasia 3

Purpose: Response rates to treatment of vulval intraepithelial neoplasia (VIN) with imiquimod and cidofovir are approximately 57% and 61%, respectively. Treatment is associated with significant side effects and, if ineffective, risk of malignant progression. Treatment response is not predicted by clinical factors. Identification of a biomarker that could predict response is an attractive prospect. This work investigated HPV DNA methylation as a potential predictive biomarker in this setting.Experimental Design: DNA from 167 cases of VIN 3 from the RT3 VIN clinical trial was assessed. HPV-positive cases were identified using Greiner PapilloCheck and HPV 16 type-specific PCR. HPV DNA methylation status was assessed in three viral regions: E2, L1/L2, and the promoter, using pyrosequencing.Results: Methylation of the HPV E2 region was associated with response to treatment. For cidofovir (n = 30), median E2 methylation was significantly higher in patients who responded (P ≤ 0.0001); E2 methylation >4% predicted response with 88.2% sensitivity and 84.6% specificity. For imiquimod (n = 33), median E2 methylation was lower in patients who responded to treatment (P = 0.03; not significant after Bonferroni correction); E2 methylation <4% predicted response with 70.6% sensitivity and 62.5% specificity.Conclusions: These data indicate that cidofovir and imiquimod may be effective in two biologically defined groups. HPV E2 DNA methylation demonstrated potential as a predictive biomarker for the treatment of VIN with cidofovir and may warrant investigation in a biomarker-guided clinical trial. Clin Cancer Res; 23(18); 5460-8. ©2017 AACR.

Functional variants of human papillomavirus type 16 demonstrate host genome integration and transcriptional alterations corresponding to their unique cancer epidemiology

BACKGROUND

Human papillomaviruses (HPVs) are a worldwide burden as they are a widespread group of tumour viruses in humans. Having a tropism for mucosal tissues, high-risk HPVs are detected in nearly all cervical cancers. HPV16 is the most common high-risk type but not all women infected with high-risk HPV develop a malignant tumour. Likely relevant, HPV genomes are polymorphic and some HPV16 single nucleotide polymorphisms (SNPs) are under evolutionary constraint instigating variable oncogenicity and immunogenicity in the infected host.

RESULTS

To investigate the tumourigenicity of two common HPV16 variants, we used our recently developed, three-dimensional organotypic model reminiscent of the natural HPV infectious cycle and conducted various "omics" and bioinformatics approaches. Based on epidemiological studies we chose to examine the HPV16 Asian-American (AA) and HPV16 European Prototype (EP) variants. They differ by three non-synonymous SNPs in the transforming and virus-encoded E6 oncogene where AAE6 is classified as a high- and EPE6 as a low-risk variant. Remarkably, the high-risk AAE6 variant genome integrated into the host DNA, while the low-risk EPE6 variant genome remained episomal as evidenced by highly sensitive Capt-HPV sequencing. RNA-seq experiments showed that the truncated form of AAE6, integrated in chromosome 5q32, produced a local gene over-expression and a large variety of viral-human fusion transcripts, including long distance spliced transcripts. In addition, differential enrichment of host cell pathways was observed between both HPV16 E6 variant-containing epithelia. Finally, in the high-risk variant, we detected a molecular signature of host chromosomal instability, a common property of cancer cells.

CONCLUSIONS

We show how naturally occurring SNPs in the HPV16 E6 oncogene cause significant changes in the outcome of HPV infections and subsequent viral and host transcriptome alterations prone to drive carcinogenesis. Host genome instability is closely linked to viral integration into the host genome of HPV-infected cells, which is a key phenomenon for malignant cellular transformation and the reason for uncontrolled E6 oncogene expression. In particular, the finding of variant-specific integration potential represents a new paradigm in HPV variant biology.

Genomic characterization of viral integration sites in HPV-related cancers

Persistent infection with carcinogenic human papillomaviruses (HPV) causes the majority of anogenital cancers and a subset of head and neck cancers. The HPV genome is frequently found integrated into the host genome of invasive cancers. The mechanisms of how it may promote disease progression are not well understood. Thoroughly characterizing integration events can provide insights into HPV carcinogenesis. Individual studies have reported limited number of integration sites in cell lines and human samples. We performed a systematic review of published integration sites in HPV-related cancers and conducted a pooled analysis to formally test for integration hotspots and genomic features enriched in integration events using data from the Encyclopedia of DNA Elements (ENCODE). Over 1,500 integration sites were reported in the literature, of which 90.8% (N = 1,407) were in human tissues. We found 10 cytobands enriched for integration events, three previously reported ones (3q28, 8q24.21 and 13q22.1) and seven additional ones (2q22.3, 3p14.2, 8q24.22, 14q24.1, 17p11.1, 17q23.1 and 17q23.2). Cervical infections with HPV18 were more likely to have breakpoints in 8q24.21 (p = 7.68 × 10(-4) ) than those with HPV16. Overall, integration sites were more likely to be in gene regions than expected by chance (p = 6.93 × 10(-9) ). They were also significantly closer to CpG regions, fragile sites, transcriptionally active regions and enhancers. Few integration events occurred within 50 Kb of known cervical cancer driver genes. This suggests that HPV integrates in accessible regions of the genome, preferentially genes and enhancers, which may affect the expression of target genes.

Surface modification of tri-calcium phosphate nanoparticles by DOPE and/or anti-E6 antibody to enhance uptake of antisense of E6 mRNA

The main aim of this study was to evaluate the uptake of E6 mRNA antisense into cervical cancer cells, induced by human papilloma virus (HPV). In this study, the carrier of the antisense was tri-calcium phosphate nanoparticles (TCP NPs) conjugated with dioleoyl phosphatidyl ethanolamine (DOPE) and/or anti-E6 antibody. At first, TCP NPs were synthesized, coated with carboxy-polyethylene glycol, and then conjugated with anti-E6 antibody and/or DOPE by carbodiimide cross-linker. Then, a single stranded DNA, which was complementary (antisense) of E6 mRNA, was attached to each one. Finally, the uptake of conjugated and unconjugated TCP NPs into HelaS3 cells was separately evaluated by Fourier transform infrared spectroscopy, optical microscopy, and fluorescent microscopy. Also, the cytotoxicity of these carriers was measured by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide (MTT) assay. Overall, 4 types of TCP NPs were used in this study, including 1) TCP NPs conjugated with DOPE (TCP NPs/DOPE), 2) TCP NPs conjugated with DOPE and antibody (TCP NPs/DOPE/Anti-E6 Ab), 3) TCP NPs conjugated with antibody (TCP NPs/Anti-E6 Ab), and 4) TCP NPs which not conjugated with DOPE and antibody (unconjugated TCP NPs). Uptake tests showed that although all types of TCP NPs could transfer antisense of E6 mRNA into HelaS3 cells, TCP NPs/DOPE and TCP NPs/DOPE/Anti-E6 Ab had more uptake than TCP NPs/Anti-E6 Ab and unconjugated TCP NPs. Moreover, MTT assay showed that TCP NPs/DOPE was more toxic than TCP NPs/DOPE/Anti-E6 Ab, TCP NPs/Anti-E6 Ab, and unconjugated TCP NPs. It can be concluded that TCP NPs/DOPE/Anti-E6 Ab is a good choice for oligonucleotide delivery, because of higher uptake and less toxicity, compared with other formulations.

Increased methylation of Human Papillomavirus type 16 DNA correlates with viral integration in Vulval Intraepithelial Neoplasia

BACKGROUND

Methylation of HPV16 DNA is a promising biomarker for triage of HPV positive cervical screening samples but the biological basis for the association between HPV-associated neoplasia and increased methylation is unclear.

OBJECTIVES

To determine whether HPV16 DNA methylation was associated with viral integration, and investigate the relationships between viral DNA methylation, integration and gene expression.

STUDY DESIGN

HPV16 DNA methylation, integration and gene expression were assessed using pyrosequencing, ligation-mediated PCR and QPCR, in biopsies from 25 patients attending a specialist vulval neoplasia clinic and in short-term clonal cell lines derived from vulval and vaginal neoplasia.

RESULTS

Increased methylation of the HPV16 L1/L2 and E2 regions was associated with integration of viral DNA into the host genome. This relationship was observed both in vivo and in vitro. Increased methylation of E2 binding sites did not appear to be associated with greater expression of viral early genes. Expression of HPV E6 and E7 did not correlate with either integration state or increased L1/L2 methylation.

CONCLUSIONS

The data suggest that increased HPV DNA methylation may be partly attributable to viral integration, and provide a biological rationale for quantification of L1/L2 methylation in triage of HPV positive cervical screening samples.