HPV16 E7 Genetic Conservation Is Critical to Carcinogenesis.

Although most cervical human papillomavirus type 16 (HPV16) infections become undetectable within 1-2 years, persistent HPV16 causes half of all cervical cancers. We used a novel HPV whole-genome sequencing technique to evaluate an exceptionally large collection of 5,570 HPV16-infected case-control samples to determine whether viral genetic variation influences risk of cervical precancer and cancer. We observed thousands of unique HPV16 genomes; very few women shared the identical HPV16 sequence, which should stimulate a careful re-evaluation of the clinical implications of HPV mutation rates, transmission, clearance, and persistence. In case-control analyses, HPV16 in the controls had significantly more amino acid changing variants throughout the genome. Strikingly, E7 was devoid of variants in precancers/cancers compared to higher levels in the controls; we confirmed this in cancers from around the world. Strict conservation of the 98 amino acids of E7, which disrupts Rb function, is critical for HPV16 carcinogenesis, presenting a highly specific target for etiologic and therapeutic research.

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

Human papillomavirus (HPV) positive oropharyngeal squamous cell carcinoma (HPV + OPSCC) is increasing in prevalence in the USA, as are cases of patients with multiple HPV + OPSCCs (mHPV + OPSCC). mHPV + OPSCCs present a unique opportunity to examine HPV + OPSCC mutation acquisition and evolution. We performed sequencing of the viral genome, somatic exome and somatic transcriptome from 8 patients each with 2 spatially distinct HPV + OPSCCs, and 37 'traditional' HPV + OPSCCs to first address if paired tumors are caused by the same viral isolate and next, if acquired alterations, and the underlying processes driving mutagenesis, are shared within pairs. All tumor pairs contained viral genomes from the same HPV type 16 sublineage and differed by 0-2 clonal single nucleotide polymorphisms (SNPs), suggesting infection with the same viral isolate. Despite this, there was significant discordance in expression profiles, mutational burden and mutational profiles between tumors in a pair, with only two pairs sharing any overlapping mutations (3/3343 variants). Within tumor pairs there was a striking discrepancy of mutational signatures, exemplified by no paired tumors sharing high APOBEC mutational burden. Here, leveraging mHPV + OPSCCs as a model system to study mutation acquisition in virally mediated tumors, in which the germline, environmental exposures, immune surveillance and tissue/organ type were internally controlled, we demonstrate that despite infection by the same viral isolate, paired mHPV + OPSCCs develop drastically different somatic alterations and even more strikingly, appear to be driven by disparate underlying mutational processes. Thus, despite a common starting point, HPV + OPSCCs evolve through variable mutational processes with resultant stochastic mutational profiles.

Association of HPV35 with cervical carcinogenesis among women of African ancestry: Evidence of viral-host interaction with implications for disease intervention.

HPV35 has been found in only ∼2% of invasive cervical cancers (ICC) worldwide but up to 10% in Sub-Saharan Africa, warranting further investigation and consideration of impact on preventive strategies. We studied HPV35 and ethnicity, in relation to the known steps in cervical carcinogenesis, using multiple large epidemiologic studies in the U.S. and internationally. Combining five U.S. studies, we measured HPV35 positivity and, in Northern California, observed HPV35 type-specific population prevalence and estimated 5-year risk of developing precancer when HPV35-positive. HPV35 genetic variation was examined for differences in carcinogenicity in 1053 HPV35+ cervical specimens from a U.S. cohort and an international collection. African-American women had more HPV35 (12.1% vs 5.1%, P < .001) and more HPV35-associated precancers (7.4% vs 2.1%, P < .001) compared to other ethnicities. Precancer risks after HPV35 infection did not vary by ethnicity (global P = .52). The HPV35 A2 sublineage showed an increased association with precancer/cancer in African-Americans (OR = 5.6 vs A1, 95% CI = 1.3-24.8) and A2 was more prevalent among ICC in Africa than other world regions (41.9% vs 10.4%, P < .01). Our analyses support a strong link between HPV35 and cervical carcinogenesis in women of African ancestry. Current HPV vaccines cover the majority of cervical precancer/cancer across all ethnic groups; additional analyses are required to determine whether the addition of HPV35 to the already highly effective nine-valent HPV vaccine would provide better protection for women in Africa or of African ancestry.

Differential microRNA regulation of HLA-C expression and its association with HIV control.

The HLA-C locus is distinct relative to the other classical HLA class I loci in that it has relatively limited polymorphism, lower expression on the cell surface, and more extensive ligand-receptor interactions with killer-cell immunoglobulin-like receptors. A single nucleotide polymorphism (SNP) 35 kb upstream of HLA-C (rs9264942; termed -35) associates with control of HIV, and with levels of HLA-C messenger RNA transcripts and cell-surface expression, but the mechanism underlying its varied expression is unknown. We proposed that the -35 SNP is not the causal variant for differential HLA-C expression, but rather is marking another polymorphism that directly affects levels of HLA-C. Here we show that variation within the 3' untranslated region (UTR) of HLA-C regulates binding of the microRNA hsa-miR-148 to its target site, resulting in relatively low surface expression of alleles that bind this microRNA and high expression of HLA-C alleles that escape post-transcriptional regulation. The 3' UTR variant associates strongly with control of HIV, potentially adding to the effects of genetic variation encoding the peptide-binding region of the HLA class I loci. Variation in HLA-C expression adds another layer of diversity to this highly polymorphic locus that must be considered when deciphering the function of these molecules in health and disease.

Genetic modifiers of neurofibromatosis type 1-associated café-au-lait macule count identified using multi-platform analysis.

Neurofibromatosis type 1 (NF1) is an autosomal dominant, monogenic disorder of dysregulated neurocutaneous tissue growth. Pleiotropy, variable expressivity and few NF1 genotype-phenotype correlates limit clinical prognostication in NF1. Phenotype complexity in NF1 is hypothesized to derive in part from genetic modifiers unlinked to the NF1 locus. In this study, we hypothesized that normal variation in germline gene expression confers risk for certain phenotypes in NF1. In a set of 79 individuals with NF1, we examined the association between gene expression in lymphoblastoid cell lines with NF1-associated phenotypes and sequenced select genes with significant phenotype/expression correlations. In a discovery cohort of 89 self-reported European-Americans with NF1 we examined the association between germline sequence variants of these genes with café-au-lait macule (CALM) count, a tractable, tumor-like phenotype in NF1. Two correlated, common SNPs (rs4660761 and rs7161) between DPH2 and ATP6V0B were significantly associated with the CALM count. Analysis with tiled regression also identified SNP rs4660761 as significantly associated with CALM count. SNP rs1800934 and 12 rare variants in the mismatch repair gene MSH6 were also associated with CALM count. Both SNPs rs7161 and rs4660761 (DPH2 and ATP6V0B) were highly significant in a mega-analysis in a combined cohort of 180 self-reported European-Americans; SNP rs1800934 (MSH6) was near-significant in a meta-analysis assuming dominant effect of the minor allele. SNP rs4660761 is predicted to regulate ATP6V0B, a gene associated with melanosome biology. Individuals with homozygous mutations in MSH6 can develop an NF1-like phenotype, including multiple CALMs. Through a multi-platform approach, we identified variants that influence NF1 CALM count.

Pemphigus is associated with KIR3DL2 expression levels and provides evidence that KIR3DL2 may bind HLA-A3 and A11 in vivo.

Although HLA-A3 and A11 have been reported to be ligands for KIR3DL2, evidence for any in vivo relevance of this interaction is still missing. To explore the functional importance of KIR3DL2 allelic variation, we analyzed the autoimmune disease pemphigus foliaceus, previously associated (lower risk) with activating KIR genes. KIR3DL2*001 was increased in patients (odds ratio (OR) = 2.04; p = 0.007). The risk was higher for the presence of both KIR3DL2*001 and HLA-A3 or A11 (OR = 3.76, p = 0.013), providing the first evidence that HLA-A3 and A11 may interact with KIR3DL2 in vivo. The nonsynonymous single nucleotide polymorphism 1190T (rs3745902) was associated with protection (OR = 0.52, p = 0.018). This SNP results in a threonine-to-methionine substitution. Individuals who have methionine in this position exhibit a lower percentage of KIR3DL2-positive natural killer (NK) cells and also lower intensity of KIR3DL2 on expressing natural killer cells; additionally, we show that the expression of KIR3DL2 is independent of other killer cell immunoglobulin-like receptors. Pemphigus foliaceus is a very unique complex disease strongly associated with immune-related genes. It is the only autoimmune disease known to be endemic, showing a strong correlation with environmental factors. Our data demonstrate that this relatively unknown autoimmune disease may facilitate understanding of the molecular mechanisms of KIR3DL2 ligand recognition.

Whole exome sequencing in families at high risk for Hodgkin lymphoma: identification of a predisposing mutation in the KDR gene.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

Genetic interplay between HLA-C and MIR148A in HIV control and Crohn disease.

Variation in the 3' untranslated region (3'UTR) of the HLA-C locus determines binding of the microRNA Hsa-miR-148a, resulting in lower cell surface expression of alleles that bind miR-148a relative to those alleles that escape its binding. The HLA-C 3'UTR variant was shown to associate with HIV control, but like the vast majority of disease associations in a region dense with causal candidates, a direct effect of HLA-C expression level on HIV control was not proven. We demonstrate that a MIR148A insertion/deletion polymorphism associates with its own expression levels, affecting the extent to which HLA-C is down-regulated, the level of HIV control, and the risk of Crohn disease only among those carrying an intact miR-148a binding site in the HLA-C 3'UTR. These data illustrate a direct effect of HLA-C expression level on HIV control that cannot be attributed to other HLA loci in linkage disequilibrium with HLA-C and highlight the rich complexity of genetic interactions in human disease.

Characterization of T gene sequence variants and germline duplications in familial and sporadic chordoma.

Chordoma is a rare bone cancer that is believed to originate from notochordal remnants. We previously identified germline T duplication as a major susceptibility mechanism in several chordoma families. Recently, a common genetic variant in T (rs2305089) was significantly associated with the risk of sporadic chordoma. We sequenced all T exons in 24 familial cases and 54 unaffected family members from eight chordoma families (three with T duplications), 103 sporadic cases, and 160 unrelated controls. We also measured T copy number variation in all sporadic cases. We confirmed the association between the previously reported variant rs2305089 and risk of familial [odds ratio (OR) = 2.6, 95% confidence interval (CI) = 0.93, 7.25, P = 0.067] and sporadic chordoma (OR = 2.85, 95% CI = 1.89, 4.29, P < 0.0001). We also identified a second common variant, rs1056048, that was strongly associated with chordoma in families (OR = 4.14, 95% CI = 1.43, 11.92, P = 0.0086). Among sporadic cases, another common variant (rs3816300) was significantly associated with risk when jointly analyzed with rs2305089. The association with rs3816300 was significantly stronger in cases with early age onset. In addition, we identified three rare variants that were only observed among sporadic chordoma cases, all of which have potential functional relevance based on in silico predictions. Finally, we did not observe T duplication in any sporadic chordoma case. Our findings further highlight the importance of the T gene in the pathogenesis of both familial and sporadic chordoma and suggest a complex susceptibility related to T.

Somatic mutations in 3929 HPV positive cervical cells associated with infection outcome and HPV type.

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.

Oropharyngeal Squamous Cell Carcinoma Morphology and Subtypes by Human Papillomavirus Type and by 16 Lineages and Sublineages.

Oropharyngeal squamous cell carcinoma (SCC) is increasing in incidence and, in Western countries, strongly associated with transcriptionally-active high-risk human papillomavirus (HPV). Within HPV-positive tumors, there is wide morphologic diversity with numerous histologic subtypes of SCC. There are also variable degrees of keratinization, anaplasia, stromal fibrosis, and maturing squamous differentiation. Unlike in the uterine cervix, where associations between HPV types and lineages/sublineages within types have been investigated with some clear correlations identified, little to no data exists for oropharyngeal SCC. In this study, for a large cohort of oropharyngeal SCC patients, we performed RTPCR for high-risk HPV. For the HPV positive patients, we sequenced the DNA of the entire HPV16 genome and determined lineages and sublineages, correlating HPV status, genotype, and HPV16 lineages/sublineages with SCC subtype and various histologic features. Of the 259 patients, 224 (86.5%) were high-risk HPV positive, of which 210/224 (93.8%) were HPV type 16 and 6/224 (2.7%) HPV type 33. Of the four HPV16 lineages, A was the most frequent (192/214 or 89.8%) and of the HPV16 A sublineages, A1 was the most frequent (112/210 or 53.3%). Patients with HPV negative tumors were more often keratinizing vs other types (23/35 or 65.7%) and thus more likely to have more maturing squamous differentiation and stromal desmoplasia. There was no significant correlation between HPV type (16 versus other), between HPV16 lineage (A versus others), or HPV16 A sublineages (A1 or A2 versus others) and morphologic type of SCC nor the various morphologic features of anaplasia/multinucleation, degree of keratinization, nor amount of stromal desmoplasia. In summary, in our cohort, there was no correlation between the type of HPV, the HPV 16 lineage or sublineage, and any of the histologic features or morphologic SCC subtypes.

Age-specific prevalence of human papillomavirus and abnormal cytology at baseline in a diverse statewide prospective cohort of individuals undergoing cervical cancer screening in Mississippi.

BACKGROUND: Mississippi (MS) has among the highest rates of cervical cancer incidence and mortality in the United States, with disproportionately higher rates among Blacks compared to Whites. Here, we evaluate the prevalence of high-risk human papillomavirus (HPV) and abnormal cytology in a representative baseline sample from a diverse statewide cohort of individuals attending cervical screening in MS from the STRIDES Study (STudying Risk to Improve DisparitiES in cervical cancer). METHODS: We included individuals aged 21-65 years undergoing screening at the University of Mississippi Medical Center (UMMC) and the Mississippi State Department of Health (MSDH) from May to November 2018. We calculated age-specific HPV prevalence, overall and by partial HPV16/18 genotyping, and abnormal cytology by race. RESULTS: A total of 6871 individuals (mean age 35.7 years) were included. HPV prevalence was 25.6% and higher in Blacks (28.0%) compared to Whites (22.4%). HPV prevalence was significantly higher in Blacks aged 21-24 years (50.2%) and 30-34 years (30.2%) compared to Whites in the same age groups (32.1% and 20.7%; p < 0.0001, respectively). The prevalence of high-grade cytologic abnormalities, a cytologic sign of cervical precancer, peaked earlier in Blacks (ages 25-29) compared to Whites (35-39). For comparison, we also analyzed HPV prevalence data from the National Health and Nutrition Examination Survey (NHANES, 2013-2016) and observed similar racial differences in HPV prevalence among women aged 21-24 years. CONCLUSIONS: Our findings suggest that Blacks undergoing cervical cancer screening in MS have higher prevalence of other high-risk 12 HPV types at younger ages and experience an earlier peak of high-grade cytologic abnormalities compared to Whites.

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.

Phylogenomic Analysis of Human Papillomavirus Type 31 and Cervical Carcinogenesis: A Study of 2093 Viral Genomes.

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.

Novel dithiolethione-modified nonsteroidal anti-inflammatory drugs in human hepatoma HepG2 and colon LS180 cells.

PURPOSE: Nonsteroidal anti-inflammatory drugs (NSAID) are promising chemopreventive agents against colon and other cancers. However, the molecular basis mediated by NSAIDs for chemoprevention has not been fully elucidated. Environmental carcinogens induce DNA mutation and cellular transformation; therefore, we examined the effect of NSAIDs on carcinogenesis mediated by the aryl hydrocarbon receptor signaling pathway. In this study, we investigated the activities of a new class of NSAIDs containing dithiolethione moieties (S-NSAID) on both arms of carcinogenesis. EXPERIMENTAL DESIGN: We investigated the effects of the S-NSAIDs, S-diclofenac and S-sulindac, on carcinogen activation and detoxification mechanisms in human hepatoma HepG2 and human colonic adenocarcinoma LS180 cells. RESULTS: We found that S-diclofenac and S-sulindac inhibited the activity and expression of the carcinogen activating enzymes, cytochromes P-450 (CYP) CYP1A1, CYP1B1, and CYP1A2. Inhibition was mediated by transcriptional regulation of the aryl hydrocarbon receptor (AhR) pathway. The S-NSAIDs down-regulated carcinogen-induced expression of CYP1A1 heterogeneous nuclear RNA, a measure of transcription rate. Both compounds blocked carcinogen-activated AhR from binding to the xenobiotic responsive element as shown by chromatin immunoprecipitation. S-diclofenac and S-sulindac inhibited carcinogen-induced CYP enzyme activity through direct inhibition as well as through decreased transcriptional activation of the AhR. S-sulindac induced expression of several carcinogen detoxification enzymes of the glutathione cycle including glutathione S-transferase A2, glutamate cysteine ligase catalytic subunit, glutamate cysteine ligase modifier subunit, and glutathione reductase. CONCLUSIONS: These results indicate that S-diclofenac and S-sulindac may serve as effective chemoprevention agents by favorably balancing the equation of carcinogen activation and detoxification mechanisms.

Genome-wide association study of café-au-lait macule number in neurofibromatosis type 1.

BACKGROUND: Neurofibromatosis type 1 (NF1) is a tumor-predisposition disorder that arises due to pathogenic variants in tumor suppressor NF1. NF1 has variable expressivity that may be due, at least in part, from heritable elements such as modifier genes; however, few genetic modifiers have been identified to date. METHODS: In this study, we performed a genome-wide association analysis of the number of café-au-lait macules (CALM) that are considered a tumor-like trait as a clinical phenotype modifying NF1. RESULTS: A borderline genome-wide significant association was identified in the discovery cohort (CALM1, N = 112) between CALM number and rs12190451 (and rs3799603, r2  = 1.0; p = 7.4 × 10-8 ) in the intronic region of RPS6KA2. Although, this association was not replicated in the second cohort (CALM2, N = 59) and a meta-analysis did not show significantly associated variants in this region, a significant corroboration score (0.72) was obtained for the RPS6KA2 signal in the discovery cohort (CALM1) using Complementary Pairs Stability Selection for Genome-Wide Association Studies (ComPaSS-GWAS) analysis, suggesting that the lack of replication may be due to heterogeneity of the cohorts rather than type I error. CONCLUSION: rs12190451 is located in a melanocyte-specific enhancer and may influence RPS6KA2 expression in melanocytes-warranting further functional studies.

The D2 and D3 Sublineages of Human Papilloma Virus 16-Positive Cervical Cancer in Guatemala Differ in Integration Rate and Age of Diagnosis.

Human papillomavirus (HPV) 16 displays substantial sequence variation; four HPV16 lineages (A, B, C, and D) have been described as well as multiple sublineages. To identify molecular events associated with HPV16 carcinogenesis, we evaluated viral variation, the integration of HPV16, and somatic mutation in 96 cervical cancer samples from Guatemala. A total of 65% (62/96) of the samples had integrated HPV16 sequences and integration was associated with an earlier age of diagnosis and premenopausal disease. HPV16 integration sites were broadly distributed in the genome, but in one tumor, HPV16 integrated into the promoter of the IFN regulatory factor 4 (IRF4) gene, which plays an important role in the regulation of the IFN response to viral infection. The HPV16 D2 and D3 sublineages were found in 23% and 30% of the tumors, respectively, and were significantly associated with adenocarcinoma. D2-positive tumors had a higher rate of integration, earlier age of diagnosis, and a lower rate of somatic mutation, whereas D3-positive tumors were less likely to integrate, had later age of diagnosis, and exhibited a higher rate of somatic mutation. In conclusion, Guatemalan cervical tumors have a high frequency of very high-risk HPV16 D2 and D3 sublineages harboring distinct histology, which may help guide future therapeutic strategies to target the tumor and reduce recurrence. SIGNIFICANCE: This study details the biological and molecular properties of the most pathogenic forms of HPV16, the cause of the majority of cervical cancers.

Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance.

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.

Human papillomavirus 16 sub-lineage dispersal and cervical cancer risk worldwide: Whole viral genome sequences from 7116 HPV16-positive women.

BACKGROUND: Human papillomavirus (HPV)16 can be separated into genetic sub-lineages (A1-4, B1-4, C1-4, D1-4) which may have differential cervical cancer risk. METHODS: A next-generation sequencing assay was used to whole-genome sequence 7116 HPV16-positive cervical samples from well-characterised international epidemiological studies, including 2076 controls, 1878 squamous cell carcinoma (SCC) and 186 adenocarcinoma/adenosquamous cell carcinoma (ADC), and to assign HPV16 sub-lineage. Logistic regression was used to estimate region-stratified country-adjusted odds ratios (OR) and 95%CI. RESULTS: A1 was the most globally widespread sub-lineage, with others showing stronger regional specificity (A3 and A4 for East Asia, B1-4 and C1-4 for Africa, D2 for the Americas, B4, C4 and D4 for North Africa). Increased cancer risks versus A1 were seen for A3, A4 and D (sub)lineages in regions where they were common: A3 in East Asia (OR=2.2, 95%CI:1.0-4.7); A4 in East Asia (6.6, 3.1-14.1) and North America (3.8, 1.7-8.3); and D in North (6.2, 4.1-9.3) and South/Central America (2.2, 0.8-5.7), where D lineages were also more frequent in ADC than SCC (3.2, 1.5-6.5; 12.1, 5.7-25.6, respectively). CONCLUSIONS: HPV16 genetic variation can strongly influence cervical cancer risk. However, burden of cervical cancer attributable to different sub-lineages worldwide is largely driven by historical HPV16 sub-lineage dispersal.