Whole-Genome Sequencing and Variant Analysis of Human Papillomavirus 16 Infections

Differential long-term bivalent HPV vaccine cross-protection by variants in the Costa Rica HPV vaccine trial

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.

The clinical and epidemiological impacts of whole genomic sequencing on bacterial and virological agents

Whole Genome Sequencing (WGS) is a molecular biology tool consisting in the sequencing of the entire genome of a given organism. Due to its ability to provide the finest available resolution of bacterial and virological genetics, it is used at several levels in the field of infectiology. On an individual scale and through application of a single technique, it enables the typological identification and characterization of strains, the characterization of plasmids, and enhanced search for resistance genes and virulence factors. On a collective scale, it enables the characterization of strains and the determination of phylogenetic links between different microorganisms during community outbreaks and healthcare-associated epidemics. The information provided by WGS enables real-time monitoring of strain-level epidemiology on a worldwide scale, and facilitates surveillance of the resistance dissemination and the introduction or emergence of pathogenic variants in humans or their environment. There are several possible approaches to completion of an entire genome. The choice of one method rather than another is essentially dictated by the matrix, either a clinical sample or a culture isolate, and the clinical objective. WGS is an advanced technology that remains costly despite a gradual decrease in its expenses, potentially hindering its implementation in certain laboratories and thus its use in routine microbiology. Even though WGS is making steady inroads as a reference method, efforts remain needed in view of so harmonizing its interpretations and decreasing the time to generation of conclusive results.

Human papillomavirus genomics: Understanding carcinogenicity

Human papillomavirus (HPV) causes virtually all cervical cancers and many cancers at other anatomical sites in both men and women. However, only 12 of 448 known HPV types are currently classified as carcinogens, and even the most carcinogenic type - HPV16 - only rarely leads to cancer. HPV is therefore necessary but insufficient for cervical cancer, with other contributing factors including host and viral genetics. Over the last decade, HPV whole genome sequencing has established that even fine-scale within-type HPV variation influences precancer/cancer risks, and that these risks vary by histology and host race/ethnicity. In this review, we place these findings in the context of the HPV life cycle and evolution at various levels of viral diversity: between-type, within-type, and within-host. We also discuss key concepts necessary for interpreting HPV genomic data, including features of the viral genome; events leading to carcinogenesis; the role of APOBEC3 in HPV infection and evolution; and methodologies that use deep (high-coverage) sequencing to characterize within-host variation, as opposed to relying on a single representative (consensus) sequence. Given the continued high burden of HPV-associated cancers, understanding HPV carcinogenicity remains important for better understanding, preventing, and treating cancers attributable to infection.

Longitudinal follow-up of HPV16 sequence after cervical infection: Low intrahost variation and no correlation with clinical evolution

Human papillomavirus (HPV) 16 exhibits different variants that may differ in their carcinogenic risk. To identify some high-risk variants, we sequenced and compared HPV16 whole genomes obtained from a longitudinal cohort of 34 HPV16-infected women who had either spontaneously cleared their infection (clearance group or "C"), or developed cervical high-grade lesions following a viral persistence (group persistence or "P"). Phylogenetic analysis of paired samples obtained at the beginning (C0 or P0) and at the end (C2 or P2) of the follow-up (median intervals between C0-C2 and between P0-P2 were 16 and 36.5 months, respectively) revealed a low genetic variability within the host compared to the genetic interhost diversity. By comparing our HPV16 sequences to a reference sequence, we observed 301 different substitutions, more often transitions (60.9%) than transversions (39.1%), that occurred throughout the viral genome, but with a low frequency in E6 and E7 oncogenes (10 and 9 substitutions), suggesting a high conservation of these genes. Deletions and insertions were mostly observed in intergenic regions of the virus. The only significant substitution found between the subgroups C2 and P2 was observed in the L2 gene (L330F), with an unclear biological relevance. Our results suggest a low longitudinal intrahost evolution of HPV16 sequences and no correlation between genetic variations and clinical evolution.

Prevalence and distribution of human papillomavirus (HPV) in Luoyang city of Henan province during 2015-2021 and the genetic variability of HPV16 and 52

Background

Persistent high-risk Human papillomavirus (HPV) subtypes infection has been implicated as a causative of cervical cancer. Distribution and genotypes of HPV infection among females and their variations would assist in the formulation of preventive strategy for cervical cancer. The purpose of the present study is to investigate the prevalence of HPV among females in central China.

Methods

The distribution and genotypes of HPV among 9943 females attending the gynecological examinations in central of China during 2015–2021 were investigated. HPV genotypes were detected using a commercial kit. Nucleotides sequences of L1, E6 and E7 genes in HPV16 or HPV52 positive samples collected in 2021 were amplified by polymerase chain reaction (PCR). Variations of L1, E6 and E7 in HPV16 and HPV52 were gained by sequencing and compared with the reference sequence. Sublineages of HPV16 and HPV52 were determined by the construction of phylogenetic tree based on L1 gene.

Results

The overall prevalence of HPV infection was 22.81%, with the infection rate of high-risk human papillomavirus (HR-HPV) was 19.02% and low-risk human papillomavirus (LR-HPV) was 6.40%. The most top five genotypes of HPV infection were HPV16 (7.49%), HPV52 (3.04%), HPV58 (2.36%), HPV18 (1.65%) and HPV51 (1.61%). Plots of the age-infection rate showed that the single HPV, multiple HPV, HR-HPV, LR-HPV infection revealed the same tendency with two peaks of HPV infection were observed among females aged ≤ 20 year-old and 60–65 year-old. The predominant sublineage of HPV16 was A1 and B2 for HPV52. For HPV16, The most prevalent mutations were T266A (27/27) and N181T (7/27) for L1, D32E for E6 and S63F for E7 in HPV16. For HPV52, all of the nucleotide changes were synonymous mutation in L1 (except L5S) and E7 genes. The K93R mutation was observed in most HPV52 E6 protein.

Conclusions

The present study provides basic information about the distribution, genotypes and variations of HPV among females population in Henan province, which would assist in the formulation of preventive strategies and improvements of diagnostic probe and vaccine for HPV in this region.

HPV16 and HPV18 type-specific APOBEC3 and integration profiles in different diagnostic categories of cervical samples

Human papillomavirus (HPV) 16 and 18 are the most predominant types in cervical cancer. Only a small fraction of HPV infections progress to cancer, indicating that additional factors and genomic events contribute to the carcinogenesis, such as minor nucleotide variation caused by APOBEC3 and chromosomal integration. We analysed intra-host minor nucleotide variants (MNVs) and integration in HPV16 and HPV18 positive cervical samples with different morphology. Samples were sequenced using an HPV whole genome sequencing protocol TaME-seq. A total of 80 HPV16 and 51 HPV18 positive samples passed the sequencing depth criteria of 300× reads, showing the following distribution: non-progressive disease (HPV16 n = 21, HPV18 n = 12); cervical intraepithelial neoplasia (CIN) grade 2 (HPV16 n = 27, HPV18 n = 9); CIN3/adenocarcinoma in situ (AIS) (HPV16 n = 27, HPV18 n = 30); cervical cancer (HPV16 n = 5). Similar numbers of MNVs in HPV16 and HPV18 samples were observed for most viral genes, with the exception of HPV18 E4 with higher numbers across clinical categories. APOBEC3 signatures were observed in HPV16 lesions, while similar mutation patterns were not detected for HPV18. The proportion of samples with integration was 13% for HPV16 and 59% for HPV18 positive samples, with a noticeable portion located within or close to cancer-related genes.

Exploring the roles of HPV16 variants in head and neck squamous cell carcinoma: current challenges and opportunities

The incidence of squamous cell carcinomas of the head and neck (HNSCC) is consistently increasing, in association with human papillomavirus (HPV) infection, especially HPV16. HPV variants show heterogeneity in the pathogenicity of cervical cancer, but little has been established about their relevance on HNSCC. This review addresses the distribution of HPV16 variants in HNSCC and their potential contribution to clinical practice. A search was performed in PubMed using the keywords HNSCC HPV16 variants. Sixty articles were identified between 2000 and 2020 and 9 articles were selected for a systematic analysis. Clinical cohorts comprised 4 to 253 patients aged between 17 and 91 years with confirmed HPV16-positive HNSCC. Samples were collected from fresh biopsies of the tumour, oral rinse or formol fixed/paraffin embedded tissue, from the oral cavity, oropharynx, hypopharynx, larynx and Waldeyer's tonsillar ring. HPV16 variants were identified using Sanger sequencing techniques. Seven studies addressed the HPV16 E6 gene, one studied E6 and E7, another studied L1 and one focused on the long control region. European variants represent 25-95%, Asian-American 5-57% and African 2-4% of the total isolates, suggesting a marked predominance of European strains. No correlations could be drawn with patient prognosis, partly because many studies relied on small patient cohorts. Additional studies are needed, particularly those employing next generation sequencing techniques (NGS), which will allow faster and accurate analysis of large numbers of samples.

Association of cervical carcinogenesis risk with HPV16 E6 and E7 variants in the Taizhou area, China

Background

Human papillomavirus (HPV) type 16 accounts for a larger share of cervical cancer and has been a major health problem worldwide for decades. The progression of initial infection to cervical cancer has been linked to viral sequence properties; however, the role of HPV16 variants in the risk of cervical carcinogenesis, especially with longitudinal follow-up, is not fully understood in China.

Methods

We aimed to investigate the genetic variability of HPV16 E6 and E7 oncogenes in isolates from cervical exfoliated cells. Between December 2012 and December 2014, a total of 310 single HPV16-positive samples were selected from women living in the Taizhou area, China. Sequences of all E6 and E7 oncogenes were analysed by PCR-sequencing assay. Detailed sequence comparison, genetic heterogeneity analyses and maximum-likelihood phylogenetic tree construction were performed with BioEdit Sequence Alignment Editor and MEGA X software. Data for cytology tests and histological diagnoses were obtained from our Taizhou Area Study with longitudinal follow-up for at least 5 years. The relationship between HPV16 variants and cervical carcinogenesis risk was analysed by the chi-square test or Fisher’s exact test.

Results

In this study, we obtained 64 distinct variation patterns with the accession GenBank numbers MT681266-MT681329. Phylogenetic analysis revealed that 98.3% of HPV16 variants belong to lineage A, in which the A4 (Asian) sublineage was dominant (64.8%), followed by A2 (12.1%), A1 (11.4%), and A3 (10.0%). The A4 (Asian) sublineage had a higher risk of CIN2+ than the A1–3 (European) sublineages (OR = 2.69, 95% CI = 1.04–6.97, P < 0.05). Furthermore, nucleotide variation in HPV16 E6 T178G is associated with the development of cervical cancer.

Conclusion

These data could provide novel insights into the role of HPV16 variants in cervical carcinogenesis risk in China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08531-y.

High Levels of Within-Host Variations of Human Papillomavirus 16 E1/E2 Genes in Invasive Cervical Cancer

Human papillomavirus type 16 (HPV16) is the most common HPV genotype found in invasive cervical cancer (ICC). Recent comprehensive genomics studies of HPV16 have revealed that a large number of minor nucleotide variations in the viral genome are present in each infected woman; however, it remains unclear whether such within-host variations of HPV16 are linked to cervical carcinogenesis. Here, by employing next-generation sequencing approaches, we explored the mutational profiles of the HPV16 genome within individual clinical specimens from ICC (n = 31) and normal cervix (n = 21) in greater detail. A total of 367 minor nucleotide variations (167 from ICC and 200 from the normal cervix) were detected throughout the viral genome in both groups, while nucleotide variations at high frequencies (>10% abundance in relative read counts in a single sample) were more prevalent in ICC (10 in ICC versus 1 in normal). Among the high-level variations found in ICC, six were located in the E1/E2 genes, and all of them were non-synonymous substitutions (Q142K, M207I, and L262V for E1; D153Y, R302T, and T357A for E2). In vitro functional analyses of these E1/E2 variants revealed that E1/M207I, E2/D153Y, and E2/R302T had reduced abilities to support viral replication, and that E2/D153Y and E2/R302T failed to suppress the viral early promoter. These results imply that some within-host variations of E1/E2 present at high levels in ICC may be positively selected for and contribute to cervical cancer development through dysfunction or de-stabilization of viral replication/transcription proteins.

Results of phase 2 trials exploring the safety and efficacy of omiganan in patients with human papillomavirus-induced genital lesions

AIMS

To assess safety and tolerability and explore pharmacodynamics and efficacy of omiganan in external anogenital warts (AGW) and vulvar high-grade squamous intraepithelial lesions (HSIL).

METHODS

Two randomized controlled trials in patients with external AGW and vulvar HSIL were conducted. Patients received topical omiganan 2.5% or placebo gel once daily for 12 weeks with a follow-up of 12 weeks. Safety and tolerability were monitored and pharmacodynamics and clinical efficacy of omiganan were assessed by analysing lesion count, size and viral load. Self-reported pain, itch and quality of life were assessed by an electronic diary and questionnaire.

RESULTS

Twenty-four AGW and 12 vulvar HSIL patients were enrolled. All patients had a high treatment adherence (99%). No serious adverse events occurred and all adverse events (n = 27) were mild, transient and self-limiting. The treatment groups were not different in terms of safety and tolerability, lesion count and size, and patient-reported outcomes pain, itch and quality of life. Human papillomavirus load significantly reduced after 12 weeks of treatment with omiganan compared to placebo (-96.6%; 95% confidence interval -99.9 to -7.4%; P = .045) in AGW patients only.

CONCLUSION

Topical omiganan appears to be safe in patients with AGW and vulvar HSIL and reduced human papillomavirus load after 12 weeks of treatment in AGW patients.

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.

Identification by high-throughput sequencing of HPV variants and quasispecies that are untypeable by linear reverse blotting assay in cervical specimens

The linear reverse blotting assays are valid methods for accurate human papillomavirus (HPV) typing required to manage women at risk of developing cervical cancer. However, some samples showed a positive signal in HPV lines but failed to display a positive signal in subsequent typing lines (designated as HPV-X), which indicate that certain types were not available on the respective typing blots. The aim of this study is to elucidate the types or variants of HPV through the high-throughput sequencing (HTS) of 54 ASCUS cervical samples in which the viruses remained untypeable with INNO LiPA HPV® assays. Low-risk (LR)-HPV types (HPV6, 30, 42, 62, 67, 72, 74, 81, 83, 84, 87, 89, 90 and 114), high-risk (HR)-HPV35 and possibly (p)HR-HPV73 were detected among HPV-X. Individual multiple infections (two to seven types) were detected in 40.7% of samples. Twenty-two specimens contained variants characterised by 2-10 changes. HPV30 reached the maximal number of 17 variants with relative abundance inferior or equal to 2.7%. The presence of L1 quasispecies explains why linear reverse blotting assays fail when variants compete or do not match the specific probes. Further studies are needed to measure the LR-HPV quasispecies dynamics and its role during persistent infection.

Genetic variability and functional implication of HPV16 from cervical intraepithelial neoplasia in Shanghai women

Human papillomavirus (HPV)16 gene mutation is usually associated with persistent HPV infection and cervical intraepithelial neoplasia (CIN). However, the functional implications of HPV16 mutations remain poorly understood.145 LCR/E6/E7 of the HPV16 isolates were amplified and sequenced, and HPV16 integration status was detected. In total, 89 SNPs (68 in the LCR, 13 in E6, 8 in E7) were discovered, 11 of which were nonsynonymous mutations (8 in E6, 3 in E7). The H85Y and E120D variants in E6 were significantly reduced in the high-grade squamous intraepithelial lesion (HSIL) group compared to the <HSIL group (P = .046 and .005), conversely the N29S in E7(P = .01). Amino acid substitutions (D32N/E, E36Q, H85Y, and E120D in E6 and N29H/S and R77C in E7) were predicted to have an effect on conserved structural and functional residues, and five amino acid substitutions (H85Y, E36Q, I34L, and D32E in E6; R77C in E7) would potentially change the secondary structure. "6329G>T," a potential binding site for TATA-binding protein, is the most common in LCR variants. A4 (Asian) was associated with an increased risk of HSIL compared to A1-3(P = .009). The H85/E120 in E6 and N29 in HPV16 E7 might play a critical role in carcinogenesis by disrupting p53 and Rb degradation due to affecting their interaction, respectively. In a word, the findings in this study provide preventative and therapeutic interventions of HPV16 -related cervical lesions/cancer.

HPV16 whole genome minority variants in persistent infections from young Dutch women

BACKGROUND

Chronic infections by one of the oncogenic human papillomaviruses (HPVs) are responsible for near 5% of the global cancer burden and HPV16 is the type most often found in cancers. HPV genomes display unexpected levels of variation when deep-sequenced. Minor nucleotide variations (MNVs) may reveal HPV genomic instability and HPV-related carcinogenic transformation of host cells.

OBJECTIVES

The objective of this study was to investigate HPV16 genome variation at the minor variant level on persisting HPV16 cervical infections from a population of young Dutch women.

STUDY DESIGN

15 HPV16 infections were sequenced using a whole-HPV genome deep sequencing protocol (TaME-seq). One infection was followed over a three-year period, eight were followed over a two-year period, three were followed over a one-year period and three infections had a single sampling point.

RESULTS AND CONCLUSIONS

Using a 1% variant frequency cutoff, we find on average 48 MNVs per HPV16 genome and 1717 MNVs in total when sequencing coverage was >100 × . We find the transition mutation T > C to be the most common, in contrast to other studies detecting APOBEC-related C > T mutation profiles in pre-cancerous and cancer samples. Our results suggest that the relative mutagenic footprint of HPV16 genomes may differ between the infections in this study and transforming lesions. In addition, we identify a number of MNVs that have previously been associated with higher incidence of high-grade lesions (CIN3+) in a population study. These findings may provide a starting point for future studies exploring causality between emerging HPV minor genomic variants and cancer development.

Abundance of HPV L1 Intra-Genotype Variants With Capsid Epitopic Modifications Found Within Low- and High-Grade Pap Smears With Potential Implications for Vaccinology

Background: The aim of this study was to explore the Human Papillomavirus (HPV) genotype composition and intra-genotype variants within individual samples of low- and high-grade cervical cytology by deep sequencing. Clinical, cytological, sequencing, and functional/structural data were forged into an integrated variant profiling pipeline for the detection of potentially vaccine-resistant genotypes or variants.

Methods: Low- and high-grade intraepithelial lesion (LSIL and HSIL) cytology samples with +HPV were subjected to amplicon (L1 gene fragment) sequencing by dideoxy (Sanger) and deep methods. Taxonomic, abundance, diversity, and phylogenetic analyses were conducted to determine HPV genotypes/sub-lineages, relative abundance, species diversity and phylogenetic distances within and between samples. Variant detection and functional analysis of translated L1 amino acid sequences determined structural variations of interest.

Results: Pure and mixed HPV infections were common among LSIL (n = 6) and HSIL (n = 6) samples. Taxonomic profiling revealed loss of species richness and gain of dominance by carcinogenic genotypes in HSIL samples. Phylogenetic analysis showed excellent correlation between HPV-type specific genetic distances and carcinogenic potential. For combined LSIL/HSIL samples (n = 12), 11 HPV genotypes and 417 mutations were detected: 375 single-nucleotide variants (SNV), 29 insertion/deletion (indel), 12 multi-nucleotide variants (MNV), and 1 replacement variant. The proportion of nonsynonymous mutations was lower for HSIL (0.38) than for LSIL samples (0.51) (p < 0.05). HPV variant analysis pinpointed nucleotide-level mutations and amino acid-level structural modifications.

Conclusion: HPV L1 intra-host and intra-genotype variants are abundant in LSIL and HSIL samples with potential functional/structural consequences. An integrated multi-omics approach to variant analysis may provide a sensitive and practical means of detecting changes in HPV evolution and dynamics within individuals or populations.

Whole-Genome Analysis of Human Papillomavirus Type 16 Prevalent in Japanese Women with or without Cervical Lesions

Recent large-scale genomics studies of human papillomaviruses (HPVs) have shown a high level of genomic variability of HPV16, the most prevalent genotype in HPV-associated malignancies, and provided new insights into the biological and clinical relevance of its genetic variations in cervical cancer development. Here, we performed deep sequencing analyses of the viral genome to explore genetic variations of HPV16 that are prevalent in Japan. A total of 100 complete genome sequences of HPV16 were determined from cervical specimens collected from Japanese women with cervical intraepithelial neoplasia and invasive cervical cancer, or without cervical malignancies. Phylogenetic analyses revealed the variant distribution in the Japanese HPV16 isolates; overall, lineage A was the most prevalent (94.0%), in which sublineage A4 was dominant (52.0%), followed by sublineage A1 (21.0%). The relative risk of sublineage A4 for cervical cancer development was significantly higher compared to sublineages A1/A2/A3 (odds ratio = 6.72, 95% confidence interval = 1.78–28.9). Interestingly, a novel cluster of variants that branched from A1/A2/A3 was observed for the Japanese HPV16 isolates, indicating that unique HPV16 variants are prevalent among Japanese women.

HPV16 variant analysis in primary and recurrent CIN2/3 lesions demonstrates presence of the same consensus variant

INTRODUCTION

Recurrent cervical intraepithelial lesions (rCIN2/3) after treatment of CIN2/3 occur in 5-15% of cases. rCIN2/3 can result from incomplete resection of CIN2/3, where the same HPV type and variant remains present. rCIN2/3 could also occur following a new infection with a different HPV variant of the same HPV type as the initial lesion. This study investigates HPV16 consensus variants in paired HPV16 positive scrapes from baseline CIN2/3 and rCIN2/3 lesions.

METHODS

Paired HPV16 positive cervical scrapes of women with CIN2/3 at baseline and rCIN2/3 6 or 12 months after treatment were selected for whole-genome amplification and Illumina sequencing. Sequences were compared and nucleotide changes over time were characterized.

RESULTS

From 14 paired samples, 10 had identical consensus variants in baseline CIN2/3 and rCIN2/3. Four paired samples showed one to three nucleotide variations at recurrent disease compared to baseline.

CONCLUSION

Identical or nearly identical HPV16 consensus variants were found in scrapes of paired HPV16 positive baseline CIN2/3 and rCIN2/3 lesions after treatment, suggesting no need for HPV variant analysis when the same HPV type is found in both lesions. These results argue for either incomplete excision of baseline CIN2/3 or inability of clearance of the original HPV infection.

Disease Manifestation and Viral Sequences in a Bonobo More Than 30 Years after Papillomavirus Infection

Pan paniscus Papillomavirus 1 (PpPV1) causes focal epithelial hyperplasia (FEH) in infected animals. Here, we analyzed the present disease manifestation and PpPV1 genomic sequence of an animal that was afflicted by an FEH epizootic outbreak in 1987 for which the sequence of the responsible PpPV1 was determined. The animal displayed FEH more than 30 years after the initial diagnosis, indicating persistence or recurrence of the disease, and evidence for active PpPV1 infection was obtained. Moreover, the sequences of the viral genomes present in the late 1980s and in 2018 differed at 23 nucleotide positions, resulting in 11 amino acid exchanges within coding regions. These findings suggest that PpPV1-induced FEH might not undergo complete and/or permanent remission in a subset of afflicted animals.

Identification of human papillomavirus type 16 variants circulating in the Calabria region by sequencing and phylogenetic analysis of HPV16 from cervical smears

Sequence analysis of HPV16 isolates reveals the presence of genome variants with characteristic mutations. The HPV16 variants have different geographical distribution and diverge into four phylogenetic lineages (A, B, C and D) and 16 sub-lineages: A1, A2, A3 (previously known as European variants), A4 (Asian variant), B1, B2, B3, B4, C1, C2, C3, and C4 (African variants), D1 (North-American variant), D2, D3 (Asian-American variants) and D4. Population studies showed that infections with viruses belonging to specific HPV16 sublineages confer different risks of viral persistence and cancer. In this study, 39 HPV16-positive cervical smears from European women living in Calabria (Italy) were analyzed for the presence of HPV16 variants. Cervical DNA extracts were processed by PCR to amplify L1, the Long Control Region (LCR), E6 and E7, which were sequenced. The sequences were concatenated and the 3169 nucleotides long fragments were characterized by BLAST and phylogenetic analysis. A total of 96 Single Nucleotide Polymorphism (SNPs) were detected, 29 of which mapping in the L1, 45 in the LCR, 15 in the E6 and 7 in the E7. The most common SNP was the T350G (29/39 samples, 74.4%), causing the L83 V amino acid change in the E6. Most of the HPV16 isolates (89.7%) had 99% of nucleotide (nt) identity to members of the A1 and A2 sublineages, while 4 isolates had 99% nt identity to members of the B2, B4, C1 and D4 sublineages. In conclusion, viruses belonging to the A1, A2, B2, B4, C1 and D4 HPV16 sublineages were found to circulate in the Calabria region.

Prevalence and concordance of oral and genital HPV in women positive for cervical HPV infection and in their sexual stable partners: An Italian screening study

OBJECTIVES

This cross-sectional study aimed to evaluate the prevalence and type of oral HPV-infection in women with a cervical HPV-lesion and in the oral and genital mucosa of their male partners.

METHODS

The study group comprised 44 sexually-active women, 20-45 years with abnormal PAP smear, not more than 6 months prior to referral together with the male partners cohabiting in stable partnerships. A detailed questionnaire was administered concerning the HPV-related risk factors. Oral swabs, oral rinses, cervical swabs and urine samples were collected. HPV DNA was detected using two different polymerase chain reactions (PCRs): MY09-11 and FAP59-64. Positive samples were genotyped by Sanger sequencing and the INNO-LiPA HPV Genotyping Extra II probe assay. The association with risk factors was assessed by fitting a generalized model, using the General Linear Model function in the R-software; correlations were calculated between all data.

RESULTS

HPV was detected in 84% of Cervical Samples, in 24.3% of oral samples and in one urine sample. Only 27% of the HPV-positive results were identical with both PCR DNA assays. 8 male had oral HPV-positive samples different from women cervical samples. In one couple the urine-male sample had the same HPV present in the female-cervical sample. A significant association resulted between women/oral sex practices and men/n. of partners.

CONCLUSIONS

This study reports that women (20.4%) with a diagnosis of cervical-HPV and their male partners (30,7%) are at high risk for subclinical oral HPV infection.

Within-Host Variations of Human Papillomavirus Reveal APOBEC Signature Mutagenesis in the Viral Genome

Persistent infection with oncogenic human papillomaviruses (HPVs) causes cervical cancer, accompanied by the accumulation of somatic mutations into the host genome. There are concomitant genetic changes in the HPV genome during viral infection; however, their relevance to cervical carcinogenesis is poorly understood. Here, we explored within-host genetic diversity of HPV by performing deep-sequencing analyses of viral whole-genome sequences in clinical specimens. The whole genomes of HPV types 16, 52, and 58 were amplified by type-specific PCR from total cellular DNA of cervical exfoliated cells collected from patients with cervical intraepithelial neoplasia (CIN) and invasive cervical cancer (ICC) and were deep sequenced. After constructing a reference viral genome sequence for each specimen, nucleotide positions showing changes with >0.5% frequencies compared to the reference sequence were determined for individual samples. In total, 1,052 positions of nucleotide variations were detected in HPV genomes from 151 samples (CIN1, n = 56; CIN2/3, n = 68; ICC, n = 27), with various numbers per sample. Overall, C-to-T and C-to-A substitutions were the dominant changes observed across all histological grades. While C-to-T transitions were predominantly detected in CIN1, their prevalence was decreased in CIN2/3 and fell below that of C-to-A transversions in ICC. Analysis of the trinucleotide context encompassing substituted bases revealed that TpCpN, a preferred target sequence for cellular APOBEC cytosine deaminases, was a primary site for C-to-T substitutions in the HPV genome. These results strongly imply that the APOBEC proteins are drivers of HPV genome mutation, particularly in CIN1 lesions.IMPORTANCE HPVs exhibit surprisingly high levels of genetic diversity, including a large repertoire of minor genomic variants in each viral genotype. Here, by conducting deep-sequencing analyses, we show for the first time a comprehensive snapshot of the within-host genetic diversity of high-risk HPVs during cervical carcinogenesis. Quasispecies harboring minor nucleotide variations in viral whole-genome sequences were extensively observed across different grades of CIN and cervical cancer. Among the within-host variations, C-to-T transitions, a characteristic change mediated by cellular APOBEC cytosine deaminases, were predominantly detected throughout the whole viral genome, most strikingly in low-grade CIN lesions. The results strongly suggest that within-host variations of the HPV genome are primarily generated through the interaction with host cell DNA-editing enzymes and that such within-host variability is an evolutionary source of the genetic diversity of HPVs.

Human papillomavirus genome variants and head and neck cancers: a perspective

Human papillomaviruses (HPV) cause infections that are responsible for diverse clinical manifestations from benign conditions to invasive cancer. As different HPV types are associated with variable pathogenic potential, minor genetic variations within a given high-risk HPV type might also be associated with distinct oncogenic capacities, through variable ability of persistence or risk of progression to precancer/cancer. Most recent HPV variant studies in the cervix using latest sequencing technology confirmed that minor changes in the HPV genome can have a major influence on carcinogenesis and have revealed key data that help better understand the carcinogenicity of HPV at a molecular level. Here we review the limited number of studies on HPV genome variants in head and neck cancers (HNC) and discuss their implications for cancer research in the light of accumulated knowledge for the cervix. Challenges in transposing HPV variant studies from the lower anogenital to the upper aerodigestive tract are also discussed, highlighting the main gaps of knowledge in the field of HPV-induced HNC. Specifically in the head and neck region, the lack of characterisation of precancerous lesions and the difficulty in sampling normal tissue will challenge the development of accurate studies. Although there is so far no indication that HPV variant research in HNC could directly translate into clinical application, such research is expected to be useful to disentangle unanswered questions in the pathogenesis of HNC. Yet, history of HPV variant research suggests that, to be successful, studies will require large international collaborative efforts.

High Whole-Genome Sequence Diversity of Human Papillomavirus Type 18 Isolates

Background: The most commonly found human papillomavirus (HPV) types in cervical cancer are HPV16 and HPV18. Genome variants of these types have been associated with differential carcinogenic potential. To date, only a handful of studies have described HPV18 whole genome sequencing results. Here we describe HPV18 variant diversity and conservation of persistent infections in a longitudinal retrospective cohort study. Methods: Cervical self-samples were obtained annually over four years and genotyped on the SPF10-DEIA-LiPA₂₅ platform. Clearing and persistent HPV18 positive infections were selected, amplified in two overlapping fragments, and sequenced using 32 sequence primers. Results: Complete viral genomes were obtained from 25 participants with persistent and 26 participants with clearing HPV18 infections, resulting in 52 unique HPV18 genomes. Sublineage A3 was predominant in this population. The consensus viral genome was completely conserved over time in persistent infections, with one exception, where different HPV18 variants were identified in follow-up samples. Conclusions: This study identified a diverse set of HPV18 variants. In persistent infections, the consensus viral genome is conserved. The identification of only one HPV18 infection with different major variants in follow-up implies that this is a potentially rare event. This dataset adds 52 HPV18 genome variants to Genbank, more than doubling the currently available HPV18 information resource, and all but one variant are unique additions.

Why Human Papillomavirus Acute Infections Matter

Most infections by human papillomaviruses (HPVs) are `acute', that is non-persistent. Yet, for HPVs, as for many other oncoviruses, there is a striking gap between our detailed understanding of chronic infections and our limited data on the early stages of infection. Here we argue that studying HPV acute infections is necessary and timely. Focusing on early interactions will help explain why certain infections are cleared while others become chronic or latent. From a molecular perspective, descriptions of immune effectors and pro-inflammatory pathways during the initial stages of infections have the potential to lead to novel treatments or to improved handling algorithms. From a dynamical perspective, adopting concepts from spatial ecology, such as meta-populations or meta-communities, can help explain why HPV acute infections sometimes last for years. Furthermore, cervical cancer screening and vaccines impose novel iatrogenic pressures on HPVs, implying that anticipating any viral evolutionary response remains essential. Finally, hints at the associations between HPV acute infections and fertility deserve further investigation given their high, worldwide prevalence. Overall, understanding asymptomatic and benign infections may be instrumental in reducing HPV virulence.