C33-A cells transfected with E6*I or E6*II the short forms of HPV-16 E6, displayed opposite effects on cisplatin-induced apoptosis

Long-read sequencing reveals oncogenic mechanism of HPV-human fusion transcripts in cervical cancer

Integration of high-risk human papillomavirus (HPV) into the host genome is a crucial event for the development of cervical cancer, however, the underlying mechanism of HPV integration-driven carcinogenesis remains unknown. Here, we performed long-read RNA sequencing on 12 high-grade squamous intraepithelial lesions (HSIL) and cervical cancer patients, including 3 pairs of cervical cancer and corresponding para-cancerous tissue samples to investigate the full-length landscape of cross-species genome integrations. In addition to massive unannotated isoforms, transcriptional regulatory events, and gene chimerism, more importantly, we found that HPV-human fusion events were prevalent in HPV-associated cervical cancers. Combined with the genome data, we revealed the existence of a universal transcription pattern in these fusion events, whereby structurally similar fusion transcripts were generated by specific splicing in E6 and a canonical splicing donor site in E1 linking to various human splicing acceptors. Highly expressed HPV-human fusion transcripts, eg, HPV16 E6*I-E7-E1_SD880-human gene, were the key driver of cervical carcinogenesis, which could trigger overexpression of E6*I and E7, and destroy the transcription of tumor suppressor genes CMAHP, TP63 and P3H2. Finally, evidence from in vitro and in vivo experiments demonstrates that the novel read-through fusion gene mRNA, E1-CMAHP (E1C, formed by the integration of HPV58 E1 with CMAHP), existed in the fusion transcript can promote malignant transformation of cervical epithelial cells via regulating downstream oncogenes to participate in various biological processes. Taken together, we reveal a previously unknown mechanism of HPV integration-driven carcinogenesis and provide a novel target for the diagnosis and treatment of cervical cancer.

Longitudinal follow-up of HPV16 sequence after cervical infection: Low intrahost variation and no correlation with clinical evolution

Human papillomavirus (HPV) 16 exhibits different variants that may differ in their carcinogenic risk. To identify some high-risk variants, we sequenced and compared HPV16 whole genomes obtained from a longitudinal cohort of 34 HPV16-infected women who had either spontaneously cleared their infection (clearance group or "C"), or developed cervical high-grade lesions following a viral persistence (group persistence or "P"). Phylogenetic analysis of paired samples obtained at the beginning (C0 or P0) and at the end (C2 or P2) of the follow-up (median intervals between C0-C2 and between P0-P2 were 16 and 36.5 months, respectively) revealed a low genetic variability within the host compared to the genetic interhost diversity. By comparing our HPV16 sequences to a reference sequence, we observed 301 different substitutions, more often transitions (60.9%) than transversions (39.1%), that occurred throughout the viral genome, but with a low frequency in E6 and E7 oncogenes (10 and 9 substitutions), suggesting a high conservation of these genes. Deletions and insertions were mostly observed in intergenic regions of the virus. The only significant substitution found between the subgroups C2 and P2 was observed in the L2 gene (L330F), with an unclear biological relevance. Our results suggest a low longitudinal intrahost evolution of HPV16 sequences and no correlation between genetic variations and clinical evolution.

New insights into the interactions of HPV-16 E6*I and E6*II with p53 isoforms and induction of apoptosis in cancer-derived cell lines

An important characteristic of cancers associated with high-risk human papillomaviruses (HR-HPV) is the inability of p53 to activate apoptosis due to the effect of the oncoprotein E6. However, the effect of HPV-16 E6 splice variant isoforms (namely E6*I and E6*II), their interaction with the existing p53 isoforms, and their influence on apoptosis is unclear. Here, we report the outcome of ectopic expression of HPV-16 E6, E6*I, and E6*II on the relative levels of p53 and p53 isoforms Δ40p53 and Δ133p53 and their interactions with these proteins. Additionally, we evaluated the effect of ectopic expression of p53, Δ40p53, and Δ133p53 on apoptosis in a p53 null pulmonary cell line (H1299) co-transfected with E6 isoforms and p53^+/+ cell lines with HR-HPV (SiHa and HeLa), transfected with p53 isoforms and treated with cisplatin, a conventional drug used to treat cervical cancer. Our results show that E6 and E6*II induced a significant decrease in p53, but only E6 triggered a Δ40p53 decrease and that E6*II interacts with p53 but not with Δ40p53 and Δ133p53. On the other hand, E6*I did not show any effect or interaction with the p53 isoforms. We found that apoptosis was elevated in H1299 cells transfected with p53 (p = 0.0001) and Δ40p53 (p = 0.0001). A weak apoptotic effect was observed when Δ133p53 was ectopically expressed (p = 0.0195). We observed that both p53 (p = 0.0006) and Δ40p53 (p = 0.0014) induced apoptosis in cisplatin-treated SiHa cells; however in cisplatin-treated HeLa cells, only p53 induced apoptosis (p = 0.0029). No significant differences in apoptosis were observed upon ectopic expression of p53, Δ40p53, and Δ133p53 in SiHa and HeLa cells. Our findings suggest a possible therapeutic application for the combining of p53 or Δ40p53 with cisplatin to induce an increased apoptosis of cancer cells expressing E6 isoforms from HPV-16.

Repurposing the Hybrid Capture 2 (HC2) screening test for whole-genome sequencing of human papillomaviruses

Simple and standardized approaches for genome analysis of human papillomavirus (HPV) by next-generation sequencing are needed. The aim of the study was to develop a protocol for direct deep sequencing of high-risk (hr) HPV strains, based on the widely used commercial Hybrid Capture 2 (QIAGEN) test, without any additional probe design. This protocol was applied to 15 HPV-positive and two HPV-negative cervical samples or cell lines and validated at the genotype level by comparing the sequencing results to those obtained using a commercial genotyping kit. The performance of our protocol, presented in this proof-of-principle study, supports its use for accurate characterization of genetic variants of hrHPV.

Extrachromosomal DNA in HPV-Mediated Oropharyngeal Cancer Drives Diverse Oncogene Transcription

PURPOSE

Human papillomavirus (HPV) plays a major role in oncogenesis and circular extrachromosomal DNA (ecDNA) is found in many cancers. However, the relationship between HPV and circular ecDNA in human cancer is not understood.

EXPERIMENTAL DESIGN

Forty-four primary tumor tissue samples were obtained from a cohort of patients with HPV-positive oropharynx squamous cell carcinoma (OPSCC). Twenty-eight additional HPV oropharyngeal cancer (HPVOPC) tumors from The Cancer Genome Atlas (TCGA) project were analyzed as a separate validation cohort. Genomic, transcriptomic, proteomic, computational, and functional analyses of HPVOPC were applied to these datasets.

RESULTS

Our analysis revealed circular, oncogenic DNA in nearly all HPVOPC, with circular human and human-viral hybrid ecDNA present in over a third of HPVOPC and viral circular DNA in remaining tumors. Hybrid ecDNA highly express fusion transcripts from HPV promoters and HPV oncogenes linked to downstream human transcripts that drive oncogenic transformation and immune evasion, and splice multiple, diverse human acceptors to a canonical SA880 viral donor site. HPVOPC have high E6*I expression with specific viral oncogene expression pattern related to viral or hybrid ecDNA composition.

CONCLUSIONS

Nonchromosomal circular oncogenic DNA is a dominant feature of HPVOPC, revealing an unanticipated link between HPV and ecDNA that leverages the power of extrachromosomal inheritance to drive HPV and somatic oncogene expression.

E6/E7 and E6* From HPV16 and HPV18 Upregulate IL-6 Expression Independently of p53 in Keratinocytes

Keratinocyte infection with high-risk human papillomavirus genotypes has been linked to cancer development. In cervix, the alpha HPV16 and HPV18 have been reported as the mayor causative agents of cervical cancer. Oncogenic progression and chronic inflammation are closely related processes, with IL-6 as one of the main pro-inflammatory cytokines involved. However, there are limited studies about the regulation of IL-6 by low and high risk HPVs and the HPV proteins implicated in this modulation. In this work, we report the overexpression of IL-6 in HPV infected cervical cancer derived cell lines (HeLa and SiHa) compared to non-tumorigenic keratinocytes (HaCaT), and in Cervical Intraepithelial Neoplasia grade 1 HPV16 and HPV18 positive cervical samples compared to HPV negative samples without lesions. Moreover, we generated HaCaT keratinocytes that express E5, E6, and E7 from high risk (16 or 18) or low risk (62 and 84) HPVs. E5 proteins do not modify IL-6 expression, while E7 modestly increase it. Interestingly, E6 proteins in HaCaT cells upregulate IL-6 mRNA expression and protein secretion. Indeed, in HaCaT cells that express high risk HPV16E6 or HPV18E6 proteins, only the truncated E6^* isoforms were expressed, showing the stronger IL-6 overexpression, while in HaCaT cells that express low risk HPV62 and HPV84 full length E6 proteins, IL-6 was also upregulated but not so drastically. Since HaCaT cells have a mutated p53 form that is not degraded by the introduction of E6 or E6/E7, it seems that E6/E7 regulate IL-6 by an additional mechanism independent of p53. In addition, basal keratinocytes showed a strong expression of IL-6R using immunohistochemistry, suggesting an autocrine mechanism over proliferative cells. Altogether, IL-6 cytokine expression in keratinocytes is upregulated by E6 and E7 proteins from HPVs 16, 18, 62, and 84, especially by high risk HPV16 and HPV18 E6^*, which may contribute to promote a pro-inflammatory and highly proliferative microenvironment that can persist over time and lead to cervical tumorigenesis.