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

Identification of HPV16 Lineages in South African and Mozambican Women with Normal and Abnormal Cervical Cytology

BACKGROUND

Human papillomavirus 16 (HPV16) is an oncogenic virus responsible for the majority of invasive cervical cancer cases worldwide. Due to genetic modifications, some variants are more oncogenic than others. We analysed the HPV16 phylogeny in HPV16-positive cervical Desoxyribonucleic Acid (DNA) samples collected from South African and Mozambican women to detect the circulating lineages.

METHODS

Polymerase chain reaction (PCR) amplification of the long control region (LCR) and 300 nucleotides of the E6 region was performed using HPV16-specific primers on HPV16-positive cervical samples collected in women from South Africa and Mozambique. HPV16 sequences were obtained through Next Generation Sequencing (NGS) methods. Geneious prime and MEGA 11 software were used to align the sequences to 16 HPV16 reference sequences, gathering the A, B, C, and D lineages and generating the phylogenetic tree. Single nucleotide polymorphisms (SNPs) in the LCR and E6 regions were analysed and the phylogenetic tree was generated using Geneious Prime software.

RESULTS

Fifty-eight sequences were analysed. Of these sequences, 79% (46/58) were from women who had abnormal cervical cytology. Fifteen SNPs in the LCR and eight in the E6 region were found to be the most common in all sequences. The phylogenetic analysis determined that 45% of the isolates belonged to the A1 sublineage (European variant), 34% belonged to the C1 sublineage (African 1 variant), 16% belonged to the B1 and B2 sublineage (African 2 variant), two isolates belonged to the D1-3 sublineages (Asian-American variant), and one to the North American variant.

CONCLUSIONS

The African and European HPV16 variants were the most common circulating lineages in South African and Mozambican women. A high-grade squamous intraepithelial lesion (HSIL) was the most common cervical abnormality observed and linked to European and African lineages. These findings may contribute to understanding molecular HPV16 epidemiology in South Africa and Mozambique.

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

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

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.

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

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

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Human papillomavirus as a driver of head and neck cancers

The human papillomavirus (HPV) family includes more than 170 different types of virus that infect stratified epithelium. High-risk HPV is well established as the primary cause of cervical cancer, but in recent years, a clear role for this virus in other malignancies is also emerging. Indeed, HPV plays a pathogenic role in a subset of head and neck cancers-mostly cancers of the oropharynx-with distinct epidemiological, clinical and molecular characteristics compared with head and neck cancers not caused by HPV. This review summarises our current understanding of HPV in these cancers, specifically detailing HPV infection in head and neck cancers within different racial/ethnic subpopulations, and the differences in various aspects of these diseases between women and men. Finally, we provide an outlook for this disease, in terms of clinical management, and consider the issues of 'diagnostic biomarkers' and targeted therapies.

Analysis of Human Papillomavirus (HPV) 16 Variants Associated with Cervical Infection in Italian Women

This study aims to evaluate HPV16 variants distribution in a population of Italian women living in two different regions (Lombardy and Sardinia) by sequence analyses of HPV16-positive cervical samples, in order to reconstruct the phylogenetic relationship among variants to identify the currently circulating lineages. Analyses were conducted starting from DNA isolated from 67 HPV16-positive cervical samples collected from two different Italian centres (31 from Lombardy and 36 from Sardinia) of women with normal and abnormal cervical cytology. The entire long control region (LCR) and 300 nt of the E6 gene was sequenced to identify intra-type variants. Sequence comparison and phylogenetic analysis were made using a distance-based neighbour joining method (NJ) and Kimura two-parameter model. Data obtained reported that Italian sequences mainly belonged to the European lineage, in particular sublineage A2. Only five sequences clustered in non-European branches: two in North American lineage (sublineage D1), two in African-1 (sublineage B1) and one in African-2. A new 27 nucleotide duplication in the central segment of the LCR region was found in a sequence obtained from a sample isolated in Sardinia. A predominance of European variants was detected, with some degree of variability among the studied HPV16 strains. This study contributes to the implementation of data regarding the molecular epidemiology of HPV16 variants.

Genetic characterization of variants of HPV‑16, HPV‑18 and HPV‑52 circulating in Italy among general and high‑risk populations

Viral factors, such as high‑risk human papillomavirus variants, can increase the risk of viral persistence and influence the progression to cancer. In the present study, the long control region (LCR) of human papillomavirus (HPV)‑16 and HPV‑52, and the L1 region of HPV‑16 and HPV‑18, identified from subjects belonging to both general and high‑risk populations (migrants, HIV+ subjects and adolescent/young people) residing in Italy, were characterized using molecular and phylogenetic techniques. To the best of our knowledge, this is the first Italian study to analyze a large number of sequences (n=458) and report phylogenetic data on the HPV‑52 variants. The phylogenetic analysis showed that 90% of the LCR variants of HPV‑16 and HPV‑52 clustered within lineage A (European lineage) and only sequences identified from subjects belonging to high‑risk populations fell into the non‑European lineages. Analysis of the LCRs revealed a high genomic diversity with a large number of changes. Several mutations in the binding sites for viral and cellular transcription factors characterized the HPV‑16 LCR variants belonging to the African lineages B and C, were observed in subjects with cytological abnormalities (high squamous intraepithelial lesions). The HPV‑16 and HPV‑18 L1 molecular characterization identified 30% of changes in the immune‑dominant epitope loops. These data give a clear picture of the situation in Italy, and a starting point for understanding the molecular pathogenesis and developing molecular diagnostics for HPV, vaccines and other therapeutic approaches in order to control and/or eliminate virus‑induced diseases.

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.

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.

Genetic characterization and clinical implications of human papillomavirus type 16 (HPV16) variants from northeastern Argentina

BACKGROUND

Human papillomavirus type 16 (HPV16) plays a central role in the development of cervical cancer. Worldwide studies indicate the existence of HPV16 variants that show different geographic distributions and oncogenic potential.

OBJECTIVE

Our goal was to describe the genetic variation of HPV16 isolates identified in urban women with different grades of cervical lesions living in northeastern Argentina.

STUDY DESIGN

We analyzed 116 HPV16-positive cervical samples (16 NLIM, 62 L-SIL, 16 H-SIL and 22 cervical cancer) from patients attending health centers in Misiones (Argentina) during 2006-13. HPV16 isolates were genetically characterized through PCR amplification and direct sequencing of 364 bp within the long control region, and the resulting sequences classified into variants based on phylogenetic analysis (lineages A, B, C and D). A potential association between HPV16 variants and lesion grade was evaluated through an odds ratio (OR) test. A temporal framework for the origin of HPV16 variants was assessed through coalescence analysis (BEAST v 1.7.5).

RESULTS

Phylogenetic analysis of HPV16 sequences showed that 92.1% of the samples clustered with lineage A, and 6.9% to lineage D. HPV16 variants from lineage D were more frequently associated with high-grade lesions and cancer (HSIL+) than lineage A variants at an OR of 13.8 (1.6-117.0). The time to most common recent ancestor (tMCRA) of all variants was 119,103 years before present (HPD 95%=48,486-197,239), a date consistent with the time frame for modern human evolution.

CONCLUSION

Our results suggest that HPV16 variants from lineage D may represent an additional risk factor for the development of cervical cancer in women living in northeastern Argentina. This study provides new information about viral isolates present in Argentina that will contribute to the monitoring of HPV16 infection in the vaccine era.