Human papillomavirus type-6 (HPV6) is the major etiological agent of anogenital warts both men and women. The present study aimed to characterize the genetic diversity among HPV6 in Southwest
China, and to investigate the origin of, selective pressure experienced by, and impact of the resultantly identified genetic variants on the HPV6 secondary structure.
Phylogenetic trees were constructed by Maximum-likelihood and the Kimura 2-parameters methods by Molecular Evolutionary Genetics Analysis version 6.0. The diversity of secondary structure was analyzed by PSIPred software. The selection pressures acting on the E6/E7 genes were estimated by Phylogenetic Analyses by Maximum Likelihood version 4.8 software.
HPV6 was the most prevalent low risk HPV type in southwest
China. In total, 143 E6 and E7 gene sequences of HPV6 isolated from patients were sequenced and compared to GenBank HPV6 reference sequence X00203. The results of these analyses revealed that both the HPV6 E6 and E7 were highly conserved within the analyzed patient samples, and comprised only 3 types of variant sequence, respectively. Furthermore, the analysis of HPV6 E6 and E7 sequences revealed seven/five single-nucleotide mutations, two/four and five/one of which were non-synonymous and synonymous, respectively. The phylogenetic analyses of the E6 and E7 sequences indicated that they belonged to sub-lineage A1 and sub-lineage B1, whereas the selective pressure analyses showed that only the E7 mutation sites 4R, 34E, and 52F were positive selection.
HPV6 (detection rate = 13.10%) was very prevalent in southwest
China, both the HPV6 E6 and E7 sequences were highly conserved within the analyzed patient samples in southwest
China, indicating that the low risk HPV6 can adapt to the environment well without much evolution.