The present study explores nucleotide variability, phylogeny and association with cervical neoplasia in high risk HPV types 16, 18, 31, 33 and 45 collected from Greek women. Of the 1894 women undergoing routine cervical cytology examination, 160 samples test positive for single infections of HPV type 16 (n = 104), HPV 31 (n = 40), HPV 33 (n = 7), HPV 18 (n = 5), and HPV 45 (n = 4) were typed by microarrays method, amplified by PCR then sequenced and phylogenetically analyzed. For HPV 16, 9 variants with nucleotide variations were included into the study. For HPV 31, 33, 18 and 45, nucleotide variations were identified in 6, 4, 2 and 3 variants, respectively. The Bayesian inference and Maximum Parsimony methods were used in order to construct the phylogenetic trees. When types were analyzed independently HPV 16 (European and non-European) and HPV 18 (African and non-African) formed distinct clades. The genomic characterization of HPV variants will be important for illuminating the geographical relatedness and biological differences and for the determination of their risk.