Human papillomavirus 33 worldwide genetic variation and associated risk of cervical cancer

Human papillomavirus (HPV) 33, a member of the HPV16-related alpha-9 species group, is found in approximately 5% of cervical cancers worldwide. The current study aimed to characterize the genetic diversity of HPV33 and to explore the association of HPV33 variants with the risk for cervical cancer. Taking advantage of the International Agency for Research on Cancer biobank, we sequenced the entire E6 and E7 open reading frames of 213 HPV33-positive cervical samples from 30 countries. We identified 28 HPV33 variants that formed 5 phylogenetic groups: the previously identified A1, A2, and B (sub)lineages and the novel A3 and C (sub)lineages. The A1 sublineage was strongly over-represented in cervical cases compared to controls in both Africa and Europe. In conclusion, we provide a classification system for HPV33 variants based on the sequence of E6 and E7 and suggest that the association of HPV33 with cervical cancer may differ by variant (sub)lineage.