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

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.

Identification and Functional Implications of the E5 Oncogene Polymorphisms of Human Papillomavirus Type 16

The persistence of the human papillomavirus type 16 (HPV16) infection on the cervical epithelium contributes to the progression of cervical cancer. Studies have demonstrated that HPV16 genetic variants may be associated with different risks of developing cervical cancer. However, the E5 oncoprotein of HPV16, which is related to several cellular mechanisms in the initial phases of the infection and thus contributes to carcinogenesis, is still little studied. Here we investigate the HPV16 E5 oncogene variants to assess the effects of different mutations on the biological function of the E5 protein. We detected and analyzed the HPV16 E5 oncogene polymorphisms and their phylogenetic relationships. After that, we proposed a tertiary structure analysis of the protein variants, preferential codon usage, and functional activity of the HPV16 E5 protein. Intra-type variants were grouped in the lineages A and D using in silico analysis. The mutations in E5 were located in the T-cell epitopes region. We therefore analyzed the interference of the HPV16 E5 protein in the NF-kB pathway. Our results showed that the variants HPV16E5_49PE and HPV16E5_85PE did not increase the potential of the pathway activation capacity. This study provides additional knowledge about the mechanisms of dispersion of the HPV16 E5 variants, providing evidence that these variants may be relevant to the modulation of the NF-κB signaling pathway.

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.

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.

Intra-Patient Genomic Variations of Human Papillomavirus Type 31 in Cervical Cancer and Precancer

Human papillomavirus type 31 (HPV31) is detected less frequently in cervical cancer than two major causative types, HPV16 and HPV18. Here, we report a comprehensive analysis of HPV31 genome sequences in cervical lesions collected from Japanese women. Of 52 HPV31-positive cervical specimens analyzed by deep sequencing, 43 samples yielded complete genome sequences of around 7900 base pairs and 9 samples yielded partially deleted genome sequences. Phylogenetic analysis showed that HPV31 variant distribution was lineage A in 19 samples (36.5%), lineage B in 28 samples (53.8%), and lineage C in 5 samples (9.6%), indicating that lineage B variants are dominant among HPV31 infections in Japan. Deletions in the viral genome were found in the region from the E1 to L1 genes, but all the deleted genomes retained the E6/E7 genes. Among intra-patient nucleotide variations relative to a consensus genome sequence in each sample, C-to-T substitutions were most frequently detected, followed by T-to-C and C-to-A substitutions. High-frequency, intra-patient mutations (>10%) in cervical cancer samples were found in the E1, E2, and E7 genes, and all of them were nonsynonymous substitutions. The enrichment of high-frequency nonsynonymous substitutions strongly suggests that these intra-patient mutations are positively selected during the development of cervical cancer/precancer.

Genetic variation within the human papillomavirus type 16 genome is associated with oropharyngeal cancer prognosis

PURPOSE

A significant barrier to adoption of de-escalated treatment protocols for human papillomavirus-driven oropharyngeal cancer (HPV-OPC) is that few predictors of poor prognosis exist. We conducted the first large whole-genome sequencing (WGS) study to characterize the genetic variation of the HPV type 16 (HPV16) genome and to evaluate its association with HPV-OPC patient survival.

PATIENTS AND METHODS

A total of 460 OPC tumor specimens from two large United States medical centers (1980-2017) underwent HPV16 whole-genome sequencing. Site-specific variable positions [single nucleotide polymorphisms (SNPs)] across the HPV16 genome were identified. Cox proportional hazards model estimated hazard ratios (HRs) and 95% confidence intervals (CIs) for overall survival by HPV16 SNPs. Harrell C-index and time-dependent positive predictive value (PPV) curves and areas under the PPV curves were used to evaluate the predictive accuracy of HPV16 SNPs for overall survival.

RESULTS

A total of 384 OPC tumor specimens (83.48%) passed quality control filters with sufficient depth and coverage of HPV16 genome sequencing to be analyzed. Some 284 HPV16 SNPs with a minor allele frequency ≥1% were identified. Eight HPV16 SNPs were significantly associated with worse survival after false discovery rate correction (individual prevalence: 1.0%-5.5%; combined prevalence: 15.10%); E1 gene position 1053 [HR for overall survival (HRos): 3.75, 95% CI 1.77-7.95; Pfdr = 0.0099]; L2 gene positions 4410 (HRos: 5.32, 95% CI 1.91-14.81; Pfdr = 0.0120), 4539 (HRos: 6.54, 95% CI 2.03-21.08; Pfdr = 0.0117); 5050 (HRos: 6.53, 95% CI 2.34-18.24; Pfdr = 0.0030), and 5254 (HRos: 7.76, 95% CI 2.41-24.98; Pfdr = 0.0030); and L1 gene positions 5962 (HRos: 4.40, 95% CI 1.88-10.31; Pfdr = 0.0110) and 6025 (HRos: 5.71, 95% CI 2.43-13.41; Pfdr = 0.0008) and position 7173 within the upstream regulatory region (HRos: 9.90, 95% CI 3.05-32.12; Pfdr = 0.0007). Median survival time for patients with ≥1 high-risk HPV16 SNPs was 3.96 years compared with 18.67 years for patients without a high-risk SNP; log-rank test P < 0.001. HPV16 SNPs significantly improved the predictive accuracy for overall survival above traditional factors (age, smoking, stage, treatment); increase in C-index was 0.069 (95% CI 0.019-0.119, P < 0.001); increase in area under the PPV curve for predicting 5-year survival was 0.068 (95% CI 0.015-0.111, P = 0.008).

CONCLUSIONS

HPV16 genetic variation is associated with HPV-OPC prognosis and can improve prognostic accuracy.

Investigation of the diversity of human papillomavirus 16 variants and L1 antigenic regions relevant for the prevention of human papillomavirus-related oropharyngeal cancer in Japan

OBJECTIVE

This study aimed to investigate the distribution of human papillomavirus 16 (HPV16) variants that contribute to the development of HPV-related oropharyngeal carcinoma (HPV-OPC) in the Japanese population and to evaluate genetic variations in the sequence encoding the L1 antigen region of the viral outer shell that is targeted by existing vaccines and is relevant for designing a prevention strategy to combat the exponential increase in HPV-OPC cases in Japan.

METHODS

Seventy Japanese HPV-OPC patients treated at Tohoku University Hospital were included in the study. DNA was extracted from formalin-fixed, paraffin-embedded tissue samples. Polymerase chain reaction and direct nucleotide sequencing were performed to determine the nucleotide polymorphisms necessary for the classification of HPV16 variants and to assess genetic diversity in the HPV16 L1 antigen region, including the BC, DE, EF, FG, and HI loops.

RESULTS

The most common variant of HPV16 was the A4 sublineage (88.6%), conventionally called the Asian type, followed by the A1/2/3 (10.0%) sublineage, classified as the European type. The only nonsynonymous substitution detected in the L1 antigen loop region was p.N181T in the EF loop, which was found in 28/70 (40%) cases. In contrast, no nonsynonymous substitutions were observed in the DE, FG, and HI loops, which are particularly important regions in the antigen loop targeted by existing HPV vaccines.

CONCLUSION

The most common HPV16 variant in Japanese HPV-OPC patients was the A4 subtype. The L1 antigen region is highly conserved, suggesting sufficient efficacy of existing HPV vaccines. These findings provide important information that will aid in the design of an HPV16 infection control strategy using existing HPV vaccines to prevent the spread of HPV-OPC in Japan.

Ancient Evolutionary History of Human Papillomavirus Type 16, 18 and 58 Variants Prevalent Exclusively in Japan

Human papillomavirus (HPV) is a sexually transmitted virus with an approximately 8-kilo base DNA genome, which establishes long-term persistent infection in anogenital tissues. High levels of genetic variations, including viral genotypes and intra-type variants, have been described for HPV genomes, together with geographical differences in the distribution of genotypes and variants. Here, by employing a maximum likelihood method, we performed phylogenetic analyses of the complete genome sequences of HPV16, HPV18 and HPV58 available from GenBank (n = 627, 146 and 157, respectively). We found several characteristic clusters that exclusively contain HPV genomes from Japan: two for HPV16 (sublineages A4 and A5), one for HPV18 (sublineage A1) and two for HPV58 (sublineages A1 and A2). Bayesian phylogenetic analyses of concatenated viral gene sequences showed that divergence of the most recent common ancestor of these Japan-specific clades was estimated to have occurred ~98,000 years before present (YBP) for HPV16 A4, ~39,000 YBP for HPV16 A5, ~38,000 YBP for HPV18 A1, ~26,000 for HPV58 A1 and ~25,000 YBP for HPV58 A2. This estimated timeframe for the divergence of the Japan-specific clades suggests that the introduction of these HPV variants into the Japanese archipelago dates back to at least ~25,000 YBP and provides a scenario of virus co-migration with ancestral Japanese populations from continental Asia during the Upper Paleolithic period.

Prognostic significance of human papillomavirus 16 viral load level in patients with oropharyngeal cancer

Human papillomavirus (HPV) infection in patients with oropharyngeal squamous cell carcinoma (OPSCC) is a major determinant for better prognosis. However, there remain HPV‐positive patients who have poor outcomes. The stratification strategy for detecting high‐risk patients among those with HPV‐positive OPSCC has not been well delineated, especially for Asian patients. We undertook a retrospective cohort study on the survival rate of 89 Japanese patients diagnosed with primary OPSCC. The tumors were concurrently analyzed for the presence of HPV E6 DNA/mRNA, viral DNA load, p16 expression, viral physical status, and viral variant lineage. Human papillomavirus 16 viral DNA was found in 45 (51%) OPSCCs. Human papillomavirus 16 DNA‐positive OPSCCs with higher viral load (classified as HPV16 DNA‐medium/high OPSCCs) showed significantly favorable overall survival and progression‐free survival compared with HPV16 DNA‐positive OPSCCs with lower viral load (<10 copies/cell; HPV16 DNA‐low OPSCCs) and HPV16 DNA‐negative OPSCCs. E6 mRNA expression was observed in all HPV16 DNA‐medium/high OPSCCs but not in HPV16 DNA‐low OPSCCs. Notably, p16‐positive and HPV16 DNA‐negative/low OPSCCs showed significantly worse survival than p16‐positive and HPV16 DNA‐medium/high OPSCCs and resembled HPV‐unrelated OPSCCs with regard to survival and risk factor profile. Although not significant, a trend toward shorter survival was observed for HPV16‐integrated OPSCCs. Phylogenetic analysis revealed two major types of HPV16 variants termed Asian (A4) and European (A1/A2/A3) variants, but no difference in survival between these variants was observed. Altogether, these findings suggest that HPV viral load is a potentially informative factor for more accurate risk stratification of patients with OPSCC.

Human Papillomavirus (HPV) Genotyping Assay Suitable for Monitoring the Impact of the 9-Valent HPV Vaccine

In Japan, a bivalent human papillomavirus (HPV) vaccine against carcinogenic HPV16/18 was licensed in 2009, and a quadrivalent vaccines against HPV16/18 and non-carcinogenic HPV6/11 was licensed in 2011. Recently, the next-generation 9-valent vaccine targeting HPV6/11/16/18/31/33/45/52/58 has been approved. Accurate HPV genotyping is essential for HPV vaccine research and surveillance. The Roche Linear Array (LA) has long been a standard assay for HPV genotyping, but its recent product discontinuation notice has urged us to introduce an alternative assay with comparable performance. In the present study, an in-house HPV genotyping assay that employs PCR with PGMY09/11 primers and reverse blotting hybridization (PGMY-CHUV) was compared with LA to assess genotype-specific agreement. A total of 100 cervical precancer specimens were subjected to both PGMY-CHUV and LA. For detection of genotypes included in the 9-valent vaccine, PGMY-CHUV completely agreed with LA for detection of HPV6, HPV11, HPV16, HPV18, HPV33 and HPV45, and showed near-complete agreement for HPV31 and HPV58 (98% and 99%, respectively). Moreover, PGMY-CHUV detected a significantly higher prevalence of HPV52 than LA (22% vs. 14%, P = 0.008 by McNemar's exact test), with 92.0% overall agreement, 63.6% positive agreement and a kappa value of 0.73. Most (87.5%) of HPV52 discordant cases involved mixed infections with HPV35 or HPV58. In conclusion, while the two assays present equivalent data for assessing the effectiveness of the bivalent and quadrivalent vaccines, PGMY-CHUV is more suitable for evaluating the impact of the current 9-valent vaccine because of its superior detection of HPV52 in co-infection cases.

Metabolic Reprogramming in Cancer: Role of HPV 16 Variants

Metabolic reprogramming is considered one of the hallmarks in cancer and is characterized by increased glycolysis and lactate production, even in the presence of oxygen, which leads the cancer cells to a process called “aerobic glycolysis” or “Warburg effect”. The E6 and E7 oncoproteins of human papillomavirus 16 (HPV 16) favor the Warburg effect through their interaction with a molecule that regulates cellular metabolism, such as p53, retinoblastoma protein (pRb), c-Myc, and hypoxia inducible factor 1α (HIF-1α). Besides, the impact of the E6 and E7 variants of HPV 16 on metabolic reprogramming through proteins such as HIF-1α may be related to their oncogenicity by favoring cellular metabolism modifications to satisfy the energy demands necessary for viral persistence and cancer development. This review will discuss the role of HPV 16 E6 and E7 variants in metabolic reprogramming and their contribution to developing and preserving the malignant phenotype of cancers associated with HPV 16 infection.