It has been suggested that DNA sequence variants of HPV16 contribute to differences in the behavior of individual cervical lesions. To address this question, we have analyzed the association of HPV16 variants with diagnostic severity in 354 HPV16-positive Oklahoman women. HPV16 variant status was determined by PCR amplification and DNA sequencing of the E6 open reading frame. European sequences were identified in 86% of samples and 14% were non-European. Of the 51 non-European cases, 61% were Asian-American, 23% African and 16% were Native American variants. European prototype and related variants were present in comparable numbers (43% each) but the relative proportion of each differed with diagnostic category. In general, the proportion of European variants and non-European variants increased with diagnostic severity while the European prototype decreased. When adjusted for age and race (white, black or Hispanic), the increased risk for carcinoma/severe dysplasia for non-European variants was statistically significant with an odds ratio of 3.8 (1.3-10.7). However, the analogous comparison for the European variants, although also showing increased association with carcinoma/severe dysplasia, did not reach statistical significance (OR = 1.6 (95% CI 0.7-3.6). Overall, HPV16 European sequences (both prototype and related variants), were predominant in Oklahoman women including those with cancers. This suggests that while there appear to be differences among the HPV16-variant categories in risk for progression to invasive cancer, all variant categories are associated with the development of invasive cancer.