Genomic diversity of human papillomavirus type 6 from patients with condyloma acuminatum in Eastern China

Abstract Objective Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts both men and women. However, there is limited data on its genomic characterization in mainland China. The aim of this study was to understand the complete genomic diversity of HPV6 from patients with condyloma acuminatum (CA) and to explore the prevalence of different variant lineages/sublineages in eastern China. Methods CA samples were collected in 3 hospitals in Shandong Province, China from January 2020 to March 2021. DNA extraction, PCR amplification, Sanger sequencing and sequence assembly were performed on HPV6-positive samples. The complete genomes obtained in this study were analyzed phylogenetically with global HPV6 sequences in GenBank database using MEGA 11. Results A total of 55 complete genomic sequences of HPV6 were obtained in this study. They were classified as HPV6 variant lineage A (n = 20), sublineage B1 (n = 34) and sublineage B3 (n = 1) by phylogenetic analysis. Sequence alignment showed E1, E5A, E5B, L1, L2, LCR were relatively highly variable regions for sublineage B1 whereas E1, E5A, L2 for lineage A. Both phylogenetic trees of lineage A and sublineage B1 composed of two main branches. Chinese sequences of lineage A segregated into the major branch while those in sublineage B1 belonged to both branches. Genomic divergence between sequences from China and other countries was 0.00% - 0.33% in lineage A and 0.00% - 0.40% in sublineage B1. Conclusions This is the first study on HPV variant lineages circulating in mainland China. The results revealed that lineage A and sublineage B1 were prevalent and they had different highly variable regions. Further surveillance is needed to understand the dynamic change of different variants in the population.