Next Article in Journal
Peptide Platform as a Powerful Tool in the Fight against COVID-19
Next Article in Special Issue
Phylogenomic Analysis of Human Papillomavirus Type 31 and Cervical Carcinogenesis: A Study of 2093 Viral Genomes
Previous Article in Journal
Antiviral Activity of Umifenovir In Vitro against a Broad Spectrum of Coronaviruses, Including the Novel SARS-CoV-2 Virus
Previous Article in Special Issue
Human Papillomavirus Detection by Whole-Genome Next-Generation Sequencing: Importance of Validation and Quality Assurance Procedures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Otolaryngology-Head and Neck Surgery, Massachusetts Eye and Ear, Boston, MA 02118, USA
2
Massachusetts General Hospital, Boston, MA 02118, USA
3
Harvard Medical School, Boston, MA 02115, USA
4
Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
5
Department of Epidemiology, University of Kentucky, Lexington, KY 40502, USA
6
Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD 20847, USA
7
Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD 21701, USA
8
Department of Pathology, Microbiology, and Immunology Vanderbilt University Medical Center, Nashville, TN 37011, USA
9
Department of Otolaryngology, University of Pittsburgh, Pittsburgh, PA 15106, USA
10
Department of Immunology, University of Pittsburgh, Pittsburgh, PA 15106, USA
11
UPMC Hillman Cancer Center, Pittsburgh, PA 15106, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally.
Submission received: 15 July 2021 / Revised: 18 August 2021 / Accepted: 20 August 2021 / Published: 23 August 2021
(This article belongs to the Special Issue HPV Genomics)

Abstract

:
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.

1. Introduction

The apolipoprotein-B mRNA-editing catalytic polypeptide-like (APOBEC) 3 family of cytidine deaminases is a major mutagenic source in human papillomavirus (HPV)-mediated cancers, including cervical and oropharyngeal squamous cell carcinoma (HPV + OPSCC) [1,2,3,4,5]. APOBEC-induced mutations constitute a high proportion of the somatic mutations in HPV + OPSCCs and can result in driver mutations, such as activating mutations in PIK3CA [3,6]. The underlying mechanisms driving APOBEC dysregulation and mutagenesis in HPV + OPSCC, and other HPV-mediated cancers, are unknown. Since APOBECs are viral restriction agents, it has been proposed that host genome mutations may actually be “collateral damage” following APOBEC activation by viral infecton [7,8].
Following HPV infection, APOBEC is upregulated by HPV oncoproteins and by IFN as part of the innate immune response, and this can inhibit HPV both directly, with the cytidine deamination of viral DNA leading to degradation, and indirectly, by decreasing HPV virion infectivity [9,10,11,12]. APOBEC-mediated mutational signatures have been identified in both viral and cancer genomes. Viral editing is identifiable and reproducible in vitro through APOBEC induction and in vivo in cervical pre-cancers and cancers [13,14,15,16]. APOBEC mutations are present in nearly all human cancers but are particularly prominent in HPV + OPSCC, where they are associated with multiple measures of immune upregulation [1,2,4,5,17,18,19,20].
We and others have shown that HPV lineages/sublineages and individual genetic variants are associated with important differences in cervical carcinogenicity [21,22,23,24,25,26,27,28,29,30,31,32,33,34,35]. In particular, our group has published large case–control studies characterizing mutations across cervical HPV16 genomes, utilizing an HPV whole genome sequencing (WGS) approach. These studies demonstrated that benign or clearing HPV16 infections have more genetic variants in specific viral regions than cervical precancers/cancers [34,36]. APOBEC-associated mutational signatures were also identified across the HPV16 genome and associated with a significant reduction in the carcinogenicity of HPV16 [36]. Thus, evidence supports the presence of APOBEC-induced mutations in both the genomes of HPV-mediated cancers and in HPV itself; however, we do not know whether APOBEC mutational activity in host and virus is directly linked. Here, using 79 HPV + OPSCC samples with paired sequencing of host somatic exomes, transcriptomes, and viral whole genomes, we apply computational approaches to characterize APOBEC activity in both host and viral genomes.

2. Materials and Methods

2.1. Samples

A total of 79 HPV + OPSCC tissue and paired blood samples were identified through existing databases at Vanderbilt University Medical Center and the University of Pittsburgh. HPV status was determined by p16 staining and subsequently by WGS for the detection of viral genomes. After the pathology review, DNA and RNA were extracted from FFPE tissue blocks using Zymo Quick DNA and RNA FFPE kits. Of these, 18 samples had been previously contributed to the Cancer Genome Atlas (TCGA) from the same institutions and had pre-existing whole exome sequencing (WES) and RNA-Seq. For these cases, DNA was extracted from FFPE blocks for viral WGS alone.

2.2. Sequencing

Somatic DNA underwent WES using the KAPA HyperPrep Library Preparation Kit, followed by hybrid capture with the Illumina Rapid Capture Exome Enrichment Kit. All sample pairs were validated with Fluidigm fingerprint data to confirm sample identity and fidelity. The average (range) coverage was 218× (70–402×). Somatic RNA underwent RNA-Seq using the Illumina TruSeq RNA Exome Kit. The average (range) number of reads per samples was 151 M (17.5 M–333 M). HPV16 DNA underwent WGS as previously described by custom AmpliSeqTM panels, followed by sequencing on an Ion Torrent platform with average (range) coverage of 23,000× (1500×–56,000×) [37]. A total of 74 WES-seq, 68 RNA-seq, and 72 viral WGS samples passed quality control. Of these samples, 63 had data available from all three platforms.

2.3. Informatics

Somatic datasets: Somatic SNVs were called using MuTect v1. Mutations were screened against Panel of Normal (PoN) databases to remove common SNPs and recurrent sequencing artifacts. SNP-counting and copy number segmentation were carried out using FACETS. SNP counts were generated with a minimum mapping quality of 15, minimum base quality of 20, pseudo-spacing of 100 and a minimum read count of 25. Copy number data were segmented using a window size of 1000, Variant Allele Fraction (VAF) threshold of 0.3 and cval of 300. Allelic amplifications/deletions were defined as regions of any size with rounded integer total copy number states above or below 2, respectively. Somatic mutation lists were combined with mutations lists from the entire TCGA cohort, and mutational signatures were deconvolved using NMF (K = 8), as described previously [1,38]. The IFN-γ score was defined as the mean expression of six genes (IFNG, IDO1, CXCL9, CXCL10, HLA-DRA and STAT1) as described previously [39]. The analysis of the PI3K/AKT/mTOR pathway was conducted using a six major effector gene signature (PI3K, AKT1, AKT2, AKT3, PTEN, and MTOR) as described previously [40], and further, hotspot loci were examined in IGV. Significantly mutated genes were identified using MutSig2CV, and significantly copy-number-altered genes were identified using GISTIC 2.0. MutSig and GISTIC p values were combined using Fisher’s method to calculate a unified p value per gene. Pathway enrichment scores were then calculated by combining unified p values for genes in a pathway, also using Fisher’s method.
Viral datasets: All HPV16 WGS data were aligned to a consensus HPV16 sequence NC_001526.4 from PaVE [41]. Position-wise base counts were obtained, and a major haplotype was created by taking the most prevalent nucleotide at each position. These haplotypes were fed into the multi-alignment tool MUSCLE, along with multiple reference builds for HPV16 (AF402678, AF472509, AF534061, AF536179, AF536180, AY686579, HQ644236, HQ644257, HQ644298, K02718), and reference sequences for HPV35, HPV33, and HPV18. The resulting alignment was fed into PhyML (run using TOPALi v2.5), assuming a transversion model and allowing for invariable sites and specifying a gamma distribution. The resulting phylogenetic tree allowed for the sublineage assignment of each sample by nearest distance to a reference sequence. All samples were realigned with respect to their relative reference. Putative mutation sites were nominated by having a minimum total coverage of 100, a minimum alternate allele fraction of 0.02, and minimum alternate reads of 2. Putative mutations were anonymized and manually reviewed in IGV alongside other samples aligned to the same reference to remove low-evidence mutations and recurrent sequencing artifacts.
Viral APOBEC-induced mutations were identified as C->T and C->G mutations at TCW motifs (W is A or T), as described previously [36]. To count the number of APOBEC targetable sites, we first counted TCW motifs across the HPV16 genome for each reference build and then the total number across all samples’ HPV16 genomes on their respective reference builds. Since there are three possible changes at each nucleotide position, APOBEC-targetable sites were counted as one-third of the total number of TCW motifs. Mutation rates were calculated as the total number of APOBEC mutations, divided by the total number of APOBEC-targetable sites at the gene level and genome level, stratified by mutation type (synonymous (S) and non-synonymous (NS)) and variant allele fraction (VAF) (low-VAF (VAF ≤ 0.5) and high-VAF (VAF > 0.5)).

3. Results

A total of 79 HPV + OPSCC primary tumor samples were included in the analysis. The average patient age was 55. A total of 82% of the cohort were male and 60% had a history of tobacco exposure of >10 pack years. Somatic exomes were pooled with TCGA data and mutational signatures extracted. As expected, HPV + OPSCC samples predominately clustered with other tumors dominated by the APOBEC mutational signature, with variable APOBEC mutational burden between samples (Figure 1A,B, Supplemental Figure S1). PIK3CA hotspot mutations, which have previously been described to be caused by APOBEC activity in HPV + OPSCC, were exclusively associated with APOBEC enriched samples (Figure 1B). The PI3K/AKT/mTOR pathway was disproportionately mutated across the cohort (p < 3.22 × 10−8) and amplifications in chromosomes arm 3q, which contains PIK3CA, were present in 41 samples. Both of these alterations have previously been described to occur frequently in HPV + OPSCC5. APOBEC mutations, A3A expression, and APOBEC RNA editing, as measured by hotspot DDOST558C > U mutations, were weakly correlated, consistent with prior reports, likely due to the episodic nature of APOBEC activity (Figure 1C) [42].
APOBEC mutations were annotated in the HPV16 genomes after first assigning each viral genome to a sublineage. A1 (n = 45) was the predominant sublineage, followed by A2 (n = 19), D3 (n = 4), A4 (n = 3) and C1 (n = 1). There was an average of five (range 0–29) mutations per viral genome. The mean APOBEC mutations in the viral genomes was one (range 0–5), with an APOBEC to non-APOBEC mutation ratio mean of 1:5. Variants were most heavily clustered in the E2 C-terminal DNA-binding domain (Q349E, S364C, S348C, E344Q), including in a position previously described [43] to impact transcription and the replication of the viral genome (D338N) (Figure 2A,B and Supplemental Table S1). There was no clear strand bias for HPV16 APOBEC mutations, suggesting that cytosines on both strands of the viral DNA are vulnerable to attack by APOBEC, likely during viral DNA replication, when both strands of viral DNA experience a single-stranded state. This is similar to APOBEC mutations in human DNA, which are also observed on both strands, but with a detectable bias toward the lagging-strand template, not observed here, possibly due to a lack of power owing to low overall viral mutation counts [44]. APOBEC mutations were found predominantly at very low and very high VAFs (Figure 2C,D). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-VAF mutations (occurring newly within the host), while non-APOBEC mutations were more likely to be high-VAF (existing prior to infection of the current host), suggesting that APOBEC mutagenesis actively occurrs in viral genomes during infection of the current host. Viral APOBEC mutations had a higher nonsynonymous/synonymous ratio, while non-APOBEC mutations showed the depletion of non-synonymous variants at high VAFs compared to expectation, suggesting the influence of purifying selection (Figure 2E, Supplemental Figure S2).
To assess the relationship between APOBEC mutations in the viral and host genomes, samples were divided into APOBEC-low and -high mutational burden groupings based on host somatic mutations (Figure 3A). APOBEC-high samples had higher viral APOBEC mutation rates (p = 0.017; Figure 3B), and APOBEC mutagenesis in host and viral genomes was positively correlated (ρ = 0.24, p = 0.056; Figure 3C). We further examined additional measures of APOBEC activity in both tumor and virus and found similar trends for APOBEC RNA editing at the DDOST hairpin hotspot (p = 0.008), A3A mRNA expression (p = 0.14), and IFN-γ scores (p = 0.053), which are tied to APOBEC activity (Figure 3D–F).

4. Discussion

We conducted the first study to-date evaluating APOBEC-induced mutations in both host and viral genomes, providing a missing link connecting APOBEC mutagenesis in host and virus and supporting a common mechanism driving APOBEC dysregulation. Our group previously characterized HPV16 APOBEC-induced mutations in over 5000 cervical benign, pre-cancer and cancer samples, and identified low VAF mutations linked to benign or clearing infections [36]. We [36] and others [16] have suggested that the low VAF viral mutations were likely induced during infection. Here, we characterized HPV16 APOBEC-induced mutations in HPV + OPSCCs using deep NGS and utilized the paired sequencing of HPV16 and host somatic genomes to determine whether APOBEC activity was concordant in host and viral genomes. We showed that OPSCC HPV16 APOBEC-induced mutations occurred predominantly at low VAF, similar to what we observed in cervical HPV16 genomes, and importantly, correlated them with host somatic mutations. These data show that APOBEC mutates both cancer and viral genomes during an infection.
We identified APOBEC-signature mutations in OPSCC HPV16 genomes, which existed at both very low (new within host mutations) and very high (old mutations carried from prior host infections) VAFs. The OPSCC HPV16 APOBEC mutation patterns observed were similar to those previously reported in cervical HPV16 genomes [36]. Viral APOBEC mutations were most prominent at low VAFs, which is of interest as HPV genomes have long been considered stable during persistent infection, yet these data suggest that APOBEC editing of viral genomes is actively occurs during infection. This is despite the fact that HPVs have evolved to have fewer APOBEC target sequence sites, likely in an effort to avoid APOBEC restriction and thus viral damage and clearance [45]. This evolutionary pattern also results in the enrichment of potential non-synonymous APOBEC mutations compared to non-APOBEC mutations. In accordance with this, we observed a high nonsynonymous-to-synonymous ratio for viral APOBEC mutations overall. The densest cluster of APOBEC mutations occurred in E2 and included D338N, an HPV16 SNP that was previously reported to reduce p53 binding [46]. Interestingly, D338N was also prevalent in cervical cancer samples from our previous study. We also found a depletion of nonsynonymous high-VAF mutations for non-APOBEC mutations, but did not observe this for APOBEC mutations. This could suggest that some nonsynonymous viral APOBEC-induced mutations may be beneficial to the virus and contribute to the evasion of host immunity by altering viral antigens [7]. It is also possible that deleterious APOBEC mutations that result in viral damage and clearance are simply not detected, as those viruses have already been cleared during initial infection stages, leaving only viral genomes with a selective growth advantage present in the cancer.
Importantly, we found that APOBEC-enriched somatic samples (greatest APOBEC mutational burden and highest level of RNA editing) had higher viral APOBEC mutation rates, providing direct evidence that APOBEC mutagenesis in host and viral genomes is linked and occurs during infection. This suggests that in HPV-mediated cancers, the underlying processes driving APOBEC dysregulation and resultant TCW context mutations may be the same for both viral and somatic mutations. Mechanisms that drive APOBEC mutagenesis in HPV + OPSCC, however, remain unknown. While APOBEC is known to be activated by HPV oncoproteins, and through IFN signaling as part of innate immunity/viral nucleic acid sensing, APOBEC mutations are prevalent in non-virally mediated tumors, such as bladder, breast, and lung cancer, as well as non-HPV mediated head and neck cancers, albeit to a lesser degree [1,2,20]. Thus, it is possible that both viral/innate immunity and immune-independent processes, such as ssDNA targeting of the lagging strand during replication stress or after double strand breaks, drive APOBEC expression in HPV + OPSCC, and the high rates of APOBEC mutations in HPV + OPSCCs are a resultant cumulative effect. This concept is supported by work from our group and others showing multiple sources of APOBEC upregulation likely contribute to APOBEC mutagenesis in HPV + OPSCC, including APOBEC germline polymorphisms and immune upregulation in response to mutation-induced neoantigens [6].
In summary, we report the first study evaluating APOBEC mutagenesis in paired human somatic and viral genomes, identifying the presence of APOBEC mutations in OPSCC HPV genomes and concordance between APOBEC mutational burden within virus and host during infection. These data provide a missing link connecting APOBEC mutagenesis in host and virus and support a common mechanism driving APOBEC dysregulation.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/v13081666/s1, Figure S1: Mutational signatures ordered by APOBEC contribution, Figure S2: Expected vs. actual nonsynonymous/synonymous ratio for APOBEC and non-APOBEC mutations by viral gene. Table S1: HPV E2 APOBEC Mutations.

Author Contributions

D.L.F., K.A.L.K., R.L.F., L.M. conceptualized the study and wrote the manuscript. M.L., A.L., M.S.L. performed data analysis and contributed to manuscript writing, J.S.L., M.P., M.Y., M.C., J.F.B., M.S. and S.B. generated data and edits the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by DLF (sundry).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Dana Farber Cancer Institute Review Board, #18-653.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alexandrov, L.B.; Nik-Zainal, S.; Wedge, D.C.; Aparicio, S.A.; Behjati, S.; Biankin, A.V.; Bignell, G.R.; Bolli, N.; Borg, A.; Borresen-Dale, A.L.; et al. Signatures of mutational processes in human cancer. Nature 2013, 500, 415–421. [Google Scholar] [CrossRef] [Green Version]
  2. Roberts, S.A.; Lawrence, M.S.; Klimczak, L.J.; Grimm, S.A.; Fargo, D.; Stojanov, P.; Kiezun, A.; Kryukov, G.V.; Carter, S.L.; Saksena, G.; et al. An APOBEC cytidine deaminase mutagenesis pattern is widespread in human cancers. Nat. Genet. 2013, 45, 970–976. [Google Scholar] [CrossRef]
  3. Henderson, S.; Chakravarthy, A.; Su, X.; Boshoff, C.; Fenton, T.R. APOBEC-mediated cytosine deamination links PIK3CA helical domain mutations to human papillomavirus-driven tumor development. Cell Rep. 2014, 7, 1833–1841. [Google Scholar] [CrossRef] [Green Version]
  4. Faden, D.L.; Thomas, S.; Cantalupo, P.G.; Agrawal, N.; Myers, J.; DeRisi, J. Multi-modality analysis supports APOBEC as a major source of mutations in head and neck squamous cell carcinoma. Oral Oncol. 2017, 74, 8–14. [Google Scholar] [CrossRef]
  5. Faden, D.L.; Langenbucher, A.; Kuhs, K.; Lewis, J.S.; Mirabello, L.; Yeager, M.; Boland, J.F.; Bass, S.; Steinberg, M.; Cullen, M.; et al. HPV+ Oropharyngeal Squamous Cell Carcinomas from patients with two tumors display synchrony of viral genomes yet discordant mutational profiles and signatures. Carcinogenesis 2020, 42, 14–20. [Google Scholar] [CrossRef]
  6. Faden, D.L.; Ding, F.; Lin, Y.; Zhai, S.; Kuo, F.; Chan, T.A.; Morris, L.G.; Ferris, R.L. APOBEC mutagenesis is tightly linked to the immune landscape and immunotherapy biomarkers in head and neck squamous cell carcinoma. Oral Oncol. 2019, 96, 140–147. [Google Scholar] [CrossRef]
  7. Smith, N.J.; Fenton, T.R. The APOBEC3 genes and their role in cancer: Insights from human papillomavirus. J. Mol. Endocrinol. 2019, 62, R269–R287. [Google Scholar] [CrossRef] [Green Version]
  8. Mori, S.; Takeuchi, T.; Ishii, Y.; Kukimoto, I. Identification of APOBEC3B promoter elements responsible for activation by human papillomavirus type 16 E6. Biochem. Biophys. Res. Commun. 2015, 460, 555–560. [Google Scholar] [CrossRef]
  9. Vieira, V.C.; Leonard, B.; White, E.A.; Starrett, G.J.; Temiz, N.A.; Lorenz, L.D.; Lee, D.; Soares, M.A.; Lambert, P.F.; Howley, P.M.; et al. Human papillomavirus E6 triggers upregulation of the antiviral and cancer genomic DNA deaminase APOBEC3B. MBio 2014, 5. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, Z.; Wakae, K.; Kitamura, K.; Aoyama, S.; Liu, G.; Koura, M.; Monjurul, A.M.; Kukimoto, I.; Muramatsu, M. APOBEC3 deaminases induce hypermutation in human papillomavirus 16 DNA upon beta interferon stimulation. J. Virol. 2014, 88, 1308–1317. [Google Scholar] [CrossRef] [Green Version]
  11. Warren, C.J.; Xu, T.; Guo, K.; Griffin, L.M.; Westrich, J.A.; Lee, D.; Lambert, P.F.; Santiago, M.L.; Pyeon, D. APOBEC3A functions as a restriction factor of human papillomavirus. J. Virol. 2015, 89, 688–702. [Google Scholar] [CrossRef] [Green Version]
  12. Westrich, J.A.; Warren, C.J.; Klausner, M.J.; Guo, K.; Liu, C.W.; Santiago, M.L.; Pyeon, D. Human Papillomavirus 16 E7 Stabilizes APOBEC3A Protein by Inhibiting Cullin 2-Dependent Protein Degradation. J. Virol. 2018, 92. [Google Scholar] [CrossRef] [Green Version]
  13. Vartanian, J.P.; Guetard, D.; Henry, M.; Wain-Hobson, S. Evidence for editing of human papillomavirus DNA by APOBEC3 in benign and precancerous lesions. Science 2008, 320, 230–233. [Google Scholar] [CrossRef] [Green Version]
  14. Wakae, K.; Aoyama, S.; Wang, Z.; Kitamura, K.; Liu, G.; Monjurul, A.M.; Koura, M.; Imayasu, M.; Sakamoto, N.; Nakamura, M.; et al. Detection of hypermutated human papillomavirus type 16 genome by Next-Generation Sequencing. Virology 2015, 485, 460–466. [Google Scholar] [CrossRef] [Green Version]
  15. Kukimoto, I.; Mori, S.; Aoyama, S.; Wakae, K.; Muramatsu, M.; Kondo, K. Hypermutation in the E2 gene of human papillomavirus type 16 in cervical intraepithelial neoplasia. J. Med. Virol. 2015, 87, 1754–1760. [Google Scholar] [CrossRef]
  16. Hirose, Y.; Onuki, M.; Tenjimbayashi, Y.; Mori, S.; Ishii, Y.; Takeuchi, T.; Tasaka, N.; Satoh, T.; Morisada, T.; Iwata, T.; et al. Within-Host Variations of Human Papillomavirus Reveal APOBEC Signature Mutagenesis in the Viral Genome. J. Virol. 2018, 92. [Google Scholar] [CrossRef] [Green Version]
  17. Gillison, M.L.; Akagi, K.; Xiao, W.; Jiang, B.; Pickard, R.K.L.; Li, J.; Swanson, B.J.; Agrawal, A.D.; Zucker, M.; Stache-Crain, B.; et al. Human papillomavirus and the landscape of secondary genetic alterations in oral cancers. Genome Res. 2019, 29, 1–17. [Google Scholar] [CrossRef] [Green Version]
  18. Seplyarskiy, V.B.; Soldatov, R.A.; Popadin, K.Y.; Antonarakis, S.E.; Bazykin, G.A.; Nikolaev, S.I. APOBEC-induced mutations in human cancers are strongly enriched on the lagging DNA strand during replication. Genome Res. 2016, 26, 174–182. [Google Scholar] [CrossRef] [Green Version]
  19. Swanton, C.; McGranahan, N.; Starrett, G.J.; Harris, R.S. APOBEC Enzymes: Mutagenic Fuel for Cancer Evolution and Heterogeneity. Cancer Discov. 2015, 5, 704–712. [Google Scholar] [CrossRef] [Green Version]
  20. Burns, M.B.; Temiz, N.A.; Harris, R.S. Evidence for APOBEC3B mutagenesis in multiple human cancers. Nat. Genet. 2013, 45, 977–983. [Google Scholar] [CrossRef] [Green Version]
  21. Hildesheim, A.; Schiffman, M.; Bromley, C.; Wacholder, S.; Herrero, R.; Rodriguez, A.C.; Bratti, M.C.; Sherman, M.E.; Scarpidis, U.; Lin, Q.-Q.; et al. Human Papillomavirus Type 16 Variants and Risk of Cervical Cancer. J. Natl. Cancer Inst. 2001, 93, 315–318. [Google Scholar] [CrossRef] [Green Version]
  22. Pientong, C.; Wongwarissara, P.; Ekalaksananan, T.; Swangphon, P.; Kleebkaow, P.; Kongyingyoes, B.; Siriaunkgul, S.; Tungsinmunkong, K.; Suthipintawong, C. Association of human papillomavirus type 16 long control region mutation and cervical cancer. Virol. J. 2013, 10, 30. [Google Scholar]
  23. Xi, L.F.; Koutsky, L.A.; Hildesheim, A.; Galloway, D.A.; Wheeler, C.M.; Winer, R.L.; Ho, J.; Kiviat, N.B. Risk for High-Grade Cervical Intraepithelial Neoplasia Associated with Variants of Human Papillomavirus Types 16 and 18. Cancer Epidemiol. Biomark. Prev. 2007, 16, 4–10. [Google Scholar] [CrossRef] [Green Version]
  24. Schiffman, M.; Rodriguez, A.C.; Chen, Z.; Wacholder, S.; Herrero, R.; Hildesheim, A.; Desalle, R.; Befano, B.; Yu, K.; Safaeian, M.; et al. A Population-Based Prospective Study of Carcinogenic Human Papillomavirus Variant Lineages, Viral Persistence, and Cervical Neoplasia. Cancer Res. 2010, 70, 3159–3169. [Google Scholar] [CrossRef] [Green Version]
  25. Cornet, I.; Gheit, T.; Iannacone, M.R.; Vignat, J.; Sylla, B.S.; Del Mistro, A.; Franceschi, S.; Tommasino, M.; Clifford, G.M. HPV16 genetic variation and the development of cervical cancer worldwide. Br. J. Cancer 2013, 108, 240–244. [Google Scholar] [CrossRef] [Green Version]
  26. Gheit, T.; Cornet, I.; Clifford, G.M.; Iftner, T.; Munk, C.; Tommasino, M.; Kjaer, S.K. Risks for Persistence and Progression by Human Papillomavirus Type 16 Variant Lineages Among a Population-Based Sample of Danish Women. Cancer Epidemiol. Biomark. Prev. 2011, 20, 1315–1321. [Google Scholar] [CrossRef] [Green Version]
  27. Zehbe, I.; Voglino, G.; Delius, H.; Wilander, E.; Tommasino, M. Risk of cervical cancer and geographical variations of human papillomavirus 16 E6 polymorphisms. Lancet 1998, 352, 1441–1442. [Google Scholar] [CrossRef]
  28. Zuna, R.E.; Moore, W.E.; Shanesmith, R.P.; Dunn, S.T.; Wang, S.S.; Schiffman, M.; Blakey, G.L.; Teel, T. Association of HPV16 E6 variants with diagnostic severity in cervical cytology samples of 354 women in a US population. Int. J. Cancer 2009, 125, 2609–2613. [Google Scholar] [CrossRef] [Green Version]
  29. Sichero, L.; Ferreira, S.; Trottier, H.; Duarte-Franco, E.; Ferenczy, A.; Franco, E.L.; Villa, L.L. High grade cervical lesions are caused preferentially by non-European variants of HPVs 16 and 18. Int. J. Cancer 2007, 120, 1763–1768. [Google Scholar] [CrossRef]
  30. Berumen, J.; Ordoñez, R.M.; Lazcano, E.; Salmeron, J.; Galvan, S.C.; Estrada, R.A.; Yunes, E.; Garcia-Carranca, A.; Gonzalez-Lira, G.; Madrigal-de la Campa, A. Asian-American Variants of Human Papillomavirus 16 and Risk for Cervical Cancer: A Case–Control Study. J. Natl. Cancer Inst. 2001, 93, 1325–1330. [Google Scholar] [CrossRef]
  31. Freitas, L.B.; Chen, Z.; Muqui, E.F.; Boldrini, N.A.T.; Miranda, A.E.; Spano, L.C.; Burk, R.D. Human Papillomavirus 16 Non-European Variants Are Preferentially Associated with High-Grade Cervical Lesions. PLoS ONE 2014, 9, e100746. [Google Scholar] [CrossRef] [Green Version]
  32. Mirabello, L.; Yeager, M.; Cullen, M.; Boland, J.F.; Chen, Z.; Wentzensen, N.; Zhang, X.; Yu, K.; Yang, Q.; Mitchell, J.; et al. HPV16 Sublineage Associations With Histology-Specific Cancer Risk Using HPV Whole-Genome Sequences in 3200 Women. J. Natl. Cancer Inst. 2016, 108, djw100. [Google Scholar]
  33. Mirabello, L.; Clarke, M.A.; Nelson, C.W.; Dean, M.; Wentzensen, N.; Yeager, M.; Cullen, M.; Boland, J.F.; Workshop, N.H.; Schiffman, M.; et al. The Intersection of HPV Epidemiology, Genomics and Mechanistic Studies of HPV-Mediated Carcinogenesis. Viruses 2018, 10, 80. [Google Scholar] [CrossRef] [Green Version]
  34. Mirabello, L.; Yeager, M.; Yu, K.; Clifford, G.M.; Xiao, Y.; Zhu, B.; Cullen, M.; Boland, J.F.; Wentzensen, N.; Nelson, C.W.; et al. HPV16 E7 Genetic Conservation Is Critical to Carcinogenesis. Cell 2017, 170, 1164–1174.e1166. [Google Scholar] [CrossRef] [Green Version]
  35. Pinheiro, M.; Gage, J.C.; Clifford, G.M.; Demarco, M.; Cheung, L.C.; Chen, Z.; Yeager, M.; Cullen, M.; Boland, J.F.; Chen, X.; et al. Association of HPV35 with cervical carcinogenesis among women of African ancestry: Evidence of viral-host interaction with implications for disease intervention. Int. J. Cancer 2020, 147, 2677–2686. [Google Scholar] [CrossRef]
  36. Zhu, B.; Xiao, Y.; Yeager, M.; Clifford, G.; Wentzensen, N.; Cullen, M.; Boland, J.F.; Bass, S.; Steinberg, M.K.; Raine-Bennett, T.; et al. Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance. Nat. Commun. 2020, 11, 886. [Google Scholar] [CrossRef] [Green Version]
  37. Cullen, M.; Boland, J.F.; Schiffman, M.; Zhang, X.; Wentzensen, N.; Yang, Q.; Chen, Z.; Yu, K.; Mitchell, J.; Roberson, D.; et al. Deep sequencing of HPV16 genomes: A new high-throughput tool for exploring the carcinogenicity and natural history of HPV16 infection. Papillomavirus Res. 2015, 1, 3–11. [Google Scholar] [CrossRef] [Green Version]
  38. TCGA. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 2015, 517, 576–582. [Google Scholar] [CrossRef] [Green Version]
  39. Ayers, M.; Lunceford, J.; Nebozhyn, M.; Murphy, E.; Loboda, A.; Kaufman, D.R.; Albright, A.; Cheng, J.D.; Kang, S.P.; Shankaran, V.; et al. IFN-gamma-related mRNA profile predicts clinical response to PD-1 blockade. J. Clin. Investig. 2017, 127, 2930–2940. [Google Scholar] [CrossRef]
  40. Fang, W.L.; Huang, K.H.; Lan, Y.T.; Lin, C.H.; Chang, S.C.; Chen, M.H.; Chao, Y.; Lin, W.C.; Lo, S.S.; Li, A.F.; et al. Mutations in PI3K/AKT pathway genes and amplifications of PIK3CA are associated with patterns of recurrence in gastric cancers. Oncotarget 2016, 7, 6201–6220. [Google Scholar] [CrossRef] [Green Version]
  41. Van Doorslaer, K.; Li, Z.; Xirasagar, S.; Maes, P.; Kaminsky, D.; Liou, D.; Sun, Q.; Kaur, R.; Huyen, Y.; McBride, A.A. The Papillomavirus Episteme: A major update to the papillomavirus sequence database. Nucleic Acids Res. 2016, 45, D499–D506. [Google Scholar] [CrossRef]
  42. Jalili, P.; Bowen, D.; Langenbucher, A.; Park, S.; Aguirre, K.; Corcoran, R.B.; Fleischman, A.G.; Lawrence, M.S.; Zou, L.; Buisson, R. Quantification of ongoing APOBEC3A activity in tumor cells by monitoring RNA editing at hotspots. Nat. Commun. 2020, 11, 2971. [Google Scholar] [CrossRef]
  43. Brown, C.; Kowalczyk, A.M.; Taylor, E.R.; Morgan, I.M.; Gaston, K. P53 represses human papillomavirus type 16 DNA replication via the viral E2 protein. Virol. J. 2008, 5, 5. [Google Scholar] [CrossRef] [Green Version]
  44. Haradhvala, N.J.; Polak, P.; Stojanov, P.; Covington, K.R.; Shinbrot, E.; Hess, J.M.; Rheinbay, E.; Kim, J.; Maruvka, Y.E.; Braunstein, L.Z.; et al. Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair. Cell 2016, 164, 538–549. [Google Scholar] [CrossRef] [Green Version]
  45. Warren, C.J.; Van Doorslaer, K.; Pandey, A.; Espinosa, J.M.; Pyeon, D. Role of the host restriction factor APOBEC3 on papillomavirus evolution. Virus Evol. 2015, 1. [Google Scholar] [CrossRef]
  46. McBride, A.A. The papillomavirus E2 proteins. Virology 2013, 445, 57–79. [Google Scholar] [CrossRef] [Green Version]
Figure 1. APOBEC is the dominant mutational signature in HPV + OPSCC. (A) Somatic mutation lists were combined with all cancers in TCGA and mutational signatures were deciphered using NMF. Mutational patterns were then projected into two-dimensional space with t-SNE. Head and neck squamous cell carcinoma tumors were colored in blue dots and HPV + OPSCCs were highlighted with blue circles. (B) Mutational signatures ordered by APOBEC contribution (magenta) in HPV+ OPSCC samples. PIK3CA hotspot mutations, E542K and E545K, were exclusively detected within APOBEC-enriched (>25%; dashed line box) tumors. (C) Weak correlations between APOBEC mutation burden and APOBEC3A (A3A) expression (left) and A3A-associated RNA editing activity (right).
Figure 1. APOBEC is the dominant mutational signature in HPV + OPSCC. (A) Somatic mutation lists were combined with all cancers in TCGA and mutational signatures were deciphered using NMF. Mutational patterns were then projected into two-dimensional space with t-SNE. Head and neck squamous cell carcinoma tumors were colored in blue dots and HPV + OPSCCs were highlighted with blue circles. (B) Mutational signatures ordered by APOBEC contribution (magenta) in HPV+ OPSCC samples. PIK3CA hotspot mutations, E542K and E545K, were exclusively detected within APOBEC-enriched (>25%; dashed line box) tumors. (C) Weak correlations between APOBEC mutation burden and APOBEC3A (A3A) expression (left) and A3A-associated RNA editing activity (right).
Viruses 13 01666 g001
Figure 2. Mutation distributions in the viral and tumor genomes. (A) All viral mutations distributed across the genome with APOBEC mutations colored in red and non-APOBEC mutations in blue. (B) Distribution of viral APOBEC mutations, with positive strand mutations in purple and negative strands in green, highlighting enrichment of variants in E2. Viral (C) and tumor (D) APOBEC and non-APOBEC mutations by variant allele fraction demonstrating predominance of very low and very high VAF mutations in the viral genome. (E) Expected (blue and red lines) vs. actual (black dots) nonsynonymous/synonymous ratio for APOBEC (red) and non-APOBEC (blue) mutations demonstrating higher non-synonymous/synonymous ratios for APOBEC mutations and depletion of non-synonymous variants at high VAFs compared to expected in non-APOBEC mutations, suggesting purifying selection.
Figure 2. Mutation distributions in the viral and tumor genomes. (A) All viral mutations distributed across the genome with APOBEC mutations colored in red and non-APOBEC mutations in blue. (B) Distribution of viral APOBEC mutations, with positive strand mutations in purple and negative strands in green, highlighting enrichment of variants in E2. Viral (C) and tumor (D) APOBEC and non-APOBEC mutations by variant allele fraction demonstrating predominance of very low and very high VAF mutations in the viral genome. (E) Expected (blue and red lines) vs. actual (black dots) nonsynonymous/synonymous ratio for APOBEC (red) and non-APOBEC (blue) mutations demonstrating higher non-synonymous/synonymous ratios for APOBEC mutations and depletion of non-synonymous variants at high VAFs compared to expected in non-APOBEC mutations, suggesting purifying selection.
Viruses 13 01666 g002
Figure 3. APOBEC mutations are concordant between tumor and viral genomes. APOBEC (A) mutations in tumor genomes, (B) mutations in the viral genome. (C) correlation between host and viral genomes, (D) RNA editing and (E) 3A gene expression, demonstrating linkage between APOBEC activity in virus and host. (F) APOBEC-enriched subjects also had higher IFN-γ scores, suggesting immune upregulation, possibly in response to viral infection.
Figure 3. APOBEC mutations are concordant between tumor and viral genomes. APOBEC (A) mutations in tumor genomes, (B) mutations in the viral genome. (C) correlation between host and viral genomes, (D) RNA editing and (E) 3A gene expression, demonstrating linkage between APOBEC activity in virus and host. (F) APOBEC-enriched subjects also had higher IFN-γ scores, suggesting immune upregulation, possibly in response to viral infection.
Viruses 13 01666 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Faden, D.L.; Kuhs, K.A.L.; Lin, M.; Langenbucher, A.; Pinheiro, M.; Yeager, M.; Cullen, M.; Boland, J.F.; Steinberg, M.; Bass, S.; et al. APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma. Viruses 2021, 13, 1666. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081666

AMA Style

Faden DL, Kuhs KAL, Lin M, Langenbucher A, Pinheiro M, Yeager M, Cullen M, Boland JF, Steinberg M, Bass S, et al. APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma. Viruses. 2021; 13(8):1666. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081666

Chicago/Turabian Style

Faden, Daniel L., Krystle A. Lang Kuhs, Maoxuan Lin, Adam Langenbucher, Maisa Pinheiro, Meredith Yeager, Michael Cullen, Joseph F. Boland, Mia Steinberg, Sara Bass, and et al. 2021. "APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma" Viruses 13, no. 8: 1666. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop