Genetic variation within the human papillomavirus type 16 genome is associated with oropharyngeal cancer prognosis

Purpose: A significant barrier to adoption of de-escalated treatment protocols for human papillomavirus-driven oropharyngeal cancer (HPV-OPC) is that few predictors of poor prognosis exist. We conducted the first large whole-genome sequencing (WGS) study to characterize the genetic variation of the HPV16 genome and to evaluate its association with HPV-OPC patient survival. Patients and Methods: 460 OPCs from 2 large US medical centers (1980-2017) underwent HPV16 WGS. Site-specific variable positions (SNPs) across the HPV16 genome were identified. Cox proportional hazards models estimated hazard ratios (HRs) and 95% confidence intervals (CIs) for overall survival by HPV16 SNPs. Harrell C-index and time-dependent positive predictive value (PPV) curves and areas under the PPV curves were used to evaluate the predictive accuracy of HPV16 SNPs for overall survival. Results: 384 OPCs (83.48%) passed quality control filters with sufficient depth and coverage of HPV16 genome sequencing to be analyzed. 284 HPV16 SNPs with a minor allele frequency >1% were identified. Eight HPV16 SNPs were significantly associated with worse survival after false discovery rate (FDR) correction (individual prevalence:1.0%-5.5%; combined prevalence: 15.10%); E1 gene position 1053 (HR for overall survival [HRos]:3.75,95%CI:1.77-7.95;Pfdr=0.0099); L2 gene positions 4410 (HRos:5.32,95%CI:1.91-14.81;Pfdr=0.0120), 4539 (HRos:6.54,95%CI:2.03-21.08;Pfdr=0.0117); 5050 (HRos:6.53,95%CI:2.34-18.24;Pfdr=0.0030) and 5254 (HRos:7.76,95%CI:2.41-24.98;Pfdr=0.0030); and L1 gene positions 5962 (HRos:4.40,95%CI:1.88-10.31;Pfdr=0.0110) and 6025 (HRos:5.71,95%CI:2.43-13.41;Pfdr=0.0008) and position 7173 within the upstream regulatory region (HRos:9.90,95%CI:3.05-32.12;Pfdr=0.0007). Median survival time for patients with ≥1 high-risk HPV16 SNPs was 3.96 years compared to 18.67 years for patients without a high-risk SNP; log-rank test P<0.001. HPV16 SNPs significantly improved the predictive accuracy for overall survival above traditional factors (age, smoking, stage, treatment): increase in C-index was 0.069 (95% CI: 0.019-0.119, P <0.001); increase in area under the PPV curve for predicting 5-year survival was 0.068 (95%CI: 0.015-0.111, P =0.008). Conclusions: HPV16 genetic variation is associated with HPV-OPC prognosis and can improve prognostic accuracy.

CTCF association with episomal HPV16 genomes regulates viral oncogene transcription and splicing

Human papillomavirus 16 (HPV16) is a high-risk alphapapillomavirus that is associated with cancers of mucosal epithelia. The virus genome exists in cells as an episome but can integrate and overexpress the E6 and E7 viral oncogenes. In related high-risk family members HPV18 and HPV31, host proteins including CTCF, an insulator, and SMC1A, a component of the cohesion complex, are known to interact with the viral genome and alter transcriptional activity, splicing patterns and episome amplification. However, the roles of these two proteins during HPV16 infection has not yet been fully examined. Here, we show during differentiation of the episomal HPV16-containing W12 cell line that CTCF association increases with the virus genome at the known E2 binding site, whilst additional CTCF binding now occurs at the putative L2 binding site, with SMC1A association occurring unchanged here. While expression of virus late transcripts (E4^L1, L2, L1) is stimulated, early transcript levels decrease by 48 hours, with the exception of the E6*IV spliced transcript. Conversely, in undifferentiated, monolayer W12 cells, CTCF knockdown increases the level of all early transcripts, whereas E6*IV level increases. Additionally, CTCF ablation as well as SMC1A knockdown results in decreases to HPV16 genome copy number. Taken together, this supports the model that while CTCF and SMC1A have a role in HPV16 genome maintenance, CTCF plays a greater part in regulating HPV16 oncogene splicing and expression during the natural lifecycle of the virus, and may be involved in a reduced risk of cancer development during episomal HPV16 infections.

Human Papillomavirus 16 E6 and E7 Synergistically Repress Innate Immune Gene Transcription

ABSTRACT Human papillomaviruses (HPV) are causative agents in 5% of all cancers, including the majority of anogenital and oropharyngeal cancers. Downregulation of innate immune genes (IIGs) by HPV to promote the viral life cycle is well documented; E6 and E7 are known repressors of these genes. More recently, we demonstrated that E2 could also repress IIGs. These studies have been carried out in cells overexpressing the viral proteins, and to further investigate the role of individual viral proteins in this repression, we introduced stop codons into E6 and/or E7 in the entire HPV16 genome and generated N/Tert-1 cells stably maintaining the HPV16 genomes. We demonstrate that E6 or E7 individually is not sufficient to repress IIG expression in the context of the entire HPV16 genome; both are required for a synergistic repression. The DNA damage response (DDR) is activated by HPV16 irrespective of E6 and E7 expression, presumably due to viral replication; E1 is a known activator of the DDR. In addition, replication stress was apparent in HPV16-positive cells lacking E6 and E7, manifested by attenuated cellular growth and activation of replication stress genes. These studies led us to the following model. Viral replication per se can activate the DDR following infection, and this activation is a known inducer of IIG expression, which may induce cellular senescence. To combat this, E6 and E7 synergistically combine to manipulate the DDR and actively repress innate immune gene expression promoting cellular growth; neither protein by itself is able to do this. IMPORTANCE The role of human papillomavirus 16 (HPV16) in human cancers is well established; however, to date there are no antiviral therapeutics that are available for combatting these cancers. To identify such targets, we must enhance the understanding of the viral life cycle. Innate immune genes (IIGs) are repressed by HPV16, and we have reported that this repression persists through to cancer. Reversal of this repression would boost the immune response to HPV16-positive tumors, an area that is becoming more important given the advances in immunotherapy. This report demonstrates that E6 and E7 synergistically repress IIG expression in the context of the entire HPV16 genome. Removal of either protein activates the expression of IIGs by HPV16. Therefore, gaining a precise understanding of how the viral oncogenes repress IIG expression represents an opportunity to reverse this repression and boost the immune response to HPV16 infections for therapeutic gain.

Human papillomavirus 16 E6 and E7 synergistically repress innate immune gene transcription

Human papillomaviruses are causative agents in 5% of all cancers, including the majority of anogenital and oropharyngeal cancers. Downregulation of innate immune genes (IIGs) by HPV to promote the viral life cycle is well documented; E6 and E7 are known repressors of these genes. More recently we demonstrated that E2 could also repress IIGs. These studies have been carried out in cells over-expressing the viral proteins and to further investigate the role of individual viral proteins in this repression we introduced stop codons into E6 and/or E7 in the entire HPV16 genome and generated N/Tert-1 cells stably maintaining the HPV16 genomes. We demonstrate that E6 or E7 individually are not sufficient to repress IIG expression in the context of the entire HPV16 genome, both are required for a synergistic repression. The DNA damage response (DDR) is activated by HPV16 irrespective of E6 and E7 expression, presumably due to viral replication; E1 is a known activator of the DDR. In addition, replication stress was apparent in the HPV16 positive cells lacking E6 and E7, manifested by attenuated cellular growth and activation of replication stress genes. These studies lead us to the following model. Viral replication per se can activate the DDR following infection, and this activation is a known inducer of IIG expression which could induce cellular senescence. To combat this, E6 and E7 synergistically combine to manipulate the DDR and actively repress innate immune gene expression promoting cellular growth; neither protein by itself is able to do this.ImportanceThe role of HPV16 in human cancers is well established; however, to date there are no anti-viral therapeutics that are available for combatting these cancers. To identify such targets, we must enhance understanding of the viral life cycle. Innate immune genes (IIGs) are repressed by HPV16, and we have reported that this repression persists through to cancer. Reversal of this repression would boost the immune response to HPV16 positive tumors, an area that is becoming more important given the advances in immunotherapy. This report demonstrates that E6 and E7 synergistically repress IIG expression in the context of the entire HPV16 genome. Removal of either protein activates the expression of IIGs by HPV16. Therefore, gaining a precise understanding of how the viral oncogenes repress IIG expression represents an opportunity to reverse this repression and boost the immune response to HPV16 infections for therapeutic gain.

Activation of NF-κB by Human Papillomavirus 16 E1 Limits E1-Dependent Viral Replication through Degradation of E1

ABSTRACTNF-κB is a family of transcription factors that regulate gene expression involved in many processes, such as the inflammatory response and cancer progression. Little is known about associations of NF-κB with the human papillomavirus (HPV) life cycle. We have developed a tissue culture system to conditionally induce E1-dependent replication of the human papillomavirus 16 (HPV16) genome in human cervical keratinocytes and found that expression of HPV16 E1, a viral helicase, results in reduction of IκBα and subsequent activation of NF-κB in a manner dependent on helicase activity. Exogenous expression of a degradation-resistant mutant of IκBα, which inhibits the activation of NF-κB, enhanced E1-dependent replication of the viral genome. Wortmannin, a broad inhibitor of phosphoinositide 3-kinases (PI3Ks), and, to a lesser extent, VE-822, an ATR kinase inhibitor, but not KU55933, an ATM kinase inhibitor, suppressed the activation of NF-κB and augmented E1-dependent replication of the HPV16 genome. Interestingly, the enhancement of E1-dependent replication of the viral genome was associated with increased stability of E1 in the presence of wortmannin as well as the IκBα mutant. Collectively, we propose that expression of E1 induces NF-κB activation at least in part through the ATR-dependent DNA damage response and that NF-κB in turn limits E1-dependent replication of HPV16 through degradation of E1, so that E1 and NF-κB may constitute a negative feedback loop.IMPORTANCEA major risk factor in human papillomavirus (HPV)-associated cancers is persistent infection with high-risk HPVs. To eradicate viruses from infected tissue, it is important to understand molecular mechanisms underlying the establishment and maintenance of persistent infection. In this study, we obtained evidence that human papillomavirus 16 (HPV16) E1, a viral DNA helicase essential for amplification of the viral genomes, induces NF-κB activation and that this limits E1-dependent genome replication of HPV16. These results suggest that NF-κB mediates a negative feedback loop to regulate HPV replication and that this feedback loop could be associated with control of the viral copy numbers. We could thus show for the first time that NF-κB activity is involved in the establishment and maintenance of persistent HPV infection.