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

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.

Characterization of the Human Papillomavirus 16 Oncogenes in K14HPV16 Mice: Sublineage A1 Drives Multi-Organ Carcinogenesis

The study of human papillomavirus (HPV)-induced carcinogenesis uses multiple in vivo mouse models, one of which relies on the cytokeratin 14 gene promoter to drive the expression of all HPV early oncogenes. This study aimed to determine the HPV16 variant and sublineage present in the K14HPV16 mouse model. This information can be considered of great importance to further enhance this K14HPV16 model as an essential research tool and optimize its use for basic and translational studies. Our study evaluated HPV DNA from 17 samples isolated from 4 animals, both wild-type (n = 2) and HPV16-transgenic mice (n = 2). Total DNA was extracted from tissues and the detection of HPV16 was performed using a qPCR multiplex. HPV16-positive samples were subsequently whole-genome sequenced by next-generation sequencing techniques. The phylogenetic positioning clearly shows K14HPV16 samples clustering together in the sub-lineage A1 (NC001526.4). A comparative genome analysis of K14HPV16 samples revealed three mutations to the human papillomaviruses type 16 sublineage A1 representative strain. Knowledge of the HPV 16 variant is fundamental, and these findings will allow the rational use of this animal model to explore the role of the A1 sublineage in HPV-driven cancer.

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

Background

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

Methods

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

Results

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

Conclusion

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

Supplementary Information

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

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

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