HPV+ oropharyngeal squamous cell carcinomas from patients with two tumors display synchrony of viral genomes yet discordant mutational profiles and signatures

Double trouble: Synchronous and metachronous primaries confound ctHPVDNA monitoring

BACKGROUND

Human papillomavirus-associated head and neck squamous cell carcinoma (HPV + HNSCC) occurs in the oropharynx (HPV + OPSCC), sinonasal cavity (HPV + SNSCC), and nasopharynx (HPV + NPC). Circulating tumor HPV DNA (ctHPVDNA) is an accurate tool for diagnosis, treatment monitoring, and recurrence detection. An emerging challenge with ctHPVDNA is that ~7.4% of HPV + HNSCC patients develop synchronous or metachronous HPV+ primaries, which could confound ctHPVDNA monitoring.

METHODS

We describe a 65-year-old patient with T2N1M0 HPV16 + OPSCC and a 55-year-old patient with T2N2M0 HPV16 + OPSCC. Both patients were enrolled in our prospective observational ctHPVDNA study with longitudinal blood collections throughout treatment. Both patients developed multiple HPV+ primaries.

RESULTS

Detailed discussion of the patients' treatment courses, the subsequent diagnoses of their second HPV+ primaries, and their ctHPVDNA monitoring is presented.

CONCLUSIONS

As ctHPVDNA use becomes more prevalent, it is important to recognize that an increase in ctHPVDNA can come not only from the primary tumor or metastatic clones, but also from synchronous or metachronous second primaries.

Single-cell transcriptomic profiling for inferring tumor origin and mechanisms of therapeutic resistance

Head and Neck Squamous Cell Carcinoma (HNSCC) is an aggressive epithelial cancer with poor overall response rates to checkpoint inhibitor therapy (CPI) despite CPI being the recommended treatment for recurrent or metastatic HNSCC. Mechanisms of resistance to CPI in HNSCC are poorly understood. To identify drivers of response and resistance to CPI in a unique patient who was believed to have developed three separate HNSCCs, we performed single-cell RNA-seq (scRNA-seq) profiling of two responding lesions and one progressive lesion that developed during CPI. Our results not only suggest interferon-induced APOBEC3-mediated acquired resistance as a mechanism of CPI resistance in the progressing lesion but further, that the lesion in question was actually a metastasis as opposed to a new primary tumor, highlighting the immense power of scRNA-seq as a clinical tool for inferring tumor origin and mechanisms of therapeutic resistance.

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 HPV type 16 (HPV16) genome and to evaluate its association with HPV-OPC patient survival.

PATIENTS AND METHODS

A total of 460 OPC tumor specimens from two large United States medical centers (1980-2017) underwent HPV16 whole-genome sequencing. Site-specific variable positions [single nucleotide polymorphisms (SNPs)] across the HPV16 genome were identified. Cox proportional hazards model 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

A total of 384 OPC tumor specimens (83.48%) passed quality control filters with sufficient depth and coverage of HPV16 genome sequencing to be analyzed. Some 284 HPV16 SNPs with a minor allele frequency ≥1% were identified. Eight HPV16 SNPs were significantly associated with worse survival after false discovery rate 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 with 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.

APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma

APOBEC is a mutagenic source in human papillomavirus (HPV)-mediated malignancies, including HPV+ oropharyngeal squamous cell carcinoma (HPV + OPSCC), and in HPV genomes. It is unknown why APOBEC mutations predominate in HPV + OPSCC, or if the APOBEC-induced mutations observed in both human cancers and HPV genomes are directly linked. We performed sequencing of host somatic exomes, transcriptomes, and HPV16 genomes from 79 HPV + OPSCC samples, quantifying APOBEC mutational burden and activity in both host and virus. APOBEC was the dominant mutational signature in somatic exomes. In viral genomes, there was a mean of five (range 0-29) mutations per genome. The mean of APOBEC mutations in viral genomes was one (range 0-5). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-variant allele fraction mutations, suggesting that APOBEC mutagenesis actively occurrs in viral genomes during infection. HPV16 APOBEC-induced mutation patterns in OPSCC were similar to those previously observed in cervical samples. Paired host and viral analyses revealed that APOBEC-enriched tumor samples had higher viral APOBEC mutation rates (p = 0.028), and APOBEC-associated RNA editing (p = 0.008), supporting the concept that APOBEC mutagenesis in host and viral genomes is directly linked and occurrs during infection. Using paired sequencing of host somatic exomes, transcriptomes, and viral genomes, we demonstrated for the first-time definitive evidence of concordance between tumor and viral APOBEC mutagenesis. This finding provides a missing link connecting APOBEC mutagenesis in host and virus and supports a common mechanism driving APOBEC dysregulation.