Direct Comparison of HPV16 Viral Genomic Integration, Copy Loss, and Structural Variants in Oropharyngeal and Uterine Cervical Cancers Reveal Distinct Relationships to E2 Disruption and Somatic Alteration

Squamous cell carcinoma of the oropharynx caused by HPV type 16 (HPV16+ OPSCC) is the most common HPV-associated malignancy in the USA and has many molecular differences from uterine cervical squamous cell carcinoma (UCSCC). Our understanding of HPV oncogenesis relied on studies of UCSCC revealing a consensus model reliant on HPV integration with a loss of E2. Here, we compare patterns of HPV integration in UCSCC and OPSCC by analysis of affinity capture sequencing of the HPV16 genome in 104 OPSCC and 44 UCSCC tumors. These cohorts were contemporaneously sequenced using an identical strategy. Integration was identified using discordant read pair clustering and assembly-based approaches. Viral integration sites, structural variants, and copy losses were examined. While large-scale deep losses of HPV16 genes were common in UCSCC and were associated with E2 loss, deep copy losses of the HPV16 genome were infrequent in HPV16+ OPSCC. Similarly, structural variants within HPV16 favored E2 loss in UCSCC but not OPSCC. HPV16 integration sites were non-random, with recurrent integration hot-spots identified. OPSCC tumors had many more integration sites per tumor when compared to UCSCC and had more integration sites in genomic regions with high gene density. These data show that viral integration and E2 disruption are distinct in UCSCC and OPSCC. Our findings also add to growing literature suggesting that HPV tumorigenesis in OPSCC does not follow the model developed based on UCSCC

Systematic analysis of SARS-CoV-2 infection of an ACE2-negative human airway cell

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) variants govern transmissibility, responsiveness to vaccination, and disease severity. In a screen for new models of SARS-CoV-2 infection, we identify human H522 lung adenocarcinoma cells as naturally permissive to SARS-CoV-2 infection despite complete absence of angiotensin-converting enzyme 2 (ACE2) expression. Remarkably, H522 infection requires the E484D S variant; viruses expressing wild-type S are not infectious. Anti-S monoclonal antibodies differentially neutralize SARS-CoV-2 E484D S in H522 cells as compared to ACE2-expressing cells. Sera from vaccinated individuals block this alternative entry mechanism, whereas convalescent sera are less effective. Although the H522 receptor remains unknown, depletion of surface heparan sulfates block H522 infection. Temporally resolved transcriptomic and proteomic profiling reveal alterations in cell cycle and the antiviral host cell response, including MDA5-dependent activation of type I interferon signaling. These findings establish an alternative SARS-CoV-2 host cell receptor for the E484D SARS-CoV-2 variant, which may impact tropism of SARS-CoV-2 and consequently human disease pathogenesis

Hallmark signatures of mutant tumors.

Prior to the next generation sequencing and characterization of the tumor genome landscape, mutations in the SWI/SNF chromatin remodeling complex and the KEAP1-NRF2 signaling pathway were underappreciated. While these two classes of mutations appeared to independently contribute to tumor development, recent reports have demonstrated a mechanistic link between these two regulatory mechanisms in specific cancer types and cell models. In this work, we expand upon these data by exploring the relationship between mutations in BAF and PBAF subunits of the SWI/SNF complex and activation of NRF2 signal transduction across many cancer types. ARID1A/B mutations were strongly associated with NRF2 transcriptional activity in head and neck squamous carcinomas (HNSC). Many additional tumor types showed significant association between NRF2 signaling and mutation of specific components of the SWI/SNF complex. Different effects of BAF and PBAF mutations on the polarity of NRF2 signaling were observed. Overall, our results support a context-dependent functional link between SWI/SNF and NRF2 mutations across human cancers and implicate ARID1A inactivation in HPV-negative HNSC in promoting tumor progression and survival through activation of the KEAP1-NRF2 signaling pathway. The tumor-specific effects of these mutations open a new area of study for how mutations in the KEAP1-NRF2 pathway and the SWI/SNF complex contribute to cancer.</div

SWI/SNF and KEAP1-NRF2 are frequently mutated in human cancers.

A. Frequency of coincident alterations of SWI/SNF subunits and KEAP1-NRF2-CUL3 in TCGA tumor types. B. Frequent alterations in SWI/SNF subunits in tumors with ≥ frequencies of KEAP1-NRF2 alterations. C. Heatmap demonstrating that ARID1A is the most frequently altered SWI/SNF subunit across TCGA tumors.</p

NRF2 gene signatures used in this study.

Prior to the next generation sequencing and characterization of the tumor genome landscape, mutations in the SWI/SNF chromatin remodeling complex and the KEAP1-NRF2 signaling pathway were underappreciated. While these two classes of mutations appeared to independently contribute to tumor development, recent reports have demonstrated a mechanistic link between these two regulatory mechanisms in specific cancer types and cell models. In this work, we expand upon these data by exploring the relationship between mutations in BAF and PBAF subunits of the SWI/SNF complex and activation of NRF2 signal transduction across many cancer types. ARID1A/B mutations were strongly associated with NRF2 transcriptional activity in head and neck squamous carcinomas (HNSC). Many additional tumor types showed significant association between NRF2 signaling and mutation of specific components of the SWI/SNF complex. Different effects of BAF and PBAF mutations on the polarity of NRF2 signaling were observed. Overall, our results support a context-dependent functional link between SWI/SNF and NRF2 mutations across human cancers and implicate ARID1A inactivation in HPV-negative HNSC in promoting tumor progression and survival through activation of the KEAP1-NRF2 signaling pathway. The tumor-specific effects of these mutations open a new area of study for how mutations in the KEAP1-NRF2 pathway and the SWI/SNF complex contribute to cancer.</div

TCGA bioinformatic pipeline.

TCGA Firehose was used to download CNV, mutation, and gene expression data. CNV and mutation data were used to identify mutant tumors in genes for each group. Confounding mutations were then filtered from the treatment group. Gene expression data were used to calculate signal-to-noise ratios within each group. GSEA was performed using signal-to-noise ratios between control and treated groups for NRF2 gene signatures and Oncogene and Hallmark signatures.</p

ARID1A and ARID1B mutant tumors alter KEAP1-NRF2 signatures.

GSEA of five NRF2 signatures was used. Normalized Enrichment Scores (NES) are shown in the heat map with Tumors (S2 Table) on the x-axis and Signatures (S3 Table) on the y-axis. A. NRF2 signature GSEA results for KEAP1-NRF2 mutant tumors (positive control). B. NRF2 signature GSEA results for ARID1A mutant tumors C. NRF2 signature GSEA results for ARID1B mutant tumors. D. NRF2 signature GSEA results for HNSC tumors comparing the CPTAC and TCGA datasets.</p