Genetic diversity of E6, E7 and the long control region in human papillomavirus type 16 variants in Beijing, China

Background

High-risk human papillomavirus type 16 (HPV16) is a risk factor for cervical cancer. The progression from initial infection to cervical cancer has been linked to properties of the viral sequences. However, the distribution of HPV16 variants among Chinese women has not been extensively addressed and the role of HPV16 variants in the risk of cervical carcinogenesis remains poorly understood.

Methods

HPV16 positive cervical exfoliated cell samples were collected from 249 women living in Beijing, China. PCR products from two fragments of E6-E7 and LCR of HPV16 in these samples were sequenced and analyzed.

Results

Lineage A was found in the subjects, including A1, A2, A3 and A4 sublineages. Based on the HPV16 reference sequences, 26 nucleotide mutations of A4 sublineage and 39 nucleotide mutations of A1-3 sublineages were found in the E6, E7 and LCR of HPV16 isolates. Point mutations T843C, A7287C and A7872G of A4 sublineage were significantly associated with high-grade cervical lesions. The high-frequency sites in HPV16 LCR located at regions that can bind to multiple transcription factors.

Conclusions

This study contributes to the identification of unique variants in the E6, E7 and LCR of HPV16 isolates infected in Chinese women. Mutations of T843C, A7287C and A7872G in A4 sublineages were significantly associated with high-grade cervical lesions, suggesting that mutations in the E7 and LCR region have potential effects on viral replication and progression of cervical cancer.

•

Lineage A is the predominant HPV16 variant in Chinese women.

•

Based on the HPV16 prototype, 25 nucleotide mutations in E6 and E7 were found and 41 nucleotide mutations in LCR were found.

•

The mutations of C843T, A7287C and A7872G were significantly correlated with the cervical lesions.

Genetic variations in the long control region of human papillomavirus type 16 isolates from India: implications for cervical carcinogenesis

Introduction. Infection with high-risk human papillomavirus (HPV) types, specifically HPV type 16 (HPV16), is considered to be the most important risk factor in the development of cervical intraepithelial neoplasia and cancer. The long control region (LCR) is a noncoding region that comprises approximately 10 % of the HPV genome and contains regulatory elements for viral transcription and replication. Sequence variations in LCR may impact on the replication efficiency and oncogenic potential of the virus.Gap statement. Studies documenting variations in LCR of HPV16 isolates pertaining to cervical neoplastic status in India are limited.Aim. The present study was designed to characterize variations in the LCR of Indian isolates of HPV16 and study their association with cervical disease grades.Methodology. The LCR was amplified and sequenced from HPV16 positive cervical samples belonging to different cervical disease grades. Sequences were aligned to identify variations and potential transcription factor binding sites (TFbs) were predicted using the JASPAR database in addition to phylogenetic studies.Results. Among the 163 HPV16 isolates analysed, 47 different nucleotide variations were detected in the LCR, of which 25 are reported for first time in Indian isolates. Point mutations were detected in 35/54 (64.8 %) samples with normal cervical status, 44/50 (88 %) samples with low-grade cervical disease and 53/59 (89.8 %) samples with high-grade cervical disease. Variations T6586C, G6657A and T6850G were significantly associated with high-grade cervical status. Thirteen LCR variations were detected in the binding sites for CEBPB, ETS1, JUN, MYB, NFIL3, PHOX2A and SOX9 transcription factors.Conclusion. The present study helped to identify unique variations in the LCRs of HPV16 Indian isolates. The variations in the A4 sub-lineage were significantly associated with high-grade disease status. The isolates belonging to the A4 and D3 sub-lineages harboured mutations in putative TFbs, implying a potential impact on viral replication and progression to cervical cancer.

Genetic variability of the HPV16 early genes and LCR. Present and future perspectives

Human papillomavirus 16 (HPV16) infection is the aetiologic factor for the development of cervical dysplasia and is regarded as highly carcinogen, because it is implicated in more than 50% of cervical cancer cases, worldwide. The tumourigenic potential of HPV16 has triggered the extensive sequence analysis of viral genome in order to identify nucleotide variations and amino acid substitutions that influence viral oncogenicity and subsequently the initiation and progression of cervical cancer. Nowadays, specific mutations of HPV16 DNA have been associated with an increased risk of high-grade squamous intraepithelial lesions and invasive cervical cancer (ICC) development, including E6: Q14H, H78Y, L83V, Ε7: N29S, S63F, E2: H35Q, P219S, T310K, E5: I65V, whereas highly conserved regions of viral DNA have been extensively characterised. In addition, numerous novel HPV16 mutations are observed among the studied populations from various geographic regions, hence advocating that different HPV16 strains seem to emerge with different tumourigenic capacities. The present review focuses on the variability of the early genes and the long control region, emphasising on the association of specific mutations with the development of severe dysplasia. Finally, it evaluates whether specific regions of HPV16 DNA are able to serve as valuable biomarkers for cervical cancer risk.

Genetic variation in the E6 and E7 genes of human papillomavirus type 16 in northeastern Argentina

The province of Misiones is considered a region with a high mortality rate due to cervical cancer (CC). To gain insight into this problem, we explored the association between genetic variation in the E6 and E7 oncogenes of HPV16 and the risk of CC. We studied 160 women with cytological diagnoses of negative for intraepithelial lesion or malignity, low-grade squamous intraepithelial lesion, and high-grade squamous intraepithelial lesion/CC and a positive test for HPV16 infection. The genetic characterization of E6 and E7 genes was undertaken through PCR amplification and direct Sanger sequencing. Phylogenetic classification was conducted using Bayesian methods. To estimate the odds ratio (OR) for an association between genetic variants in the E6 and E7 genes and the risk of CC, we used ordinal logistic regression adjusted by age. The final data set comprised 112 samples. Diagnostic single-nucleotide polymorphisms (SNPs) and phylogenetic trees confirmed the presence of Lineage A (95.5%) and D (4.5%) in the samples. For the E6 gene, we identified eleven different sequences, with the most common ones being Lineage A E6 350G (58.9%) and E6 350T (37.5%). The E6 350G was associated with progression to HSIL/CC, with an OR of 19.41 (4.95-76.10). The E7 gene was more conserved than E6, probably due to the functional constraints of this small protein. Our results confirmed the association of the E6 350G SNP with a higher risk of developing CC. These data will contribute to understanding the biological bases of CC incidence in this region.

Association of Human Papillomavirus Type 16 Long Control Region Variations with Cervical Cancer in a Han Chinese Population

Objective: High-risk human papillomavirus (HPV) E6 and E7 proteins are the major oncoproteins involved in the tumorigenesis of cervical cancer. The long control region (LCR) in HPV plays an important role in regulating the expression of the E6 and E7 oncogenes. In the current study, we investigated the association of HPV16 LCR variations with cervical cancer. Methods: A total of 139 HPV16-positive cervical cancer patients (case group) and 116 HPV16-positive asymptomatic individuals (control group) were enrolled in the current study. Then, the HPV16 LCR was sequenced to determine the association between LCR variations and cervical cancer. Results: In the current study, HPV16 A1-A3 (19.4%), A4 (78.4%) and D3 (2.2%) variants were found in the case group. However, only A1-A3 (34.5%) and A4 variants (65.5%) were found in the control group. The distribution of the HPV16 variants between the case and control groups was significantly different (P=0.009). Moreover, a total of eleven variations (A7167G, A7173C, C7176T, C7200T, T7269C, C7286A, C7729A, C7763T, A7841G, G7867A and T24C) were significantly different between the case and control groups (P<0.05). For the sub-lineage analysis, only C7873G variations were significantly different between the case and control groups in the A4 (As) variant (P=0.039). Conclusion: Our results showed that specific variations in the HPV16 LCR were associated with cervical cancer. Our study will provide a good reference for further understanding of the relationship between HPV16 LCR variation and cervical cancer.

Mutational landscape and intra-host diversity of human papillomavirus type 16 long control region and E6 variants in cervical samples

Human papillomavirus genotype 16 (HPV16) is the most frequent high-risk HPV (HR-HPV) identified in cervical precursor lesions and cervical cancer (CC) worldwide. The oncogenic potential of HPV16 is partly dependent on the lineage involved in the infection and the presence of clinically relevant mutations. In this report, we present the distribution of HR-HPV and the mutational profile and intra-host variability of HPV16 lineages, based on analysis of the long control region (LCR) and the E6 gene in samples with normal cytology (n = 39), squamous intraepithelial lesions (n = 25), and CC (n = 39). HR-HPV genotyping was performed using multiplex real-time PCR. HPV16 lineage assignments and mutation frequencies were determined by conventional PCR and Sanger DNA sequencing, and intra-patient viral populations were analyzed using next-generation sequencing (NGS). The most frequent HR-HPV type was HPV16, followed by HPV31 and HPV18. The frequency of HPV16 sublineages was A1/A2 > D2 > D3 and B1. Moreover, the most frequent mutations, both in samples from this study and in the available sequences from Mexican isolates in the GenBank database were LCR-G7518A, which is involved in carcinogenesis, and E6-T350G (producing L83V), associated with persistence of infection. Otherwise, deep sequencing revealed high conservation of viral lineages and mutations, independently of the stages studied. In conclusion, the high frequency and stability of these molecular markers, as well as the circulating viral lineages, could be related to the incidence of CC associated with HPV16. Hence, they deserve a broader analysis to determine the risk of specific populations for progression of the disease.

Distribution of HPV 16 E6 gene variants in screening women and its associations with cervical lesions progression

The aim of this study is to investigate distribution of human papillomavirus (HPV) 16 variants in screening healthy women and the potential association between HPV 16 variants and progression of cervical lesions. For this study a total of 2000 healthy women in Tianjin urban area and 212 patients who were HPV 16 positive and underwent colposcopy were analyzed for HPV 16 variants by pyrosequencing. The results show that the HPV 16 was the most prevalent genotype in Tianjin healthy women and five HPV 16 variant types were detected. The HPV 16 variants were determined by sequencing parital E6 region and the detected variants were European prototype E-T350 (E-p), E-G350, E-C109 G, Asian (As) and Asian-American (AA), among which the E-p variant was the most prevalent (82.76%) followed by As variant. Interestingly, in patients with suspected cervical lesions the most prevalent variant was As variant (54.9%) by increasing significance with severity of cervical diseases (OR 4.337; 95% CI 1.248-15.067; P = 0.021), and followed by HPV 16 E-p variant while E-G350 variant only appeared in HSIL and cervical cancer. Our results show that HPV 16 E-p variant was more prevalent than As in Tianjin healthy screening women while As variant was the most frequently type in HSIL and cervical cancer. It is suggested that the mutation of HPV 16 Asian variants, comparing with HPV 16 E-p variants, might contribute to the transformation from HPV 16 persistent infection to cervical cancer.

Analysis of Nucleotide Alterations in the E6 Genomic Region of Human Papillomavirus Types 6 and 11 in Condyloma Acuminatum Samples from Brazil

Condyloma acuminata (CA), or genital warts, are benign proliferative epidermal or mucous lesions that are caused by infection with human papillomavirus (HPV), mainly the low-risk types 6 and 11. HPV variants are defined as viral sequences that share identity in the nucleotide sequence of the L1 gene greater than 98%. Based on this criterion, HPV6 and 11 variant lineages have been studied, and there are ongoing attempts to correlate these genetic variants with different clinical findings of infection. Therefore, the aims of this study were to detect variants and nucleotide alterations present in the E6 regions of HPV types 6 and 11 found in CA samples, to correlate the HPV presence with the clinical-pathological data of the patients, and to determine phylogenetic relationships with variants from other places in the world. The E6 regions of 25 HPV6 samples and 7 HPV11 samples from CA were amplified using PCR with specific primers. The products were ligated to a cloning vector and five colonies of each sample were sequenced to observe the nucleotide alterations. Twelve samples were identified as the HPV6B3 variant, presenting the mutation (guanine) G474A (adenine), and one of them also showed the mutation (thymine) T369G. The other 13 patients were positive for HPV6B1 without nucleotide alterations. In the analysis of the HPV11 samples, all patients showed the mutations T137C and (cytosine) C380T. One patient also presented the nucleotide alteration T410C. None of the mutations found in the 32 analyzed samples resulted in amino acid changes. Patient age, local occurrence, and HIV infection did not show significant association with HPV infection. Besides, the data found in this study did not show a relationship with the geographical region of isolation when compared to other data from different regions of the world. In this way, despite the nucleotide alterations found, it was not possible to observe amino acid changes and variants grouping according to geographical region.

Genetic variations in E6, E7 and the long control region of human papillomavirus type 16 among patients with cervical lesions in Xinjiang, China

Background

Xinjiang is one of the areas with the highest incidence of cervical cancer in China. Genetic variation in Human papillomavirus type 16 (HPV16) may increase the ability of the virus to mediate carcinogenesis and immune escape, which are risk factors for the progression of cervical cancer. We investigated polymorphism in HPV16 and the distribution of its sub-lineages in the region by analyzing the E6, E7 and long control region (LCR) gene sequences from women with HPV16-positive cervical samples in Xinjiang.

Methods

A total of 138 cases of cervical lesions and squamous cell carcinoma with infection of HPV16 virus were collected. The E6 and E7 genes and LCR of HPV16 virus were sequenced and compared with the HPV16 European prototype reference and other HPV16 mutants for single nucleotide polymorphisms. Neighbor-joining phylogenetic trees were constructed using E6, E7 and LCR sequences.

Results

Fourteen missense mutations were found in the E6 gene; the loci with the highest mutation frequency were T350G (36/75, 48%) and T178G (19/75, 25.3%). In the E7 gene, the locus with the highest mutation frequency was A647G (18/75, 24%). A total of 33 polymorphic sites were found in the LCR, of which T7447C (39/95, 40.1%) was the most frequent.

Conclusion

HPV16 in Xinjiang is mainly of the European variant, followed by the Asian variant type; no Africa 1, 2 or Asia–America variant types were found.

The polymorphisms of LCR, E6, and E7 of HPV-58 isolates in Yunnan, Southwest China

Backgroud

Variations in HPV LCR/E6/E7 have been shown to be associated with the viral persistence and cervical cancer development. So far, there are few reports about the polymorphisms of the HPV-58 LCR/E6/E7 sequences in Southwest China. This study aims to characterize the gene polymorphisms of the HPV-58 LCR/E6/E7 sequences in women of Southwest China, and assess the effects of variations on the immune recognition of viral E6 and E7 antigens.

Methods

Twelve LCR/E6/E7 of the HPV-58 isolates were amplified and sequenced. A neighbor-joining phylogenetic tree was constructed by MEGA 7.0, followed by the secondary structure prediction of the related proteins using PSIPRED v3.3. The selection pressure acting on the HPV-58 E6 and E7 coding regions was estimated by Bayes empirical Bayes analysis of PAML 4.8. Meanwhile, the MHC class-I and II binding peptides were predicted by the ProPred-I server and ProPred server. The transcription factor binding sites in the HPV-58 LCR were analyzed using the JASPAR database.

Results

Twenty nine SNPs (20 in the LCR, 3 in the E6, 6 in the E7) were identified at 27 nucleotide sites across the HPV-58 LCR/E6/E7. From the most variable to the least variable, the nucleotide variations were LCR > E7 > E6. The combinations of all the SNPs resulted in 11 unique sequences, which were clustered into the A lineage (7 belong to A1, 2 belong to A2, and 2 belong to A3). An insertion (TGTCAGTTTCCT) was found between the nucleotide sites 7280 and 7281 in 2 variants, and a deletion (TTTAT) was found between 7429 and 7433 in 1 variant. The most common non-synonymous substitution V77A in the E7 was observed in the sequences encoding the α-helix. 63G in the E7 was determined to be the only one positively selected site in the HPV-58 E6/E7 sequences. Six non-synonymous amino acid substitutions (including S71F and K93 N in the E6, and T20I, G41R, G63S/D, and V77A in the E7) were affecting multiple putative epitopes for both CD4⁺ and CD8⁺ T-cells. In the LCR, C7265G and C7266T were the most variable sites and were the potential binding sites for the transcription factor SOX10.

Conclusion

These results provide an insight into the intrinsic geographical relatedness and biological differences of the HPV-58 variants, and contribute to further research on the HPV-58 epidemiology, carcinogenesis, and therapeutic vaccine development.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-0986-7) contains supplementary material, which is available to authorized users.

Genetic variability and functional implication of the long control region in HPV-16 variants in Southwest China

HPV-16 long control region (LCR) has been shown to be the most variable region of the HPV-16 genome and may play important roles in viral persistence and the development of cervical cancer. This study aimed to assess the risk of HPV-16 LCR variants for cervical cancer in women of Southwest China. 2146 cervical scrapings of volunteer outpatients and 74 cervical cancer tissues were screened.14 entire HPV-16 LCRs from asymptomatic carriers and 34 entire HPV-16 LCRs from cervical cancer patients were successfully amplified and sequenced to align to others described. 58 different point mutations were detected in 54 nucleotide sites of HPV-16 LCR. G7193T and G7521A variants, accounting for 100% of the infections, were predicted to locate at the binding site for FOXA1 and SOX9, respectively. A7730C variant which showed a high mutation frequency in cervical cancer was predicted to be a binding site for the cellular transcription factor PHOX2A. In addition, phylogenetic analysis displayed a high prevalence of A lineage in HPV-16 LCR in this Southwest China population. This study may help understanding of the intrinsic geographical relatedness and the correlations between LCR mutations and the development of carcinogenic lesions in Southwest China population. And it provides useful data for the further study of the biological function of HPV-16 LCR variants.

The role of interleukin 10 in human papilloma virus infection and progression to cervical carcinoma

Although Human Papillomavirus (HPV) exerts a vital influence on cervical carcinogenesis, other factors influence the development of a squamous intraepithelial lesion (SIL) that may or not progress to cervical cancer. Among several cytokines, Interleukin 10 (IL-10) stands out as an important anti-inflammatory factor, leading to immune system evasion through an immunosuppressive state. In the cervical microenvironment, during different stages of HPV infection, IL-10 production can be induced and maintained by different cell sources, including infected keratinocytes, some subsets of dendritic cells (DC), tumor associated macrophages (TAM), T regulatory cells (Treg) and tumor cells. Further, a wide range of effects can be exerted by IL-10 on different cell populations, such as inhibiting proinflammatory cytokine production, DCs differentiation, antigen presenting function and T-helper 1 (Th1) polarization. IL-10 is one of several cytokines involved in cancer development and sustenance, although its role in cancer is still controversial and poorly understood. However, cervical IL-10 levels tend to increase in parallel to SIL development and are even higher within cervical tumors. Accumulating data have shown that after HPV infection, IL-10 levels are enhanced as a result of HPV E2, E6 and E7 proteins action over IL-10 gene transcription, while IL-10 stimulates HPV E6 and E7 expression. Therefore, this interplay between HPV and IL-10 creates a vicious cycle that could favor an immunosuppressive microenvironment in the cervix, facilitating the progression of a simple HPV infection to SIL or cervical cancer.

p53, Cyclin D1, p21 (WAF1) and Ki-67 (MIB1) Expression at Invasive Tumour Fronts of Oral Squamous Cell Carcinomas and Development of Local Recurrence

BACKGROUND

Expression of p53, cyclin D1, p21 (WAF1) and Ki-67 (MIB1) was evaluated in oral squamous cell carcinoma (OSCC) to test whether levels of these markers at invasive tumour fronts (ITFs) could predict the development of local recurrence.

MATERIALS AND METHODS

Archived paraffin-embedded specimens from 51 patients with T1/T2 tumours were stained immunohistochemically and analysed quantitatively. Local recurrence-free survival was tested with Kaplan-Meier survival plots (log-rank test) using median values to define low and high expression groups and with a Cox's proportional hazards model in which the expression scores were entered as continuous variables.

RESULTS

The assessment of expression of all markers was highly reliable, univariate analysis showing that patients with clear surgical margins, with low cyclin D1 and high p21 expression at the ITF had the best local recurrence-free survival. Multivariate analysis showed that these three parameters were independent prognostic factors but that neither p53 nor MIB1 expression were of prognostic value.

CONCLUSIONS

Assessment of p53, cyclin D1, p21 (WAF1), and Ki-67 (MIB1) at the ITF could help to predict local recurrence in early stage oral squamous cell carcinoma cases.

Variation of HPV Subtypes with Focus on HPV-Infection and Cancer in the Head and Neck Region

The human papillomavirus (HPV) comprises a heterogeneous group of double-strand DNA viruses with variable potential to infect human epithelial cells and trigger neoplastic transformation. Its 8 kb genome encodes proteins required for virus replication and self-organized formation of infectious particles but also for early proteins E6 and E7 able to trigger neoplastic transformation. E6 and E7 of high-risk (HR) HPV subtypes can bind to p53 or release E2F and abrogate replication control. Due to variable amino acid sequence (AAS) in the binding sites of E6 and E7 particular HR-HPV variants within subtypes are essentially heterogeneous in efficacy triggering neoplastic transformation and cancer development. This could explain differences in the clinical course of HPV-driven head and neck cancer.

Molecular analysis of human Papillomavirus detected among women positive for cervical lesions by visual inspection with acetic acid/Lugol's iodine (VIA/VILI) in Libreville, Gabon

Background

The human papillomavirus (HPV) is the causative agent of cervical cancer, which is the leading cancer-related cause of death for women in Sub-Saharan Africa. In 2013, the Gabonese Ministry of Health and the Sylvia Bongo Ondimba Foundation implemented cervical cancer screening programs based on the detection of cancerous lesions by visual inspection with acetic acid and/or Lugol’s iodine (VIA/VILI). This pilot study was set up to determine the HPV profile and analyze the nucleotide sequence variation of HPV16 circulating in patients with cervical abnormalities detected by VIA/VILI testing.

Methods

The cervical abnormalities observed upon VIA/VILI were confirmed by liquid-based cytology for all tested women. Nested PCR using the MY09/11 and GP5+/6+ primer sets was used to detect HPVs present in the extracted DNA. HPV genotypes were determined after sequencing of amplicons based on a high-throughput sequencing approach. For isolates of the HPV16 genotype, the E6 gene and the long control region (LCR) were directly sequenced using Sanger method.

Results

The study included 87 women who showed a positive VIA/VILI result. Cervical abnormalities were found in 40.23 % of women and 40 % were classified as high-grade lesions. The HPV detection rate was 82.9 % among women with abnormal cytology. Among all the identified high-risk HPV genotypes, HPV16, 18 and 33 were the most frequent. Multiple HPV infections were observed in 42.31 % of HPV-infected women. Analysis of the HPV16 sequence variation in the E6 gene and in the LCR showed that 85.3 and 14.7 % belonged to the African and European lineages, respectively. Among the African branch variants, Af2 was the most frequently identified in this study.

Conclusion

This study offers the first report of the HPV detection rate and molecular epidemiology among Gabonese women with a positive result in a VIA/VILI screening test. Moreover, data on the HPV16 sequence variation confirm the predominance of African variants in high-grade lesions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13027-016-0098-1) contains supplementary material, which is available to authorized users.

Characterization of Intra-Type Variants of Oncogenic Human Papillomaviruses by Next-Generation Deep Sequencing of the E6/E7 Region

Different human papillomavirus (HPV) types are characterized by differences in tissue tropism and ability to promote cell proliferation and transformation. In addition, clinical and experimental studies have shown that some genetic variants/lineages of high-risk HPV (HR-HPV) types are characterized by increased oncogenic activity and probability to induce cancer. In this study, we designed and validated a new method based on multiplex PCR-deep sequencing of the E6/E7 region of HR-HPV types to characterize HPV intra-type variants in clinical specimens. Validation experiments demonstrated that this method allowed reliable identification of the different lineages of oncogenic HPV types. Advantages of this method over other published methods were represented by its ability to detect variants of all HR-HPV types in a single reaction, to detect variants of HR-HPV types in clinical specimens with multiple infections, and, being based on sequencing of the full E6/E7 region, to detect amino acid changes in these oncogenes potentially associated with increased transforming activity.

Changes in global gene expression profiles induced by HPV 16 E6 oncoprotein variants in cervical carcinoma C33-A cells

We analyzed the effects of the expression of HPV 16 E6 oncoprotein variants (AA-a, AA-c, E-A176/G350, E-C188/G350, E-G350), and the E-Prototype in global gene expression profiles in an in vitro model. E6 gene was cloned into an expression vector fused to GFP and was transfected in C33-A cells. Affymetrix GeneChip Human Transcriptome Array 2.0 platform was used to analyze the expression of over 245,000 coding transcripts. We found that HPV16 E6 variants altered the expression of 387 different genes in comparison with E-Prototype. The altered genes are involved in cellular processes related to the development of cervical carcinoma, such as adhesion, angiogenesis, apoptosis, differentiation, cell cycle, proliferation, transcription and protein translation. Our results show that polymorphic changes in HPV16 E6 natural variants are sufficient to alter the overall gene expression profile in C33-A cells, explaining in part the observed differences in oncogenic potential of HPV16 variants.