Abstract
The infection with high-risk human papillomavirus (HR-HPV) is the most important risk factor for development of cervical cancer. The intra-type variations of HPV have different biological and pathological consequences with respect to disease progression. In the present study, six major Indian variants were experimentally identified in E6 gene of HPV-16 and showed their impact on immunogenicity by in silico methods. Four different phylogenetic lineages were observed in sequences including European (E) prototype, European variant, Asian and American Asian variant classes and complete absence of African phylogenetic lineages. On the prediction of B- and T-cell epitopes, 18 and 23 potent epitopes for MHC-II alleles, 10 potent MHC-I and 15 B-cell epitopes in each reference and variant sequence were identified. Interestingly, the presence of variation H78Y and L83V result in creation of four new epitopes for the HLA-DQA1*0101/DQB1*0501. Out of 15 B-cell predicted epitopes, three most potent epitopes were identified in both reference and variant sequence. Notably the amino acid stretch from amino acid 16-60 and 76-94 are very important for the immunological properties of E6 protein because these regions contain majority of the predicted epitopes. In future, this could control the cervical cancer by targeting these amino acid stretches for the development of HPV-16 vaccine.