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

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

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

Impact of HPV-16 Lineages Infection in Response to Radio-Chemotherapy in Cervical Cancer

BACKGROUND

HPV is strongly related to cervical cancer. HPV lineages can contribute to a response to cervical cancer therapy. The aim of this research was to estimate the frequency of human papillomavirus (HPV)-16 lineages in specimens of cervical cancer, relate the pathological factors in these variants, and assess their response to treatment with radical chemoradiotherapy.

METHODS

Samples of cervical cancer were collected from women who were referred to a reference cancer hospital to test the presence of human papillomavirus-type DNA. The standard protocol of this service consisted of cisplatin-based chemotherapy of 40 mg/m², plus conventional pelvic irradiation in doses of 45-50.4 Gy and high dose-rate brachytherapy of 28-30 Gy to Point A. The response to chemotherapy was evaluated after three months in patients with the HPV-16 lineage.

RESULTS

HPV DNA was detected in 104 (88.1%) of the 118 patients. HPV-16 was present in 63 patients (53%). Lineages of HPV-16 were identified in 57 patients and comprised 33 instances of (57.8%) lineage A, 2 instances of lineage B (3.5%), 2 instances of lineage C (3.5%), and 20 instances of (35.0%) lineage D. The median age of the patients was 48.4 years (range 25-85 years). Squamous cell carcinoma was detected 48 times (84.2%). Adenocarcinoma was more likely to occur in lineage D, as three of the four cases occurred in this lineage. A total of 11 patients with the HPV-16 variant were treated with chemoradiotherapy. After three months, it was observed that nine of the eleven patients (81.8%) achieved a complete response, five with the lineage A type, two with the lineage C type, and two with the lineage D type. The two cases of partial response and disease progression, one of each, occurred in lineage A.

CONCLUSIONS

In addition to the small number of patients and HPV variants, we noticed a better response in patients with the HPV-16 lineage A. Increasing the sample size could be helpful to better assess the impact of HPV variants on cervical cancer treatment.

Transcriptional activity of the long control region in human papillomavirus type 33 intratype variants

BACKGROUND

High-risk human papillomaviruses (HPVs) are responsible for the development of cervical and other anogenital cancers. Intratype sequence variants of certain high-risk HPV types (e.g. 16, 18 and 31) are thought to have different oncogenic potential, partly due to nucleotide sequence variation in the viral long control region (LCR). The LCR has an important role in the regulation of viral replication and transcription. The purpose of this study was to explore sequence variation in the LCR of HPV 33 intratype variants in Hungary and to see whether there are differences in the transcriptional activities of the variants.

METHODS

The complete HPV 33 LCR was amplified from HPV 33 positive cervical samples. After sequencing the LCR variants, multiple sequence alignment and phylogenetic analyses were carried out. Representative HPV 33 LCR sequence variants were selected for cloning and functional analysis. After transient transfection of HeLa cells, luciferase reporter assays were used to analyse the transcriptional activities of different LCR variants.

RESULTS

Altogether 10 different variants were identified by sequence analysis of the HPV 33 LCR. The results of phylogenetic analysis showed that 3 variants belonged to sublineage A1, while the other 7 variants clustered with sublineage A2. Variants belonging to sublineage A2 had significantly lower transcriptional activities than variants belonging to sublineage A1. Within sublineage A2, the two variants analysed had significantly different transcriptional activities, which was shown to be caused by the A7879G variation.

CONCLUSIONS

Nucleotide variation in the HPV 33 LCR can result in altered transcriptional activity of the intratype variants. Our results can help to understand the correlation between LCR polymorphism and the oncogenic potential of HPV 33 variants.

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.

Characterization of Human Papillomavirus 16 from Kinshasa (Democratic Republic of the Congo)-Implications for Pathogenicity and Vaccine Effectiveness

Human Papillomavirus (HPV) type 16 is the main etiological agent of cervical cancer worldwide. Mutations within the virus genome may lead to an increased risk of cancer development and decreased vaccine response, but there is a lack of information about strains circulating in Sub-Saharan Africa. Endocervical cytology samples were collected from 480 women attending a voluntary cervical cancer screening program at Monkole Hospital and four outpatient centers in Kinshasa, Democratic Republic of the Congo (DRC). The prevalence of HPV infection was 18.8% and the most prevalent high-risk types were HPV16 (12.2%) followed by HPV52 (8.8%) and HPV33/HPV35 (7.8% each). HPV16 strains were characterized: 57.1% were classified as C lineage; two samples (28.6%) as A1 and one sample belonged to B1 lineage. HPV33, HPV35, HPV16, and HPV58 were the most frequent types associated with low-grade intraepithelial lesion while high-grade squamous intraepithelial lesions were predominantly associated with HPV16. Several L1 mutations (T266A, S282P, T353P, and N181T) were common in Kinshasa, and their potential effect on vaccine-induced neutralization, especially the presence of S282P, should be further investigated. Long control region (LCR) variability was high with frequent mutations like G7193T, G7521A, and G145T that could promote malignancy of these HPV16 strains. This study provides a helpful basis for understanding HPV16 variants circulating in Kinshasa and the potential association between mutations of LCR region and malignancy and of L1 and vaccine activity.

Human papillomavirus integration perspective in small cell cervical carcinoma

Small cell cervical carcinoma (SCCC) is a rare but aggressive malignancy. Here, we report human papillomavirus features and genomic landscape in SCCC via high-throughput HPV captured sequencing, whole-genome sequencing, whole-transcriptome sequencing, and OncoScan microarrays. HPV18 infections and integrations are commonly detected. Besides MYC family genes (37.9%), we identify SOX (8.4%), NR4A (6.3%), ANKRD (7.4%), and CEA (3.2%) family genes as HPV-integrated hotspots. We construct the genomic local haplotype around HPV-integrated sites, and find tandem duplications and amplified HPV long control regions (LCR). We propose three prominent HPV integration patterns: duplicating oncogenes (MYCN, MYC, and NR4A2), forming fusions (FGFR3-TACC3 and ANKRD12-NDUFV2), and activating genes (MYC) via the cis-regulations of viral LCRs. Moreover, focal CNA amplification peaks harbor canonical cancer genes including the HPV-integrated hotspots within MYC family, SOX2, and others. Our findings may provide potential molecular criteria for the accurate diagnosis and efficacious therapies for this lethal disease.

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.

Epidemiology and Molecular Biology of HPV Variants in Cervical Cancer: The State of the Art in Mexico

Cervical cancer (CC) continues to be a major public health problem in Mexico, ranking second among cancers in women. A persistent infection with human papillomaviruses (HPV) is the main risk factor for CC development. In addition, a significant fraction of other cancers including those of the anus, oropharynx, and penis are also related to HPV infection. In CC, HPV-16 is the most prevalent high-risk HPV type, followed by HPV-18, both being responsible for 70% of cases. HPV intratype variant lineages differ in nucleotide sequences by 1–10%, while sublineages differ by 0.5–1%. Several studies have postulated that the nucleotide changes that occur between HPV intratype variants are reflected in functional differences and in pathogenicity. Moreover, it has been demonstrated that HPV-16 and -18 intratype variants differentially affect molecular processes in infected cells, changing their biological behavior that finally impacts in the clinical outcome of patients. Mexico has participated in providing knowledge on the geographical distribution of intratype variants of the most prevalent HPVs in premalignant lesions of the cervix and cervical cancer, as well as in other HPV-related tumors. In addition, functional studies have been carried out to assess the cellular effects of intratype variations in HPV proteins. This review addresses the state of the art on the epidemiology of HPV-16 and HPV-18 intratype variants in the Mexican population, as well as their association with persistence, precancer and cervical cancer, and functional aspects related to their biological behavior.

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.

Genetic Diversity of HPV 16 and HPV 18 Based on Partial Long Control Region in Iranian Women

Background

Human papillomavirus (HPV) 16 and HPV 18 account for 75% of all cervical cancers. The L1 gene, encoding the major surface protein (MSP), is used to classify HPV types (lineages and sublineages), genotypes, and intratypic variants. Therefore, this study aimed to investigate the lineages, sublineages, genetic variabilities, and mutation effects on transcription factor binding sites by using partial sequences of the HPV 16 and HPV 18 long control regions (LCRs) in these samples.

Materials and Methods

After DNA isolation from 56 positive samples, the LCR of HPV 16 and HPV 18 were amplified using specific primers, and phylogenetic trees were drawn through MEGA X. Compared to the reference sequences, single nucleotide polymorphisms (SNPs) were identified. The transcription binding sites were also evaluated using the online PROMO database.

Results

The LCRs of 52 samples were successfully sequenced. Overall, 81.58% of all HPV 16 variants belonged to the D1 sublineage, followed by A4 (13.16%), A1 (2.63%), and C1 (2.63%) sublineages. All HPV 18 isolates belonged to A sublineage, 92.85% to A3 sublineage, and 7.15% to A4 sublineage. Out of 27 SNPs in the HPV 16 LCR, A7382T, T7384G, C7387T, C7393G, A7431G, T7448C, and C7783A were HPV 16-specific. Also, among 14 SNPs in the HPV 18 LCR, C7577A and A7943T were not previously reported. An insertion (C) between 7432 and 7433 positions was identified in all studied HPV 16 variants. Besides, most of the HPV 16 mutations were embedded in the YY1, TFIID, Oct-2, and NF-1 binding sites, while c-Fos and MBF1, as the most common binding sites, were affected by HPV 18 LCR mutations.

Conclusion

The present results showed that D1 and A3 were the dominant sublineages of HPV 16 and HPV 18, respectively. Therefore, women infected with these variants need to be examined in further longitudinal studies to obtain more information about the oncogenic potential of these dominant variants in Iran. Besides, YY1, TFIID, Oct-2, NF-1, c-Fos, and MBF1 were the most frequent binding sites, which were influenced by the mutations.

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

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

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

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

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.

The frequency of high-risk human papillomavirus types, HPV16 lineages, and their relationship with p16INK4a and NF-κB expression in head and neck squamous cell carcinomas in Southwestern Iran

High-risk human papillomaviruses (hr-HPVs) are the key risk factors implicated in the development of a significant proportion of head and neck squamous cell carcinomas (HNSCCs). We aimed to investigate the distribution of hr-HPV types and HPV16 lineages in a sample of patients with HNSCC and the possible association between HPV status and the expression of P16INK4A and NF-κB in Iranian HNSCC patients. We examined 108 formalin-fixed, paraffin-embedded (FFPE) histologically confirmed primary SCC tissue specimens of different head and neck anatomical sites. HPV types and HPV16 lineages were determined by nested PCR and overlapping nested PCR assays, respectively, followed by gene sequencing and phylogenetic analysis. The expression of p16INK4a and NF-κB was evaluated by immunohistochemistry. Twenty-five (23.1%) HNSCC tissue specimens were tested positive for HPV infection. The most prevalent HPV type was HPV-16, followed by HPV18 and HPV11. HPV16 variants belonged to the lineage A and lineage D which were further sorted into sublineages A1, A2, and D2. A significant association between HPV status and p16INK4a immunoreactivity was observed in more than 76% of the HPV-related HNSCCs (P < 0.0001). The overexpression of p16INK4a and cytoplasmic NF-κB was more common in low-grade HNSCC tumors. Our data highlights that HPV16, in particular the A2 sublineage, followed by A1 and D2 sublineages are the major agents associated with HNSCCs in Iran. Based on HPV16 predominance and its lineage distribution pattern, it seems that the prophylactic vaccines developed for cervical cancer prevention could also be applicable for the prevention of HPV-related HNSCCs in our population.

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.

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

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

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

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

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.

Redetection of human papillomavirus type 16 infections of the cervix in mid-adult life

PURPOSE

To assess whether HPV 16 originally detected in adolescent women can be redetected in adulthood.

METHODS

A convenience sample of 27 adult women with known HPV 16 detection during adolescence was assessed for HPV 16 redetection. A comparison of the long control region (LCR) DNA sequences was performed on some of the original and redetected HPV 16 isolates.

RESULTS

Median age at reenrollment was 27.5 years (interquartile range of 26.7-29.6). Reenrollment occurred six years on average after the original HPV 16 detection. Eleven of 27 women had HPV 16 redetected. Some of these HPV 16 infections had apparently cleared during adolescence. LCR sequencing was successful in paired isolates from 6 women; in 5 of 6 cases the redetected HPV 16 isolates were identical to those detected during adolescence, CONCLUSIONS: HPV 16 may be episodically detected in young women, even over long time periods. HPV 16 redetection with identical LCR sequences suggests low-level persistent infection rather than true clearance, although newly acquired infection with an identical HPV 16 isolate cannot be excluded. However, this study suggests that a new HPV 16-positive test in a clinical setting may not indicate a new infection.

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

Backgroud

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

Methods

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

Results

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

Conclusion

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

Electronic supplementary material

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

Expression of Yin Yang 1 in cervical cancer and its correlation with E-cadherin expression and HPV16 E6

The molecular mechanisms of normal cervical squamous epithelium advancing to cervical intraepithelial neoplasia (CIN) and eventually to cervical squamous cell carcinoma (CSCC) are largely unknown. This study explored abnormal expression of Yin Yang 1 (YY1) in cervical cancer and its correlation with the expression of E-cadherin and human papillomavirus (HPV) 16 E6. YY1, E-cadherin and HPV16 E6 expression were detected by immunohistochemistry in 90 cervical tissue specimens collected from 30 patients with hysteromyoma, 15 patients with CIN I, 15 patients with CIN II-III, and 30 patients with CSCC. The H-score method was employed to measure the expression of YY1, E-cadherin and HPV16 E6. Increased expression of YY1 and HPV16 E6, and the decreased expression levels of E-cadherin were strongly associated with malignant transformation of the cervical epithelium and the histological progression of CSCC. The expression of YY1 in cervical tissues was inversely correlated with E-cadherin expression, and positively correlated with HPV16 E6 expression. Expression of YY1 in CSCC tissues was not significantly correlated with tumor differentiation, but was significantly correlated with an advanced clinical stage of CSCC. These results suggest that up-regulation of YY1 is closely associated with the progression of CSCC, and YY1 may play an important role in the pathogenesis of cervical cancer by modulating the expression of E-cadherin and HPV16 E6.

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.

Novel insights into epigenetic drivers of oropharyngeal squamous cell carcinoma: role of HPV and lifestyle factors

In the last years, the explosion of high throughput sequencing technologies has enabled epigenome-wide analyses, allowing a more comprehensive overview of the oropharyngeal squamous cell carcinoma (OPSCC) epigenetic landscape. In this setting, the cellular pathways contributing to the neoplastic phenotype, including cell cycle regulation, cell signaling, DNA repair, and apoptosis have been demonstrated to be potential targets of epigenetic alterations in OPSCC. Of note, it has becoming increasingly clear that HPV infection and OPSCC lifestyle risk factors differently drive the epigenetic machinery in cancer cells. Epigenetic changes, including DNA methylation, histone modifications, and non-coding RNA expression, can be used as powerful and reliable tools for early diagnosis of OPSCC patients and improve prognostication. Since epigenetic changes are dynamic and reversible, epigenetic enzymes may also represent suitable targets for the development of more effective OPSCC therapeutic strategies. Thus, this review will focus on the main known epigenetic modifications that can occur in OPSCC and their exploitation as potential biomarkers and therapeutic targets. Furthermore, we will address epigenetic alterations to OPSCC risk factors, with a particular focus on HPV infection, tobacco exposure, and heavy alcohol consumption.

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

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

Study of association and molecular analysis of human papillomavirus in breast cancer of Indian patients: Clinical and prognostic implication

OBJECTIVES

Human papillomavirus (HPV) causes tumors primarily Cervical cancer. Recently, inconsistent reports came up in Breast cancer (BC) too. In India, despite treatment 70,218 BC patients die each year. So, we explored the association of HPV, if any, with BC prognosis in Indian pre-therapeutic (PT) and Neo-adjuvant chemotherapy (NACT) patients with subsequent analysis of HPV profile.

METHODS

HPV prevalence was checked and analysis of physical status, copy number, genome variation, promoter methylation and expression (mRNA and protein) of the prevalent subtype was done.

RESULTS

High prevalence of HPV was observed in both PT (64.0%) and NACT (71.0%) cases with significant association with younger (20-45 yrs) PT patients. Interestingly, HPV infection was significantly increased from adjacent normal breast (9.5%, 2/21), fibro adenomas (30%, 3/10) to tumors (64.8%, 203/313) samples. In both PT and NACT cases, HPV16 was the most prevalent subtype (69.0%) followed by HPV18 and HPV33. Survival analysis illustrated hrHPV infected PT patients had worst prognosis. So, detailed analysis of HPV16 profile was done which showed Europian-G350 as the most frequent HPV16 variant along with high rate of integration. Moreover, low copy number and hyper-methylation of P97 early promoter were concordant with low HPV16 E6 and E7 mRNA and protein expression. Notably, four novel variations (KT020838, KT020840, KT020841 and KT020839) in the LCR region and two (KT020836 and KT020837) in the E6 region were identified for the first time along with two novel E6^E7*I (KU199314) and E6^E7*II (KU199315) fusion transcript variants.

CONCLUSION

Thus, significant association of hrHPV with prognosis of Indian BC patients led to additional investigation of HPV16 profile. Outcomes indicated a plausible role of HPV in Indian BC patients.

Different Association of Human Papillomavirus 16 Variants with Early and Late Presentation of Cervical Cancer

The median age of cervical cancer (CC) presentation coincides with the mean age of menopause presentation (49 years) in Mexico. Here, we investigated the association between different HPV16 variants and early (≤ 49 years) or delayed (≥ 50 years) CC presentation. We conducted a case-case study that included 462 CCs, 386 squamous cell carcinomas (SCC), 63 adenocarcinomas (ACC), and 13 additional cell types. Variants were identified by PCR and DNA sequencing. The risk conferred by each variant for developing CC earlier than 50 years was analyzed using a univariate logistic regression model considering old-aged patients (≥ 50 years) and non-HPV16 cases as the reference variables. Overall, the frequency of HPV16 was 50.9%, and the only identified variants were the European A1/2 (31.2%) and the Asian-American D2 (10.8%), and D3 (8.9%). D2 was mainly associated with ≤ 49-year-old patients (15.9%); A1/2 was uniformly distributed between the two age groups (~31%), whereas D3 increased with age to a frequency of 11.8% in the older group. Only the D2 variant conferred a 3.3-fold increase in the risk of developing CC before 50 years of age (OR = 3.3, 95% CI = 1.7-6.6, p < 0.001) in relation with non-HPV16 cases. Remarkably, this risk was higher for ACC (OR = 6.0, 95% CI = 1.1-33, p < 0.05) than for SCC (OR = 2.8, 95% CI = 1.3-5.9, p < 0.01). Interestingly, when analyzing only the HPV16-positive CC, D2 increases (OR = 2.5, 95% CI = 1.2-5, p < 0.05) and D3 decreases (OR = 0.45, 95% CI 0.2-0.9, p < 0.05) the risk to develop CC before 50 years old in relation with A1/2 variant. These results indicated that D2 variant is associated with early and D3 with delayed CC presentation, whereas A1/2 variant was uniformly distributed between the two age groups.

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

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.

Phylogenetic analysis and predicted functional effect of protein mutations of E6 and E7 HPV16 strains isolated in Indonesia

Background: E6 and E7 are oncoproteins of HPV16. Natural amino acid variation in HPV16 E6 can alter its carcinogenic potential. The aim of this study was to analyze phylogenetically E6 and E7 genes and proteins of HPV16 from Indonesia and predict the effects of single amino acid substitution on protein function. This analysis could be used to reduce time, effort, and research cost as initial screening in selection of protein or isolates to be tested in vitro or in vivo.Methods: In this study, E6 and E7 gene sequences were obtained from 12 samples of  Indonesian isolates, which  were compared with HPV16R (prototype) and 6 standard isolates in the category of European (E), Asian (As), Asian-American (AA), African-1 (Af-1), African-2 (Af-2), and North American (NA) branch from Genbank. Bioedit v.7.0.0 was used to analyze the composition and substitution of single amino acids. Phylogenetic analysis of E6 and E7 genes and proteins was performed using Clustal X (1.81) and NJPLOT softwares. Effects of single amino acid substitutions on protein function of E6 and E7 were analysed by SNAP.Results: Java variants and isolate ui66* belonged to European branch, while the others belonged to Asian and African branches. Twelve changes of amino acids were found in E6 and one in E7 proteins. SNAP analysis showed two non neutral mutations, i.e. R10I and C63G in E6 proteins. R10I mutations were found in Af-2 genotype (AF472509) and Indonesian isolates (Af2*), while C63G mutation was found only in Af2*.Conclusion: E6 proteins of HPV16 variants were more variable than E7. SNAP analysis showed that only E6 protein of African-2 branch had functional differences compared to HPV16R.

Deep sequencing of HPV16 genomes: A new high-throughput tool for exploring the carcinogenicity and natural history of HPV16 infection

For unknown reasons, there is huge variability in risk conferred by different HPV types and, remarkably, strong differences even between closely related variant lineages within each type. HPV16 is a uniquely powerful carcinogenic type, causing approximately half of cervical cancer and most other HPV-related cancers. To permit the large-scale study of HPV genome variability and precancer/cancer, starting with HPV16 and cervical cancer, we developed a high-throughput next-generation sequencing (NGS) whole-genome method. We designed a custom HPV16 AmpliSeq™ panel that generated 47 overlapping amplicons covering 99% of the genome sequenced on the Ion Torrent Proton platform. After validating with Sanger, the current "gold standard" of sequencing, in 89 specimens with concordance of 99.9%, we used our NGS method and custom annotation pipeline to sequence 796 HPV16-positive exfoliated cervical cell specimens. The median completion rate per sample was 98.0%. Our method enabled us to discover novel SNPs, large contiguous deletions suggestive of viral integration (OR of 27.3, 95% CI 3.3-222, P=0.002), and the sensitive detection of variant lineage coinfections. This method represents an innovative high-throughput, ultra-deep coverage technique for HPV genomic sequencing, which, in turn, enables the investigation of the role of genetic variation in HPV epidemiology and carcinogenesis.

Polymorphisms and functional analysis of the intact human papillomavirus16 e2 gene

High risk human papillomavirus (HR-HPV) E2 proteins play roles in transcriptional regulation and are commonly functionally disrupted when the HPV genome integrates into host chromosomes. Some 15-40% of cancer cases, however, contain an intact E2 gene or episomal HPV. In these cases, polymorphism of the E2 gene might be involved. This study aimed to determine polymorphisms of the E2 gene in episomal HPV16 detected in high grade squamous intraepithelial lesions and squamous cell carcinomas and altered functions compared to the E2 prototype. The E2 gene was amplified and sequenced. Two expression vectors containing E2 gene polymorphisms were constructed and transfected in SiHa and C33A cells, then E6 gene as well as Il- 10 and TNF-α expression was determined by quantitative RT-PCR. Expression vectors and reporter vectors containing the HPV16 long control region (LCR) were co-transfected and transcriptional activity was determined. The results showed that a total of 32 nucleotides and 23 amino acids were changed in all 20 cases of study, found in the transactivation (TA) domain, hinge (H) region and DNA binding (DB) domain with 14, 5 and 13 nucleotide positions. They mostly caused amino acid change. The expressing vectors containing different E2 gene polymorphisms showed E6 mRNA suppression, TNF-α mRNA suppression and IL-10 induction but no statistically significant differences when compared to the E2 prototype. Moreover, promoter activity in HPV16 LCR was not affected by E2 protein with different gene polymorphisms, in contrast to nucleotide variations in LCR that showed an effect on transcription activity. These results demonstrated that E2 gene polymorphisms of episomal HPV16 did not affect transcriptional regulation and suggested that nucleotide variation as well as epigenetic modification of the LCR might play a role in inducing malignant transformation of cells containing episomal HPV16.

Activity of Andrographolide and Its Derivatives on HPV16 Pseudovirus Infection and Viral Oncogene Expression in Cervical Carcinoma Cells

Andrographolide (Androg) has been reported to contain antiviral and antitumor activities, but the effects of Androg on human papillomavirus (HPV) infection and cervical cancer have not been elucidated. This study investigated the effects of Androg and its derivatives, namely, 14-deoxy-11,12-didehydroandrographolide (14-DDA) and 3,19-isopropylidene andrographolide (IPAD), on HPV16 pseudovirus (HPV16PsV) infectivity, HPV16 E6 oncogene expression and cervical cancer cell apoptosis. The result demonstrated that all compounds inhibited HPV16PsV infection and that 14-DDA showed the highest potency. Only Androg suppressed long control region (LCR) transcription activity of HPV16 in transiently transfected C33A cells and significantly inhibited E6 oncogene expression in SiHa cells in a dose-dependent manner. A twofold subcytotoxic concentration of IPAD exhibited an inhibitory effect on E6 oncogene expression at 48-h posttreatment. Interestingly, p53 protein was restored in a downstream process and was detected earlier by IPAD treatment than by Androg treatment. This result corresponded to the level of cell apoptosis and cell cycle arrest at the G2/M phase. E6 oncogene expression was also suppressed in CaSki cells treated with Androg and IPAD leading to cell apoptosis. These findings imply that Androg and its derivatives have different activities and may be effective agents for HPV prevention and cervical cancer treatment.

Sequence variation in the E2-binding domain of HPV16 and biological function evaluation in Tunisian cervical cancers

HPV16 E2 variants have different effects on the transcriptional activity of the LCR. In this study, we examined the nucleotide and amino acid sequence variation within the HPV16 E2 gene and to correlate with disease progression. E2 gene disruption was detected by PCR amplification of the entire E2 gene using a single set of primers. Nucleotide variations were analyzed by bidirectional sequencing. mRNA expression patterns of E6 and E7 gene transcripts were evaluated by a reverse transcriptase-PCR method (RT-PCR). The detection of intact E2 genes was significantly higher among controls than cases (81.8% versus 37.5%, resp., P < 0.05). Among the E subgroup, variation at position 3684 C>A results in the amino acid substitution T310K and was more common among the E2 undisrupted cases (7/9; 77.7%), compared to controls (2/9; 22.2%). In addition, specific sequence variations identified in the E2 ORF at positions 3684 C>A were associated with increased viral oncogenes E6-E7 production. Besides HPV16 E2 disruption, the 3684 C>A variation within undisrupted E2 genes could be involved in an alternative mechanism for deregulating the expression of the HPV16 E6 and E7 oncogenes and appears to be a major factor contributing to the development of cervical cancer in Tunisian women.

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.

The modulation of apoptosis by oncogenic viruses

Transforming viruses can change a normal cell into a cancer cell during their normal life cycle. Persistent infections with these viruses have been recognized to cause some types of cancer. These viruses have been implicated in the modulation of various biological processes, such as proliferation, differentiation and apoptosis. The study of infections caused by oncogenic viruses had helped in our understanding of several mechanisms that regulate cell growth, as well as the molecular alterations leading to cancer. Therefore, transforming viruses provide models of study that have enabled the advances in cancer research. Viruses with transforming abilities, include different members of the Human Papillomavirus (HPV) family, Hepatitis C virus (HCV), Human T-cell Leukemia virus (HTLV-1), Epstein Barr virus (EBV) and Kaposi’s Sarcoma Herpesvirus (KSHV).

Apoptosis, or programmed cell death, is a tightly regulated process that plays an important role in development and homeostasis. Additionally, it functions as an antiviral defense mechanism. The deregulation of apoptosis has been implicated in the etiology of diverse diseases, including cancer. Oncogenic viruses employ different mechanisms to inhibit the apoptotic process, allowing the propagation of infected and damaged cells. During this process, some viral proteins are able to evade the immune system, while others can directly interact with the caspases involved in apoptotic signaling. In some instances, viral proteins can also promote apoptosis, which may be necessary for an accurate regulation of the initial stages of infection.