Viral whole genome sequencing reveals high variations in APOBEC3 editing between HPV risk categories

High-risk human papillomavirus (HPV) infections are responsible for cervical cancer. However, little is known about the differences between HPV types and risk categories regarding their genetic diversity and particularly APOBEC3-induced mutations – which contribute to the innate immune response to HPV. Using a capture-based next-generation sequencing, 156 HPV whole genome sequences covering 43 HPV types were generated from paired cervical and anal swabs of 30 Togolese female sex workers (FSWs) sampled in 2017. Genetic diversity and APOBEC3-induced mutations were assessed at the viral whole genome and gene levels. Thirty-four pairwise sequence comparisons covering 24 HPV types in cervical and anal swabs revealed identical infections in the two anatomical sites. Differences in genetic diversity among HPV types was observed between patients. The E6 gene was significantly less conserved in low-risk HPVs (lrHPVs) compared to high-risk HPVs (hrHPVs) (p = 0.009). APOBEC3-induced mutations were found to be more common in lrHPVs than in hrHPVs (p = 0.005), supported by our data and by using large HPV sequence collections from the GenBank database. Focusing on the most common lrHPVs 6 and 11 and hrHPVs 16 and 18, APOBEC3-induced mutations were predominantly found in the E4 and E6 genes in lrHPVs, but were almost absent in these genes in hrHPVs. The variable APOBEC3 mutational signatures could contribute to the different oncogenic potentials between HPVs. Further studies are needed to conclusively determine whether APOBEC3 editing levels are associated to the carcinogenic potential of HPVs at the type and sublineage scales.