Explore the vast range of titles available.

HPV16 causes half of cervical cancers worldwide; for unknown reasons, most infections resolve within two years. Here, we analyze the viral genomes of 5,328 HPV16-positive case-control samples to investigate mutational signatures and the role of human APOBEC3-induced mutations in viral clearance and cervical carcinogenesis. We identify four de novo mutational signatures, one of which matches the COSMIC APOBEC-associated signature 2. The viral genomes of the precancer/cancer cases are less likely to contain within-host somatic HPV16 APOBEC3-induced mutations (Fisher’s exact test, P = 6.2 x 10 −14 ), and have a 30% lower nonsynonymous APOBEC3 mutation burden compared to controls. We replicate the low prevalence of HPV16 APOBEC3-induced mutations in 1,749 additional cases. APOBEC3 mutations also historically contribute to the evolution of HPV16 lineages. We demonstrate that cervical infections with a greater burden of somatic HPV16 APOBEC3-induced mutations are more likely to be benign or subsequently clear, suggesting they may reduce persistence, and thus progression, within the host. The APOBEC mutational signature is prevalent in different tumour types. Here, using HPV16- positive cervical samples, the authors show that the signature is more prevalent in the viral genome of benign or clearing HPV16 infections compared to the viral genomes of the more advanced precancerous lesions or cervical cancer.

Invasive cervical cancers (ICC), caused by HPV infections, have a heterogeneous molecular landscape. We investigate the detection, timing, and HPV type specificity of somatic mutations in 3929 HPV-positive exfoliated cervical cell samples from individuals undergoing cervical screening in the U.S. using deep targeted sequencing in ICC cases, precancers, and HPV-positive controls. We discover a subset of hotspot mutations rare in controls (2.6%) but significantly more prevalent in precancers, particularly glandular precancer lesions (10.2%), and cancers (25.7%), supporting their involvement in ICC carcinogenesis. Hotspot mutations differ by HPV type, and HPV18/45-positive ICC are more likely to have multiple hotspot mutations compared to HPV16-positive ICC. The proportion of cells containing hotspot mutations is higher (i.e., higher variant allele fraction) in ICC and mutations are detectable up to 6 years prior to cancer diagnosis. Our findings demonstrate the feasibility of using exfoliated cervical cells for detection of somatic mutations as potential diagnostic biomarkers. Invasive cervical cancer is caused by HPV infection, but the disease itself is highly variable. Here, the authors use deep targeted sequencing to identify hotspot mutations in routine screening samples prior to diagnosis, which differed depending on HPV type.

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.

Human papillomavirus (HPV) type 31 (HPV31) is closely related to the most carcinogenic type, HPV16, but only accounts for 4% of cervical cancer cases worldwide. Viral genetic and epigenetic variations have been associated with carcinogenesis for other high-risk HPV types, but little is known about HPV31. We sequenced 2093 HPV31 viral whole genomes from two large studies, one from the U.S. and one international. In addition, we investigated CpG methylation in a subset of 175 samples. We evaluated the association of HPV31 lineages/sublineages, single nucleotide polymorphisms (SNPs) and viral methylation with cervical carcinogenesis. HPV31 A/B clade was >1.8-fold more associated with cervical intraepithelial neoplasia grade 3 and cancer (CIN3+) compared to the most common C lineage. Lineage/sublineage distribution varied by race/ethnicity and geographic region. A viral genome-wide association analysis identified SNPs within the A/B clade associated with CIN3+, including H23Y (C626T) (odds ratio = 1.60, confidence intervals = 1.17–2.19) located in the pRb CR2 binding-site within the E7 oncogene. Viral CpG methylation was higher in lineage B, compared to the other lineages, and was most elevated in CIN3+. In conclusion, these data support the increased oncogenicity of the A/B lineages and suggest variation of E7 as a contributing risk factor.

Multiple primary tumors (MPT) have been described in carriers of inherited cancer predisposition genes. However, the genetic etiology of a large proportion of MPT cases remains unclear. We reviewed 267 patients with hereditary cancer predisposition syndromes (HCPS) that underwent genetic counseling and selected 22 patients with MPT to perform genomic analysis (CytoScan HD Array, Affymetrix) aiming to identify new alterations related to a high risk of developing MPT. Twenty patients had a positive family history of cancer and 11 met phenotypic criteria for HCPS. Genetic testing for each of the genes associated with these syndromes revealed negative results for pathogenic mutations. Seventeen rare germline copy number variations (CNVs) covering 40 genes were identified in 11 patients, including an EPCAM / MSH2 deletion in one Lynch syndrome patient. An enrichment analysis revealed a significant number of genes (where the CNVs are mapped) associated with carcinogenesis and/or related to functions implicated with tumor development, such as proliferation and cell survival. An interaction network analysis highlighted the importance of TP53 pathway in cancer emergence. A high number of germline copy-neutral loss of heterozygosity (cnLOH) was identified in nine cases, particularly in two patients. Eighteen genes were covered by both rare CNVs and cnLOH, including 14 related to tumorigenesis and seven genes ( ABCC1 , KDM4C , KIAA0430 , MYH11, NDE1 , PIWIL2 , and ULK2 ) specifically associated with cellular growth and proliferation. Overall, we identified 14 cases with rare CNVs and/or cnLOH that may contribute to the risk of MPT development. Key message CNVs may explain the risk of hereditary cancer syndromes in MPT patients. CNVs affecting genes related to cancer are candidates to be involved in MPT risk. EPCAM / MSH2 deletions should be investigated in patients suspected to have LS. Gene enrichment related to the TP53 network is associated with MPT development. cnLOH and CNVs contribute to the risk of MPT development.

The AS04-adjuvanted human papillomavirus (HPV)16/18 vaccine, an L1-based vaccine, provides strong vaccine efficacy (VE) against vaccine-targeted type infections, and partial cross-protection to phylogenetically-related types, which may be affected by variant-level heterogeneity. We compared VE against incident HPV31, 33, 35, and 45 detections between lineages and SNPs in the L1 region among 2846 HPV-vaccinated and 5465 HPV-unvaccinated women through 11-years of follow-up in the Costa Rica HPV Vaccine Trial. VE was lower against HPV31-lineage-B (VE=60.7%;95%CI = 23.4%,82.8%) compared to HPV31-lineage-A (VE=94.3%;95%CI = 83.7%,100.0%) (VE-ratio = 0.64;95%CI = 0.25,0.90). Differential VE was observed at several lineage-associated HPV31-L1-SNPs, including a nonsynonymous substitution at position 6372 on the FG-loop, an important neutralization domain. For HPV35, the only SNP-level difference was at position 5939 on the DE-loop, with significant VE against nucleotide-G (VE=65.0%;95%CI = 28.0,87.8) but not for more the common nucleotide-A (VE=7.4%;95%CI = −34.1,36.7). Because of the known heterogeneity in precancer/cancer risk across cross-protected HPV genotype variants by race and region, our results of differential variant-level AS04-adjuvanted HPV16/18 vaccine efficacy has global health implications.

Abstract High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant's allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

Abstract Background Pediatric cancers are the leading cause of death by disease in children despite improved survival rates overall. The contribution of germline genetic susceptibility to pediatric cancer survivors has not been extensively characterized. We assessed the frequency of pathogenic or likely pathogenic (P/LP) variants in 5451 long-term pediatric cancer survivors from the Childhood Cancer Survivor Study. Methods Exome sequencing was conducted on germline DNA from 5451 pediatric cancer survivors (cases who survived ≥5 years from diagnosis; n = 5105 European) and 597 European cancer-free adults (controls). Analyses focused on comparing the frequency of rare P/LP variants in 237 cancer-susceptibility genes and a subset of 60 autosomal dominant high-to-moderate penetrance genes, for both case-case and case-control comparisons. Results Of European cases, 4.1% harbored a P/LP variant in high-to-moderate penetrance autosomal dominant genes compared with 1.3% in controls (2-sided P = 3 × 10-4). The highest frequency of P/LP variants was in genes typically associated with adult onset rather than pediatric cancers, including BRCA1/2, FH, PALB2, PMS2, and CDKN2A. A statistically significant excess of P/LP variants, after correction for multiple tests, was detected in patients with central nervous system cancers (NF1, SUFU, TSC1, PTCH2), Wilms tumor (WT1, REST), non-Hodgkin lymphoma (PMS2), and soft tissue sarcomas (SDHB, DICER1, TP53, ERCC4, FGFR3) compared with other pediatric cancers. Conclusion In long-term pediatric cancer survivors, we identified P/LP variants in cancer-susceptibility genes not previously associated with pediatric cancer as well as confirmed known associations. Further characterization of variants in these genes in pediatric cancer will be important to provide optimal genetic counseling for patients and their families.

Multiple primary thyroid cancer (TC) and breast cancer (BC) are commonly diagnosed, and the lifetime risk for these cancers is increased in patients with a positive family history of both TC and BC. Although this phenotype is partially explained by TP53 or PTEN mutations, a significant number of patients are negative for these alterations. We judiciously recruited patients diagnosed with BC and/or TC having a family history of these tumors and assessed their whole-exome sequencing. After variant prioritization, we selected MUS81 c.1292G>A (p.R431H) for further investigation. This variant was genotyped in a healthy population and sporadic BC/TC tissues and investigated at the protein level and cellular models. MUS81 c.1292G>A was the most frequent variant (25%) and the strongest candidate due to its function of double-strand break repair. This variant was confirmed in four relatives from two families. MUS81 p.R431H protein exhibited lower expression levels in tumors from patients positive for the germline variant, compared with wild-type BC, and normal breast and thyroid tissues. Using cell line models, we showed that c.1292G>A induced protein instability and affected DNA damage response. We suggest that MUS81 is a novel candidate involved in familial BC/TC based on its low frequency in healthy individuals and proven effect in protein stability.

Abstract APOBEC (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like) is a major mutagenic source in human papillomavirus positive oropharyngeal squamous cell carcinoma (HPV+ OPSCC). Why APOBEC mutations predominate in HPV+OPSCC remains an area of active investigation. Prevailing theories focus on APOBECs role as a viral restriction agent. APOBEC-induced mutations have been identified in both human cancers and HPV genomes, but whether they are directly linked in HPV+OPSCCs remains unknown. We performed sequencing of host somatic exomes, transcriptomes and HPV16 genomes from 79 HPV+ OPSCC samples, quantifying APOBEC mutational burden and activity in both the host and virus. APOBEC was the dominant mutational signature in somatic exomes. APOBEC vulnerable PIK3CA hotspot mutations were exclusively present in APOBEC enriched samples. In viral genomes, there was a mean (range) of 5 (0-29) mutations per genome. Mean (range) of APOBEC mutations in the viral genomes was 1 (0-5). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-variant allele frequency mutations, suggesting that APOBEC mutagenesis is actively occurring in viral genomes during infection. 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) suggesting that APOBEC mutagenesis in host and viral genomes are directly linked. Using paired sequencing of host somatic exomes, transcriptomes, and viral genomes from HPV+OPSCC samples, here, we show concordance between tumor and viral APOBEC mutagenesis, suggesting that APOBEC-mediated viral restriction results in off-target host-genome mutations. These data provide a missing link connecting APOBEC mutagenesis in host and virus and support a common mechanism driving APOBEC dysregulation. Competing Interest Statement The authors have declared no competing interest. Footnotes * The authors declare no potential conflicts of interest.