Prevalence of oncogenic Human papillomavirus and genetic diversity in the L1 gene of HPV16 HPV 18 HPV31 and HPV33 found in women from Vojvodina Province Serbia

The main purpose of this paper is to estimate the pre-vaccination prevalence of 12 hrHPV types among 564 women from Vojvodina province (Serbia). The corrected contingency coefficient (Ccorr) was used to estimate the importance of association of examined HPV types and cytological diagnosis. The highest association with the abnormal cytology was observed for HPV 16 (Ccorr = 0.493) in all age groups of participants. The effect of HPV 16 was especially clear within the group of women older than 35 years (Ccorr = 0.691), compared with women younger than 35 (Ccorr = 0.333). The molecular characterization at the level of L1 gene of HPV 16, 18, 31 and 33 variants was for the first time assessed in our region. Nearly all HPV 16 isolates cluster with variant lineage A (96.4%) the remaining isolates clustering with variant lineage D. All of HPV 18 and HPV 33 isolates are clustering within the lineage A while isolates of HPV 31 group with lineages A and C. This contributes to understanding of intrinsic geographical and biological differences of examined HPV types and could be useful for development of cervical cancer screening strategies in Vojvodina (Serbia) and diagnosis of HPV related cervical cancer in general.