Genetic variations of E6 and long control region of human papillomavirus type 16 from patients with cervical lesion in Liaoning, China

Background High-risk human papillomavirus type 16 (HPV16) is a risk factor for cervical cancer. Previous studies suggest that polymorphisms in the E6 gene or the long control region(LCR)of HPV16 may alter the oncogenic potential of the virus. The aims of this study were to investigate the genetic variations of HPV16 E6 gene and LCR in isolates from Chinese population and correlation of the E6 and LCR polymorphisms with disease status of infected patients. Methods HPV16 positive endocervical specimens were collected from 304 women living in Northeast of China. Sequences of E6 gene and LCR were analyzed by PCR-sequencing. Results Two lineages were found in the populations, including EUR lineage and As lineage. Based on the HPV16 prototype, the most frequent variation in the E6 gene was T178A/G (48.7%), followed by mutations of G94A (12.2%) and T350G (9.9%). The rank orders of incidence of E6 variations in amino acid were as follows: D25E (46.3%), L83V (9.9%) and H78Y (4.3%). Nucleotide variations in LCR were found in all the 304 isolates from HPV16 positive cervical samples. The most commonly observed LCR variations were the transition replacement G7193T, 7434CIns, G7521A and 7863ADel (100%). The As lineage was associated with HPV persistent infections and with disease status of ≥CIN2,3. The EUR lineage variants showed a negative trend of association with the severity of ≥CIN2,3. Among 41 variations found in LCR, 25 (61.0%) were located at the binding sites for transcription factors. Occurrence of ≥CIN2,3 was significantly associated with the mutations of R10G/L83V in E6 and the C7294T co-variation in LCR, after adjusting for ages of infected patients. Conclusions Associations between As lineage and HPV persistent infections, and with disease status of ≥CIN2,3, and an association between the EUR lineage and negative trend of association with the severity of ≥CIN2,3 were found in this study. An association between a co-variation of R10G/L83V in E6 and C7294T in LCR and an increased risk for developing CIN-2,3 was found in a HPV16 infected population of Chinese women. These findings indicate that HPV16 polymorphism influences development of CIN-2,3.