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.

Insights into the mechanisms and structure of breakage-fusion-bridge cycles in cervical cancer using long-read sequencing.

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cancer cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and our study applies long-read sequencing to this important chromosomal rearrangement type. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes, and only one BFB breakpoint showed chromothripsis. Five cell lines have a chromosomal region 11q BFB event, with YAP1-BIRC3-BIRC2 amplification. Indeed, YAP1 amplification is associated with a 10-year-earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that individuals with cervical cancer and YAP1-BIRC3-BIRC2 amplification, especially those of African ancestry, might benefit from targeted therapy. In summary, we uncovered valuable insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

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.

An optimized protocol for retina single-cell RNA sequencing.

Purpose: Single-cell RNA sequencing (scRNA-seq) is a powerful technique used to explore gene expression at the single cell level. However, appropriate preparation of samples is essential to obtain the most information out of this transformative technology. Generating high-quality single-cell suspensions from the retina is critical to preserve the native expression profile that will ensure meaningful transcriptome data analysis. Methods: We modified the conditions for rapid and optimal dissociation of retina sample preparation. We also included additional filtering steps in data analysis for retinal scRNA-seq. Results: We report a gentle method for dissociation of the mouse retina that minimizes cell death and preserves cell morphology. This protocol also results in detection of higher transcriptional complexity. In addition, the modified computational pipeline leads to better-quality single-cell RNA-sequencing data in retina samples. We also demonstrate the advantages and limitations of using fresh versus frozen retinas to prepare cell or nuclei suspensions for scRNA-seq. Conclusions: We provide a simple yet robust and reproducible protocol for retinal scRNA-seq analysis, especially for comparative studies.

Correction to: Low-cost HPV testing and the prevalence of cervical infection in asymptomatic populations in Guatemala.

Following publication of the original article [1], an error was reported in the tagging of Joël Fokom Domgue in the author group. The tagging in this correction article has been fixed.

Low-cost HPV testing and the prevalence of cervical infection in asymptomatic populations in Guatemala.

BACKGROUND: A low cost and accurate method for detecting high-risk (HR) human papillomavirus (HPV) is important to permit HPV testing for cervical cancer prevention. We used a commercially available HPV method (H13, Hybribio) which was documented to function accurately in a reduced volume of cervical specimen to determine the most prevalent HPV types and the distribution of HPV infections in over 1795 cancer-free women in Guatemala undergoing primary screening for cervical cancer by cytology. METHODS: HR-HPV detection was attempted in cervical samples from 1795 cancer-free women receiving Pap smears using the Hybribio™ real-time PCR assay of 13 HR types. The test includes a globin gene internal control. HPV positive samples were sequenced to determine viral type. Age-specific prevalence of HPV was also assessed in the study population. RESULTS: A total of 13% (226/1717) of women tested HPV+, with 78 samples (4.3%) failing to amplify the internal control. The highest prevalence was found in younger women (< 30 years, 22%) and older ones (≥60 years, 15%). The six most common HR-HPV types among the 148 HPV+ typed were HPV16 (22%), HPV18 (11%), HPV39 (11%), HPV58 (10%), HPV52 (8%), and HPV45 (8%). CONCLUSIONS: In this sample of cancer free women in Guatemala, HPV16 was the most prevalent HR type in Guatemala and the age-specific prevalence curve peaked in younger ages. Women in the 30-59-year age groups had a prevalence of HR-HPV of 8%, however, larger studies to better describe the epidemiology of HPV in Guatemala are needed.

Novel and known ribosomal causes of Diamond-Blackfan anaemia identified through comprehensive genomic characterisation.

BACKGROUND: Diamond-Blackfan anaemia (DBA) is an inherited bone marrow failure syndrome (IBMFS) characterised by erythroid hypoplasia. It is associated with congenital anomalies and a high risk of developing specific cancers. DBA is caused predominantly by autosomal dominant pathogenic variants in at least 15 genes affecting ribosomal biogenesis and function. Two X-linked recessive genes have been identified. OBJECTIVES: We aim to identify the genetic aetiology of DBA. METHODS: Of 87 families with DBA enrolled in an institutional review board-approved cohort study (ClinicalTrials.gov Identifier:NCT00027274), 61 had genetic testing information available. Thirty-five families did not have a known genetic cause and thus underwent comprehensive genomic evaluation with whole exome sequencing, deletion and CNV analyses to identify their disease-associated pathogenic variant. Controls for functional studies were healthy mutation-negative individuals enrolled in the same study. RESULTS: Our analyses uncovered heterozygous pathogenic variants in two previously undescribed genes in two families. One family had a non-synonymous variant (p.K77N) in RPL35; the second family had a non-synonymous variant (p. L51S) in RPL18. Both of these variants result in pre-rRNA processing defects. We identified heterozygous pathogenic variants in previously known DBA genes in 16 of 35 families. Seventeen families who underwent genetic analyses are yet to have a genetic cause of disease identified. CONCLUSIONS: Overall, heterozygous pathogenic variants in ribosomal genes were identified in 44 of the 61 families (72%). De novo pathogenic variants were observed in 57% of patients with DBA. Ongoing studies of DBA genomics will be important to understand this complex disorder.

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.

Whole-exome sequencing and functional studies identify RPS29 as a novel gene mutated in multicase Diamond-Blackfan anemia families.

Diamond-Blackfan anemia (DBA) is a cancer-prone inherited bone marrow failure syndrome. Approximately half of DBA patients have a germ-line mutation in a ribosomal protein gene. We used whole-exome sequencing to identify disease-causing genes in 2 large DBA families. After filtering, 1 nonsynonymous mutation (p.I31F) in the ribosomal protein S29 (RPS29[AUQ1]) gene was present in all 5 DBA-affected individuals and the obligate carrier, and absent from the unaffected noncarrier parent in 1 DBA family. A second DBA family was found to have a different nonsynonymous mutation (p.I50T) in RPS29. Both mutations are amino acid substitutions in exon 2 predicted to be deleterious and resulted in haploinsufficiency of RPS29 expression compared with wild-type RPS29 expression from an unaffected control. The DBA proband with the p.I31F RPS29 mutation had a pre-ribosomal RNA (rRNA) processing defect compared with the healthy control. We demonstrated that both RPS29 mutations failed to rescue the defective erythropoiesis in the rps29(-/-) mutant zebra fish DBA model. RPS29 is a component of the small 40S ribosomal subunit and essential for rRNA processing and ribosome biogenesis. We uncovered a novel DBA causative gene, RPS29, and showed that germ-line mutations in RPS29 can cause a defective erythropoiesis phenotype using a zebra fish model.

A comprehensive resequence-analysis of 250 kb region of 8q24.21 in men of African ancestry.

BACKGROUND: Genome-wide association studies (GWAS) have identified that a ∼1 M region centromeric to the MYC oncogene on chromosome 8q24.21 harbors at least five independent loci associated with prostate cancer risk and additional loci associated with cancers of breast, colon, bladder, and chronic lymphocytic leukemia (CLL). Because GWAS identify genetic markers that may be indirectly associated with disease, fine-mapping based on sequence analysis provides important insights into patterns of linkage disequilibrium (LD) and is critical in defining the optimal variants to nominate for biological follow-up. METHODS: To catalog variation in individuals of African ancestry, we resequenced a region (250 kb; chr8:128,050, 768­128, 300,801, hg19) containing several prostate cancer susceptibility loci as well as a locus associated with CLL. Our samples included 78 individuals from Ghana and 47 of African-Americans from Johns Hopkins University. RESULTS: After quality control metrics were applied to next-generation sequence data, 1,838 SNPs were identified. Of these, 285 were novel and not yet reported in any public database. Using genotypes derived from sequencing, we refined the LD and recombination hotspots within the region and determined a set of tag SNPs to be used in future fine-mapping studies. Based on LD, we annotated putative risk loci and their surrogates using ENCODE data, which should help guide laboratory studies. CONCLUSIONS: In comparison to the 1000 Genome Project data, we have identified additional variants that could be important in establishing priorities for future functional work designed to explain the biological basis of associations between SNPs and both prostate cancer and CLL.

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.

The new sequencer on the block: comparison of Life Technology's Proton sequencer to an Illumina HiSeq for whole-exome sequencing.

We assessed the performance of the new Life Technologies Proton sequencer by comparing whole-exome sequence data in a Centre d'Etude du Polymorphisme Humain trio (family 1463) to the Illumina HiSeq instrument. To simulate a typical user's results, we utilized the standard capture, alignment and variant calling methods specific to each platform. We restricted data analysis to include the capture region common to both methods. The Proton produced high quality data at a comparable average depth and read length, and the Ion Reporter variant caller identified 96 % of single nucleotide polymorphisms (SNPs) detected by the HiSeq and GATK pipeline. However, only 40 % of small insertion and deletion variants (indels) were identified by both methods. Usage of the trio structure and segregation of platform-specific alleles supported this result. Further comparison of the trio data with Complete Genomics sequence data and Illumina SNP microarray genotypes documented high concordance and accurate SNP genotyping of both Proton and Illumina platforms. However, our study underscored the problem of accurate detection of indels for both the Proton and HiSeq platforms.

Extrachromosomal Amplification of Human Papillomavirus Episomes Is a Mechanism of Cervical Carcinogenesis.

SIGNIFICANCE: Multimers of the HPV genome are generated in cervical tumors replicating as extrachromosomal episomes, which is associated with deletion and rearrangement of the HPV genome and provides a mechanism for oncogenesis without integration.

Insights into the Mechanisms and Structure of Breakage-Fusion-Bridge Cycles in Cervical Cancer using Long-Read Sequencing.

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and this is one of the first analyses of these events using long-read sequencing. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes and only one BFB breakpoint showed chromothripsis. Five cell lines have a Chr11q BFB event, with YAP1/BIRC2/BIRC3 gene amplification. Indeed, YAP1 amplification is associated with a 10-year earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that cervical cancer patients with YAP1/BIRC2/BIRC3-amplification, especially those of African American ancestry, might benefit from targeted therapy. In summary, we uncovered new insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

HPV16 E7 Nucleotide Variants Found in Cancer-Free Subjects Affect E7 Protein Expression and Transformation.

The human papillomavirus (HPV) type 16 E7 oncogene is critical to carcinogenesis and highly conserved. Previous studies identified a preponderance of non-synonymous E7 variants amongst HPV16-positive cancer-free controls compared to those with cervical cancer. To investigate the function of E7 variants, we constructed full-length HPV16 E7 genes and tested variants at positions H9R, D21N, N29S, E33K, T56I, D62N, S63F, S63P, T64M, E80K, D81N, P92L, and P92S (found only in controls); D14E, N29H cervical intraepithelial neoplasia (CIN2), and P6L, H51N, R77S (CIN3). We determined the steady-state level of cytoplasmic and nuclear HPV16 E7 protein. All variants from controls showed a reduced level of E7 protein, with 7/13 variants having lower protein levels. In contrast, 2/3 variants from the CIN3 precancer group had near-wild type E7 levels. We assayed the activity of representative variants in stably transfected NIH3T3 cells. The H9R, E33K, P92L, and P92S variants found in control subjects had lower transforming activity than D14E and N29H variants (CIN2), and the R77S (CIN3) had activity only slightly reduced from wild-type E7. In addition, R77S and WT E7 caused increased migration of NIH3T3 cells in a wound-healing assay compared with H9R, E33K, P92L, and P92S (controls) and D14E (CIN2). These data provide evidence that the E7 variants found in HPV16-positive cancer-free women are partially defective for transformation and cell migration, further demonstrating the importance of fully active E7 in cancer development.

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.

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.