Genetic diversity and bioinformatic analysis in the L1 gene of HPV genotypes 31, 33, and 58 circulating in women with normal cervical cytology

HPV-31, -33, and -58, along with HPV-45 and -52, account for almost 11% of HPV-associated cancers. Our previous studies showed that after HPV-16 and -51, HPV-58 was common and HPV-31 was as frequent as HPV-18 among Iranian women with normal cytology. Hence, in this study, we aimed to investigate the intra-type variations in L1 genes of HPV-58, -31, and -33 to find the predominant lineages circulating in women with normal cytology. Complete coding sequencing of the L1 gene was amplified and nucleotide and amino acid sequences were compared to those of the references. The selective pressure on L1 protein and whether the variations of the L1 genes embed in L1 loops, or N-glycosylated sites were also investigated. B1, A, and A1 (sub)lineages were common in the HPV-58, -33, and -31 samples, respectively. Ninety nucleotide mutations were observed. Twenty nine nucleotide changes corresponded to nonsynonymous substitutions in which seventeen mutations were located in L1 loops. Only one codon position in HPV-58 sequences was found as the positive selection. No difference was observed in N-glycosylation sites between reference and understudied amino acid sequences. In the current study, we reported, for the first time, the (sub) lineages, amino acid, and genetic diversity in the L1 gene of circulating HPV-58, -33, and -31, in women with normal cytology, in Iran. Such studies can not only have epidemiological values, but also aid to set vaccination programs.