AbstractCervical cancer is one of the most common gynecological cancers in the world but in India, it is the top most cancer among women. Persistent infection with high-risk human papillomaviruses (HR-HPVs) is the most important risk factor. The sequence variation(s) in the most common HR-HPV i.e. HPV type 16 leads to altered biological functions with possible clinical significance in the different geographical locations. Sixteen major variants (V1-V16) in full length L1 gene of HPV-16 were identified following analysis of 250 prospectively collected cervical cancer tissue biopsies and their effect on immunogenicity was studied. The effect of these major variations on the epitopes were predicted by in silico methods and the immunogenicity of variants and respective reference DNA vaccine constructs were evaluated by administration of prepared DNA vaccine constructs in female BALB/c mice to evaluate antibody titer. In the present study, L500F (V16) variation showed a significant ~2.7 fold (p < 0.002) increase in antibody titer, whereas T379P (V8) showed ~0.4 fold (p < 0.328) decrease after final injection. These results showed a promising roadmap for the development of DNA based vaccine and for the generation of effective response, though there is a need to study more prevalent variants of HPV in the Indian population.