Human papillomaviruses (HPVs) co-evolve slowly with the human host and each HPV genotype displays epithelial tropisms. We assessed the evolution of intra HPV genotype variants within samples, and their association to anogenital site, cervical cytology and HIV status. Variability in the L1 gene of 35 HPV genotypes was characterized phylogenetically using maximum likelihood, and portrayed by phenotype. Up to a thousand unique variants were identified within individual samples. In-depth analyses of the most prevalent genotypes, HPV16, HPV18 and HPV52, revealed that the high diversity was dominated by a few abundant variants. This suggests high intra-host mutation rates. Clades of HPV16, HPV18 and HPV52 were associated to anatomical site and HIV co-infection. Particularly, we observed that one HPV16 clade was specific to vaginal cells and one HPV52 clade was specific to anal cells. One major HPV52 clade, present in several samples, was strongly associated with cervical neoplasia. Overall, our data suggest that tissue tropism and HIV immunosuppression are strong shapers of HPV evolution. Keywords: HPV variability, HPV phylogenetics, Anogenital, Tissue tropism, HIV