Identification of novel E6-E7 sequence variants of human papillomavirus 16

The rate of evolution of the human papillomavirus 16 (HPV16) genome is low. However, the ability of the E6 oncoprotein to interact with distinct p53 variants causes selective pressure on the E6 gene. In addition, intratypic recombination events in the HPV16 E6 and E7 genes have been characterized as extraordinary phenomena during the evolutionary history of virus. In the present study, we identified two new sequence variants through nucleotide analysis of the E6-E7 region of the HPV16 genome. Maximum-likelihood and empirical Bayesian methods were used in order to identify positive selection at particular residues of the E6 and E7 genes. Using the single recombination breakpoint (SBP) method, we found evidence of recombination events in the E6 ORF. Nucleotide sequence analysis showed that the new sequence variants are phylogenetically distant from the other members of the population. Our results indicate that new evolutionary intermediates of HPV16 might be formed either though positive selective pressure or through recombination events by multiple infections with distinct HPV16 variants.