Abstract
Human papillomaviruses (HPVs) are widespread, sexually transmitted group of viruses that infect most individuals at some stage, causing genital warts and cancers. They are members of the Papillomaviridae family, which contains about 400 HPV types. China is among the high HPV burden countries with reported infections of multiple HPV types, accounting for 17.3% of global deaths and 18.2% of global new cases. Thus, understanding the genetic variation and geographic diversity characteristics of HPVs isolated in China is critical for global HPV prevention strategies. Thus, we analyzed the available HPV genome sequences isolated in China that grouped into two categories (alpha- and gamma-papillomaviruses) based on full-length genomes. The most common were HPV-16, -6, -58, and -52 respectively. In addition, four of the novel strains isolated in China, e.g. TG550, JDFY01, CH2, and L55 clustered with the HPV-mSK 159, 244, 201, and 200 respectively. Our phylogeographic network analysis indicated that the L55, TG550, and CH2 are genetically identical to the mSK 200, 046, and 201 respectively, while JDFY01 appeared separately, connected to the mSK-040 following five mutational steps. Also, we found ten recombination events among HPV-6/11 types within their E1, E2, E7, L1/L2 proteins, and Long Control Region ORFs. We achieved the consensus amino acid sequences of HPV proteins and found a conserved stretch of amino acids within E5A of all HPVs circulating in China. These findings offer valued insights into the genetic relationships, distribution, and evolution of the HPVs in China that may assist in adapting effective HPV preventive measures.