Risk of cervical cancer and geographical variations of human papillomavirus 16 E6 polymorphisms

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.

HPV Testing and its Role in Cervical Cancer Screening

The recognition that persistent infection with carcinogenic human papillomavirus (HPV) is a necessary cause of cervical precancer and cancer has led to the introduction of HPV testing into cervical cancer screening, either as a primary screening test or in conjunction with cervical cytology (i.e., co-testing). HPV testing has much higher sensitivity for detection of cervical precancer and provides greater long-term reassurance if negative compared to cytology. However, most HPV infections are transient, and do not progress to invasive cancer, thus triage tests are required to identify individuals who should be referred to colposcopy for diagnostic evaluation. This chapter begins with a description of the biology, natural history, and epidemiology of HPV as a foundation for understanding the role of HPV in cervical carcinogenesis. This section is followed by a detailed discussion regarding the introduction of HPV-based testing and triage into cervical cancer screening and management. Summarized triage tests include cervical cytology, HPV genotyping, p16/Ki-67 dual stain, and HPV and cellular methylation markers. The final section of this chapter includes an important discussion on cervical cancer disparities, particularly within the United States, followed by concluding remarks.

Human papillomavirus genomics: Understanding carcinogenicity

Human papillomavirus (HPV) causes virtually all cervical cancers and many cancers at other anatomical sites in both men and women. However, only 12 of 448 known HPV types are currently classified as carcinogens, and even the most carcinogenic type - HPV16 - only rarely leads to cancer. HPV is therefore necessary but insufficient for cervical cancer, with other contributing factors including host and viral genetics. Over the last decade, HPV whole genome sequencing has established that even fine-scale within-type HPV variation influences precancer/cancer risks, and that these risks vary by histology and host race/ethnicity. In this review, we place these findings in the context of the HPV life cycle and evolution at various levels of viral diversity: between-type, within-type, and within-host. We also discuss key concepts necessary for interpreting HPV genomic data, including features of the viral genome; events leading to carcinogenesis; the role of APOBEC3 in HPV infection and evolution; and methodologies that use deep (high-coverage) sequencing to characterize within-host variation, as opposed to relying on a single representative (consensus) sequence. Given the continued high burden of HPV-associated cancers, understanding HPV carcinogenicity remains important for better understanding, preventing, and treating cancers attributable to infection.

Genetic variation within the human papillomavirus type 16 genome is associated with oropharyngeal cancer prognosis

PURPOSE

A significant barrier to adoption of de-escalated treatment protocols for human papillomavirus-driven oropharyngeal cancer (HPV-OPC) is that few predictors of poor prognosis exist. We conducted the first large whole-genome sequencing (WGS) study to characterize the genetic variation of the HPV type 16 (HPV16) genome and to evaluate its association with HPV-OPC patient survival.

PATIENTS AND METHODS

A total of 460 OPC tumor specimens from two large United States medical centers (1980-2017) underwent HPV16 whole-genome sequencing. Site-specific variable positions [single nucleotide polymorphisms (SNPs)] across the HPV16 genome were identified. Cox proportional hazards model estimated hazard ratios (HRs) and 95% confidence intervals (CIs) for overall survival by HPV16 SNPs. Harrell C-index and time-dependent positive predictive value (PPV) curves and areas under the PPV curves were used to evaluate the predictive accuracy of HPV16 SNPs for overall survival.

RESULTS

A total of 384 OPC tumor specimens (83.48%) passed quality control filters with sufficient depth and coverage of HPV16 genome sequencing to be analyzed. Some 284 HPV16 SNPs with a minor allele frequency ≥1% were identified. Eight HPV16 SNPs were significantly associated with worse survival after false discovery rate correction (individual prevalence: 1.0%-5.5%; combined prevalence: 15.10%); E1 gene position 1053 [HR for overall survival (HRos): 3.75, 95% CI 1.77-7.95; Pfdr = 0.0099]; L2 gene positions 4410 (HRos: 5.32, 95% CI 1.91-14.81; Pfdr = 0.0120), 4539 (HRos: 6.54, 95% CI 2.03-21.08; Pfdr = 0.0117); 5050 (HRos: 6.53, 95% CI 2.34-18.24; Pfdr = 0.0030), and 5254 (HRos: 7.76, 95% CI 2.41-24.98; Pfdr = 0.0030); and L1 gene positions 5962 (HRos: 4.40, 95% CI 1.88-10.31; Pfdr = 0.0110) and 6025 (HRos: 5.71, 95% CI 2.43-13.41; Pfdr = 0.0008) and position 7173 within the upstream regulatory region (HRos: 9.90, 95% CI 3.05-32.12; Pfdr = 0.0007). Median survival time for patients with ≥1 high-risk HPV16 SNPs was 3.96 years compared with 18.67 years for patients without a high-risk SNP; log-rank test P < 0.001. HPV16 SNPs significantly improved the predictive accuracy for overall survival above traditional factors (age, smoking, stage, treatment); increase in C-index was 0.069 (95% CI 0.019-0.119, P < 0.001); increase in area under the PPV curve for predicting 5-year survival was 0.068 (95% CI 0.015-0.111, P = 0.008).

CONCLUSIONS

HPV16 genetic variation is associated with HPV-OPC prognosis and can improve prognostic accuracy.

Prevalence and distribution of human papillomavirus (HPV) in Luoyang city of Henan province during 2015-2021 and the genetic variability of HPV16 and 52

Background

Persistent high-risk Human papillomavirus (HPV) subtypes infection has been implicated as a causative of cervical cancer. Distribution and genotypes of HPV infection among females and their variations would assist in the formulation of preventive strategy for cervical cancer. The purpose of the present study is to investigate the prevalence of HPV among females in central China.

Methods

The distribution and genotypes of HPV among 9943 females attending the gynecological examinations in central of China during 2015–2021 were investigated. HPV genotypes were detected using a commercial kit. Nucleotides sequences of L1, E6 and E7 genes in HPV16 or HPV52 positive samples collected in 2021 were amplified by polymerase chain reaction (PCR). Variations of L1, E6 and E7 in HPV16 and HPV52 were gained by sequencing and compared with the reference sequence. Sublineages of HPV16 and HPV52 were determined by the construction of phylogenetic tree based on L1 gene.

Results

The overall prevalence of HPV infection was 22.81%, with the infection rate of high-risk human papillomavirus (HR-HPV) was 19.02% and low-risk human papillomavirus (LR-HPV) was 6.40%. The most top five genotypes of HPV infection were HPV16 (7.49%), HPV52 (3.04%), HPV58 (2.36%), HPV18 (1.65%) and HPV51 (1.61%). Plots of the age-infection rate showed that the single HPV, multiple HPV, HR-HPV, LR-HPV infection revealed the same tendency with two peaks of HPV infection were observed among females aged ≤ 20 year-old and 60–65 year-old. The predominant sublineage of HPV16 was A1 and B2 for HPV52. For HPV16, The most prevalent mutations were T266A (27/27) and N181T (7/27) for L1, D32E for E6 and S63F for E7 in HPV16. For HPV52, all of the nucleotide changes were synonymous mutation in L1 (except L5S) and E7 genes. The K93R mutation was observed in most HPV52 E6 protein.

Conclusions

The present study provides basic information about the distribution, genotypes and variations of HPV among females population in Henan province, which would assist in the formulation of preventive strategies and improvements of diagnostic probe and vaccine for HPV in this region.

Repurposing the Hybrid Capture 2 (HC2) screening test for whole-genome sequencing of human papillomaviruses

Simple and standardized approaches for genome analysis of human papillomavirus (HPV) by next-generation sequencing are needed. The aim of the study was to develop a protocol for direct deep sequencing of high-risk (hr) HPV strains, based on the widely used commercial Hybrid Capture 2 (QIAGEN) test, without any additional probe design. This protocol was applied to 15 HPV-positive and two HPV-negative cervical samples or cell lines and validated at the genotype level by comparing the sequencing results to those obtained using a commercial genotyping kit. The performance of our protocol, presented in this proof-of-principle study, supports its use for accurate characterization of genetic variants of hrHPV.

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.

APOBEC Mutagenesis Is Concordant between Tumor and Viral Genomes in HPV-Positive Head and Neck Squamous Cell Carcinoma

APOBEC is a mutagenic source in human papillomavirus (HPV)-mediated malignancies, including HPV+ oropharyngeal squamous cell carcinoma (HPV + OPSCC), and in HPV genomes. It is unknown why APOBEC mutations predominate in HPV + OPSCC, or if the APOBEC-induced mutations observed in both human cancers and HPV genomes are directly linked. We performed sequencing of host somatic exomes, transcriptomes, and HPV16 genomes from 79 HPV + OPSCC samples, quantifying APOBEC mutational burden and activity in both host and virus. APOBEC was the dominant mutational signature in somatic exomes. In viral genomes, there was a mean of five (range 0-29) mutations per genome. The mean of APOBEC mutations in viral genomes was one (range 0-5). Viral APOBEC mutations, compared to non-APOBEC mutations, were more likely to be low-variant allele fraction mutations, suggesting that APOBEC mutagenesis actively occurrs in viral genomes during infection. HPV16 APOBEC-induced mutation patterns in OPSCC were similar to those previously observed in cervical samples. Paired host and viral analyses revealed that APOBEC-enriched tumor samples had higher viral APOBEC mutation rates (p = 0.028), and APOBEC-associated RNA editing (p = 0.008), supporting the concept that APOBEC mutagenesis in host and viral genomes is directly linked and occurrs during infection. Using paired sequencing of host somatic exomes, transcriptomes, and viral genomes, we demonstrated for the first-time definitive evidence of concordance between tumor and viral APOBEC mutagenesis. This finding provides a missing link connecting APOBEC mutagenesis in host and virus and supports a common mechanism driving APOBEC dysregulation.

Metabolic Reprogramming in Cancer: Role of HPV 16 Variants

Metabolic reprogramming is considered one of the hallmarks in cancer and is characterized by increased glycolysis and lactate production, even in the presence of oxygen, which leads the cancer cells to a process called “aerobic glycolysis” or “Warburg effect”. The E6 and E7 oncoproteins of human papillomavirus 16 (HPV 16) favor the Warburg effect through their interaction with a molecule that regulates cellular metabolism, such as p53, retinoblastoma protein (pRb), c-Myc, and hypoxia inducible factor 1α (HIF-1α). Besides, the impact of the E6 and E7 variants of HPV 16 on metabolic reprogramming through proteins such as HIF-1α may be related to their oncogenicity by favoring cellular metabolism modifications to satisfy the energy demands necessary for viral persistence and cancer development. This review will discuss the role of HPV 16 E6 and E7 variants in metabolic reprogramming and their contribution to developing and preserving the malignant phenotype of cancers associated with HPV 16 infection.

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.

Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance

HPV16 causes half of cervical cancers worldwide; for unknown reasons, most infections resolve within two years. Here, we analyze the viral genomes of 5,328 HPV16-positive case-control samples to investigate mutational signatures and the role of human APOBEC3-induced mutations in viral clearance and cervical carcinogenesis. We identify four de novo mutational signatures, one of which matches the COSMIC APOBEC-associated signature 2. The viral genomes of the precancer/cancer cases are less likely to contain within-host somatic HPV16 APOBEC3-induced mutations (Fisher's exact test, P = 6.2 x 10-14), and have a 30% lower nonsynonymous APOBEC3 mutation burden compared to controls. We replicate the low prevalence of HPV16 APOBEC3-induced mutations in 1,749 additional cases. APOBEC3 mutations also historically contribute to the evolution of HPV16 lineages. We demonstrate that cervical infections with a greater burden of somatic HPV16 APOBEC3-induced mutations are more likely to be benign or subsequently clear, suggesting they may reduce persistence, and thus progression, within the host.

High prevalent human papillomavirus infections of the oral cavity of asymptomatic HIV-positive men

Background

Incidence of anal and oral infections with Human Papillomavirus (HPV) is increasing, particularly among Human Immunodeficiency Virus-positive (HIV+) men. HPV type 16 has exhibited the highest incidence and only limited data is available on other prevalent types, variants of HPV16, as well as associated factors. We were interested in identifying prevalent HPV types, variants of type 16, as well as factors associated with HPV16 infections in the oral cavity of HIV+ men who have sex with men (MSM).

Methods

A cross-sectional study of oral cavity samples from HIV+ MSM, that in a previous study were identified as positive for HPV16 in the anal canal. Cells from the oral cavity (102 samples, paired with 102 from the anal canal of same patient) were used to extract DNA and detect HPV infections using INNO-LiPA HPV Genotyping Extra II, and PCR. From these, 80 samples (paired, 40 anal and 40 oral) were used to identify variants of type 16 by sequencing. Statistical differences were estimated by the X² test, and p values equal to or less than 0.05 were considered significant. SPSS ver. Twenty-four statistical software (IBM Corp) was used.

Results

We found a high prevalence of High-Risk HPV (HR-HPV) and Low-Risk HPV (LR-HPV). Patients were positive in the oral cavity for HR types; 16, 39 and 18 (80.4, 61.8 and 52.9% respectively) and LR types 11 and 6 (53.9 and 34.3% respectively). Surprisingly, only European variants of type 16 were found in the oral cavity, although American Asian (22.5%) and African (2.5%) variants were identified in the anal canal. The analysis showed that CD4 counts could be the most important risk factor associated with HR-HPV infections in the oral cavity, anal canal or both anatomical regions. The risk of infection of the oral cavity with type 18 increased in men diagnosed with HIV for more than 6 years.

Conclusions

Prevalence of both HR and LR HPV’s in the oral cavity of Mexican HIV+ MSM is very high. The fact that only European variants of HPV16 were found in the oral cavity suggest a possible tropism not previously described.

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.

Human papillomavirus 16 sub-lineage dispersal and cervical cancer risk worldwide: Whole viral genome sequences from 7116 HPV16-positive women

BACKGROUND

Human papillomavirus (HPV)16 can be separated into genetic sub-lineages (A1-4, B1-4, C1-4, D1-4) which may have differential cervical cancer risk.

METHODS

A next-generation sequencing assay was used to whole-genome sequence 7116 HPV16-positive cervical samples from well-characterised international epidemiological studies, including 2076 controls, 1878 squamous cell carcinoma (SCC) and 186 adenocarcinoma/adenosquamous cell carcinoma (ADC), and to assign HPV16 sub-lineage. Logistic regression was used to estimate region-stratified country-adjusted odds ratios (OR) and 95%CI.

RESULTS

A1 was the most globally widespread sub-lineage, with others showing stronger regional specificity (A3 and A4 for East Asia, B1-4 and C1-4 for Africa, D2 for the Americas, B4, C4 and D4 for North Africa). Increased cancer risks versus A1 were seen for A3, A4 and D (sub)lineages in regions where they were common: A3 in East Asia (OR=2.2, 95%CI:1.0-4.7); A4 in East Asia (6.6, 3.1-14.1) and North America (3.8, 1.7-8.3); and D in North (6.2, 4.1-9.3) and South/Central America (2.2, 0.8-5.7), where D lineages were also more frequent in ADC than SCC (3.2, 1.5-6.5; 12.1, 5.7-25.6, respectively).

CONCLUSIONS

HPV16 genetic variation can strongly influence cervical cancer risk. However, burden of cervical cancer attributable to different sub-lineages worldwide is largely driven by historical HPV16 sub-lineage dispersal.

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.

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.

Eurogin roadmap 2017: Triage strategies for the management of HPV-positive women in cervical screening programs

Cervical cancer screening will rely, increasingly, on HPV testing as a primary screen. The requirement for triage tests which can delineate clinically significant infection is thus prescient. In this EUROGIN 2017 roadmap, justification behind the most evidenced triages is outlined, as are challenges for implementation. Cytology is the triage with the most follow-up data; the existence of an HR-HPV-positive, cytology-negative group presents a challenge and retesting intervals for this group (and choice of retest) require careful consideration. Furthermore, cytology relies on subjective skills and while adjunctive dual-staining with p16/Ki67 can mitigate inter-operator/-site disparities, clinician-taken samples are required. Comparatively, genotyping and methylation markers are objective and are applicable to self-taken samples, offering logistical advantages including in low and middle income settings. However, genotyping may have diminishing returns in immunised populations and type(s) included must balance absolute risk for disease to avoid low specificity. While viral and cellular methylation markers show promise, more prospective data are needed in addition to refinements in automation. Looking forward, systems that detect multiple targets concurrently such as next generation sequencing platforms will inform the development of triage tools. Additionally, multistep triage strategies may be beneficial provided they do not create complex, unmanageable pathways. Inevitably, the balance of risk to cost(s) will be key in decision making, although defining an acceptable risk will likely differ between settings. Finally, given the significant changes to cervical screening and the variety of triage strategies, appropriate education of both health care providers and the public is essential.

Human Papillomavirus Genotypes and HPV16 E6/E7 Variants among Patients with Genital Cancers in Vietnam

We previously reported human papillomavirus type 52 (HPV52) as the most prevalent high-risk genotype in non-cancer individuals in Vietnam. This study aimed to evaluate HPV genotypes and HPV16 E6 and E7 (E6/E7) gene variations in Vietnamese patients with genital cancers. Biopsy samples were collected from 124 Vietnamese patients with genital cancers (20 with vaginal, 50 with vulvar, and 54 with penile cancer). The HPV-DNA was amplified and genotyped, and HPV16 E6/E7 genes were compared with those previously reported for women with normal cervical cytology (N = 23). HPV-DNA was detected in 80.6% (100/124) of the cancer patients (80.0% of vaginal, 82.0% of vulvar, and 79.6% of penile), with HPV16/18 in 86.0% (86/100) and HPV52 in 7.0% (7/100) of the HPV-positive samples. The HPV-DNA prevalence and HPV genotype distribution did not significantly differ among the genital cancer patients (both P = 0.95). Significantly fewer instances of the HPV16 A4 sublineage (34.8% vs. 82.6%, P < 0.0001) and HPV16 E7 29S (36.4% vs. 87.0%, P = 0.0002) occurred in the cancer patients than in the women with normal cytology. Our results indicate that HPV16/18 accounts for more than 85% of genital cancers in Vietnam, and the HPV16 sublineage A4 containing E7 29S may be less oncogenic.

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.

The Intersection of HPV Epidemiology, Genomics and Mechanistic Studies of HPV-Mediated Carcinogenesis

Of the ~60 human papillomavirus (HPV) genotypes that infect the cervicovaginal epithelium, only 12–13 “high-risk” types are well-established as causing cervical cancer, with HPV16 accounting for over half of all cases worldwide. While HPV16 is the most important carcinogenic type, variants of HPV16 can differ in their carcinogenicity by 10-fold or more in epidemiologic studies. Strong genotype-phenotype associations embedded in the small 8-kb HPV16 genome motivate molecular studies to understand the underlying molecular mechanisms. Understanding the mechanisms of HPV genomic findings is complicated by the linkage of HPV genome variants. A panel of experts in various disciplines gathered on 21 November 2016 to discuss the interdisciplinary science of HPV oncogenesis. Here, we summarize the discussion of the complexity of the viral–host interaction and highlight important next steps for selected applied basic laboratory studies guided by epidemiological genomic findings.

Phylogenetic analysis of Human papillomavirus 16 variants isolated from Indian Breast cancer patients showed difference in genetic diversity with that of cervical cancer isolates

The genetic variations of HPV16 in Breast Cancer (BC) are not well studied unlike HPV16 in Cervical Cancer (CACX). In this study, the genetic variations of HPV16 in BC were compared with HPV16 in CACX. In sequencing analysis of LCR, E6 and E7 regions of HPV16 in BC and CACX the A lineage was seen to be 64.2% and 66.6% respectively. The other lineages showed differential frequency in BC and CACX. The mutation frequency index of the regions in BC and CACX was in the following order: LCR>E6>E7. However, the inter-patient genetic diversity in LCR and E6/E7 regions was high in BC than CACX. The LCR region showed more variations than the E6/E7 region in BC. Apart from some common variations, some unique tissue specific variants in LCR and E6/E7 region were seen in BC and in CACX. Besides the selection of some common variants in both BC and CACX, some unique variants in BC (D98Y; 395 G>T) and CACX (R48W; 245 G>T) were observed. The 7521 G>A variant of LCR showed association with Luminal B subtype of BC and progression of CACX. Whereas, 145 G>T (Q14H) and 335 C>T (H78Y) variants of E6 showed association with either early invasiveness of BC and/or poor outcome of the patients. Thus, this study indicates that there may be a difference in the genetic variation of HPV16 in BC and in CACX.

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.

Genomic diversity and phylogenetic relationships of human papillomavirus 16 (HPV16) in Nepal

OBJECTIVES/BACKGROUND

Sequence variants in HPV16 confer differences in oncogenic potential; however, to date there have not been any HPV sequence studies performed in Nepal. The objective of this study was to characterize HPV16 viral genome sequences from Nepal compared to a reference sequence in order to determine their lineages. Additionally, we sought to determine if five High-grade Squamous Intraepithelial Lesion (HSIL) subjects were genetically distinct from the non-HSIL subjects.

METHODS

DNA was isolated from exfoliated cervical cells from 17 individuals in Nepal who were previously identified to be HPV16-positive. A custom HPV16 Ion Ampliseq panel of multiplexed degenerate primers was designed that generated 47 overlapping amplicons and covered 99% of the viral genome for all known HPV16 variant lineages. All sequence data were processed through a custom quality control and analysis pipeline of sequence comparisons and phylogenetic analysis.

RESULTS

There were high similarities across the genomes, with two major indels observed in the non-coding region between E5 and L2. Compared to the PAVE reference HPV16 genome, there were up to 9, 4, 38, 27, 8, 7, 52, and 32 nucleotide variants in the E6, E7, E1, E2, E4, E5, L2, and L1 genes in the Nepalese samples, respectively. Based on sequence variation, HPV16 from Nepal falls across the A, C, and D lineages in this study. We found no evidence of genetic distinctness between HSIL and non-HSIL subjects.

CONCLUSIONS

The evolutionary and pathological characteristics of the representative HPV16 genomes from Nepal seem similar to results from other parts of the world and provide the basis for further studies.

HPV16 Sublineage Associations With Histology-Specific Cancer Risk Using HPV Whole-Genome Sequences in 3200 Women

BACKGROUND

HPV16 is a common sexually transmitted infection although few infections lead to cervical precancer/cancer; we cannot distinguish nor mechanistically explain why only certain infections progress. HPV16 can be classified into four main evolutionary-derived variant lineages (A, B, C, D) that have been previously suggested to have varying disease risks.

METHODS

We used a high-throughput HPV16 whole-genome sequencing assay to investigate variant lineage risk among 3215 HPV16-infected women. Using sublineages A1/A2 as the reference, we assessed all variant lineage associations with infection outcome over three or more years of follow-up: 1107 control subjects (

RESULTS

A4 sublineage was associated with an increased risk of cancer, specifically adenocarcinoma (OR = 9.81, 95% CI = 2.02 to 47.69, P = 4.7x10(-03)). Lineage B had a lower risk of CIN3 (OR = 0.51, 95% CI = 0. 28 to 0.91, P = 02) while lineage C showed increased risk (OR = 2.06, 95% CI = 1.09 to 3.89, P = 03). D2/D3 sublineages were strongly associated with an increased risk of CIN3 and cancer, particularly D2 (OR for cancer = 28.48, 95% CI = 9.27 to 87.55, P = 5.0x10(-09)). D2 had the strongest increased risk of glandular lesions, AIS (OR = 29.22, 95% CI = 8.94 to 95.51, P = 2.3x10(-08)), and adenocarcinomas (OR = 137.34, 95% CI = 37.21 to 506.88, P = 1.5x10(-13)). Moreover, the risk of precancer and cancer for specific variant lineages varied by a women's race/ethnicity; those women whose race/ethnicity matched that of the infecting HPV16 variant had an increased risk of CIN3 + (P < 001).

CONCLUSIONS

Specific HPV16 variant sublineages strongly influence risk of histologic types of precancer and cancer, and viral genetic variation may help explain its unique carcinogenic properties.

Collapse

Key Words Collapse

MESH Headings

• Adenocarcinoma/ethnology
• Adenocarcinoma/virology
• Adenocarcinoma in Situ/ethnology
• Adenocarcinoma in Situ/virology
• Adult
• Aftercare
• California/epidemiology
• Carcinoma/ethnology
• Carcinoma/virology
• Carcinoma, Squamous Cell/ethnology
• Carcinoma, Squamous Cell/virology
• Female
• Genome
• High-Throughput Nucleotide Sequencing
• Human papillomavirus 16/classification
• Human papillomavirus 16/genetics
• Human papillomavirus 16/pathogenicity
• Humans
• Middle Aged
• Papillomavirus Infections/ethnology
• Papillomavirus Infections/virology
• Phylogeny
• Precancerous Conditions/ethnology
• Precancerous Conditions/virology
• Risk Factors
• Uterine Cervical Neoplasms/ethnology
• Uterine Cervical Neoplasms/virology
• Young Adult
• Uterine Cervical Dysplasia/ethnology
• Uterine Cervical Dysplasia/virology

Collapse

Grants

• U01 CA078527 NCI NIH HHS

Collapse

Affiliation(s)

Lisa Mirabello Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Meredith Yeager Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Michael Cullen Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Joseph F Boland Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Zigui Chen Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Nicolas Wentzensen Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Xijun Zhang Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Kai Yu Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Qi Yang Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Jason Mitchell Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
David Roberson Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Sara Bass Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Yanzi Xiao Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Laurie Burdett Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Tina Raine-Bennett Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Thomas Lorey Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Philip E Castle Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Robert D Burk Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)
Mark Schiffman Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Rockville, MD (LM, MY, MC, JFB, NW, XZ, KY, QY, JM, DR, SB, LB, YX, MS); Cancer Genomics Research Laboratory, Leidos Biomedical Research, Inc., Frederick, MD (MY, MC, JFB, XZ, QY, JM, DR, SB, LB); Department of Microbiology, The Chinese University of Hong Kong, Hong Kong (ZC); Women's Health Research Institute, Division of Research, Kaiser Permanente Northern California, Oakland CA (TRB); Regional Laboratory, Kaiser Permanente Northern California, Oakland, CA (TL); Department of Epidemiology and Population Health, at Albert Einstein College of Medicine, Bronx, NY (PEC, RDB); Departments of Pediatrics; Microbiology & Immunology; and, Obstetrics, Gynecology and Women's Health, at Albert Einstein College of Medicine, Bronx, NY (RDB)

Collapse

Variants of human papillomaviruses 16 (HPV16) in Uigur women in Xinjiang, China

BACKGROUND

Persistent infection of high-risk human papillomaviruses 16 (HPV16) has been considered as the leading cause of cervical cancer. In this study we assessed HPV16 sequence variation and genetic diversity of HPV16 variants in cervical cancer in Uigur women in Xinjiang, China. We analyzed the nucleotide sequences of the open reading frames of E6 and E7, and part of the open reading frames of L1 of HPV16 in Uigur women.

METHODS

Biopsies of histologically confirmed HPV16 infections with cervical cancer were obtained from 43 Uigur women in Xinjiang, China. E6, E7 and L1 genes of HPV16 of all samples were amplified and sequenced; the sequences were used in phylogenetic analysis of HPV16 variants.

RESULTS

Our analysis revealed nine nucleotide changes in E6 (five changes), E7 (one change) and L1 (three changes) gene. The most frequently observed variations were T350G (79.1 %). One variation T295G (D64E) at E6 were detected in 6 cases (KT959536, KT959542, KT959546, KT959550, KT959553, KT959558). Deletion (464Asp) along with insertion (448Ser) were observed in L1 (100 %). Most variants were European lineage (97.7 %); only one belongs to Asia variants with common T178G (D25E) in E6 and A647G (N29S) in E7.

CONCLUSION

The most prevalent HPV16 variants in the Uigur women we studied were of the European lineage. Our results indicate that HPV16 European lineage may serve as a harmful factor associated with the development and progression of cervical cancer.

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

Human papillomavirus type 16 is commonly implicated in HPV-related cancers. However, only a small number of infected individuals progress to this stage. Epidemiological evidence demonstrated that oncogenic risk is population-specific and variations within the viral oncogene, E6, have been suggested to play a role in these findings. Of focus in this study is the European-T350G variant, which is characterized by an L>V amino acid substitution at residue 83 of the prototype E6 protein. To elucidate the functional effects of this polymorphism, we followed keratinocytes transduced with E-T350G E6 for over 60 passages and compared them to keratinocytes transduced, in parallel, with prototype or Asian-American (Q14H/L83V/H78Y) E6. We found that although E-T350G E6 immortalized transduced keratinocytes in the absence of E7, these cells were not fully transformed. We also found that E-T350G down-regulated E-cadherin compared to the other variants, providing a possible link between its population-based oncogenicity and host genetic variations.

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.

Tumourigenesis driven by the human papillomavirus type 16 Asian-American e6 variant in a three-dimensional keratinocyte model

Infection with a transforming human papillomavirus (HPV) such as type 16 (of species Alphapapillomavirus 9) causes ano-genital and oral tumours via viral persistence in human squamous cell epithelia. Epidemiological studies showed that the naturally occurring HPV16 Asian-American (AA) variant (sublineage D2/D3) is found more often than the European Prototype (EP) (sublineage A1) in high-grade cervical neoplasia and tumours compared to non-cancer controls. Just three amino acid changes within the early gene, E6, of HPV16 AA have been linked to this augmented tumourigenicity. The AAE6 variant's greater immortalizing and transforming potential over EPE6 has recently been confirmed in retrovirally-transduced keratinocytes expressing the E6 gene only. However, the tumourigenic role of the full-length viral genome of HPV16 has not yet been addressed with regard to these E6 variants. To investigate this process in the context of these two HPV16 E6 genotypes, an organotypic tissue culture model was used to simulate the HPV infectious life cycle. The AAE6 variant demonstrated an enhanced ability over EPE6 to drive the viral life cycle toward tumourigenesis, as evidenced phenotypically—by a more severe grade of epithelial dysplasia with higher proliferation and deregulated differentiation, and molecularly—by high viral oncogene E6 and E7 expression, but lack of productive viral life cycle markers. In contrast, EPE6 had low E6 and E7 but high E1∧E4 expression, indicative of a productive life cycle. We suggest increased viral integration into the host genome for AAE6 as one possible mechanism for these observed differences from EPE6. Additionally, we found downstream effects on immortalization and host innate immune evasion. This study highlights how minor genomic variations in transforming viruses can have a significant affect on their tumourigenic ability.

Intratypic variants of human papillomavirus type 16 and risk of cervical neoplasia in Taiwan

The associations between variants of human papillomavirus (HPV) 16 and risk of cervical neoplasia have been reported, but nucleotide variations of HPV 16 in Asian populations and their association with cervical neoplasia have not been evaluated extensively. During 1991-1992, 11,923 women from seven townships in Taiwan were enrolled. The HPV DNA in cervical cells was detected and genotyped using EasyChip HPV blot. Nucleotide variations in the long control region (LCR), E6, and E7 genes were determined using DNA sequencing for 170 HPV 16-positive cervical samples. The Asian variant was the most prevalent variant (81.8%) of HPV 16 in Taiwan, and was also associated with increased prevalence of histologically confirmed cervical intraepithelial neoplasia grade 3 or worse, showing an age-adjusted odds ratio (exact confidence limits) of 10.70 (1.62-451.05; P = 0.0049) compared to the HPV 16 European variant. Similar significant associations with cervical intraepithelial neoplasia grade 3 or worse were also observed for distinct nucleotide substitutions, including T178A/G, A647G, A7730C/G, T7781C, G7842A, and C24T/G. These results demonstrate that non-European variants (non-E) of HPV 16, predominantly Asian variants, are associated with increased risk for severe cervical neoplasia, compared with European variants. Molecular mechanisms accounting for varied cervical neoplasia risk among different HPV 16 variants warrant further investigation.

Analysis of mutations in the E6 oncogene of human papillomavirus 16 in cervical cancer isolates from Moroccan women

BACKGROUND

Worldwide, cervical cancer is the second most common cancer in women. High-risk human papillomavirus (HPV) play a crucial role in the etiology of cervical cancer and the most prevalent genotype is HPV16. HPV 16 intratypic variants have been reported to differ in their prevalence, biological and biochemical properties. The present study was designed to analyze and identify HPV type 16 E6 variants among patients with cervical cancer in Morocco.

METHODS

A total of 103 HPV16 positive samples were isolated from 129 cervical cancer cases, and variant status was subsequently determined by DNA sequencing of the E6 gene.

RESULTS

Isolates from patients were grouped into the European (E), African (Af) and North-American (NA1) phylogenetic clusters with a high prevalence of E lineage (58.3%). The Af and NA1 variants were detected in 31.1% and 11.6% of the HPV16 positive specimens, respectively, whereas, only 3% of cases were prototype E350T. No European-Asian (EA), Asian (As) or Asian-American (AA) variants were observed in our HPV16-positive specimens. At the amino acid level, the most prevalent non-synonymous variants were L83V (T350G), H78Y (C335T), E113D (A442C), Q14D (C143G/G145T) and R10I (G132T), and were observed respectively in 65%, 41.8%, 38.8%, 30.1% and 23.3% of total samples.Moreover, HPV16 European variants were mostly identified in younger women at early clinical diagnosis stages. Whereas, HPV16 Af variants were most likely associated with cervical cancer development in older women with pronounced aggressiveness.

CONCLUSION

This study suggests a predominance of E lineage strains among Moroccan HPV 16 isolates and raises the possibility that HPV16 variants have a preferential role in progression to malignancy and could be associated with the more aggressive nature of cervical cancer.

Functional similarity between E6 proteins of cutaneous human papillomaviruses and the adenovirus E1A tumor-restraining module

The adenovirus E1A C-terminal region restrains oncogenic transformation through interaction with three distinct cellular protein complexes that include the DYRK1A/1B/HAN11 complex. The E6 proteins of beta-human papillomaviruses (beta-HPVs) also interact with the DYRK1/HAN11 complex. A variant of HPV5 E6 frequently found in epidermodysplasia verruciformis skin lesions interacted less efficiently with DYRK1A/HAN11. The E6 variant and E7 of HPV5 efficiently coimmortalized primary epithelial cells, suggesting that naturally arising variants may contribute potential oncogenic activities of beta-HPV E6 proteins.

HPV-16 E6 and E7 protein T cell epitopes prediction analysis based on distributions of HLA-A loci across populations: an in silico approach

Human papillomavirus type 16 (HPV-16) is the most prevalent virus in human cervical cancers, as it is present in more than half of all cases. Many studies have found continued expression of E6 and E7 proteins in the majority of cervical cancer cases, but not in normal tissues. These results indicated that the E6 and E7 proteins could be ideal candidate therapeutic vaccines against HPV-16 infection and cervical cancer. Using the Immune Epitope Database Analysis Resource, cytotoxic T lymphocyte (CTL) epitopes of the HPV-16 E6 and E7 proteins were predicted according to worldwide frequency distributions of HLA-A alleles (HLA-A*01:01, -A*02:01, -A*02:06, -A*03:01, -A*11:01, -A*24:02, -A*26:01, -A*31:01 and -A*33:03). Our results predicted a total of 81 epitopes of HPV-16 E6 (n=59) and E7 (n=22). Epitope cluster analysis showed that among the 20 clusters of HPV-16 E6, cluster 3 contained the most epitopes (10 epitopes), which was represented by HLA-A*31:01 and -A*33:03. Of the 10 clusters of HPV-16 E7, cluster 3 contained the most epitopes (5 epitopes), which was represented by HLA-A*01:01 and -A*26:01. Our results indicated that the combination of epitopes FAFRDLCIVYR₅₂₋₆₂ of E6 (HLA-A*02:06, HLA-A*31:01, and HLA-A*33:03), PYAVCDKCLKF₆₆₋₇₆ of E6 (HLA-A*11:01 and HLA-A*24:02), HGDTPTLHEY₂₋₁₁ of E7 (HLA-A*01:01 and HLA-A*26:01), and YMLDLQPETT₁₁₋₂₀ of E7 (HLA-A*02:01) could vaccinate >50% of all individuals worldwide. Our results propose CTL epitopes or combinations of them predicted in current study for candidate therapeutic vaccines to effectively control HPV-16 infection and development of cervical cancer.

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.

Association of human papillomavirus type 16 long control region mutation and cervical cancer

BACKGROUND

The variation of human papillomavirus (HPV) genes or HPV variants demonstrates different risks of cervical cancer. Mutation in the long control region (LCR) at YY1-motifs is one of the mechanisms for enhancing viral oncogene expression during the course of cancer cell progression. In Thai women, cervical cancers are almost always associated with HPV16 variant sub-lineage Asian (HPV16As); however, the mechanism involved remains elusive. The aim of this study was to understand further the oncogenic potential of HPV16As.

METHODS

A total of 82 HPV16-positive specimens from Thai women were selected from formalin-fixed paraffin-embedded cervical tissues, and the full length E6 gene of each specimen was amplified and sequenced. LCRs of the HPV16As-positive cases were amplified and sequenced to analyze their polymorphisms. Transcriptional activities of the HPV16As LCRs were then compared with sub-lineage European (EUR), sub-lineage Asian-American 1 (AA1) and HPV16 prototype by insertion of the LCRs into the pGL3-Basic vector.

RESULTS

The HPV16 DNA sequences were classified as HPV16 prototype (18.3%), Asian (As, 61%), Asian American-1 (AA1, 8.5%), European (EUR, 7.3%), Asian African-2 (AFR2, 3.7%) and Java-135C (J135C, 1.2%). The prevalence of HPV16As was 30% in low-grade squamous intraepithelial lesion (LSIL), while that in high-grade squamous intraepithelial lesion (HSIL) and squamous cell cervical carcinoma (SCC) were 63.9% and 66.7%, respectively, which demonstrates a significant association of HPV16As with the disease severity. LCR polymorphisms from 43 HPV16As positive cases were analyzed by PCR-sequencing. Thirty-eight nucleotide variation positions spanned nucleotide positions 7157-82. Ten new mutations found in the HPV16As LCRs were located predominantly at the enhancer and proximal to the 3'-end of the early promoter. The LCRs of the common HPV16As, EUR and AA1 showed 5, 13 and 23-fold higher activity than the HPV16 prototype LCR, while those of the new nucleotide variations of As showed 19 (As-sv1) and 30 (As-sv14) -fold higher activity than the HPV16 prototype.

CONCLUSIONS

HPV16As DNA sequence variation, especially at the proximal to early promoter in the LCR, enhances transcriptional activity. This could be one of the possible mechanisms for HPV16As-associated cervical cancer development.

Codon 72 polymorphism of p53 and HPV type 16 E6 variants as risk factors for patients with squamous epithelial lesion of the uterine cervix

The Arg/Arg genotype versus Arg/Pro or Pro/Pro at codon 72 of the p53 gene in association with human papillomavirus (HPV) 16 E6 variants has been implicated as a risk marker in cervical neoplasia. However, research on this topic has produced controversial results. The association of p53 codon 72 polymorphism alone and in combination with specific HPV 16 E6 variants with risk of developing squamous intraepithelial cervical lesion has been investigated in low and high-grade squamous intraepithelial lesions and in HPV-negative controls from an Italian population. The data obtained showed statistically significant different distribution of p53 genotypes between healthy controls and precursor lesions, with the p53 arginine homozygous increased in high-grade squamous intraepithelial lesions. The T350G HPV 16 variant was the most frequent variant observed in the analyzed group of Italian women, showing a slight decreasing with the severity of the lesion. At the same time, the number of the prototype T350 slightly increased with the severity of the cytological lesions. In conclusion, p53 arginine homozygous was found to be increased in high-grade lesions, supporting the results of previous investigations indicating that HPV-positive patients with p53 Arg/Arg have an increased risk of developing pre-cancerous lesions. In addition, T350G HPV 16 variant was over-represented in p53 Arg homozygous women with cervical lesions. When p53 genotype and HPV 16 variants are considered together, no difference emerges between cases and controls so is not possible to assess that the oncogenic effect of HPV 16 T350G variant may be influenced by the p53 genotype.

Human papillomavirus infection and the risk of cervical cancer in Japan

Infection with a high-risk human papillomavirus (HPV) is established as a major risk factor for cervical cancer. In Japan, eight HPV genotypes (16, 18, 31, 33, 35, 45, 52 and 58) confer a much higher risk of cervical cancer and its immediate premalignant lesions (cervical intraepithelial neoplasia grade 3 [CIN 3]) than do other high-risk and the low-risk HPV types. However, only a small fraction of infected women develop CIN 3 or invasive cervical cancer, suggesting the involvement of additional cofactors in cervical carcinogenesis. A multicenter prospective cohort study of Japanese women with low-grade cervical abnormalities has demonstrated that tobacco smoking is significantly associated with an increased risk of persistent low-grade cervical lesions, while the presence of human leukocyte antigen Class II DRB1*1302 protects against progression to CIN 3. Understanding the epidemiology of HPV infections and cervical cancer may improve prevention and lead to novel clinical management strategies for cervical precursor lesions. Viral and host factors that affect the risk of cervical cancer are discussed in this review.

Characterization and whole genome analysis of human papillomavirus type 16 e1-1374^63nt variants

BACKGROUND

The variation of the most common Human papillomavirus (HPV) type found in cervical cancer, the HPV16, has been extensively investigated in almost all viral genes. The E1 gene variation, however, has been rarely studied. The main objective of the present investigation was to analyze the variability of the E6 and E1 genes, focusing on the recently identified E1-1374^63nt variant.

METHODOLOGY/PRINCIPAL FINDINGS

Variation within the E6 of 786 HPV16 positive cervical samples was analyzed using high-resolution melting, while the E1-1374^63nt duplication was assayed by PCR. Both techniques were supplemented with sequencing. The E1-1374^63nt duplication was linked with the E-G350 and the E-C109/G350 variants. In comparison to the referent HPV16, the E1-1374^63nt E-G350 variant was significantly associated with lower grade cervical lesions (p = 0.029), while the E1-1374^63nt E-C109/G350 variant was equally distributed between high and low grade lesions. The E1-1374^63nt variants were phylogenetically closest to E-G350 variant lineage (A2 sub-lineage based on full genome classification). The major differences between E1-1374^63nt variants were within the LCR and the E6 region. On the other hand, changes within the E1 region were the major differences from the A2 sub-lineage, which has been historically but inconclusively associated with high grade cervical disease. Thus, the shared variations cannot explain the particular association of the E1-1374^63nt variant with lower grade cervical lesions.

CONCLUSIONS/SIGNIFICANCE

The E1 region has been thus far considered to be well conserved among all HPVs and therefore uninteresting for variability studies. However, this study shows that the variations within the E1 region could possibly affect cervical disease, since the E1-1374^63nt E-G350 variant is significantly associated with lower grade cervical lesions, in comparison to the A1 and A2 sub-lineage variants. Furthermore, it appears that the silent variation 109T>C of the E-C109/G350 variant might have a significant role in the viral life cycle and warrants further study.

Human papillomavirus type 16 genetic variants: phylogeny and classification based on E6 and LCR

Naturally occurring genetic variants of human papillomavirus type 16 (HPV16) are common and have previously been classified into 4 major lineages; European-Asian (EAS), including the sublineages European (EUR) and Asian (As), African 1 (AFR1), African 2 (AFR2), and North-American/Asian-American (NA/AA). We aimed to improve the classification of HPV16 variant lineages by using a large resource of HPV16-positive cervical samples collected from geographically diverse populations in studies on HPV and/or cervical cancer undertaken by the International Agency for Research on Cancer. In total, we sequenced the entire E6 genes and long control regions (LCRs) of 953 HPV16 isolates from 27 different countries worldwide. Phylogenetic analyses confirmed previously described variant lineages and subclassifications. We characterized two new sublineages within each of the lineages AFR1 and AFR2 that are robustly classified using E6 and/or the LCR. We could differentiate previously identified AA1, AA2, and NA sublineages, although they could not be distinguished by E6 alone, requiring the LCR for correct phylogenetic classification. We thus provide a classification system for HPV16 genomes based on 13 and 32 phylogenetically distinguishing positions in E6 and the LCR, respectively, that distinguish nine HPV16 variant sublineages (EUR, As, AFR1a, AFR1b, AFR2a, AFR2b, NA, AA1, and AA2). Ninety-seven percent of all 953 samples fitted this classification perfectly. Other positions were frequently polymorphic within one or more lineages but did not define phylogenetic subgroups. Such a standardized classification of HPV16 variants is important for future epidemiological and biological studies of the carcinogenic potential of HPV16 variant lineages.

Gene expression profiles are altered in human papillomavirus-16 E6 D25E-expressing cell lines

Previously, we have reported that the human papillomavirus (HPV) type 16 E6 D25E is the most prevalent variant in Korean women at high risk for cervical cancers. Several studies have identified an association between the increased frequency of this variant and the elevated risk of cervical intraepithelial neoplasia and invasive cervical carcinoma. To investigate whether the HPV-16 E6 D25E variant might influence cervical cancer progression, we used an oligonucleotide microarray approach to identify transcriptionally altered gene expression patterns in recombinant wild-type E6 or E6 D25E variant-expressing HPV-negative cancer cells. We found that 211 genes were significantly up- or down-regulated (at least 1.5-fold, p < 0.05). We identified 14 genes, nine down-regulated and five up-regulated upon E6 D25E expression, compared with wild-type E6 expression. These results further emphasize the unique biological activity of the HPV-16 E6 D25E variant.

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.

Rare human papillomavirus 16 E6 variants reveal significant oncogenic potential

The aim of this study was to determine whether low prevalence human papillomavirus (HPV) 16 E6 variants differ from high prevalence types in their functional abilities. We evaluated functions relevant to carcinogenesis for the rarely-detected European variants R8Q, R10G and R48W as compared to the commonly detected L83V. Human immortalized keratinocytes (NIKS) stably transduced with the E6 variants were used in most functional assays. Low and high prevalence E6 variants displayed similar abilities in abrogation of growth arrest and inhibition of p53 elevation induced by actinomycin D. Differences were detected in the abilities to dysregulate stratification and differentiation of NIKS in organotypic raft cultures, modulate detachment induced apoptosis (anoikis) and hyperactivate Wnt signaling. No distinctive phenotype could be assigned to include all rare variants. Like L83V, raft cultures derived from variants R10G and R48W similarly induced hyperplasia and aberrantly expressed keratin 5 in the suprabasal compartment with significantly lower expression of keratin 10. Unlike L83V, both variants, and particularly R48W, induced increased levels of anoikis upon suspension in semisolid medium. R8Q induced a unique phenotype characterized by thin organotypic raft cultures, low expression of keratin 10, and high expression of keratins 5 and 14 throughout all raft layers. Interestingly, in a reporter based assay R8Q exhibited a higher ability to augment TCF/β-catenin transcription. The data suggests that differences in E6 variant prevalence in cervical carcinoma may not be related to the carcinogenic potential of the E6 protein.