Genetic variability analysis of human papillomavirus 58: Novel sublineage identification and persistent infection association

This study aims to characterize the genetic variability of HPV58, identify novel lineages and sublineages, and explore the association between persistent/multiple HPV58 infections and genetic variation. In this study, samples from 124 women with HPV58 infection in Eastern China were collected and 81 isolates of E6 and L1 full-length genes were successfully amplified from 55 samples. We evaluated the diversity of genetic variants and performed correlation analyses between genetic variability and pathology, vaccination, multiple infections, and persistent infections. Among the E6 and L1 gene sequences collected, the dominant prevailing sublineages were A1 (46.2%) and A2 (23.1%). In addition, we found two potential novel sublineages denoted as the A4 and A5 sublineage. A total of 50 nucleotide substitutions, including 28 synonymous substitutions and 22 nonsynonymous substitutions, were observed in the E6 and L1 genes. Among them, variants with A388C/K93N substitutions in the E6 gene correlated with persistent infection (>= 1 and >= 2 years) (p < 0.005), and C307T/C66C was associated with persistent infection (>= 2 years) (p < 0.005). Notably, two mutations above were detected in the isolate from the patient with breakthrough vaccine infection. Our study found two novel sublineages and sites of genetic variability in multiple and persistent infection variants. In addition, we identified two mutational sites associated with persistent infection. This study provides new insight into the clinical characteristics of HPV 58 genetic variations and offers new ideas for research on next-generation vaccines in Eastern China.