Chikungunya virus: genomic microevolution in Eastern India and its in-silico epitope prediction

This is the first study reporting whole genome sequences of two CHIKV strains (KJ679577 and KJ679578) isolated from Eastern Indian patients sera during 2010–2011 outbreak, both of which were of ECSA genotype, but from different subgroups: Indian Ocean outbreak and ECSA subtypes. Furthermore, viral sequences were analyzed using different in-silico approaches to identify potential genetic variations that might have functional implications on various aspects of virus replication, viral protein functionality, immunogenicity and transmission. Epitope prediction analysis revealed 70.9% increase in number of MHC Class-II interacting epitopes of KJ679578 and 25–28% increase in Class-I interacting epitopes of KJ679577 and KJ679578 compared to that of EF027141 (CHIKV of Asian genotype circulating in India during 1973, after which CHIKV infection disappeared from India for three decades). CHIKV peptides DLAKLAFKRSSKYDLECAQIPVHMKSDA and KVVLCGDPKQCGFFNMMQMKYNYNHNI were predicted to interact with maximum number of HLA Class-I (68 and 76.5%, respectively) and Class-II (47 and 100%, respectively) alleles present within Indian population with allele frequency of > 0.1 and were also recognized as predicted B-cell epitopes with BCPred score between 0.766 and 0.961 and with antigenicity ranging from 0.52 to 1.69; thus these peptides might be used to induce T- and B-cell-mediated immunity against CHIKV. Thus, the present study might help to bridge the gap between virus microevolution and its implication in host immunity by taking into account viral genetic and conformational changes. Predicted epitopes might be used as promising targets for peptide-based vaccine development and rapid diagnostics against CHIKV infection. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s13205-018-1339-3) contains supplementary material, which is available to authorized users.