HPV16 Sublineage Associations With Histology-Specific Cancer Risk Using HPV Whole-Genome Sequences in 3200 Women

DNA methylation and copy number alterations in the progression of HPV-associated high-grade vulvar intraepithelial lesion

Human papillomavirus (HPV)-associated high-grade vulvar intraepithelial lesion (HSIL) is a precursor of vulvar squamous cell carcinoma (VSCC). Because of the 8% cancer risk, many vulvar HSIL patients undergo aggressive and mutilating treatments. Characterizing HSIL by their progression risk can help individualize treatment strategies. Accordingly, copy number alterations (CNAs) and DNA methylation have been identified as biomarkers for cancer risk stratification of HSIL. Here, we assessed their potential correlation, and relation to HPV16 (sub)lineages and progression to vulvar cancer. Eighty-two vulvar formalin-fixed paraffin-embedded (FFPE) samples, including controls, HSIL, HSIL adjacent to VSCC and VSCC, with previously determined DNA methylation profiles, were analysed for CNAs using mFAST-SeqS. Genome-wide z-scores were calculated to determine overall aneuploidy (aneuploidy scores), and compared to the methylation levels and status of marker panel ZNF582/SST/miR124-2. For 52 HPV16-positive cases, HPV (sub)lineages were determined by Sanger sequencing. HPV16 lineage A was predominant (86.4%), followed equally by lineages B, C, and D. Frequent chromosomal alterations included chr1pq, chr3q, chr9q gains, and chr2q, chr4q losses. Median aneuploidy scores increased across disease categories, from 0 in controls, to 3 in HSIL, 16 in HSIL adjacent to VSCC and 29 in VSCC. A positive relationship between aneuploidy scores and DNA methylation levels was found (ρ = 0.61, Spearman's rank correlation test). Aneuploidy scores were significantly higher in methylation-positive samples (p < .001). In conclusion, we showed that DNA methylation and CNAs both rise with increasing severity of disease, indicating their prognostic value for cancer risk stratification of HSIL, while no relation to HPV16 (sub)lineages was found.

The Lineage and Sublineage Investigation of Human Papillomavirus Type 16 in Tehran, Iran, During 2022-2023: A Cross-Sectional Study

Background and Aims

Diverse HPV 16 variants have been known to vary in geographic distribution and oncogenic potential. On this point, the nucleotide changes of the E6 gene were studied to find lineages of HPV 16 in normal, premalignant, and malignant stages of uterine cervical samples.

Methods

In this study, 120 HPV 16-infected samples were investigated using PCR and sequencing.

Results

In our samples, three lineages A, C, and D were found and lineage D was predominant (79.2%) followed by lineages A (20%) and C (0.8%). Concerning the association between histopathological stages and lineages, no statistically important differences were found.

Conclusion

Our finding revealed that two lineages A and D were circulating, with the dominancy of lineage, in Tehran, Iran. This finding emphasizes the geographic diversity of distinct HPV 16 (sub)lineages is different in the world.

Exploring HPV-linked head and neck cancer in Southern Punjab, Pakistan: Insights from HPV-16 phylogenetic analysis

INTRODUCTION/OBJECTIVES

Head and neck region squamous cell carcinoma (HNSCC) is a heterogeneous disease that can be categorized into human papillomavirus (HPV)-positive (20 %) and HPV-negative (80 %) subtypes. However, the prevalence of HPV genotypes is not clear in Pakistan. This study investigated how common the HPV-16 genotype is in patients with HNSCC in the Southern Punjab region of Pakistan, and the specific molecular features of this genotype.

METHODS

For this cross-sectional study, 85 tissue samples were collected from diagnosed cases of HNSCC. Formalin-fixed paraffin-embedded tissue sections were used for genomic DNA extraction. The L1 region was amplified using GP5+ and GP6+ primers to detect HPV DNA. Real-time PCR was conducted to genotype high-risk HPV (HR-HPV). Whole genome sequencing was used for phylogenetic analysis of HPV-16 and to detect mutations/single nucleotide polymorphisms (SNPs).

RESULTS

Among the 85 samples, 7.1 % were positive for HPV, where 4.7 % were positive for HPV-16 and 2.4 % were positive for HPV-18. A significant association was found between HR-HPV positivity and histological grade (p < 0.05). The HPV-16 genome sequence obtained in this study was closely related to those from Thailand, the United States, India, China, and Europe, and 11 mutations/SNPs were detected in the sequenced genome, where four were novel.

CONCLUSION

The findings obtained in the present study demonstrate the low prevalence of HR-HPV associated HNSCC in Pakistan. Phylogenetic analysis showed that HPV-16 genome isolated and sequenced in this study had a distinct genetic structure and it also shared similarities with genomes reported from Thailand, the United States, India, China, and Europe.

Multiple DNA Viruses and HPV Integration in Inverted Papilloma and Associated Sinonasal Carcinoma

OBJECTIVES

Sinonasal inverted papilloma (IP) has a locally destructive growth pattern, can relapse, and can undergo malignant transformation (IP-associated sinonasal squamous cell carcinoma (IP-SNSCC)). Human papillomaviruses (HPV)-6 and -16 are frequently detected in IPs. To clarify the possible roles of other DNA viruses in IPs, we explored viruses not studied in this context before. With the setting of pre- and post-malignant transformation samples, we investigated HPV genomes in depth to assess the integration of HPV into the human genome and the presence of minor intratypic variants.

MATERIALS AND METHODS

We analyzed 35 IP samples representing 28 individuals, of which six had IP-SNSCC. For virus screening, we applied qPCR to detect 16 different DNA viruses in three virus families, comprising herpesviruses, parvoviruses, and polyomaviruses. In addition, targeted next generation sequencing (NGS) was used for detailed HPV analysis.

RESULTS

We detected herpes-, parvo-, and polyomaviruses in 13/28 (46%) patients, with codetections of multiple viruses in six (21%) patients. NGS revealed HPV16 DNA in 2/6 IP-SNSCC and in their respective earlier benign IP samples, as well as in a plasma sample from one of these patients. HPV6 was detected in two IP samples without subsequent malignant transformation. We identified sequence reads containing junctions of HPV6 and HPV16 and host genome suggestive of viral integration. HPV6 and HPV16 minor intratypic variants were present across pre- and post-malignant transformation, with mostly nonsynonymous mutations.

CONCLUSIONS

Multiple DNA viruses were present in IPs. HPV16 was detected only in IP-SNSCCs or in tumors that later underwent malignant transformation.

LEVEL OF EVIDENCE

3 Laryngoscope, 135:677-686, 2025.

Analysis of the progression of cervical cancer in a low-and-middle-income country: From pre-malignancy to invasive disease

To better understand cervical cancer progression, we analyzed RNA from 262 biopsies from women referred for colposcopy. We determined the HPV type and analyzed the expression of 51 genes. HPV31 was significantly more prevalent in precancer than stage 1 cancer and invasive cancer (p < 0.0001), and HPV16 increased in invasive disease (p < 0.0001). CCNE1, MELTF, and ULBP2 were significantly increased in HPV16-positive compared to HPV31 precancers, while NECTIN2 and HLA-E expression decreased. Markers of the innate immune system, DNA repair genes, and cell cycle genes are significantly increased during cancer progression (p = 0.0001). In contrast, the TP53 and RB1 tumor suppressor gene expression is significantly decreased in cancer cells. The T cell markers CD28 and FLT3LG expression decreased in cancer while FOXP3, IDO1, and ULBP2 expression increased. There is a significantly higher survival rate in individuals with increased expression of CD28 (p = 0.0005), FOXP3 (p = 0.0002), IDO1 (p = 0.038), FLT3LG (p = 0.026), APOBEC3B (p = 0.0011), and RUNX3 (p = 0.019), and a significantly lower survival rate in individuals with increased expression of ULBP2 (p = 0.035). These results will help us elucidate the molecular factors influencing the progression of cervical precancer to cancer. Understanding the risk of progression of specific HPV types and sublineages may aid in the triage of positive patients, and better knowledge of the immune response may aid in developing and applying immunotherapies.

Molecular epidemiology of human papillomavirus variants in cervical cancer in India

Background & objectives Cervical cancer (CC) has been documented as the fourth most common cancer worldwide. Persistent infections with high-risk human papillomavirus (hr-HPV) have been suggested in the development of CC. Although prophylactic vaccines are available for the prevention of prevalent hr-HPV types, intra-type variations exist within a particular HPV type that has varying oncogenic potential as well as the mechanism of pathogenicity and varying neutralization by antibodies. Therefore, we carried out a systematic review to determine the distribution of HPV intra-typic variations in different geographical locations of India and their reported implications. Methods Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) guidelines were followed to retrieve relevant articles from the standard databases using appropriate keywords. Consequently, 17 articles were included in the current review after screening based on inclusion and exclusion criteria. Results The majority of articles included in this review reported variations within the HPV16 E6 gene, followed by the L1 and E7 genes. Analysis of available data indicated the differential regional distribution of some variations. These variations have also been reported to impact the biological functions of various viral proteins. Interpretation & conclusions The distribution of lineages varied with the different genomic regions sequenced. Additionally, there were certain unique and common variations in the HPV genome with respect to geographical regions. Hence, we suggest the identification of region-specific variations for the development of diagnostic and prognostic interventions.

Why HPV16? Why, now, HPV42? How the discovery of HPV42 in rare cancers provides an opportunity to challenge our understanding about the transition between health and disease for common members of the healthy microbiota

In 2022, a bioinformatic, agnostic approach identified HPV42 as causative agent of a rare cancer, later confirmed experimentally. This unexpected association offers an opportunity to reconsider our understanding about papillomavirus infections and cancers. We have expanded our knowledge about the diversity of papillomaviruses and the diseases they cause. Yet, we still lack answers to fundamental questions, such as what makes HPV16 different from the closely related HPV31 or HPV33; or why the very divergent HPV13 and HPV32 cause focal epithelial hyperplasia, while HPV6 or HPV42 do not, despite their evolutionary relatedness. Certain members of the healthy skin microbiota are associated to rare clinical conditions. We propose that a focus on cellular phenotypes, most often transient and influenced by intrinsic and extrinsic factors, may help understand the continuum between health and disease. A conceptual switch is required towards an interpretation of biology as a diversity of states connected by transition probabilities, rather than quasi-deterministic programs. Under this perspective, papillomaviruses may only trigger malignant transformation when specific viral genotypes interact with precise cellular states. Drawing on Canguilhem's concepts of normal and pathological, we suggest that understanding the transition between fluid cellular states can illuminate how commensal-like infections transition from benign to malignant.

Human Papillomavirus (HPV) Vaccination: Progress, Challenges, and Future Directions in Global Immunization Strategies

Human papillomavirus (HPV) is a widespread viral pathogen, responsible for a significant burden of cervical and other cancers worldwide. Over the past decades, the development and widespread adoption of prophylactic HPV vaccines have dramatically reduced the incidence of HPV-related diseases. However, despite the efficacy of these vaccines, global immunization efforts still face several obstacles, including low vaccination coverage in low- and middle-income countries, vaccine hesitancy, and disparities in access to healthcare. This review aims to provide a comprehensive overview of the current state of HPV vaccines, including their mechanisms of action, safety profiles, and real-world efficacy. We will explore the impact of HPV vaccines on cancer prevention, examine the challenges related to vaccine distribution and uptake, and assess the role of public health policies in improving global vaccination rates. Additionally, the review will highlight the latest advancements in therapeutic HPV vaccines, ongoing research into next-generation vaccines, and the potential of HPV vaccination strategies in the context of personalized medicine. By examining these factors, we aim to provide insights into the future directions of HPV vaccination and its role in global public health.

Meeting report: Considerations for trial design and endpoints in licensing therapeutic HPV16/18 vaccines to prevent cervical cancer

Cervical cancer is a major cause of morbidity and mortality globally with a disproportionate impact on women in low- and middle-income countries. In 2021, the World Health Organization (WHO) called for increased vaccination, screening, and treatment to eliminate cervical cancer. However, even with widespread rollout of human papillomavirus (HPV) prophylactic vaccines, millions of women who previously acquired HPV infections will remain at risk for progression to cancer for decades to come. The development and licensing of an affordable, accessible therapeutic HPV vaccine, designed to clear or control carcinogenic HPV and/or to induce regression precancer could significantly contribute to the elimination efforts, particularly benefiting those who missed out on the prophylactic vaccine. One barrier to development of such vaccines is clarity around the regulatory pathway for licensure. In Washington, D.C. on September 12-13, 2023, a meeting was convened to provide input and guidance on trial design with associated ethical and regulatory considerations. This report summarizes the discussion and conclusions from the meeting. Expert presentation topics included the current state of research, potential regulatory challenges, WHO preferred product characteristics, modeling results of impact of vaccine implementation, epidemiology and natural history of HPV infection, immune responses related to viral clearance and/or precancer regression including potential biomarkers, and ethical considerations. Panel discussions were held to explore specific trial design recommendations to support the licensure process for two vaccine indications: (1) treatment of prevalent HPV infection or (2) treatment of cervical precancers. Discussion covered inclusion/exclusion criteria, study endpoints, sample size and power, safety, study length, and additional data needed, which are reported here. Further research of HPV natural history is needed to address identified gaps in regulatory guidance, especially for therapeutic vaccines intended to treat existing HPV infections.

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.

Comparative Analysis of HPV16 Variants in the Untranslated Regulatory Region, L1, and E6 Genes among Vaccinated and Unvaccinated Young Women: Assessing Vaccine Efficacy and Viral Diversity

HPV16 is occasionally detected in vaccinated women who received the bivalent HPV16/18 vaccine, usually at low viral loads. This study explored potential differences in HPV16 variants between vaccinated and unvaccinated women. HPV16-postive viral loads were detected in 1.9% (17/875) and 13% (162/760) of vaccinated and unvaccinated women, respectively, showcasing the vaccine's high efficacy. The L1, E6, and URR regions of HPV16 were sequenced from genital swabs from 16 vaccinated and 25 unvaccinated women in the HAVANA (HPV Among Vaccinated And Non-vaccinated Adolescents) study. The majority of HPV16 variants from vaccinated and unvaccinated women clustered similarly with sub-lineages A1 and A2. Additionally, a separate cluster within lineage A was found, with the variants sharing the L1-located SNP A753G (synonymous) and the URR-located SNP T340C, which did not occur in the other variants. Furthermore, four variants from vaccinated women had relatively long branches, but were not characterized by specific SNPs. The frequency of G712A in the URR was the only SNP observed to be marginally higher among vaccinated women than unvaccinated women. Non-synonymous SNPs T266A in the FG-loop of L1 and L83V in E6 were common among variants from vaccinated and unvaccinated women, but present in similar frequencies. In conclusion, the detection of HPV16 in vaccinated (and unvaccinated) women seemed to be the result of random circulation within this study population.

Identification of HPV16 Lineages in South African and Mozambican Women with Normal and Abnormal Cervical Cytology

BACKGROUND

Human papillomavirus 16 (HPV16) is an oncogenic virus responsible for the majority of invasive cervical cancer cases worldwide. Due to genetic modifications, some variants are more oncogenic than others. We analysed the HPV16 phylogeny in HPV16-positive cervical Desoxyribonucleic Acid (DNA) samples collected from South African and Mozambican women to detect the circulating lineages.

METHODS

Polymerase chain reaction (PCR) amplification of the long control region (LCR) and 300 nucleotides of the E6 region was performed using HPV16-specific primers on HPV16-positive cervical samples collected in women from South Africa and Mozambique. HPV16 sequences were obtained through Next Generation Sequencing (NGS) methods. Geneious prime and MEGA 11 software were used to align the sequences to 16 HPV16 reference sequences, gathering the A, B, C, and D lineages and generating the phylogenetic tree. Single nucleotide polymorphisms (SNPs) in the LCR and E6 regions were analysed and the phylogenetic tree was generated using Geneious Prime software.

RESULTS

Fifty-eight sequences were analysed. Of these sequences, 79% (46/58) were from women who had abnormal cervical cytology. Fifteen SNPs in the LCR and eight in the E6 region were found to be the most common in all sequences. The phylogenetic analysis determined that 45% of the isolates belonged to the A1 sublineage (European variant), 34% belonged to the C1 sublineage (African 1 variant), 16% belonged to the B1 and B2 sublineage (African 2 variant), two isolates belonged to the D1-3 sublineages (Asian-American variant), and one to the North American variant.

CONCLUSIONS

The African and European HPV16 variants were the most common circulating lineages in South African and Mozambican women. A high-grade squamous intraepithelial lesion (HSIL) was the most common cervical abnormality observed and linked to European and African lineages. These findings may contribute to understanding molecular HPV16 epidemiology in South Africa and Mozambique.

HPV16 Phylogenetic Variants in Anogenital and Head and Neck Cancers: State of the Art and Perspectives

HPV16 is responsible for approximately 60% and 90% of global HPV-induced cervical and oropharyngeal cancers, respectively. HPV16 intratype variants have been identified by HPV genome sequencing and classified into four phylogenetic lineages (A-D). Our understanding of HPV16 variants mostly derives from epidemiological studies on cervical cancer (CC) in which HPV16 B, C, and D lineages (previously named "non-European" variants) were mainly associated with high-grade cervical lesions and cancer. Although a predominance of HPV16 lineage A (previously named "European variants") has been observed in head and neck squamous cell carcinoma (HNSCC), epidemiological and in vitro biological studies are still limited for this tumor site. Next Generation Sequencing (NGS) of the entire HPV genome has deepened our knowledge of the prevalence and distribution of HPV variants in CC and HNSCC. Research on cervical cancer has shown that certain HPV16 sublineages, such as D2, D3, A3, and A4, are associated with an increased risk of cervical cancer, and sublineages A4, D2, and D3 are linked to a higher risk of developing adenocarcinomas. Additionally, lineage C and sublineages D2 or D3 of HPV16 show an elevated risk of developing premalignant cervical lesions. However, it is still crucial to conduct large-scale studies on HPV16 variants in different HPV-related tumor sites to deeply evaluate their association with disease development and outcomes. This review discusses the current knowledge and updates on HPV16 phylogenetic variants distribution in HPV-driven anogenital and head and neck cancers.

Distribution of human papillomavirus among Vietnamese women with cervical cancer and unusual genetic variability of HPV16

HPV16, with typical mutations that differ in geographical distribution and carcinogenic potency, has implications for cervical cancer screening, clinical diagnosis, and treatment. DNASTAR and MEGA were used to identify HPV16 variants and construct a phylogenetic tree. The most prevalent HPV genotypes were HPV16 (63.9%), HPV18 (26.7%), and other HPV (6.9%). HPV16 alterations were found in all E6, E7, and L1 genes, including 15 missense and 18 synonymous mutations. Missense mutations include R10G, Q14H, D25E, H78Y, L83V (E6); M29V, R35K, L78R, L95P (E7); H73Y, T176 N, N178T, T317P, T386S, L472F/I (L1). HPV16 sublineages include A1 (17.2%), A2 (0.9%), A3 (56.0%), A4 (19.0%), D1 (4.3%), and D3 (2.6%). Although several mutations in the oncoproteins E6, E7, and L1 have been detected, mutations known to be associated with cervical cancer risk, such as D25E and L83V, occur at a relatively low frequency. This suggests that HPV16 mutations are associated with cervical cancer through a complicated mechanism.

High-Risk Genotypes of Human Papillomavirus at Diverse Anogenital Sites among Chinese Women: Infection Features and Potential Correlation with Cervical Intraepithelial Neoplasia

BACKGROUND

Both cervical cancer and cervical intraepithelial neoplasia (CIN) are associated with human papillomavirus (HPV) infection at different anogenital sites, but the infection features of high-risk (HR) HPVs at these sites and their association with cervical lesions have not been well characterized. Given the limitation of cervical HPV 16/18 test in screening patients with high-grade CIN (CIN 2+), studies on whether non-16/18 HR-HPV subtype(s) have potential as additional indicator(s) to improve CIN 2+ screening are needed.

METHODS

The infection of 15 HR-HPVs in vulva, anus, vagina, and cervix of 499 Chinese women was analyzed, and CIN lesion-associated HR-HPV subtypes were revealed.

RESULTS

In addition to the well-known cervical-cancer-associated HPV 16, 52, and 58, HPV 51, 53, and 56 were also identified as high-frequency detected subtypes prevalently and consistently present at the anogenital sites studied, preferentially in multi-infection patterns. HPV 16, 52, 58, 56, and 53 were the top five prevalent subtypes in patients with CIN 2+. In addition, we found that cervical HPV 33/35/52/53/56/58 co-testing with HPV 16/18 might improve CIN 2+ screening performance.

CONCLUSION

This study provided a new insight into HR-HPV screening strategy based on different subtype combinations, which might be used in risk stratification clinically.

Changes in intrahost genetic diversity according to lesion severity in longitudinal HPV16 samples

Human papillomavirus type 16 (HPV16) is the most common cause of cervical cancer, but most infections are transient with lesions not progressing to cancer. There is a lack of specific biomarkers for early cancer risk stratification. This study aimed to explore the intrahost HPV16 genomic variation in longitudinal samples from HPV16-infected women with different cervical lesion severity (normal, low-grade, and high-grade). The TaME-seq deep sequencing protocol was used to generate whole genome HPV16 sequences of 102 samples collected over time from 40 individuals. Single nucleotide variants (SNVs) and intrahost SNVs (iSNVs) were identified in the viral genomes. A majority of individuals had a unique set of SNVs and these SNVs were stable over time. Overall, the number of iSNVs and APOBEC3-induced iSNVs were significantly lower in high-grade relative to normal and low-grade samples. A significant increase in the number of APOBEC3-induced iSNVs over time was observed for normal samples when compared to high-grade. Our results indicates that the lower incidence of iSNVs and APOBEC3-induced iSNVs in high-grade lesions may have implications for novel biomarkers discoveries, potentially aiding early stratification of HPV-induced cervical precancerous lesions.

Immunoinformatics Design and In Vivo Immunogenicity Evaluation of a Conserved CTL Multi-Epitope Vaccine Targeting HPV16 E5, E6, and E7 Proteins

Human papillomavirus type 16 (HPV16) infection is responsible for more than 50% of global cervical cancer cases. The development of a vaccine based on cytotoxic T-lymphocyte (CTL) epitopes is a promising strategy for eliminating pre-existing HPV infections and treating patients with cervical cancer. In this study, an immunoinformatics approach was used to predict HLA-I-restricted CTL epitopes in HPV16 E5, E6, and E7 proteins, and a set of conserved CTL epitopes co-restricted by human/murine MHCs was screened and characterized, with the set containing three E5, four E6, and four E7 epitopes. Subsequently, the immunogenicity of the epitope combination was assessed in mice, and the anti-tumor effects of the multi-epitope peptide vaccine E5E6E7pep11 and the recombinant protein vaccine CTB-Epi11E567 were evaluated in the TC-1 mouse tumor model. The results demonstrated that mixed epitope peptides could induce antigen-specific IFN-γ secretion in mice. Prophylactic immunization with E5E6E7pep11 and CTB-Epi11E567 was found to provide 100% protection against tumor growth in mice. Moreover, both types of the multi-epitope vaccine significantly inhibited tumor growth and prolonged mouse survival. In conclusion, in this study, a multi-epitope vaccine targeting HPV16 E5, E6, and E7 proteins was successfully designed and evaluated, demonstrating potential immunogenicity and anti-tumor effects and providing a promising strategy for immunotherapy against HPV-associated tumors.

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.

Global evaluation of lineage-specific human papillomavirus capsid antigenicity using antibodies elicited by natural infection

Human Papillomavirus (HPV) type variants have been classified into lineages and sublineages based upon their whole genome sequence. Here we have examined the specificity of antibodies generated following natural infection with lineage variants of oncogenic types (HPV16, 18, 31, 33, 45, 52 and 58) by testing serum samples assembled from existing archives from women residing in Africa, The Americas, Asia or Europe against representative lineage-specific pseudoviruses for each genotype. We have subjected the resulting neutralizing antibody data to antigenic clustering methods and created relational antigenic profiles for each genotype to inform the delineation of lineage-specific serotypes. For most genotypes, there was evidence of differential recognition of lineage-specific antigens and in some cases of a sufficient magnitude to suggest that some lineages should be considered antigenically distinct within their respective genotypes. These data provide compelling evidence for a degree of lineage specificity within the humoral immune response following natural infection with oncogenic HPV.

Genetic diversity of HPV35 in Chad and the Central African Republic, two landlocked countries of Central Africa: A cross-sectional study

Human Papillomavirus (HPV)-35 accounts for up 10% of cervical cancers in Sub-Saharan Africa. We herein assessed the genetic diversity of HPV35 in HIV-negative women from Chad (identified as #CHAD) and HIV-infected men having sex with men (MSM) in the Central African Republic (CAR), identified as #CAR. Ten HPV35 DNA from self-collected genital secretions (n = 5) and anal margin samples (n = 5) obtained from women and MSM, respectively, were sequenced using the ABI PRISM® BigDye Sequencing technology. All but one HPV35 strains belonged to the A2 sublineage, and only #CAR5 belonged to A1. HPV35 from #CAR had higher L1 variability compared to #CHAD (mean number of mutations: 16 versus 6). L1 of #CAR5 showed a significant variability (2.29%), suggesting a possible intra-type divergence from HPV35H. Three (BC, DE, and EF) out of the 5 capsid loops domains remained totally conserved, while FG- and HI- loops of #CAR exhibited amino acid variations. #CAR5 also showed the highest LCR variability with a 16bp insertion at binding sites of the YY1. HPV35 from #CHAD exhibited the highest variability in E2 gene (P<0.05). E6 and E7 oncoproteins remained well conserved. There is a relative maintenance of a well conserved HPV35 A2 sublineage within heterosexual women in Chad and MSM with HIV in the Central African Republic.

Human Papillomavirus 16 Lineage A Variants Associated With Persistent Genital Infections in Men: The HPV Infection in Men (HIM) Study

BACKGROUND

Human papillomavirus (HPV) 16 non-A lineage variants have higher carcinogenic potential for cervical cancer. HPV-16 variants natural history among males is not established. We evaluated HPV-16 variants prevalence and persistence in the external genitalia of men enrolled in the prospective HPV Infection in Men (HIM) Study.

METHODS

The HIM Study included men from the United States, Brazil, and Mexico. HPV-16 variants were distinguished using polymerase chain reaction sequencing. The prevalence of HPV-16 variants was assessed, and associations with infection persistence were estimated.

RESULTS

We characterized the HPV-16 variants for 1700 genital swab samples from 753 men and 22 external genital lesions in 17 men. The prevalence of HPV-16 lineages differed by country and marital status (P < .001). Overall, 90.9% of participants harbored lineage A variants. The prevalence of non-A lineages was heterogenous among countries. HPV-16 lineage A variants were associated with a 2.69-fold increased risk of long-term persistent infections compared with non-A lineages. All high-grade penile intraepithelial neoplasia harbored lineage A variants and occurred in the context of long-term persistent infections with the same variants.

CONCLUSIONS

The prevalence and persistence of HPV-16 variants observed at the male external genitalia suggest differences in the natural history of these variants between men and women, which may be associated with intrinsic differences in the infected genital epithelia.

HPV16 Intratypic Variants in Head and Neck Cancers: A North American Perspective

Human papillomavirus (HPV) is the major causative agent for cervical and many head and neck cancers (HNCs). HPVs randomly acquire single nucleotide polymorphisms (SNPs) that may become established via positive selection. Within an HPV type, viral isolates differing by <2% in the L1 region are termed "variants" and classified based on combinations of SNPs. Studies in cervical cancer demonstrate clear differences between HPV16 intratypic variants in terms of persistence of infection, tumor histology, cancer risk, and death. Much less is known about the frequency of HPV16 variants in HNC, and their effects on clinical outcomes. We combined HPV16 positive (HPV16+) HNC samples from a local Southwestern Ontario, Canada cohort with those from the Cancer Genome Atlas to create a larger North American cohort of 149 cases with clinical data and determined the distribution of intratypic variants and their impact on clinical outcomes. Most isolates were lineage A, sublineage A1, or A2, with roughly half exhibiting the T350G polymorphism in E6. Univariable analysis identified significant differences between 350T and 350G intratypic variants in clinical T, N, and O staging, as well as disease-free survival. Multivariable analysis failed to identify any clinical factor as a statistically significant covariate for disease-free survival differences between 350T and 350G. Significant differences in several measures of B-cell mediated immune response were also observed between 350T and 350G intratypic variants. We suggest that HPV genetic variation may be associated with HNC clinical characteristics and may have prognostic value.

Distribution of high- and low-risk human papillomavirus genotypes and their prophylactic vaccination coverage among West African women: systematic review

INTRODUCTION

The second most deadly gynecological cancer worldwide, cervical cancer is steadily on the rise in sub-Saharan Africa, while vaccination programs are struggling to get off the ground. This systematic review's aim was to assess the prevalence and distribution of high- and low-risk HPV genotypes in West African women.

METHODS

Original studies were retrieved from PubMed/Medline, Embase, Scopus, Google Scholar, and Science Direct. In these studies, Human papillomavirus (HPV) DNA was assessed in cervical samples by polymerase chain reaction (PCR), Hybrid capture, and sequencing. The quality of the articles was assessed and the results were extracted and reviewed.

RESULTS

Thirty-nine studies from 10 West African countries were included for the systematic review including 30 for the pooled analysis. From an overall of 17358 participants, 5126 of whom were infected with at least one HPV genotype, the systematic review showed a prevalence varying from 8.9% to 81.8% in the general population. In contrast, the pooled prevalence of infection was 28.6% (n = 3890; 95% CI 27.85-29.38), and HPV-52 (13.3%), HPV-56 (9.3%), and HPV-35 (8.2) were the most frequent. Quadrivalent and nonavalent vaccines covered 18.2% and 55.8% of identified genotypes respectively.

CONCLUSION

Faced with this growing public health challenge in West Africa, it would be necessary for all its countries to have reliable data on HPV infection and to introduce the nonavalent vaccine. A study of the genotypic distribution of HPV in high-grade precancerous lesions and cervical cancer would be very useful in West Africa.

The E6 gene polymorphism of Human papillomavirus 16 in relation to the risk of cervical cancer in Tunisian women

Human papillomavirus type 16 (HPV-16) is the most prevalent HPV type worldwide and in Tunisia and the major carcinogenic HPV type found in cervical precancers and cancers. Previous studies have reported that genetic diversity of HPV16-E6 oncoprotein might be associated with cervical intraepithelial neoplasia progression. In this study we aimed to investigate the prevalence of HPV-16 E6 variants in precancerous lesions in Tunisian population to assess potential correlation with disease severity. Positive HPV cervical samples were obtained from the Laboratory of Anatomy Pathology of Pasteur Institute of Tunis. Cytological study was performed to identify cervical precancerous lesions. HPVs were typed using Reverse Line Hybridization. Only samples with HPV-16 single infection were selected for HP16-E6 genetic diversity investigation. HPV-16 E6 gene amplification was performed by PCR using specific primers and sequenced by Sanger Sequencing. The multiple alignment of generated sequences was performed using MEGAX software. Phylogenetic tree was constructed using Maximum Likehood method. The ternary complex of E6, E6AP and p53 core domain was used to perform in silico point mutations and thermodynamic calculations to assess stability and binding affinity. Genetic analysis of Tunisian E6-HPV16 sequences showed the presence of three lineages: European (A), African (C) and Asian American (D). Interestingly, the EUR variants were identified as the dominant lineage of HPV-16 and HPV-16 E6 350 G (L83V) was the most detected mutation in precancerous lesions. Modelling data showed that African variants induced the largest destabilizing effect on E6 structure and decreasing thereby in the affinity toward E6AP. Therefore, women infected with European variants are associated with low and high intraepithelial lesions. The findings give useful information for personalized decision algorithms of intra-epithelial cervical neoplasia in Tunisian women.

Genetic variability analysis of human papillomavirus 58: Novel sublineage identification and persistent infection association

This study aims to characterize the genetic variability of HPV58, identify novel lineages and sublineages, and explore the association between persistent/multiple HPV58 infections and genetic variation. In this study, samples from 124 women with HPV58 infection in Eastern China were collected and 81 isolates of E6 and L1 full-length genes were successfully amplified from 55 samples. We evaluated the diversity of genetic variants and performed correlation analyses between genetic variability and pathology, vaccination, multiple infections, and persistent infections. Among the E6 and L1 gene sequences collected, the dominant prevailing sublineages were A1 (46.2%) and A2 (23.1%). In addition, we found two potential novel sublineages denoted as the A4 and A5 sublineage. A total of 50 nucleotide substitutions, including 28 synonymous substitutions and 22 nonsynonymous substitutions, were observed in the E6 and L1 genes. Among them, variants with A388C/K93N substitutions in the E6 gene correlated with persistent infection (≥1 and ≥2 years) (p < 0.005), and C307T/C66C was associated with persistent infection (≥2 years) (p < 0.005). Notably, two mutations above were detected in the isolate from the patient with breakthrough vaccine infection. Our study found two novel sublineages and sites of genetic variability in multiple and persistent infection variants. In addition, we identified two mutational sites associated with persistent infection. This study provides new insight into the clinical characteristics of HPV 58 genetic variations and offers new ideas for research on next-generation vaccines in Eastern China.

High risk HPV-positive women cervicovaginal microbial profiles in a Greek cohort: a retrospective analysis of the GRECOSELF study

Increasing evidence supports a role for the vaginal microbiome (VM) in the severity of HPV infection and its potential link to cervical intraepithelial neoplasia. However, a lot remains unclear regarding the precise role of certain bacteria in the context of HPV positivity and persistence of infection. Here, using next generation sequencing (NGS), we comprehensively profiled the VM in a series of 877 women who tested positive for at least one high risk HPV (hrHPV) type with the COBAS® 4,800 assay, after self-collection of a cervico-vaginal sample. Starting from gDNA, we PCR amplified the V3-V4 region of the bacterial 16S rRNA gene and applied a paired-end NGS protocol (Illumina). We report significant differences in the abundance of certain bacteria compared among different HPV-types, more particularly concerning species assigned to Lacticaseibacillus, Megasphaera and Sneathia genera. Especially for Lacticaseibacillus, we observed significant depletion in the case of HPV16, HPV18 versus hrHPVother. Overall, our results suggest that the presence or absence of specific cervicovaginal microbial genera may be linked to the observed severity in hrHPV infection, particularly in the case of HPV16, 18 types.

Genetic variation of E6 and E7 genes of human papillomavirus type 16 from central China

BACKGROUND

Persistent high-risk human papillomavirus (HR-HPV) infection is an important factor in the development of cervical cancer, and human papillomavirus type 16 (HPV-16) is the most common HR-HPV type worldwide. The oncogenic potential of HPV-16 is closely related to viral sequence variation.

METHODS

In order to clarify the variant characteristics of HPV-16 E6 and E7 genes in central China, E6 and E7 sequences of 205 HPV-16 positive samples were amplified by polymerase chain reaction. PCR products of E6 and E7 genes were further sequenced and subjected to variation analysis, phylogenetic analysis, selective pressure analysis and B-cell epitope prediction.

RESULTS

Twenty-six single nucleotide variants were observed in E6 sequence, including 21 non-synonymous and 5 synonymous variants. Twelve single nucleotide variants were identified in E7 sequence, including 6 non-synonymous and 6 synonymous variants. Four new variants were found. Furthermore, nucleotide variation A647G (N29S) in E7 was significantly related to the higher risk of HSIL and cervical cancer. Phylogenetic analysis showed that the E6 and E7 sequences were all distributed in A lineage. No positively selected site was found in HPV-16 E6 and E7 sequences. Non-conservative substitutions in E6, H31Y, D32N, D32E, I34M, L35V, E36Q, L45P, N65S and K75T, affected multiple B-cell epitopes. However, the variation of E7 gene had little impact on the corresponding B-cell epitopes (score < 0.85).

CONCLUSION

HPV-16 E6 and E7 sequences variation data may contribute to HR-HPV prevention and vaccine development in Jingzhou, central China.

HPV Testing and its Role in Cervical Cancer Screening

The recognition that persistent infection with carcinogenic human papillomavirus (HPV) is a necessary cause of cervical precancer and cancer has led to the introduction of HPV testing into cervical cancer screening, either as a primary screening test or in conjunction with cervical cytology (i.e., co-testing). HPV testing has much higher sensitivity for detection of cervical precancer and provides greater long-term reassurance if negative compared to cytology. However, most HPV infections are transient, and do not progress to invasive cancer, thus triage tests are required to identify individuals who should be referred to colposcopy for diagnostic evaluation. This chapter begins with a description of the biology, natural history, and epidemiology of HPV as a foundation for understanding the role of HPV in cervical carcinogenesis. This section is followed by a detailed discussion regarding the introduction of HPV-based testing and triage into cervical cancer screening and management. Summarized triage tests include cervical cytology, HPV genotyping, p16/Ki-67 dual stain, and HPV and cellular methylation markers. The final section of this chapter includes an important discussion on cervical cancer disparities, particularly within the United States, followed by concluding remarks.

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.

Human Papillomavirus 16 Lineage D is Associated with High Risk of Cervical Cancer in the Brazilian Northeast Region

OBJECTIVE

 Similar to Human Papillomavirus (HPV) genotypes, different lineages of a genotype also have different carcinogenic capabilities. Studies have shown that specific genotype lineages of oncogenic HPV are associated with variable risks for the development of cervical intraepithelial neoplasia (CIN2/CIN3) and cervical cancer. The present study aimed to analyze the genetic diversity of the HPV16 genotype in women with CIN2/CIN3 and cervical cancer, from the northeast region of Brazil.

METHODS

 A cross-sectional multicenter study was conducted in the northeast region of Brazil, from 2014 to 2016. This study included 196 cases of HPV16 variants (59 and 137 cases of CIN2/CIN3 and cervical cancer, respectively). The difference of proportion test was used to compare patients with CIN2/CIN3 and cervical cancer, based on the prevalent HPV16 lineage (p < 0.05).

RESULTS

 According to the histopathological diagnosis, the percentage of lineage frequencies revealed a marginal difference in the prevalence of lineage A in CIN2/CIN3, compared with that in cervical cancer (p = 0.053). For lineage D, the proportion was higher in cancer cases (32.8%), than in CIN2/CIN3 cases (16.9%), with p = 0.023.

CONCLUSION

 HPV16 lineage A was the most frequent lineage in both CIN2/CIN3 and cervical cancer samples, while lineage D was predominant in cervical cancer, suggesting a possible association between HPV16 lineage D and cervical cancer.

Distinguishing Genetic Drift from Selection in Papillomavirus Evolution

Pervasive purifying selection on non-synonymous substitutions is a hallmark of papillomavirus genome history, but the role of selection on and the drift of non-coding DNA motifs on HPV diversification is poorly understood. In this study, more than a thousand complete genomes representing Alphapapillomavirus types, lineages, and SNP variants were examined phylogenetically and interrogated for the number and position of non-coding DNA sequence motifs using Principal Components Analyses, Ancestral State Reconstructions, and Phylogenetic Independent Contrasts. For anciently diverged Alphapapillomavirus types, composition of the four nucleotides (A, C, G, T), codon usage, trimer usage, and 13 established non-coding DNA sequence motifs revealed phylogenetic clusters consistent with genetic drift. Ancestral state reconstruction and Phylogenetic Independent Contrasts revealed ancient genome alterations, particularly for the CpG and APOBEC3 motifs. Each evolutionary analytical method we performed supports the unanticipated conclusion that genetic drift and different evolutionary drivers have structured Alphapapillomavirus genomes in distinct ways during successive epochs, even extending to differences in more recently formed variant lineages.

Impact of HPV-16 Lineages Infection in Response to Radio-Chemotherapy in Cervical Cancer

BACKGROUND

HPV is strongly related to cervical cancer. HPV lineages can contribute to a response to cervical cancer therapy. The aim of this research was to estimate the frequency of human papillomavirus (HPV)-16 lineages in specimens of cervical cancer, relate the pathological factors in these variants, and assess their response to treatment with radical chemoradiotherapy.

METHODS

Samples of cervical cancer were collected from women who were referred to a reference cancer hospital to test the presence of human papillomavirus-type DNA. The standard protocol of this service consisted of cisplatin-based chemotherapy of 40 mg/m², plus conventional pelvic irradiation in doses of 45-50.4 Gy and high dose-rate brachytherapy of 28-30 Gy to Point A. The response to chemotherapy was evaluated after three months in patients with the HPV-16 lineage.

RESULTS

HPV DNA was detected in 104 (88.1%) of the 118 patients. HPV-16 was present in 63 patients (53%). Lineages of HPV-16 were identified in 57 patients and comprised 33 instances of (57.8%) lineage A, 2 instances of lineage B (3.5%), 2 instances of lineage C (3.5%), and 20 instances of (35.0%) lineage D. The median age of the patients was 48.4 years (range 25-85 years). Squamous cell carcinoma was detected 48 times (84.2%). Adenocarcinoma was more likely to occur in lineage D, as three of the four cases occurred in this lineage. A total of 11 patients with the HPV-16 variant were treated with chemoradiotherapy. After three months, it was observed that nine of the eleven patients (81.8%) achieved a complete response, five with the lineage A type, two with the lineage C type, and two with the lineage D type. The two cases of partial response and disease progression, one of each, occurred in lineage A.

CONCLUSIONS

In addition to the small number of patients and HPV variants, we noticed a better response in patients with the HPV-16 lineage A. Increasing the sample size could be helpful to better assess the impact of HPV variants on cervical cancer treatment.

The combined finding of HPV 16, 18, or 45 and cytologic Atypical Glandular Cells (AGC) indicates a greatly elevated risk of in situ and invasive cervical adenocarcinoma

BACKGROUND

Cervical screening has not effectively controlled cervical adenocarcinoma (AC). Human papillomavirus (HPV) testing is recommended for cervical screening but the optimal management of HPV-positive individuals to prevent AC remains a question. Cytology and HPV typing are two triage options to predict the risk of AC. We combined two potential biomarkers (atypical glandular cell, AGC, cytology and HPV-types 16, 18, or 45) to assess their joint effect on detecting AC.

METHODS

Kaiser Permanente Northern California (KPNC) used triennial co-testing with cytology and HPV testing (positive/negative) for routine cervical screening between 2003 and 2020. HPV typing of a sample of residual HPV test specimens was performed on a separate cohort selected from KPNC (Persistence and Progression, PaP, cohort). We compared risk of prevalent and incident histologic AC/AIS (adenocarcinoma in situ) associated with preceding combinations of cytologic results and HPV typing. Risk of squamous cell cancer (SCC)/cervical intraepithelial neoplasia grade 3 (CIN3) (SCC/CIN3) was also included for comparison.

RESULTS

Among HPV-positive individuals in PaP cohort, 99% of prevalent AC and 96% of AIS were linked to HPV-types 16, 18, or 45 (denoted HPV 16/18/45). Although rare (0.09% of screening population), the concurrent detection of HPV 16/18/45 with AGC cytology predicted a highly elevated relative risk of underlying histologic AC/AIS; the absolute risk of diagnosing AC/AIS was 12% and odds ratio (OR) was 1341 (95%CI:495-3630) compared to patients with other high-risk HPV types and normal cytology. Cumulatively (allowing non-concurrent results), approximately one-third of the AC/AIS cases ever had HPV 16/18/45 and AGC cytology (OR = 1785; 95%CI:872-3656). AGC was not as strongly associated with SCC/CIN3.

CONCLUSION

Detection of HPV 16/18/45 positivity elevates risk of adenocarcinoma, particularly if AGC cytology is also found.

Impact of HPV testing in opportunistic cervical screening: Support for primary HPV screening in the United States

Human papillomavirus (HPV) testing for cervical screening increases diagnosis of precancer and reduces the incidence of cervical cancer more than cytology alone. However, real-world evidence from diverse practice settings is lacking for the United States (U.S.) to support clinician adoption of primary HPV screening. Using a population-based registry, which captures all cervical cytology (with or without HPV testing) and all cervical biopsies, we conducted a real-world evidence study of screening in women aged 30 to 64 years across the entire state of New Mexico. Negative cytology was used to distinguish cotests from reflex HPV tests. A total of 264 198 cervical screening tests (with exclusions based on clinical history) were recorded as the first screening test between 2014 and 2017. Diagnoses of cervical intraepithelial neoplasia grades 2 or 3 or greater (CIN2+, CIN3+) from 2014 to 2019 were the main outcomes. Of cytology-negative screens, 165 595 (67.1%) were cotests and 4.8% of these led to biopsy within 2 years vs 3.2% in the cytology-only group. Among cytology-negative, HPV tested women, 347 of 398 (87.2%) CIN2+ cases were diagnosed in HPV-positive women, as were 147 of 164 (89.6%) CIN3+ cases. Only 29/921 (3.2%) CIN3+ and 67/1964 (3.4%) CIN2+ cases were diagnosed in HPV-negative, cytology-positive women with biopsies. Under U.S. opportunistic screening, across a diversity of health care delivery practices, and in a population suffering multiple disparities, we show adding HPV testing to cytology substantially increased the yield of CIN2+ and CIN3+. CIN3+ was rarely diagnosed in HPV-negative women with abnormal cytology, supporting U.S. primary HPV-only screening.

The features of high-risk human papillomavirus infection in different female genital sites and impacts on HPV-based cervical cancer screening

BACKGROUND

The causal role of high-risk Human papillomavirus (HR-HPV) in the pathogenesis of anogenital cancers is well established. In contrast, information on HR-HPV distribution of continuous anatomic sites within the female genital tract is limited, and the impact of sample type on the clinical performance in HPV-based cervical cancer screening warrants investigation.

METHODS

A total of 2,646 Chinese women were enrolled in the study from May 2006 to April 2007. We analyzed the infection features by infection status and pathological diagnoses of 489 women with complete HR-HPV type and viral load data on the cervix, upper vagina, lower vagina, and perineum samples. Additionally, we assessed the clinical performance for detecting high-grade cervical intraepithelial neoplasia of grade two or worse (≥ CIN2) among these four types of samples.

RESULTS

HR-HPV positivity rate was lower in the cervix (51.53%) and perineum (55.83%), higher in the upper (65.64%) and lower vagina (64.42%), and increased with the severity of cervical histological lesions (all P<0.001). Single infection was more dominant than multiple infections at each anatomic site of the female genital tract. The proportion of single HR-HPV infection decreased successively from the cervix (67.05%) to the perineum (50.00%) (P_trend=0.019) in cervical intraepithelial neoplasia grade 1 (CIN1) and was higher in samples of the cervix (85.11%) and perineum (72.34%) in ≥ CIN2. In addition, the highest viral load was observed in the cervix compared to the other three sites. The overall agreement of the cervical and perineum samples was 79.35% and increased continuously from normal (76.55%) to ≥ CIN2 (91.49%). As for the detection of ≥ CIN2, the sensitivity was 100.00%, 97.87%, 95.74%, and 91.49% for the cervix, upper vagina, lower vagina, and perineum samples, respectively.

CONCLUSIONS

Single HR-HPV infection predominated throughout the female genital tract, but the viral load was lower compared to multiple HR-HPV infections. Despite the decreasing viral load from cervix to perineum, the clinical performance for detecting ≥ CIN2 of the perineum sample was comparable to that of the cervix.

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.

Diversity of cervicovaginal human papillomavirus (HPV) genotypes and naturally occurring E6/E7 DNA polymorphisms of HPV-16 in Ghana

Human papillomavirus (HPV) E6 and E7 oncogene expression is essential for cervical carcinogenesis. Evidence exists that E6/E7 variants may have different transforming activities while the risk of HPV-16 variants (A/D) differs by race/ethnicity. We determined the type-specific diversity of HPV infection in women with high grade cervical disease or cervical cancer in Ghana and investigated naturally occurring E6/E7 DNA variants in this population. HPV genotyping was carried out on 207 cervical swab samples collected from women referred to a gynaecology clinic at two teaching hospitals in Ghana. HPV-16, HPV-18 and HPV-45 were detected in 41.9%, 23.3% and 16.3% of cases respectively. HPV-16 E6/E7 DNA sequencing was performed in 36 samples. Thirty samples contained E6/E7 variants of the HPV-16-B/C lineage. 21/36 samples were of the HPV-16C1 sublineage variant and all contained the E7 A647G(N29S) single nucleotide polymorphism (SNP). This study reveals the diversity of E6/E7 DNA and the dominance of HPV16 B/C variants in cervicovaginal HPV infection in Ghana. Type-specific HPV diversity analysis indicates that most Ghanaian cervical disease cases are vaccine preventable. The study provides an important baseline from which for the impact of vaccine and antivirals on clinically relevant HPV infection and associated disease can be measured.

Genetic diversity and functional implication of the long control region in human papillomavirus types 52, 58, and 16 from Central China

OBJECT

High-risk human papillomavirus (HR-HPV) is a main reason for cervical cancer. The long control region (LCR) of the genome plays a variety of roles in the transcription of the virus.

METHODS

LCR sequences were amplified by polymerase chain reaction (PCR) and confirmed by DNA sequencing. MEGA 11.0 software and NCBI blast were used to analyze the sequences and construct the Neighbor-Joining tree. In addition, the JASPAR database was used to predict the potential transcription factor binding sites (TFBS).

RESULTS

For HPV-52 LCR, 68 single nucleotide polymorphisms (SNPs), 8 deletions, and 1 insertion were found, 17 of which were novel variations. Most of the variants were clustered in B2 sub-lineage (96.22%). For HPV-58 LCR, 25.43% of samples were prototype. 49 SNPs, 2 deletions, and 1 insertion were observed in the remaining samples. A1 sub-lineage was the most frequent (64.16%). For HPV-16 LCR, 75 SNPs and 2 deletions were identified, 13 of which were newly identified. A total of 55.68% of the variants were distributed in A4 sub-lineage. The JASPAR results suggested that multiple variations occurred in TFBSs, which might affect the function of transcription factors.

CONCLUSIONS

This study provides experimental data for further studies on the epidemiology and biological function of LCR. Various LCR mutational data may prove useful for exploring the carcinogenic mechanism of HPV.

Distribution and Clinical Significance of HPV16 Variants in Head and Neck Squamous Cell Carcinomas: Data from a Portuguese Cohort and Systematic Review

INTRODUCTION

Genomic variants of the human papillomavirus type 16 (HPV16) are thought to play differential roles in the susceptibility to head and neck squamous cell carcinomas (HNSCC) and its biological behaviour. This study aimed to establish the prevalence of HPV16 variants in an HNSCC cohort and associate them with clinical pathological characteristics and patient survival.

METHODS

We retrieved samples and clinical data from 68 HNSCC patients. DNA samples were available from tumour biopsy at the time of the primary diagnosis. Targeted next-generation sequencing was used to obtain whole-genome sequences, and variants were established based on phylogenetic classification.

RESULTS

74% of samples clustered in lineage A, 5.7% in lineage B, 2.9% in lineage C, and 17.1% in lineage D. Comparative genome analysis revealed 243 single nucleotide variations. Of these, one hundred were previously reported, according to our systematic review. No significant associations with clinical pathological variables or patient survival were observed. The E6 amino acid variations E31G, L83V, and D25E and E7 N29S, associated with cervical cancer, were not observed, except for N29S in a single patient.

CONCLUSION

These results provide a comprehensive genomic map of HPV16 in HSNCC, highlighting tissue-specific characteristics which will help design tailored therapies for cancer patients.

[New recommendations for cervical cancer screening in France]

In France, recent advances in cervical cancer screening include an organized cervical cancer-screening program and the introduction of HPV testing as a first-line test for women aged 30-65 years. The HPV test, performed on a cervical smear taken by a health professional, could also be performed on a vaginal self-sample in certain indications. The detection kits used to test for HPV should target high-risk HPV, be validated for screening and meet the performance requirements for this indication. Although no longer used as a first-line test in women aged 30-65 years, cytological examination of cervical cells remains important, particularly in the triage of HPV positive women. The interest of other biological techniques, such as HPV genotyping, viral load, cellular expression of p16/Ki-67 proteins and the methylation of cellular or viral genes, still needs to be clarified, but they could help to refine the triage strategy of HPV-positive women and limit the need for colposcopy and unnecessary stress for patients.

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.

Racial and Ethnic Disparities in Cervical Cancer Incidence, Survival, and Mortality by Histologic Subtype

PURPOSE

We conducted an integrated population-based analysis of histologic subtype-specific cervical cancer incidence, survival, and incidence-based mortality by race and ethnicity, with correction for hysterectomy prevalence.

METHODS

Using the SEER 21 and 18 registries, we selected primary cases of malignant cervical cancer diagnosed among women ≥ 15 years. We evaluated age-adjusted incidence rates among cases diagnosed between 2000 and 2018 (SEER21) and incidence-based mortality rates among deaths from 2005 to 2018 (SEER18), per 100,000 person-years. Rates were stratified by histologic subtype and race/ethnicity (incidence and mortality), and stage, age at diagnosis, and county-level measures of social determinants of health (incidence only). Incidence and mortality rates were corrected for hysterectomy using data from the Behavioral Risk Factor Surveillance System. We estimated 5-year relative survival by histologic subtype and stratified by stage at diagnosis.

RESULTS

Incidence rates of cervical squamous cell carcinoma were highest in Black and Hispanic women, while incidence rates of cervical adenocarcinoma (ADC) were highest among Hispanic and White women, particularly for localized ADC. County-level income and education variables were inversely associated with squamous cell carcinoma incidence rates in all racial and ethnic groups but had less influence on ADC incidence rates. Black women had the highest overall mortality rates and lowest 5-year relative survival, irrespective of subtype and stage. Disparities in survival were particularly pronounced for Black women with regional and distant ADC, compared with other racial/ethnic groups.

CONCLUSION

Although Black women are less likely to be diagnosed with ADC compared with all other racial/ethnic groups, they experience the highest mortality rates for this subtype, likely attributed to the poor survival observed for Black women with regional and distant ADC.

Mapping HPV 16 Sub-Lineages in Anal Cancer and Implications for Disease Outcomes

The incidence of anal cancer is rising worldwide. As identified in cervical cancer management, an improvement in the early detection and management of anal pre-cancer is essential. In other cancers associated with human papillomavirus (HPV), HPV 16 sub-lineages have been shown to be associated with disease status and prognosis. However, in anal cancer, they have been under-explored. A total of 119 HPV 16-positive anal cancer lesions diagnosed between 2009 and 2018 in Scotland and 134 HPV 16-positive residual rectal swabs from asymptomatic men collected in 2016/7 were whole genome sequenced. The association of HPV 16 sub-lineages with underlying disease status (cancer vs. asymptomatic) and overall survival in anal cancer samples was assessed (comparing A1 vs non-A1 sub-lineages). A1 was the dominant sub-lineage present in the anal cancer (76.5%) and the asymptomatic (76.1%) cohorts. A2 was the second most dominant sub-lineage in both groups (16.8% and 17.2%, respectively). We did not observe significant associations of sub-lineage with demographics, clinical variables or survival (A1 vs. non-A1 sub-lineages (HR 0.83, 0.28-2.46 p = 0.743)). HPV 16 sub-lineages do to not appear to cluster with disease vs asymptomatic carriage or be independently associated with outcomes in anal cancer patients. Further international studies on anal HPV sub-lineage mapping will help to determine whether this is a consistent observation.

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.

Comprehensive Viral Genotyping Reveals Prognostic Viral Phylogenetic Groups in HPV16-Associated Squamous Cell Carcinoma of the Oropharynx

Human papillomavirus-positive (HPV+) squamous cell carcinoma of the oropharynx (OPSCC) is the most prevalent HPV-associated malignancy in the United States and is primarily caused by HPV subtype 16 (HPV16). Favorable treatment outcomes have led to increasing interest in treatment deescalation to reduce treatment-related morbidity. Prognostic biomarkers are needed to identify appropriately low-risk patients for reduced treatment intensity. Targeted DNA sequencing including all HPV16 open reading frames was performed on tumors from 104 patients with HPV16+ OPSCC treated at a single center. Genotypes closely related to the HPV16-A1 reference were associated with increased numbers of somatic copy-number variants in the human genome and poor recurrence-free survival (RFS). Genotypes divergent from HPV16-A1 were associated with favorable RFS. These findings were independent of tobacco smoke exposure. Total RNA sequencing was performed on a second independent cohort of 89 HPV16+ OPSCC cases. HPV16 genotypes divergent from HPV16-A1 were again validated in this independent cohort, to be prognostic of improved RFS in patients with moderate (less than 30 pack-years) or low (no more than 10 pack-years) of tobacco smoke exposure. In summary, we show in two independent cohorts that viral sequence divergence from the HPV16-A1 reference is correlated with improved RFS in patients with moderate or low tobacco smoke exposure.

IMPLICATIONS

HPV16 genotype is a potential biomarker that could be easily adopted to guide therapeutic decision-making related to deescalation therapy.

Genetic diversity of E6, E7 and the long control region in human papillomavirus type 16 variants in Beijing, China

Background

High-risk human papillomavirus type 16 (HPV16) is a risk factor for cervical cancer. The progression from initial infection to cervical cancer has been linked to properties of the viral sequences. However, the distribution of HPV16 variants among Chinese women has not been extensively addressed and the role of HPV16 variants in the risk of cervical carcinogenesis remains poorly understood.

Methods

HPV16 positive cervical exfoliated cell samples were collected from 249 women living in Beijing, China. PCR products from two fragments of E6-E7 and LCR of HPV16 in these samples were sequenced and analyzed.

Results

Lineage A was found in the subjects, including A1, A2, A3 and A4 sublineages. Based on the HPV16 reference sequences, 26 nucleotide mutations of A4 sublineage and 39 nucleotide mutations of A1-3 sublineages were found in the E6, E7 and LCR of HPV16 isolates. Point mutations T843C, A7287C and A7872G of A4 sublineage were significantly associated with high-grade cervical lesions. The high-frequency sites in HPV16 LCR located at regions that can bind to multiple transcription factors.

Conclusions

This study contributes to the identification of unique variants in the E6, E7 and LCR of HPV16 isolates infected in Chinese women. Mutations of T843C, A7287C and A7872G in A4 sublineages were significantly associated with high-grade cervical lesions, suggesting that mutations in the E7 and LCR region have potential effects on viral replication and progression of cervical cancer.

•

Lineage A is the predominant HPV16 variant in Chinese women.

•

Based on the HPV16 prototype, 25 nucleotide mutations in E6 and E7 were found and 41 nucleotide mutations in LCR were found.

•

The mutations of C843T, A7287C and A7872G were significantly correlated with the cervical lesions.

The T350G Variation of Human Papillomavirus 16 E6 Gene Prevails in Oropharyngeal Cancer from a Small Cohort of Greek Patients

Recent trends have shown a dramatic rise in the incidence of oropharyngeal squamous cell carcinoma strongly associated with high-risk human papillomavirus (HPV) of type 16. The genetic variability of HPV16 has been extensively studied in cervical cancer but there are very limited published data concerning the genetic variations of this HPV type in oropharyngeal cancer. In the present study, the genetic variations of HPV16 E6 gene sequences originated from a small cohort of Greek patients diagnosed with oropharyngeal cancer were assessed. The vast majority of the sequences clustered within the European variant branch. The T350G variation was found to be the predominant one. This finding may indicate the need for further studies that could explain the possible impact of this variant in the pathomechanisms of oropharyngeal cancer.

Epidemiology and Molecular Biology of HPV Variants in Cervical Cancer: The State of the Art in Mexico

Cervical cancer (CC) continues to be a major public health problem in Mexico, ranking second among cancers in women. A persistent infection with human papillomaviruses (HPV) is the main risk factor for CC development. In addition, a significant fraction of other cancers including those of the anus, oropharynx, and penis are also related to HPV infection. In CC, HPV-16 is the most prevalent high-risk HPV type, followed by HPV-18, both being responsible for 70% of cases. HPV intratype variant lineages differ in nucleotide sequences by 1–10%, while sublineages differ by 0.5–1%. Several studies have postulated that the nucleotide changes that occur between HPV intratype variants are reflected in functional differences and in pathogenicity. Moreover, it has been demonstrated that HPV-16 and -18 intratype variants differentially affect molecular processes in infected cells, changing their biological behavior that finally impacts in the clinical outcome of patients. Mexico has participated in providing knowledge on the geographical distribution of intratype variants of the most prevalent HPVs in premalignant lesions of the cervix and cervical cancer, as well as in other HPV-related tumors. In addition, functional studies have been carried out to assess the cellular effects of intratype variations in HPV proteins. This review addresses the state of the art on the epidemiology of HPV-16 and HPV-18 intratype variants in the Mexican population, as well as their association with persistence, precancer and cervical cancer, and functional aspects related to their biological behavior.

A Straightforward HPV16 Lineage Classification Based on Machine Learning

Human Papillomavirus (HPV) is the causal agent of 5% of cancers worldwide and the main cause of cervical cancer and it is also associated with a significant percentage of oropharyngeal and anogenital cancers. More than 60% of cervical cancers are caused by HPV16 genotype, which has been classified into lineages (A, B, C, and D). Lineages are related to the progression of cervical cancer and the current method to assess lineages is by building a Maximum Likelihood Tree (MLT); which is slow, it cannot assess poor sequenced samples, and annotation is done manually. In this study, we have developed a new model to assess HPV16 lineage using machine learning tools. A total of 645 HPV16 genomes were analyzed using Genome-Wide Association Study (GWAS), which identified 56 lineage-specific Single Nucleotide Polymorphisms (SNPs). From the SNPs found, training-test models were constructed using different algorithms such as Random Forest (RF), Support Vector Machine (SVM), and K-nearest neighbor (KNN). A distinct set of HPV16 sequences (n = 1,028), whose lineage was previously determined by MLT, was used for validation. The RF-based model allowed a precise assignment of HPV16 lineage, showing an accuracy of 99.5% in the known lineage samples. Moreover, the RF model could assess lineage to 273 samples that MLT could not determine. In terms of computer consuming time, the RF-based model was almost 40 times faster than MLT. Having a fast and efficient method for assigning HPV16 lineages, could facilitate the implementation of lineage classification as a triage or prognostic marker in the clinical setting.