Uneven distribution of HPV 16 E6 prototype and variant (L83V) oncoprotein in cervical neoplastic lesions

The difference of transcriptome of HPV-infected patients contributes more to the occurrence of cervical cancer than the mutations of E6 and E7 genes in HPV16

Human papillomavirus (HPV) E6 and E7 genes are biomarkers and drivers of the progression of cervical cancer (CxCa). The aim of this study was to investigate the relationship between HPV16 E6, E7 gene mutations and the occurrence and development of CxCa. Cervical exfoliated cells and clinical data of patients with cervical diseases were collected. Sample DNA was extracted, the E6 and E7 gene fragments were amplified by PCR, and the mutations were detected by Sanger sequencing and compared with standard sequences. Microarray was used to sequence the transcriptome of cells. Data of transcriptome analyzed and visualized using R software and its packages. Analysis of clinical characteristics demonstrated the association of HPV16 infection with CxCa (P < .05). Sanger sequencing results showed that the mutation sites of E6 gene included T178G/A, T350G, A131C, and T241G; among these, A131C and T241G were synonymous mutations. The mutation sites of E7 gene included A647G, T846C, G666A, T843C, and T760C, and all of them were synonymous mutations except A647G. There was no significant difference in the distribution of HPV16 E6, E7 mutations among CxCa, cervical intraepithelial neoplasia, and infection groups (P > .05). Compared with the non- CxCa group, gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis of differentially expressed genes (DEGs) showed more significant enrichment of DEGs in the biological processes, pathways, and diseases closely related to cancer. Compared with the non-mutation group, the DEGs in the E6, E7 gene mutation group were significantly enriched in the events related to infection and immunity. To summarize, HPV16 may be associated with the occurrence and development of CxCa, but HPV16 E6 and E7 gene mutations have little effect on the occurrence and development of CxCa. Individual differences may have a greater effect on the progression of CxCa.

Genetic variation of E6 and E7 genes of human papillomavirus type 16 from central China

BACKGROUND

Persistent high-risk human papillomavirus (HR-HPV) infection is an important factor in the development of cervical cancer, and human papillomavirus type 16 (HPV-16) is the most common HR-HPV type worldwide. The oncogenic potential of HPV-16 is closely related to viral sequence variation.

METHODS

In order to clarify the variant characteristics of HPV-16 E6 and E7 genes in central China, E6 and E7 sequences of 205 HPV-16 positive samples were amplified by polymerase chain reaction. PCR products of E6 and E7 genes were further sequenced and subjected to variation analysis, phylogenetic analysis, selective pressure analysis and B-cell epitope prediction.

RESULTS

Twenty-six single nucleotide variants were observed in E6 sequence, including 21 non-synonymous and 5 synonymous variants. Twelve single nucleotide variants were identified in E7 sequence, including 6 non-synonymous and 6 synonymous variants. Four new variants were found. Furthermore, nucleotide variation A647G (N29S) in E7 was significantly related to the higher risk of HSIL and cervical cancer. Phylogenetic analysis showed that the E6 and E7 sequences were all distributed in A lineage. No positively selected site was found in HPV-16 E6 and E7 sequences. Non-conservative substitutions in E6, H31Y, D32N, D32E, I34M, L35V, E36Q, L45P, N65S and K75T, affected multiple B-cell epitopes. However, the variation of E7 gene had little impact on the corresponding B-cell epitopes (score < 0.85).

CONCLUSION

HPV-16 E6 and E7 sequences variation data may contribute to HR-HPV prevention and vaccine development in Jingzhou, central China.

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.

Association of cervical carcinogenesis risk with HPV16 E6 and E7 variants in the Taizhou area, China

Background

Human papillomavirus (HPV) type 16 accounts for a larger share of cervical cancer and has been a major health problem worldwide for decades. The progression of initial infection to cervical cancer has been linked to viral sequence properties; however, the role of HPV16 variants in the risk of cervical carcinogenesis, especially with longitudinal follow-up, is not fully understood in China.

Methods

We aimed to investigate the genetic variability of HPV16 E6 and E7 oncogenes in isolates from cervical exfoliated cells. Between December 2012 and December 2014, a total of 310 single HPV16-positive samples were selected from women living in the Taizhou area, China. Sequences of all E6 and E7 oncogenes were analysed by PCR-sequencing assay. Detailed sequence comparison, genetic heterogeneity analyses and maximum-likelihood phylogenetic tree construction were performed with BioEdit Sequence Alignment Editor and MEGA X software. Data for cytology tests and histological diagnoses were obtained from our Taizhou Area Study with longitudinal follow-up for at least 5 years. The relationship between HPV16 variants and cervical carcinogenesis risk was analysed by the chi-square test or Fisher’s exact test.

Results

In this study, we obtained 64 distinct variation patterns with the accession GenBank numbers MT681266-MT681329. Phylogenetic analysis revealed that 98.3% of HPV16 variants belong to lineage A, in which the A4 (Asian) sublineage was dominant (64.8%), followed by A2 (12.1%), A1 (11.4%), and A3 (10.0%). The A4 (Asian) sublineage had a higher risk of CIN2+ than the A1–3 (European) sublineages (OR = 2.69, 95% CI = 1.04–6.97, P < 0.05). Furthermore, nucleotide variation in HPV16 E6 T178G is associated with the development of cervical cancer.

Conclusion

These data could provide novel insights into the role of HPV16 variants in cervical carcinogenesis risk in China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08531-y.

HPV16 E6 regulates the proliferation, invasion, and apoptosis of cervical cancer cells by downregulating miR-504

BACKGROUND

Human papillomavirus (HPV) 16 infection is a necessary condition for the pathogenesis and development of cervical cancer. The E6 protein is expressed by the HPV16 E6 gene and promotes malignant phenotype transformation, which is an important mechanism for the occurrence and development of cervical cancer. MicroRNA-504 (miR-504) has been reported as an oncogene or tumor suppressor gene; the expression of miR-504 in cervical cancer has been found to be negatively correlated with HPV infection. However, the relationship between HPV16 E6 and miR-504 and the role of miR-504 in cervical cancer are not clear. In the current study, we observed the effect of HPV16 E6 on the expression of miR-504 in cervical cancer cells, and analyzed whether HPV16 E6 affects proliferation, invasion, and apoptosis in cervical cancer cells by regulating the expression of miR-504.

METHODS

Cervical cancer cells (SiHa) were divided into four groups: the empty vector group, E6 overexpression group, E6 overexpression + miR-NC group, and E6 overexpression+miR-504 group. The expressions levels of HPV16 E6 mRNA and miR-504 were detected by real-time polymerase chain reaction (PCR), and the expression level of HPV16 E6 protein was detected by Western blot. Cell proliferation, invasion, and apoptosis were detected by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay, Tastelessly, and flow cytometry, respectively.

RESULTS

The expression level of miR-504 was significantly decreased in E6 overexpression cells compared to the control cells (P<0.05); the overexpression of miR-504 with miR-504 mimic significantly reversed the downregulation of miR-504 in E6 overexpression SiHa cells (P<0.05). MTT and Transwell assays showed that the overexpression of E6 significantly increased proliferation and invasion of SiHa cells (P<0.05). The overexpression of miR-504 reversed the role of HPV16 E6 on the proliferation and invasion in E6 overexpression SiHa cells, and the difference was statistically significant (P<0.05). Further analysis showed that the overexpression of E6 significantly reduced apoptosis of SiHa cells (P<0.05). The overexpression of miR-504 reversed the role of HPV16 E6 on apoptosis in E6 overexpression SiHa cells, and the difference was statistically significant (P<0.05).

CONCLUSIONS

HPV16 E6 may promote the proliferation and invasion, and inhibit the apoptosis, of cervical cancer SiHa cells by downregulating miR-504 expression.

Bioinformatics analysis of HPV-68 E6 and E7 oncoproteins for designing a therapeutic epitope vaccine against HPV infection

The incidence and mortality of cervical cancer, which mainly results from the infection of human papillomavirus (HPV) is significantly increasing in Xinjiang. According to the previous research, the incidence of HPV-68 in cervical cancer patients in Xinjiang is significantly higher than in other parts of China. HPV E6 and E7 oncoproteins play a crucial role in cervical cancer, and can be used as ideal targets for therapeutic vaccines. Therefore, we analyzed and identified the possible T-cell and B-cell dominant epitopes and various aspects of HPV-68 E6 and E7 oncoproteins, including the physicochemical properties, secondary and tertiary structures using a bioinformatic approach, which provided a basis for designing an effective HPV infection therapeutic vaccine. The results showed that E6 oncoproteins was an unstable and hydrophilic protein, while E7 oncoproteins was unstable and hydrophilic protein. The secondary structure of the E6 oncoproteins consisted of 45.57% alpha helixes, 14.56% extended strands, 4.43% beta turns and 35.44% random coils. The secondary structure of E7 oncoproteins consisted of 35.45% alpha helixes, 17.27% extended strands, 0.91% beta turns and 46.36% random coils. Moreover, our results identified 5 dominant T-cell epitopes and 6 dominant B-cell epitopes in the E6 oncoproteins structure and 5 dominant T-cell epitopes and 3 dominant B-cell epitopes in E7 oncoproteins. In conclusion, this study provides comprehensive biological information about the HPV-68 E6 and E7 oncoproteins, which will lay a theoretical foundation for multi-epitope vaccines against HPV infection.

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.

Association of Human Papillomavirus Genotype 16 Viral Variant and Viral Load with Cervical High-grade Intraepithelial Lesions

Human papillomavirus genotype 16 (HPV16) is by far the genotype most strongly associated with cervical cancer; viral variant and/or viral load of HPV16 could modulate this association. The objective was to determine the association between the viral variant and viral load of HPV16 and the presence of cervical high-grade lesions. This cross-sectional study included all women in whom HPV infection was found by cervical smear during routine gynecologic health checks. Women with single or multiple HPV16 infections (n = 176) were selected for viral variant and viral load analysis. Smear results were classified using the Bethesda system. HPV types were classified according to the International Agency for Research on Cancer. Odds ratios (OR) with their 95% confidence intervals (CI) were estimated by logistic regression, adjusted for age, immigrant status, and coinfection with other high-risk genotypes. No statistically significant associations were found regarding the detected viral variants. A viral load above the median (>1,367.79 copies/cell) was associated with a significant risk of high-grade epithelial lesion or carcinoma, after adjusting for age, immigrant status, coinfections, and viral variant: (adjusted OR 7.89; 95% CI: 2.75-22.68). This relationship showed a statistically significant dose-response pattern after categorizing by viral load tertiles: adjusted OR for a viral load greater than the third tertile was 17.23 (95% CI: 4.20-70.65), with adjusted linear P _trend = 0.001. In patients infected with HPV16, viral load is associated with high-grade intraepithelial lesions or cervical carcinoma. This could be useful as prognostic biomarker of neoplastic progression and as screening for cervical cancer.

Identification of variants and therapeutic epitopes in HPV-33/HPV-58 E6 and E7 in Southwest China

BACKGROUND

Human papillomavirus (HPV) E6 and E7 oncoproteins play a crucial role in HPV-related diseases, such as cervical cancer, and can be used as ideal targets for therapeutic vaccines. Human leukocyte antigen (HLA) participates in the immune response to block HPV infection and invasion by its target/recognition function. HPV-33 and HPV-58 are highly prevalent among Chinese women. Therefore, it is of great significance to study the E6 and E7 region-specific gene polymorphisms of HPV-33 and HPV-58 in Southwest China and to identify ideal epitopes for vaccine design. Both HPV-33 and HPV-58 belong to α-9 genus HPV and are highly homologous, so their correlations are included in our research.

METHODS

To study the E6 and E7 variations and polymorphisms of HPV-33 and HPV-58 in Southwest China, we collected samples, extracted and sequenced DNA, and identified variants. Nucleotide sequences were translated into amino acids by Mega 6.0 software. The physical/chemical properties, amino acid-conserved sequences and secondary structure of protein sequences were analysed by the Protparam server, ConSurf server and PSIPRED software. The T and B cell epitopes of the E6/E7 reference and variant sequences in HPV-33 and HPV-58 were predicted by the Immune Epitope Database (IEDB) analysis server and the ABCpred server, respectively.

RESULTS

Five and seven optimal HLA-I restricted T cell epitopes were selected from HPV-33 and HPV-58 E6, respectively, and these optimal epitopes are mainly located in _41-58EVYDFAFADLTVVYREGN of HPV-33 E6 and _40-60SEVYDFVFADLRIVYRDGNPF of HPV-58 E6. Six optimal HLA-I-restricted T cell epitopes were selected from HPV-33 and HPV-58 E7, and these epitopes are mainly located in _77-90RTIQQLLMGTVNIV of HPV-33 E7 and _78-91RTLQQLLMGTCTIV of HPV-58 E7.

CONCLUSIONS

HPV-33/HPV-58 E6/E7 gene polymorphisms and T/B cell epitopes of their reference and variant sequences were studied, and candidate epitopes were selected by bioinformatics techniques for therapeutic vaccine design for people in Southwest China. This study was the first to investigate the correlation of epitopes between HPV-33 and HPV-58. After experimental validation, these selected epitopes will be employed to induce a wide range of immune responses in heterogeneous HLA populations.

Human papillomavirus type 16 E6 and E7 gene variations associated with cervical cancer in a Han Chinese population

BACKGROUD

Human papillomavirus type 16 (HPV16) is a high-risk HPV subtype and a potent carcinogen. The HPV16 E6 and E7 genes are considered oncogenes that play a core role in the development of cervical cancer.

METHODS

In the current study, we enrolled 97 HPV16-positive cervical cancer patients (case group) and 136 HPV16-positive asymptomatic individuals (control group) in a study to analyse the association between HPV16 E6 and E7 gene variations and cervical cancer.

RESULTS

Our results showed that three HPV16 sub-lineages (A1-A3, A4 and D3) were present; the distribution of these variants between the case and control group was not significantly different (P = 0.178). When the distribution of the HPV16 E6 and E7 gene variations was compared, the distribution of only A131C (R10R) in the E6 gene showed a different trend between the case and control groups and C749T (S63F) in the E7 gene was significantly different between the case and control groups (P = 0.071 and P = 4.861 × 10^-10, respectively). Regarding the sub-lineages, no variations in the E6 gene were significantly different between the case and control group for the A4 (As) and A1-A3 (EUR) sub-lineages. However, the distribution of C749T (S63F) in the E7 gene was significantly different between the case and control groups for the A4 (As) and A1-A3 (EUR) sub-lineages (P = 1.815 × 10^-8 and P = 0.008). In the current study, we found that the C749T (S63F) variation in the HPV16 E7 gene was associated with cervical cancer not only in the A4 (As) sub-lineage but also in the A1-A3 (EUR) sub-lineage.

CONCLUSION

Our study will provide a good reference for further functional studies of the relationship between cervical cancer carcinogenesis and the HPV16 E6 and E7 genes.

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.

Molecular approaches for HPV genotyping and HPV-DNA physical status

Persistent infection with high-risk human papillomavirus (HPV) genotypes is the leading cause of cervical cancer development. To this end several studies have focused on designing molecular assays for HPV genotyping, which are considered as the gold standard for the early diagnosis of HPV infection. Moreover, the tendency of HPV DNA to be integrated into the host chromosome is a determining event for cervical oncogenesis. Thus, the establishment of molecular techniques was promoted in order to investigate the physical status of the HPV DNA and the locus of viral insertion into the host chromosome. The molecular approaches that have been developed recently facilitate the collection of a wide spectrum of valuable information specific to each individual patient and therefore can significantly contribute to the establishment of a personalised prognosis, diagnosis and treatment of HPV-positive patients. The present review focuses on state of the art molecular assays for HPV detection and genotyping for intra-lesion analyses, it examines molecular approaches for the determination of HPV-DNA physical status and it discusses the criteria for selecting the most appropriate regions of viral DNA to be incorporated in HPV genotyping and in the determination of HPV-DNA physical status.

Identification of the impact on T- and B- cell epitopes of human papillomavirus type-16 E6 and E7 variant in Southwest China

Cervical cancers almost are infected by human papillmavirus (HPV), encoding E6 and E7 oncoproteins which are regard as ideal targets on the mechanism of this disease and development of vaccines. HLA (human leukocyte antigen) participates in the local immune response to prevent tumor invasion and progression. But due to highly polymorphism of HLA, prediction shows its importance in this study. More effective immunoinformatics was used for predicting epitopes from HPV-16 E6 and E7, including T- and B-cell epitopes. Eight substitutions are detected. Specifically speaking, for HLA-I, HLA-A*33:03 (26), HLA-B*13:01 (14), HLA-C*03:02 (5) for E6 and HLA-A*02:01 (6), HLA-B*40:01 (5), HLA-C*03:04 (4) for E7 are most frequency. Epitope _41-48EVYDFAFR for HLA-A*33:03 (0.1) for E6 has best binding affinity, as well as HLA*02:01 and HLA-B*40:01 (0.2) for E7. The mutations of D25E and L83V of E6 and N29S of E7 produce new epitopes, and the percentile values change with them. For HLA-II, seventeen epitopes in the reference at percentile value from 0.22 to 4.76, while in variant from 0.22 to 4.96. For the B-cell epitopes, three most potent epitopes for E6 were listed, and N29S lead the growth of score from 0.81 to 0.83. In summary, E6_40-55REVYDFAFRDLCIVYR and E7_11-22YMLDLQPETTDL are the important regions, containing the majority of predicted epitopes. E6 72-83 for HLA-A*02:01 and E6 74-84 for HLA-B*15:02 maybe are the new direct for therapeutic vaccine aimed at L83V variants. HLA-DRB1*15:02 is better binder with T cell in our HLA class II. It is a systematic, detail recognition for T- and B-cell epitopes of HPV-16 E6 and E7 from Southwest China, which may be helpful to design vaccines specifically for women in Southwest China and testing methods specifically for this region. The results of our study may contribute to future researches on vaccines improvement, or screening methods for a particular population.

CD8+ T cell response to human papillomavirus 16 E7 is able to predict survival outcome in oropharyngeal cancer

INTRODUCTION

Immunological response to human papillomavirus (HPV) in the development and progression of HPV16+ oropharyngeal squamous cell carcinoma (OPSCC) (accounting for the majority of viral associated cases) is largely unknown and may provide important insights for new therapeutic strategies.

METHODS

In this prospective clinical trial (UKCRN11945), we examined cell-mediated immune responses to HPV16 E2, E6 and E7 in peripheral blood using IFN-γ enzyme-linked immunosorbent spot assay. CD56⁺, CD4⁺, CD8⁺ and regulatory T cell frequencies were also discerned by flow cytometry. Fifty-one study participants with oropharyngeal carcinoma were recruited. Control subjects were those undergoing tonsillectomy for benign disease. All patients were treated with curative intent by radiotherapy ± chemotherapy. Disease-specific survival was investigated by multivariate analysis.

RESULTS

HPV16 DNA was detected in 41/51 of the OPSCC participants. T cell responses against HPV16 E6 or E7 peptides were detected in 33/51 evaluable patients, respectively and correlated with HPV status. Matched pre- and post-treatment T cell responses were available for 39/51 OPSCC cases. Within the whole cohort, elevated post-treatment CD8⁺ response to HPV16 E7 correlated with longer disease free survival (multivariate DFS p < 0.03). Within the HPV + OPSCC cohort, a significant increase in regulatory T cells (p < 0.02) was noted after treatment.

CONCLUSIONS

This is the first study to provide survival data in OPSCC stratified by cell-mediated immune response to HPV16 peptides. Within the HPV16+ OPSCC cohort, enhanced immunoreactivity to antigen E7 was linked to improved survival. An increase in regulatory T cell frequencies after treatment may suggest that immunosuppression can contribute to a reduced HPV-specific cell-mediated response.

Association of human papillomavirus type 16 and its genetic variants with cervical lesion in Korea

Persistent human papillomavirus type 16 (HPV16) is the major risk factor for cervical cancer. HPV16 intratypic variants differ in their geographical distribution and oncogenic potential. This study aimed to analyze the distribution of HPV16 variants and their association with cervical lesion histopathology in Korean women. In total, 133 HPV16-positive cervical samples from women admitted to Seoul National University Boramae Hospital were analyzed by sequencing E6, E7, and L1 genes and the long control region (LCR), and the variant distribution according to cervical lesion grade was determined. Isolates were grouped into a phylogenetic lineage, and A1-3, A4, C, and D sublineages were detected in 54.1, 37.8, 0.7, and 7.4% of samples, respectively. The most commonly observed LCR variations were 7521G>A (91.5%), 7730A>C (59.6%), and 7842G>A (59.6%). Furthermore, A4 or D sublineage-positive women had a higher risk for cervical cancer than women who were positive for A1-3. Among HPV phylogenetic clusters, A1-3 was the predominant sublineage, and within A1-3, the 350G polymorphism was highly frequent. These results differed from those of previous studies in Korea and other Asian countries. The findings suggest that cervical neoplasia incidence in HPV16-infected patients could be affected by the distribution of HPV16 variants in the population.

Genetic variability and lineage phylogeny of human papillomavirus type-16 and -53 based on the E6, E7, and L1 genes in Southwest China

Human papillomaviruses (HPVs) are circular double-stranded DNA viruses that are highly prevalent in the general population, and account for the cervical cancer burden in women worldwide. In this study, we analyzed HPV-16, the most prevalent type worldwide, and HPV-53, a possible high-risk type from infected women in Southwest China. To characterize mutations, intratypic variants, and genetic variability in the E6, E7, and L1 genes of HPV-16 (n=97) and HPV-53 (n=15), these genes were sequenced and submitted to GenBank. Phylogenetic trees were constructed using Bayesian trees, followed by secondary structure analysis and B-cell epitope prediction. Moreover, the selection pressures of the E6, E7, and L1 genes were estimated. In total, 27 novel variants of HPV-16 and 11 novel variants of HPV-53 were identified. In the HPV-16 E6-E7-L1 sequences, 73 nucleotide changes were observed with 40/73 being non-synonymous mutations (two in the alpha helix and five in the beta sheet) and 33/73 being synonymous. In the HPV-53 E6-E7-L1 sequences, 64 nucleotide changes were observed with 26/64 being non-synonymous mutations (three in the alpha helix and one in the beta sheet) and 38/64 being synonymous. Selective pressure analysis showed that most of these mutations did not reflect positive selection. The maximal divergence between any two variants within each gene of these two HPV types ranging from 0.94%(HPV-16 L1 gene)to 2.80%(HPV-53 E6 gene). Identifying new variants of HPV-16 and -53 from women in Southwest China may be helpful to design vaccines specifically for women in Southwest China and testing methods specifically for this region. The results of our study may contribute to future researches in diagnostic probes, vaccines improvement, or screening methods for a particular population.

Genetic variations of human papillomavirus type 16: implications for cervical carcinogenesis

Human papillomaviruses (HPVs) are the causative agent of cervical cancer, and among approximately 15 high-risk genotypes, HPV16 accounts for more than half the cases of cervical cancer worldwide. Recent progress in determining HPV genomic sequences from clinical samples has revealed a wide variety in HPV16 genome sequences, and has allowed for comprehensive classification of intratype HPV16 variants. These consist of four variant lineages containing nucleotide variations in 1.0%-10.0% of the complete viral genome sequence. Epidemiological data suggest that the non-European-Asian lineages of HPV16 entail a higher risk of progression to invasive cervical cancer than the European-Asian lineage. Deep sequencing analysis has recently demonstrated that HPV16 genome sequences are highly homogeneous in individual clinical specimens compared with those of RNA viruses. However, an extremely sensitive PCR method, differential DNA denaturation PCR, has detected hypermutations from C to T or G to A in the E2 gene and the long control region of the HPV16 genome, which suggests the involvement of cellular apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) proteins in this hypermutation. The quasispecies status of the HPV16 genome in the infected cervix may affect the development of cervical cancer and warrants further investigation.

Human papillomavirus 16 non-European variants are preferentially associated with high-grade cervical lesions

HPV16 accounts for 50–70% of cervical cancer cases worldwide. Characterization of HPV16 variants previously indicated that they differ in risks for viral persistence, progression to cervical precancer and malignant cancer. The aim of this study was to examine the association of severity of disease with HPV16 variants identified in specimens (n = 281) obtained from a Cervical Pathology and Colposcopy outpatient clinic in the University Hospital of Espírito Santo State, Southeastern Brazil, from April 2010 to November 2011. All cytologic and histologic diagnoses were determined prior to definitive treatment. The DNA was isolated using QIAamp DNA Mini Kit and HPV was detected by amplification with PGMY09/11 primers and positive samples were genotyped by RFLP analyses and reverse line blot. The genomes of the HPV16 positive samples were sequenced, from which variant lineages were determined. Chi² statistics was performed to test the association of HPV16 variants between case and control groups. The prevalence of HR-HPV types in <CIN1, CIN2 and CIN3+ were 33.7%, 84.4% and 91.6%, respectively. Thirty-eight of 49 (78%) HPV16 positive samples yielded HPV16 sequence information; of which, 32 complete genomes were sequenced and an additional 6 samples were partially sequenced. Phylogenetic analysis and patterns of variations identified 65.8% (n = 25) as HPV16 European (E) and 34.2% (n = 13) as non-European (NE) variants. Classification of disease into CIN3+ vs. <CIN3 indicated that NE types were associated with high-grade disease with an OR = 4.6 (1.07–20.2, p = 0.05). The association of HPV16 NE variants with an increased risk of CIN3+ is consistent with an HPV16 genetically determined enhanced oncogenicity. The prevalence of genetic variants of HPV16 is distributed across different geographical areas and with recent population admixture, only empiric data will provide information on the highest risk HPV16 variants within a given population.

Human papillomavirus type 16 long control region and E6 variants stratified by cervical disease stage

•

We sequenced HPV16 LCR–E6 variants in cervical disease samples from England.

•

98% of variants were of the EUR lineage.

•

Site-specific entropy identified several variable sites in the LCR and E6.

•

No single or combination of sites were associated with disease, including E6 T350G.

Objective

Certain intra-type variants of HPV16 have been shown to be associated with an increased risk of developing high grade cervical disease, but their potential association is confounded by apparent geographic and phylogenetic lineage dependency. The objective of this study was to evaluate the relationship between HPV16 sequence variants and cervical disease stage in monospecific infection samples from a single lineage (European, EUR) in England.

Methods

One hundred and twelve women singly infected with HPV16 and displaying normal and abnormal cytology grades were selected. An 1187 bp fragment encompassing the entire LCR and a portion of the E6 open reading frame was sequenced to identify intra-type variants. Intra-type diversity was estimated using Shannon entropy.

Results

Almost all samples (110/112; 98%) were assigned to the EUR lineage, one sample was classified as European-Asian (EAS) and another African (Afr1a). The mean pairwise distance of the EUR sequences in this study was low (0.29%; 95%CI 0.13–0.45%) but there were nevertheless several sites in the LCR (n = 5) and E6 (n = 2) that exhibited a high degree of entropy. None of these sites, however, including the T350G non-synonymous (L83V) substitution in E6, alone or in combination, were found to be associated with cervical disease stage.

Conclusions

Despite using single infection samples and samples from a single variant lineage, intra-type variants of HPV16 were not differentially associated with cervical disease. Monitoring intra-lineage, site-specific variants, such as T350G, is unlikely to be of diagnostic value.

Molecular and evolutionary analysis of HPV16 E6 and E7 genes in Greek women

Human papillomavirus type 16 (HPV16) non-European variants have been associated with persistent infection and cervical cancer development, while the L83V variant of the E6 gene has been correlated with the progression of cervical malignancy. The present study investigated the presence of the HPV16 L83V variant in Greek women. Molecular evolutionary analysis of the HPV16 E6 and E7 oncogenes was conducted in order to estimate the evolution of the HPV16 genome in the Greek population. The E6 L83V variant was found in 78.2 % of high- and 64.28 % of low-grade specimens. Moreover, the prototype and E6 L83V variants were both prevalent in high- and low-grade malignancies in Greek women. Selective pressure analysis of the individual amino acid residues of HPV16 sequences from the Greek population indicates that codon 83 of the E6 protein, as well as codon 85 of the E7 protein, are undergoing positive selection. Novel sequence variations were recorded within the E6 and E7 genes in cervical samples, characterized as (T350G) European variants. However, no signal of intratypic recombination event was identified within the E6–E7 region. Molecular and evolutionary analyses of HPV16 genomes from distinct geographical locations might provide valuable information about viral evolution and oncogenecity.

Degradation of p53 by natural variants of the E6 protein of human papillomavirus type 16

The degradation of p53 by high-risk human papillomavirus (HR-HPV) E6 proteins is recognized as necessary for the immortalization of mammary epithelial cells and the progression of cancer. The HR-HPV type 16 E6 proteins exhibit numerous variants associated with different risk factors for the development of cervical cancer. Two variants of E6 proteins, D25E and L83V, are common in cervical carcinomas among Asian and European populations. In the present study, we compared the effect of two E6 variants on p53 degradation by a prototype E6 protein. We demonstrate that both the D25E and L83V variants downregulate p53 through a ubiquitin-proteasome pathway, and that the effect is very similar to that of the prototype E6 protein. The reduction in the p53 protein levels was induced through the ubiquitin-proteasome pathway via interaction with E6 proteins. The expression of p21 CIP1/WAF1, a downstream molecule of p53, was similarly reduced in both prototype and variant E6 protein-expressing cell lines, leading to aberrant G1/S cell cycle arrest. These results suggest that the natural variants, E6 D25E and L83V, similar to the prototype E6 protein, contribute to tumorigenesis by degrading p53.

Whole genome analysis of human papillomavirus type 16 multiple infection in cervical cancer patients

The characterization of the whole genome of human papillomavirus type 16 (HPV16) from cervical cancer specimens with multiple infections in comparison with single infection samples as the oncogenic potential of the virus may differ. Cervical carcinoma specimens positive for HPV16 by PCR and INNO-LiPA were randomly selected for whole genome characterization. Two HPV16 single infection and six HPV16 multiple infection specimens were subjected to whole genome analysis by using conserved primers and subsequent sequencing. All HPV16 whole genomes from single infection samples clustered in the European (E) lineage while all multiple infection specimens belonged to the non-European lineage. The variations in nucleotide sequences in E6, E7, E2, L1 and Long control region (LCR) were evaluated. In the E6 region, amino acid changes at L83V were related to increased cancer progression. An amino acid variation N29S within the E7 oncoprotein significantly associated with severity of lesion was also discovered. In all three domains of the E2 gene non synonymous mutations were found. The L1 region showed various mutations which may be related to conformation changes of viral epitopes. Some transcription factor binding sites in the LCR region correlated to virulence were shown on GRE/1, TEF- 1, YY14 and Oct-1. HPV16 European variant prone to single infection may harbor a major variation at L83V which significantly increases the risk for developing cervical carcinoma. HPV16 non-European variants prone to multiple infections may require many polymorphisms to enhance the risk of cervical cancer development.

Whole genome sequencing and evolutionary analysis of human papillomavirus type 16 in central China

Human papillomavirus type 16 plays a critical role in the neoplastic transformation of cervical cancers. Molecular variants of HPV16 existing in different ethnic groups have shown substantial phenotypic differences in pathogenicity, immunogenicity and tumorigenicity. In this study, we sequenced the entire HPV16 genome of 76 isolates originated from Anyang, central China. Phylogenetic analysis of these sequences identified two major variants of HPV16 in the Anyang area, namely the European prototype (E(p)) and the European Asian type (E(As)). These two variants show a high degree of divergence between groups, and the E(p) comprised higher genetic diversity than the E(As). Analysis with two measurements of genetic diversity indicated that viral population size was relatively stable in this area in the past. Codon based likelihood models revealed strong statistical support for adaptive evolution acting on the E6 gene. Bayesian analysis identified several important amino acid positions that may be driving adaptive selection in the HPV 16 population, including R10G, D25E, L83V, and E113D in the E6 gene. We hypothesize that the positive selection at these codons might be a contributing factor responsible for the phenotypic differences in carcinogenesis and immunogenicity among cervical cancers in China based on the potential roles of these molecular variants reported in other studies.

Human papillomavirus type 16 variant analysis of E6, E7, and L1 [corrected] genes and long control region in [corrected] cervical carcinomas in patients in northeast China

Human papillomavirus type 16 (HPV 16) plays a cardinal role in the pathogenesis of cervical cancer. HPV 16 has intratypic variants which show different geographical distributions and different oncogenic potentials. To analyze the presence of sequence variations of HPV 16 variants in northeast China, 71 cervical carcinomas were identified by HPV typing. HPV 16-positive specimens were analyzed by PCR-directed sequencing in the E6, E7, and L1 genes and the LCR (long control region). The variation data were compared with those of neighboring districts. In this hospital-based study, HPV 16 was the most common type (73.24%). In HPV 16-positive specimens, 67.31% belonged to the European (E) lineage, while 32.69% were Asian (As) variants. The Asian-American (AA), African-1 (Af-1), African-2 (Af-2), and northern American (NA) lineages were not detected. The most frequently observed variation sites were T178G (32.69%) in E6; A647G (34.62%), G666A (38.46%), and T846C (32.69%) in E7; C6826T (36.17%) and G7060A (61.70%) in L1; and G7521A (98.08%) in the LCR. The most prevalent amino acid variations were D25E in E6 and N29S in E7. In addition, 28 novel variations of HPV 16 were reported. Some covariations between different genes were obtained. In this study, HPV 16 variants belonged to the European lineage and the Asian lineage. Compared with neighboring districts, the distribution of HPV 16 variants in northeast China had a typical pattern. As the first report on HPV 16 variants in northeast China, it should be helpful for designing a HPV vaccine and HPV vaccination program in China.

Genetic Diversity of HPV-16E6,E7, andL1Genes in Women With Cervical Lesions in Liaoning Province, China

Introduction

High-risk human papillomaviruses (HPVs) play a cardinal role in the etiology of cervical cancer. The most prevalent type, HPV-16, shows intratypic sequence variants that are known to differ in oncogenic potential and geographic distribution. Intratype variations in oncogenic E6/E7 and capsid L1 proteins of HPV-16 are associated with risk of viral persistence and progression.

Methods

This study was designed to analyze sequence variations inE6,E7, andL1genes of HPV-16 in patients with cervical lesion to identify the most prevalent and novel HPV-16 variants in northern China.

Results

Our results showed that HPV-16 variants with respect to E6 and E7 were high prevalence of the Asian lineage: 48.3% and 51.4%, respectively. Sequences of theE6gene revealed 4 amino acid changes of variants D25E and L83V, with 48.3% (69/143) and 11.2% (16/143), respectively, and variants H78Y and E113D in this study. The results also showed the prevalence of 4 hot spots of E7 nucleotide variations leading to N29H, N29S, and 2 silent variations, nucleotide G666A and nucleotide T846C, with 4.2% (6/142), 43% (61/142), 32.4% (46/142), and 43% (61/142), respectively. The following L1 variations were found in this study: L103F, P104K, P104Y, P104S, D105G, P106S, N108P, F109V, C172S, H228D, and T292A. It was also found that 448S was inserted and 465D was deleted in the L1 amino acid sequences of all the samples. No significant relationship between HPV-16 variants and high-grade lesions was found.

Conclusions

The study provides some new data on the genetic diversity of HPV-16, which may help to understand the oncogenic potential of the virus and design the diagnosis reagents and vaccine of HPV in China. Furthermore, in-depth studies are needed to determine the clinical and biological effects of these variants.

Human papillomavirus (HPV) genotypes and HPV16 variants and risk of adenocarcinoma and squamous cell carcinoma of the cervix

OBJECTIVE

Human papillomavirus (HPV) genotypes have been extensively studied in uterine cervix squamous cell carcinoma and HPV16 variants have been found to be associated with increased cancer risk, but few reports have been published on genotype distribution and HPV16 variant prevalence in adenocarcinoma tumors. The objective of this study was to analyze viral genotypes and HPV16 intratypic variants in cervical adenocarcinoma and squamous cell carcinoma of Italian women.

METHODS

A total of 39 invasive adenocarcinoma and 132 squamous cell carcinoma were reviewed and classified according to the modified WHO classification. HPV sequences were detected by nested PCR, using the broad spectrum consensus-primer pairs MY09/MY11 and the GP5+/GP6+ system, and genotyped by nucleotide sequence analysis. The HPV16-positive cases were amplified with E6-specific oligonucleotides and amplimers subjected to direct nucleotide sequence for variant identification.

RESULTS

The prevalence rate of any HPV infection was 72% in adenocarcinoma, and 85% in cervical squamous cell carcinoma. Among the 140 HPV-positive cancer cases, a total of nine mucosal HPV genotypes (HPV16, 18, 31, 33, 35, 39, 45, 58, 82) epidemiologically classified as carcinogenic or probably carcinogenic viruses were identified. The HPV type 16 was the most common viral type representing 64% and 73% of all infections in adenocarcinoma and squamous cell carcinoma, respectively. The E6 nucleotide sequence analysis of HPV16 isolates allowed the identification of Asian American (AA) variants in 33% of adenocarcinoma and in 20% of squamous cell carcinoma suggesting their stronger association with cancer of glandular origin.

CONCLUSION

These results suggest that HPV16 has a high prevalence in both invasive adenocarcinoma and squamous cell carcinoma from Italian patients. Moreover this study confirms previous observations, summarized in a systematic review of the literature, on the increased cancer risk of HPV16 AA class in adenoglandular cancer, possibly related to their more oncogenic behavior compared to HPV16 European variants.

E6 molecular variants of human papillomavirus (HPV) type 16: an updated and unified criterion for clustering and nomenclature

Reports on taxonomic identification of E6/HPV 16 variants, don't have a worldwide, updated and unified criterion for clustering and nomenclature. Our aim was to update the existing criterion and propose a new one for clustering and nomenclature for E6/HPV 16 molecular variants based on the descriptive and comparative analyses of nucleotide sequences. A systematic search of the publications between 1991 and 2010 was carried out in PUBMED and manually. 240 E6/HPV 16 variants were identified. 157 were classified as European (E), 24 as Asian (As), 14 as Asian American (AA), 11 as North American 1 (NA 1), 19 as African 1 (Af 1) and 15 as African 2 (Af 2). Three classes were determined for the E, 3 each for the As, Af 2 and AA branches, 4 classes for the NA 1 and 6 for the Af 1 branch. This study reports our results and proposes an updated criterion for clustering and nomenclature that will be useful for E6 variant identification.

[Human papillomavirus type 16 variant analysis of upstream regulatory region and E6, E7 oncogene from cervical cancer patients in Beijing]

To investigate distributional characteristics of mutations of HPV16 upstream regulatory region (URR) and E6 and E7 oncogene in the patients with cervical cancer in Beijing and to explore the potential association between oncogenesis of cervical cancer and HPV variants in this region, cervical cancer tissue from 31 cases with positive HPV16 were subjected to regular DNA extraction procedure. The corresponding primers were then designed to amplify the target sequence of URR, E6 and E7. The PCR products were sequenced and blast analysis against GenBank was carried out to evaluate the gene mutation and identify the phylogenetic branches. Among all the cases studied, URR was found to be the most frequent mutation fragments, followed by E7, and E6 was the most conservative sequence. A total of 8 hot mutation spot was identi-fied, which were URR G7521A (100%), C7435G (96.77%), C24T (45.16%), A7729C (45.16%), G7839A (45.16%), E6 T178G (41.94%), E7 A647G (45.16%), and T846C (45.16%). The most frequent HPV 16 branch was type As (54.84%), followed by type E (45.16%). Our results suggested that the mutations of G7521A, A7729C, G7839A, T178G, T350G, A647G, and G658A were likely to be associated with the enhanced oncogenic potential of HPV16 and oncogenesis of cer-vical cancer. In Beijing area, two major branches of HPV16 were type As and E. This finding provides valuable information for HPV vaccine development and infection treatment. Type As and E variants had different distributions among various ages and clinic stage groups. It might lead to a new explanation for the getting younger trend of cervical cancer.

Association of HPV16 E6 variants with diagnostic severity in cervical cytology samples of 354 women in a US population

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.

Entire genome characterization of human papillomavirus type 16 from infected Thai women with different cytological findings

Global prevalence of human papillomavirus type 16 (HPV16) exceeds that of other types. This project has been aimed at attaining basic molecular knowledge of HPV16 by sequencing the whole genome of HPV16 isolated from Thai women at various clinical stages of disease progression. Our group analyzed seven samples of HPV16 in infected women ranging from normal to cervical cancer and discovered two critical non-synonymous changes within the coding region converting the E2-219P prototype to E2-219T in cervical cancer and the L2-269S prototype to L2-269D in CIN III, respectively. Phylogenetic analysis based on the whole genome with special emphasis on the genes E2, E6, L1, and L2 showed the Thai samples to be more closely related to the European than the non-European strains. The vaccine strain's L1 polypeptides showed close relationship to our samples. The results provide basic data for future research on cervical cancer pathogenesis and representative data of HPV16 genome in Southeast Asia.

Human papillomavirus type 16 variant analysis of E6, E7, and L1 genes and long control region in biopsy samples from cervical cancer patients in north India

High-risk human papillomaviruses (HPVs), particularly HPV types 16 and 18 (HPV-16 and HPV-18, respectively), play a cardinal role in the etiology of cervical cancer. The most prevalent type, HPV-16, shows intratypic sequence variants that are known to differ in oncogenic potential and geographic distribution. This study was designed to analyze sequence variations in E6, E7, and L1 genes and the LCR (for long control region) of HPV-16 in cervical cancer patients to identify the most prevalent and novel HPV-16 variants and to correlate them with the severity of the disease. Cervical biopsies from 60 HPV-16-positive cancer cases were analyzed by PCR and DNA sequencing. The most frequently observed variations were T350G (100%) in E6, T789C (87.5%) in E7, A6695C (54.5%) in L1, and G7521A (91.1%) in the LCR. In addition, only one novel variant (T527A) in E6 and four new variants each in L1 (A6667C, A6691G, C6906T, and A6924C) and in the LCR (C13T, A7636C, C7678T, and G7799A) were identified. While E7 was found to be highly conserved, the variant 350G of E6 was the most prevalent in all of the histopathological grades. The majority of LCR variants were found at the YY1 transcription factor binding sites. Interestingly, a complete absence of the Asian lineage and a high prevalence of European lineages in E6, E7, L1, and the LCR (85, 86.7, 67.7, and 63.3%, respectively) indicate a possible epidemiological linkage between Europe and India with regard to the dissemination of HPV-16 infections in India.

An oligoarray for the detection of human papillomavirus type 16 variants

On the basis of human papillomavirus (HPV) E6 gene mutations, there are more than five variants of HPV 16. We applied a sensitive and specific stacking hybridization assay using an oligoarray for the detection of Asian–American (AA) and European (E) (E350G) HPV 16 variants. A simple glass slide was coated with capture probes consisting of short oligonucleotide DNA sequences (7–9 mers) specific for AA and E variants. Two different regions of the E6 HPV 16 gene were amplified with a set of two primers, which were used as target DNA. These targets were preannealed with auxiliary labeled oligonucleotides and hybridized to the oligoarray in the presence of specific and complementary capture probes. Our designed array based on shorter capture probes successfully discriminated between HPV 16 AA and E variants. The present DNA oligoarray system could be useful as a reliable technique for HPV 16 detection and does not require specialized equipment; nevertheless, further intra- and interlaboratory studies are needed.

Risk for high-grade cervical intraepithelial neoplasia associated with variants of human papillomavirus types 16 and 18

BACKGROUND

Although the variant lineages of human papillomavirus (HPV) types 16 and 18 are well established, their individual associations with high-grade cervical intraepithelial neoplasia (CIN) have not been extensively evaluated.

METHODS

Study subjects were women participating in the Atypical Squamous Cells of Undetermined Significance/Low-Grade Squamous Intraepithelial Lesion Triage Study who were positive for HPV16 or HPV18 at enrollment. These women were followed every 6 months for 2 years. Viral isolates from enrollment samples were characterized by DNA sequencing and classified as variant lineages.

RESULTS

Over a 2-year study period, CIN3 was histologically diagnosed in 291 of the 779 HPV16-positive women and 47 of the 275 HPV18-positive women. Among women without CIN2-3 at enrollment, the risk of subsequent CIN3 was 2.7-fold greater for those with HPV16 African-2 [95% confidence interval (95% CI), 1.0-7.0] and 3.1-fold greater for those with HPV16 Asian American (95% CI, 1.6-6.0), compared with European variants. Relative to infection with HPV18 African variants, the risk associating subsequent CIN3 was 3.8 (95% CI, 0.9-17.2) for infection with HPV18 European variants and 4.8 (95% CI, 1.0-23.6) for infection with HPV18 Asian American variants. Similar associations were observed when the 2-year prevalence of CIN3 was used as the end point. Further, for those with HPV16 European variants, the 2-year prevalence of CIN3 was higher in White women than in African American women (P = 0.01); this trend was reversed for those with HPV16 African-1 variants (P = 0.22). A similar pattern was present for infections with HPV18 European versus African variants.

CONCLUSIONS

The lineages of HPV16 and HPV18 variants are associated with differing risks for high-grade CIN.

Molecular variants of human papillomavirus type 16 and risk for cervical neoplasia in South Africa

Non-European variants of human papillomavirus (HPV) type 16 are generally associated with a greater risk of cervical neoplasia than European prototype variants. We investigated whether this association would persist in a population in which non-European HPV 16 variants were more common. We sequenced HPV 16 isolates in cervical samples collected from 93 Black South African women enrolled in a cervical cancer screening study and examined associations between cervical neoplasia identified though colposcopy with cervical biopsy and the specific HPV 16 variant identified. The European prototype variant (EP) was the most commonly identified variant in this population (47% of all isolates), but African variants (Af-1 and Af-2) were also quite common (41% of all isolates). In contrast to previous studies, we found no evidence that non-European variants were associated with an increased risk of neoplasia. Rather, most of the HPV 16-associated cancers were found in association with EP (71% of 14 cases). In this setting where African HPV 16 variants were common, no increased risk for cervical neoplasia was found among women with these variants compared with other HPV 16 variants.

Analysis of mutations in the E6/E7 oncogenes and L1 gene of human papillomavirus 16 cervical cancer isolates from China

Human papillomavirus type 16 (HPV16) has a number of intratypic variants; each has a different geographical distribution and some are associated with enhanced oncogenic potential. Cervical samples were collected from 223 cervical cancer patients and from 196 age-matched control subjects in China. DNA samples were amplified by using primers specific for the E6, E7 and partial L1 regions. Products were sequenced and analysed. It was found by using a PCR-sequence-based typing method that HPV infection rates in China were 92.8 % in cervical cancer patients and 15.8 % in healthy controls. HPV16 was detected in 70.4 % of cervical cancer patients and in 6.1 % of controls. In HPV16-positive cervical cancers, 23.6 % belonged to the prototype, 65.5 % were of the Asian variant, 5.5 % were of African type 1 and 3.6 % were European variants, whilst only one was a new variant that differed from any variant published so far. Prevalences of HPV16 E6 D25E and E113D variants were 67.3 and 9 %, respectively. In addition to D25E and E113D, the following E6 variations were found in this study: R129K, E89Q, S138C, H78Y, L83V and F69L. The results also showed that the prevalences of three hot spots of E7 nucleotide variation, N29S, S63F and a silent variation, nt T846C, were 70.2 % (33/47), 51.1 % (24/47) and 61.7 % (29/47), respectively. The following L1 variations were found in this study: S377A, K387E, E378D, K382E and T379P. It was also found that the average age of Asian variant-positive cervical cancer patients (42.98+/-10.43 years) was 7.56 years lower than that of prototype-positive patients (50.54+/-10.91). It is suggested that the high frequency of HPV16 Asian variants might contribute to the high incidence of cervical cancer in China.

Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses

Human papillomavirus type 16 (HPV16) is the primary etiological agent of cervical cancer, the second most common cancer in women worldwide. Complete genomes of 12 isolates representing the major lineages of HPV16 were cloned and sequenced from cervicovaginal cells. The sequence variations within the open reading frames (ORFs) and noncoding regions were identified and compared with the HPV16R reference sequence. This whole-genome approach gives us unprecedented precision in detailing sequence-level changes that are under selection on a whole-viral-genome scale. Of 7,908 base pair nucleotide positions, 313 (4.0%) were variable. Within the 2,452 amino acids (aa) comprising 8 ORFs, 243 (9.9%) amino acid positions were variable. In order to investigate the molecular evolution of HPV16 variants, maximum likelihood models of codon substitution were used to identify lineages and amino acid sites under selective pressure. Five codon sites in the E5 (aa 48, 65) and E6 (aa 10, 14, 83) ORFs were demonstrated to be under diversifying selective pressure. The E5 ORF had the overall highest nonsynonymous/synonymous substitution rate (omega) ratio (M3 = 0.7965). The E2 gene had the next-highest omega ratio (M3 = 0.5611); however, no specific codons were under positive selection. These data indicate that the E6 and E5 ORFs are evolving under positive Darwinian selection and have done so in a relatively short time period. Whether response to selective pressure upon the E5 and E6 ORFs contributes to the biological success of HPV16, its specific biological niche, and/or its oncogenic potential remains to be established.

Polymorphism in the E6 gene of human papillomavirus type 16 in the cervical tissues of Korean women

The aim of this study was to identify sequence variants in the HPV 16 E6 gene in Korean women and to examine the possible association between these sequence variants and cervical cancer development. We examined the HPV 16 DNA of 215 patients with no cervical disease (NCD) (n = 105) or with cervical neoplasia (n = 110) [cervical intraepithelial neoplasia (CIN), n = 61; invasive cervical carcinoma (ICC), n = 49] using the nested polymerase chain reaction (PCR) and PCR-directed sequencing methods. Fifty-four (NCD, n = 10; CIN, n = 17; ICC, n = 27) of the 215 samples contained HPV 16 E6 DNA, but only two (7.4%) of 27 ICC samples had prototype sequences. The most frequently found variation was D25E (in NCD, n = 8, 80%; in CIN, n = 9, 52.9%; in ICC, n = 23, 85.2%). This is a rare variation in western countries. No significance difference was found between the frequencies of D25E variation in cancerous and non-cancerous lesions. Among the 11 kinds of variants identified, four variants were novel and have been registered with GenBank. This study demonstrates that the D25 variant is the most prevalent E6 genomic variant type in Korean population. However, it was not found to be associated with an increased risk of ICC.

Human papillomavirus type 16 E6 amino acid 83 variants enhance E6-mediated MAPK signaling and differentially regulate tumorigenesis by notch signaling and oncogenic Ras

Oncogenically high-risk human papillomaviruses (HPVs) are causally associated with the progression of major human neoplasia-like cancers of the cervix. Several studies have defined functions of the key E6 and E7 oncoproteins in epithelial cell immortalization. The roles of these oncogenes in the progression of immortalized epithelial cells to invasive tumors are still poorly understood. Here, we establish a novel link between the E6 oncoprotein and activation of mitogen-activated protein kinase (MAPK) signaling and show that this signaling involves Rap1. We find that activated MAPK signaling cooperates with deregulated Notch1 signaling to recreate features of HPV-driven invasive cervical carcinomas. We extend our analysis to evaluate an E6 (amino acid [aa] 83) variant that has been linked to invasive tumors. The variant enhances MAPK signaling and cooperative transformation with deregulated Notch1 signaling. Unlike E6, this variant surprisingly inhibits oncogenic Ras-mediated transformation. Our data reveal that the quantitative differences in activation of MAPK signaling by E6 and its variant correlate with differences in cooperative transformation with other signaling pathways, thus suggesting that thresholds of MAPK activation may define permissive conditions for other signaling pathways in tumorigenesis. Epidemiological studies have suggested the importance of E6 aa 83 variants in invasive carcinomas; our data support a key deterministic role for this variant in human cervical tumorigenesis. These observations, along with our recent data showing that deregulated Notch signaling activates phosphatidylinositol 3-kinase signaling, strengthen the possibility of the existence of Ras-independent mechanisms to recreate signaling through classical Ras effector pathways.

p53-Dependent Apoptotic Mechanism of a New Designer Bimetallic Compound Tri-phenyl Tin Benzimidazolethiol Copper Chloride (TPT-CuCl2): In Vivo Studies in Wistar Rats as Well as in Vitro Studies in Human Cervical Cancer Cells

We have studied the effect of tri-phenyl tin benzimadazolethiolcopper chloride (TPT-CuCl(2)), a novel bimetallic compound, on the regulation of apoptosis in HeLa cells, MCF-7 cells, and in vivo Wistar rat model. TPT-CuCl(2) induces significant apoptosis in HeLa cell line characterized by DNA fragmentation and chromosome condensation. Comet assay revealed that TPT-CuCl(2) targets and causes severe damage to the DNA. Treatment of HeLa cells with TPT-CuCl(2) rescues the accumulation of p53 from the suppression of human papilloma virus E6, resulting in a dramatic up-regulation of Bax and Bak and down-regulation of the antiapoptotic factor Survivin. Apoptotic induction by TPT-CuCl(2) was shown to mediate in a p53-depedent manner; loss of p53 impairs the release of cytochrome c and Smac/DIABLO from mitochondria to cytosol. Moreover, we have shown that TPT-CuCl(2) induced-apoptosis was through an intrinsic mitochondrial pathway, which was inhibited by viral oncoprotein E1B19K. Caspase-3 was found to be indispensable in TPT-CuCl(2)-triggered apoptosis signaling pathway, because caspase-3 deficient cell line MCF-7 was resistant to TPT-CuCl(2). Furthermore, in vivo studies using C6 glioblastoma xenograft rat model revealed that TPT-CuCl(2) exhibits significant antiproliferative activity against tumor development with minimal cytotoxicity toward normal physiological function of the experimental rats. These findings imply the attractiveness of TPT-CuCl(2) as a drug candidate for further development.

Distribution of human papillomavirus type 16 variants in human immunodeficiency virus type 1-positive and -negative women

The prevalence of human papillomavirus type 16 E6 variant lineages was characterized in a cross-sectional study of 24 human immunodeficiency virus type 1 (HIV)-positive and 33 HIV-negative women in New Orleans. The European prototype was the predominant variant in the HIV-negative women (39·4 %), while in the HIV-positive women the European 350G variant was predominant (29·1 %). In exact logistic regression models, HIV-positive women were significantly more likely to harbour any variant with a nucleotide G-350 mutation compared with HIV-negative women [58·3 % vs 21·1 %; adjusted odds ratio (AOR)=6·28, 95 % confidence interval (CI)=1·19–46·54]. Models also revealed a trend towards increased prevalence of Asian–American lineage in HIV-positive women compared with HIV-negative women (25·0 % vs 6·0 %; AOR=6·35, 95 % CI=0·77–84·97). No association was observed between any variant and cytology or CD4 cell counts or HIV-1 viral loads. These observations reflect a difference in the distribution of HPV-16 variants among HIV-positive and -negative women, indicating that HIV-positive status may lead to increased prevalence of a subset of variants.

Human papillomavirus type 16 E6 variants and HLA class II alleles among Japanese women with cervical cancer

The enhanced oncogenicity of particular human papillomavirus type 16 (HPV16) E6 variants is population-dependent, implying the involvement of additional genetic cofactors. This study was designed to investigate the association between E6 variants and human leukocyte antigen (HLA) polymorphism within a Japanese population. Fifty-seven women with HPV16-positive cervical cancer were analyzed for E6 sequence variation and its relationship to HLA class II alleles. Compared with local controls (n = 138) and published controls (n = 916), DRB1*1501 and DQB1*0602 frequencies were significantly increased among patients with HPV16 E6 prototype (n = 11). Additionally, DRB1*1502 was positively associated with a particular E6 variant designated D25E (n = 25), although we could not find a significant association between HLA class II alleles and L83V variants (n = 16). Our observations suggest that a specific match between E6 variant proteins and HLA types may contribute to HPV16-related cervical carcinogenesis.

Human papillomavirus type-16 variants in Quechua aboriginals from Argentina

Cervical carcinoma is the leading cause of cancer death in Quechua indians from Jujuy (northwestern Argentina). To determine the prevalence of HPV-16 variants, 106 HPV-16 positive cervical samples were studied, including 33 low-grade squamous intraepithelial lesions (LSIL), 28 high-grade squamous intraepithelial lesions (HSIL), 9 invasive cervical cancer (ICC), and 36 samples from women with normal colposcopy and cytology. HPV genome variability was examined in the L1 and E6 genes by PCR-hybridization. In a subset of 20 samples, a LCR fragment was also analyzed by PCR-sequencing. Most variants belonged to the European branch with subtle differences that depended on the viral gene fragment studied. Only about 10% of the specimens had non-European variants, including eight Asian-American, two Asian, and one North-American-1. E6 gene analysis revealed that 43% of the samples were identical to HPV-16 prototype, while 57% corresponded to variants. Interestingly, the majority (87%) of normal smears had HPV-16 prototype, whereas variants were detected mainly in SIL and ICC. LCR sequencing yielded 80% of variants, including 69% of European, 19% Asian-American, and 12% Asian. We identified a new variant, the Argentine Quechua-51 (AQ-51), similar to B-14 plus two additional changes: G7842-->A and A7837-->C; phylogenetic inference allocated it in the Asian-American branch. The high proportion of European variants may reflect Spanish colonial influence on these native Inca descendants. The predominance of HPV-16 variants in pathologic samples when compared to normal controls could have implications for the natural history of cervical lesions.

Human papillomavirus type 16 E6 and E7 gene variations in Indian cervical cancer

OBJECTIVES

Human papillomavirus type 16 is a causative factor for development of cervical cancer. The E6 and E7 genes of HPV 16 are critical to the process of immortalization and transformation of host cells. Recent reports suggest that variants of these two genes may contribute to the risk of malignant progression of cancer in the uterine cervix. However, no data exist on sequence variations of HPV 16 E6 and E7 genes that may exist in India. Therefore, we examined intratype variations in the E6 and E7 viral genes in DNA isolated from HPV 16-positive cervical scrapes and biopsies.

METHODS

The open reading frames of the E6 and E7 genes were amplified by PCR and then directly sequenced by the fluorescent dye dideoxy termination method.Results. In addition to the prototype E6 gene sequence, five sets of mutations of the E6 gene were identified. The European prototype (350T) was detected in 9.1% of the study group while the European variant (350G) was seen in 28% of patients. The remaining variants (a combination of the 350G mutation with 335T, 145T, or 419G) were significantly associated with cases compared to controls. The 350G + 145T variant was found at much higher incidence in cases in younger women, suggesting that this variant may be associated with aggressive tumor behavior. Interestingly the 350G + 419G combination was found only in controls. There was no significant association between the four genotypes of E7 and any stage of tumor progression or age.

CONCLUSIONS

The results indicate that specific mutations in the E6 gene are found in young Indian women with high-grade squamous intraepithelial lesions and invasive cancer, suggesting that these mutations represent more oncogenically active HPV 16. Whether this increased oncogenecity is due to differences in p53 inactivation, ineffective keratinocyte differentiation, and/or altered response to the immune system by these oncogenic E6 mutants remains to be clarified.

Host and viral genetics and risk of cervical cancer: a review

Infection with human papillomaviruses (HPV) is known to play a central role in the development of cervical cancer. Both host and viral genetic factors have been postulated to be important determinants of risk of HPV progression to neoplasia among infected individuals. In this report, we review epidemiological studies that have evaluated the role in cervical cancer pathogenesis of genetic variation in human leukocyte antigen (HLA) genes and in the HPV genome itself. A protective effect of HLA Class II DRB1*13/DBQ1*0603 alleles is the most consistent HLA finding in the published literature. A consistent association between HPV16 non-European variants and risk of disease is also evident from published work. These findings are discussed. Gaps in our understanding and future research needs are also discussed.