The human papillomavirus 16 European-T350G E6 variant can immortalize but not transform keratinocytes in the absence of E7

HPV16 Phylogenetic Variants in Anogenital and Head and Neck Cancers: State of the Art and Perspectives

HPV16 is responsible for approximately 60% and 90% of global HPV-induced cervical and oropharyngeal cancers, respectively. HPV16 intratype variants have been identified by HPV genome sequencing and classified into four phylogenetic lineages (A-D). Our understanding of HPV16 variants mostly derives from epidemiological studies on cervical cancer (CC) in which HPV16 B, C, and D lineages (previously named "non-European" variants) were mainly associated with high-grade cervical lesions and cancer. Although a predominance of HPV16 lineage A (previously named "European variants") has been observed in head and neck squamous cell carcinoma (HNSCC), epidemiological and in vitro biological studies are still limited for this tumor site. Next Generation Sequencing (NGS) of the entire HPV genome has deepened our knowledge of the prevalence and distribution of HPV variants in CC and HNSCC. Research on cervical cancer has shown that certain HPV16 sublineages, such as D2, D3, A3, and A4, are associated with an increased risk of cervical cancer, and sublineages A4, D2, and D3 are linked to a higher risk of developing adenocarcinomas. Additionally, lineage C and sublineages D2 or D3 of HPV16 show an elevated risk of developing premalignant cervical lesions. However, it is still crucial to conduct large-scale studies on HPV16 variants in different HPV-related tumor sites to deeply evaluate their association with disease development and outcomes. This review discusses the current knowledge and updates on HPV16 phylogenetic variants distribution in HPV-driven anogenital and head and neck cancers.

HPV16 Intratypic Variants in Head and Neck Cancers: A North American Perspective

Human papillomavirus (HPV) is the major causative agent for cervical and many head and neck cancers (HNCs). HPVs randomly acquire single nucleotide polymorphisms (SNPs) that may become established via positive selection. Within an HPV type, viral isolates differing by <2% in the L1 region are termed "variants" and classified based on combinations of SNPs. Studies in cervical cancer demonstrate clear differences between HPV16 intratypic variants in terms of persistence of infection, tumor histology, cancer risk, and death. Much less is known about the frequency of HPV16 variants in HNC, and their effects on clinical outcomes. We combined HPV16 positive (HPV16+) HNC samples from a local Southwestern Ontario, Canada cohort with those from the Cancer Genome Atlas to create a larger North American cohort of 149 cases with clinical data and determined the distribution of intratypic variants and their impact on clinical outcomes. Most isolates were lineage A, sublineage A1, or A2, with roughly half exhibiting the T350G polymorphism in E6. Univariable analysis identified significant differences between 350T and 350G intratypic variants in clinical T, N, and O staging, as well as disease-free survival. Multivariable analysis failed to identify any clinical factor as a statistically significant covariate for disease-free survival differences between 350T and 350G. Significant differences in several measures of B-cell mediated immune response were also observed between 350T and 350G intratypic variants. We suggest that HPV genetic variation may be associated with HNC clinical characteristics and may have prognostic value.

The analysis of HPV integration sites based on nanopore sequencing and the profiling changes along the course of photodynamic therapy

OBJECTIVE

To detect the HPV genotype and integration sites in patients with high-risk HPV infection at different stages of photodynamic therapy using nanopore technology and to evaluate the treatment effect.

METHODS

Four patients with HPV infection were selected and subjected to photodynamic therapy, and cervical exfoliated cell was sampled at before treatment, after three courses of treatment and six courses of treatment, their viral abundance and insertion sites were analyzed by nanopore technology, and pathological examinations were performed before and after treatment. In this study, we developed a novel assay that combined viral sequence enrichment and Nanopore sequencing for identification of HPV genotype and integration sites at once. The assay has obvious advantages over qPCR or NGS-based methods, as it has better sensitivity after viral sequences enrichment and can generate long-reads (kb to Mb) for better detection rate of structure variations, moreover, fast turn-around time for real-time viral sequencing and analysis.

RESULTS

The pathological grade was reduced in all four patients after photodynamic therapy. Virus has been cleared in two cases after treatment, the virus amount reduced after treatment but not completely cleared in one case, and two type viruses were cleared and one type virus persisted after treatment in the last patient with multiple infection. Viral abundance and the number of integration sites were positively correlated. Gene enrichment analysis showed complete viral clearance in 1 patient and 3 patients required follow-up.

CONCLUSION

Nanopore sequencing can effectively monitor the abundance of HPV viruses and integration sites to show the presence status of viruses, and combined with the results of gene enrichment analysis, the treatment effect can be dynamically assessed.

The T350G Variation of Human Papillomavirus 16 E6 Gene Prevails in Oropharyngeal Cancer from a Small Cohort of Greek Patients

Recent trends have shown a dramatic rise in the incidence of oropharyngeal squamous cell carcinoma strongly associated with high-risk human papillomavirus (HPV) of type 16. The genetic variability of HPV16 has been extensively studied in cervical cancer but there are very limited published data concerning the genetic variations of this HPV type in oropharyngeal cancer. In the present study, the genetic variations of HPV16 E6 gene sequences originated from a small cohort of Greek patients diagnosed with oropharyngeal cancer were assessed. The vast majority of the sequences clustered within the European variant branch. The T350G variation was found to be the predominant one. This finding may indicate the need for further studies that could explain the possible impact of this variant in the pathomechanisms of oropharyngeal 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.

AutoVEM2: A flexible automated tool to analyze candidate key mutations and epidemic trends for virus

In our previous work, we developed an automated tool, AutoVEM, for real-time monitoring the candidate key mutations and epidemic trends of SARS-CoV-2. In this research, we further developed AutoVEM into AutoVEM2. AutoVEM2 is composed of three modules, including call module, analysis module, and plot module, which can be used modularly or as a whole for any virus, as long as the corresponding reference genome is provided. Therefore, it's much more flexible than AutoVEM. Here, we analyzed three existing viruses by AutoVEM2, including SARS-CoV-2, HBV and HPV-16, to show the functions, effectiveness and flexibility of AutoVEM2. We found that the N501Y locus was almost completely linked to the other 16 loci in SARS-CoV-2 genomes from the UK and Europe. Among the 17 loci, 5 loci were on the S protein and all of the five mutations cause amino acid changes, which may influence the epidemic traits of SARS-CoV-2. And some candidate key mutations of HBV and HPV-16, including T350G of HPV-16 and C659T of HBV, were detected. In brief, we developed a flexible automated tool to analyze candidate key mutations and epidemic trends for any virus, which would become a standard process for virus analysis based on genome sequences in the future.

Virus-Host Protein-Protein Interactions between Human Papillomavirus 16 E6 A1 and D2/D3 Sub-Lineages: Variances and Similarities

High-risk strains of human papillomavirus are causative agents for cervical and other mucosal cancers, with type 16 being the most frequent. Compared to the European Prototype (EP; A1), the Asian-American (AA; D2/D3) sub-lineage seems to have increased abilities to promote carcinogenesis. Here, we studied protein–protein interactions (PPIs) between host proteins and sub-lineages of the key transforming E6 protein. We transduced human keratinocyte with EP or AA E6 genes and co-immunoprecipitated E6 proteins along with interacting cellular proteins to detect virus–host binding partners. AAE6 and EPE6 may have unique PPIs with host cellular proteins, conferring gain or loss of function and resulting in varied abilities to promote carcinogenesis. Using liquid chromatography-mass spectrometry and stringent interactor selection criteria based on the number of peptides, we identified 25 candidates: 6 unique to AAE6 and EPE6, along with 13 E6 targets common to both. A novel approach based on pathway selection discovered 171 target proteins: 90 unique AAE6 and 61 unique EPE6 along with 20 common E6 targets. Interpretations were made using databases, such as UniProt, BioGRID, and Reactome. Detected E6 targets were differentially implicated in important hallmarks of cancer: deregulating Notch signaling, energetics and hypoxia, DNA replication and repair, and immune response.

High prevalence of human papillomavirus and European variants of HPV 16 infecting concomitantly to cervix and oral cavity in HIV positive women

OBJECTIVE

Identify the prevalence of HPV infections in the uterine cervix and oral cavity and HPV16 variants in HIV+ women.

METHODS

A total of 174 HIV+ women attended an HIV+ specialized clinic in Mexico City. Cells were obtained from the oral cavity and cervix to extract DNA. Polymerase chain reaction (PCR) was used to amplify the HPV sequence with generic primers. We detected specific HPV types using the INNO-LiPA HPV Genotyping Extra II Kit (INNOGENETICS). The identification of variants was studied by sequencing the E6 gene with a Big Dye Terminator Kit and an Applied Biosystems 3500/3500xL genetic analyzer.

RESULTS

HPV infection was very high in the uterine cervix (168/174, 96.6%) and oral cavity (161/174, 92.5%). The prevalence of HPV concurrent infections in the cervix and oral cavity was 155/174 (89.1%). We found hrHPVs to be more prevalent than low-risk HPVs (lrHPVs) in the oral cavity (90.2% versus 45.4%) and that infections simultaneously affected the cervix (94.3% versus 36.2%) and oral cavity (85.1% versus 20.1%). Surprisingly, only European variants of HPV type 16 were found in the uterine cervix of women and the oral cavity of all tested samples (52 oral cavity samples and 52 uterine cervix samples).

CONCLUSIONS

The high prevalence of HPV, multiple infections and presence of the EP350G intravariant in both anatomical regions are strongly related to the persistence of the virus, which is fundamental for the development of cancer. Therefore, it is very important to control and monitor this high-risk population as well as implement programs for the early detection of HPV and vaccination.

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.

An epithelial organoid model with Langerhans cells for assessing virus-host interactions

Persistent infection with oncogenic human papillomavirus (HPV) may lead to cancer in mucosal and skin tissue. Consequently, HPV must have developed strategies to escape host immune surveillance. Nevertheless, most HPV infections are cleared by the infected host. Our laboratory investigates Langerhans cells (LCs), acting at the interface between innate and adaptive immunity. We hypothesize that this first line of defence is vital for potential HPV elimination. As an alternative to animal models, we use smaller-scale epithelial organoids grown from human primary keratinocytes derived from various anatomical sites. This approach is amenable to large sample sizes-an essential aspect for scientific rigour and statistical power. To evaluate LCs phenotypically and molecularly during the viral life cycle and onset of carcinogenesis, we have included an engineered myeloid cell line with the ability to acquire an LC phenotype. This model is accurately tailored for the crucial time-window of early virus elimination in a complex organism and will shed more light on our long-standing research question of how naturally occurring HPV variants influence disease development. It may also be applied to other microorganism-host interaction research or enquiries of epithelium immunobiology. Finally, our continuously updated pathogen-host analysis tool enables state-of-the-art bioinformatics analyses of next-generation sequencing data. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Oncogenicitiy Comparison of Human Papillomavirus Type 52 E6 Variants

Human papillomavirus (HPV) infection contributes to virtually all cases of cervical cancer, the fourth most common cancer affecting women worldwide. The oncogenicity of HPV is mainly attributable to the E6 and E7 oncoproteins. HPV-52 is the seventh most common HPV type globally, but it has a remarkably high prevalence in East Asia. In previous studies it has been speculated that the oncogenicity might vary among different HPV-52 variants. In the present study, we compared the oncogenicity of E6 derived from the HPV-52 prototype and three commonly found variants, V1 (K93R), V2 (E14D/V92L) and V3 (K93R/N122K), through molecular and phenotypic approaches. We demonstrated that cells containing V1 achieved higher colony formation and showed greater cell migration ability when compared to other variants, but no difference in cell immortalization ability was observed. At the molecular level, the three variants formed complexes with E6-associated protein (E6AP) and p53 as efficiently as the prototype. They degraded p53 and PSD95/Dlg/ZO-1(PDZ) proteins, including MAGI-1c and Dlg, to a similar extent. They also exhibited a similar subcellular localization, and shared a half-life of approximately 45 min. Our findings provide a clearer picture of HPV-52 E6 variant oncogenicity, which is important for further studies aiming to understand the unusually high prevalence of HPV-52 among cervical cancers in East Asia.

Identification of human papillomavirus type 16 variants circulating in the Calabria region by sequencing and phylogenetic analysis of HPV16 from cervical smears

Sequence analysis of HPV16 isolates reveals the presence of genome variants with characteristic mutations. The HPV16 variants have different geographical distribution and diverge into four phylogenetic lineages (A, B, C and D) and 16 sub-lineages: A1, A2, A3 (previously known as European variants), A4 (Asian variant), B1, B2, B3, B4, C1, C2, C3, and C4 (African variants), D1 (North-American variant), D2, D3 (Asian-American variants) and D4. Population studies showed that infections with viruses belonging to specific HPV16 sublineages confer different risks of viral persistence and cancer. In this study, 39 HPV16-positive cervical smears from European women living in Calabria (Italy) were analyzed for the presence of HPV16 variants. Cervical DNA extracts were processed by PCR to amplify L1, the Long Control Region (LCR), E6 and E7, which were sequenced. The sequences were concatenated and the 3169 nucleotides long fragments were characterized by BLAST and phylogenetic analysis. A total of 96 Single Nucleotide Polymorphism (SNPs) were detected, 29 of which mapping in the L1, 45 in the LCR, 15 in the E6 and 7 in the E7. The most common SNP was the T350G (29/39 samples, 74.4%), causing the L83 V amino acid change in the E6. Most of the HPV16 isolates (89.7%) had 99% of nucleotide (nt) identity to members of the A1 and A2 sublineages, while 4 isolates had 99% nt identity to members of the B2, B4, C1 and D4 sublineages. In conclusion, viruses belonging to the A1, A2, B2, B4, C1 and D4 HPV16 sublineages were found to circulate in the Calabria region.

Human papillomavirus genome variants and head and neck cancers: a perspective

Human papillomaviruses (HPV) cause infections that are responsible for diverse clinical manifestations from benign conditions to invasive cancer. As different HPV types are associated with variable pathogenic potential, minor genetic variations within a given high-risk HPV type might also be associated with distinct oncogenic capacities, through variable ability of persistence or risk of progression to precancer/cancer. Most recent HPV variant studies in the cervix using latest sequencing technology confirmed that minor changes in the HPV genome can have a major influence on carcinogenesis and have revealed key data that help better understand the carcinogenicity of HPV at a molecular level. Here we review the limited number of studies on HPV genome variants in head and neck cancers (HNC) and discuss their implications for cancer research in the light of accumulated knowledge for the cervix. Challenges in transposing HPV variant studies from the lower anogenital to the upper aerodigestive tract are also discussed, highlighting the main gaps of knowledge in the field of HPV-induced HNC. Specifically in the head and neck region, the lack of characterisation of precancerous lesions and the difficulty in sampling normal tissue will challenge the development of accurate studies. Although there is so far no indication that HPV variant research in HNC could directly translate into clinical application, such research is expected to be useful to disentangle unanswered questions in the pathogenesis of HNC. Yet, history of HPV variant research suggests that, to be successful, studies will require large international collaborative efforts.

Synthetic siRNA targeting human papillomavirus 16 E6: a perspective on in vitro nanotherapeutic approaches

High-risk human papillomaviruses infect skin and mucosa, causing approximately 5% of cancers worldwide. In the search for targeted nanotherapeutic approaches, siRNAs against the viral E6 transcript have been molecules of interest but have not yet seen successful translation into the clinic. By reviewing the past approximately 15 years of in vitro literature, we identify the need for siRNA validation protocols which concurrently evaluate ranges of key treatment parameters as well as characterize downstream process restoration in a methodical, quantitative manner and demonstrate their implementation using our own data. We also reflect on the future need for more appropriate cell culture models to represent patient lesions as well as the application of personalized approaches to identify optimal treatment strategies.

Association of codon 72 polymorphism of p53 with the severity of cervical dysplasia, E6-T350G and HPV16 variant lineages in HPV16-infected women

Purpose. Polymorphic variability in the tumour-suppressor protein p53 at codon 72 has a considerable impact on cervical cancer development. The present study clarified the association between p53 codon 72 genotypes and the risk of cervical disease in Greek patients. We also examined whether the presence of specific p53 genotypes in combination with HPV16 variants or E6 T350G sequence variation can modify an individual's susceptibility to cervical disease.Methodology. The analysis of p53 genotypes was performed through PCR-RFLP. Sequence and phylogenetic tree analyses of the HPV16 E6 gene were also performed in order to identify HPV16 variants and T350G sequence variation.Results/Key findings. The outcomes of the present analysis revealed that women who are homozygous for the arg genotype are at a 4.17-fold higher risk of developing HPV16-associated HSIL+ (OR=4.17, 95 % CI:1.48-4.9, P=0.0049). Moreover, p53 arg/arg patients infected by an HPV16 prototype strain were associated with an increased risk of more severe lesions, while a significant relationship between the p53 arg/arg genotype in patients with T350G sequence variation and the risk of high-grade squamous intraepithelial lesions (HSILs) was revealed.Conclusion. The oncogenic potential of the virus is increased by the presence of the p53 arg/arg genotype in the Greek population in such a way that the specific protein interaction E6 (L83V)-p53 (Arg-72) can modify an individual's susceptibility to cervical disease.

Two common variants of human papillomavirus type 16 E6 differentially deregulate sugar metabolism and hypoxia signalling in permissive human keratinocytes

Human papillomavirus type 16 (HPV16) is responsible for most cancers attributable to HPV infection and naturally occurring variants of the HPV16 E6 oncoprotein predispose individuals to varying risk for developing cancer. Population studies by us and others have demonstrated that the common Asian-American E6 (AAE6) variant is a higher risk factor for cervical cancer than the E6 of another common variant, the European prototype (EPE6). However, a complete understanding of the molecular processes fundamental to these epidemiological findings is still lacking. Our previously published functional studies of these two E6 variants showed that AAE6 had a higher immortalization and transformation potential than EPE6. Proteomic analysis revealed markedly different protein patterns between these variants, especially with respect to key cellular metabolic enzymes. Here, we tested the Warburg effect and hypoxia signalling (hallmarks of cancer development) as plausible mechanisms underlying these observations. Lactate and glucose production were enhanced in AAE6-transduced keratinocytes, likely due to raised levels of metabolic enzymes, but independent of hypoxia-inducible factor 1 alpha (HIF-1α) activity. The HIF-1α protein level and activity were elevated by AAE6 in hypoxic conditions, leading to a hypoxia-tolerant phenotype with enhanced migratory potential. The deregulation of HIF-1α was caused by the AAE6 variant's ability to augment mitogen-activated protein kinase/extracellular related kinase signalling. The present study reveals prominent underlying mechanisms of the AAE6's enhanced oncogenic potential.

Human papillomavirus type 16 lineage analysis based on E6 region in cervical samples of Iranian women

It is suggested that distinct HPV 16 variants differ in oncogenic potential and geographic distribution. As such, understanding the regional variants of HPV 16 would be of great importance for evolutionary, epidemiological and biological analysis. In this regard, the sequence variations of E6 gene were investigated to characterize more common variants of HPV 16 in normal cells, premalignant and malignant lesions of the cervix. In total, 106 isolates of HPV 16 were analyzed by PCR and sequencing. Overall, two different lineages (A and D) were identified. Lineage D comprised 70.7% of samples and the remaining 29.3% belonged to lineage A. Regarding to cytology/histology, lineage D was dominant in both normal+CIN I-II and CIN III+ICC groups as it was detected in 80% and 66.2% of cases, respectively. The comparison of the lineages between different groups (35 normal+CIN I-II samples and 71 CIN III+ICC samples) revealed that lineage A is more prevalent in cervical cancer cases (7 (20%) vs. 24 (33.8%)) although the difference observed did not reach statistical significance (p=0.07). In conclusion, our findings confirm that HPV lineages A and D are more prevalent in Iran, with the lineage D as the most dominant in all studied groups.

Functional variants of human papillomavirus type 16 demonstrate host genome integration and transcriptional alterations corresponding to their unique cancer epidemiology

BACKGROUND

Human papillomaviruses (HPVs) are a worldwide burden as they are a widespread group of tumour viruses in humans. Having a tropism for mucosal tissues, high-risk HPVs are detected in nearly all cervical cancers. HPV16 is the most common high-risk type but not all women infected with high-risk HPV develop a malignant tumour. Likely relevant, HPV genomes are polymorphic and some HPV16 single nucleotide polymorphisms (SNPs) are under evolutionary constraint instigating variable oncogenicity and immunogenicity in the infected host.

RESULTS

To investigate the tumourigenicity of two common HPV16 variants, we used our recently developed, three-dimensional organotypic model reminiscent of the natural HPV infectious cycle and conducted various "omics" and bioinformatics approaches. Based on epidemiological studies we chose to examine the HPV16 Asian-American (AA) and HPV16 European Prototype (EP) variants. They differ by three non-synonymous SNPs in the transforming and virus-encoded E6 oncogene where AAE6 is classified as a high- and EPE6 as a low-risk variant. Remarkably, the high-risk AAE6 variant genome integrated into the host DNA, while the low-risk EPE6 variant genome remained episomal as evidenced by highly sensitive Capt-HPV sequencing. RNA-seq experiments showed that the truncated form of AAE6, integrated in chromosome 5q32, produced a local gene over-expression and a large variety of viral-human fusion transcripts, including long distance spliced transcripts. In addition, differential enrichment of host cell pathways was observed between both HPV16 E6 variant-containing epithelia. Finally, in the high-risk variant, we detected a molecular signature of host chromosomal instability, a common property of cancer cells.

CONCLUSIONS

We show how naturally occurring SNPs in the HPV16 E6 oncogene cause significant changes in the outcome of HPV infections and subsequent viral and host transcriptome alterations prone to drive carcinogenesis. Host genome instability is closely linked to viral integration into the host genome of HPV-infected cells, which is a key phenomenon for malignant cellular transformation and the reason for uncontrolled E6 oncogene expression. In particular, the finding of variant-specific integration potential represents a new paradigm in HPV variant biology.

Analysis of human papillomavirus 16 variants and risk for cervical cancer in Chinese population

HPV16 is the most carcinogenic HPV type, but only a minority of HPV16 infections progress to cancer. Intratype genetic variants of HPV16 have been suggested to confer differential carcinogenicity. To investigate risk implications of HPV16 variants among Chinese women, a case-control study was conducted with 298 cervical cancer patients and 85 controls (all HPV16-positive). HPV16 isolates were predominantly of the A variant lineage, and variants of A4 (previously named "Asian") sublineage were common. A4/Asian variants were significantly associated with increased risk of cervical cancer compared to A1-3 (OR=1.72, 95% CI=1.04-2.85). Furthermore, a meta-analysis including 703 cases and 323 controls from East Asia confirmed the association (OR=2.82, 95% CI=1.44-5.52). In conclusion, A4 variants appear to predict higher risk of cervical cancer among HPV16-positive women, which may provide clues to the genetic basis of differences in the carcinogenicity of HPV16 variants.

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.