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

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

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

Identification of novel human papillomavirus lineages and sublineages in HIV/HPV-coinfected pregnant women by next-generation sequencing

Infection by human papillomavirus (HPV) is a necessary condition for development of cervical cancer, and has also been associated with malignancies of other body anatomical sites. Specific HPV types have been associated with premalignant lesions and invasive carcinoma, but mounting evidence suggests that within-type lineages and sublineages also display distinct biological characteristics associated with persistent infections and evolution to cervical cancer. In the present study, we have assessed HPV multiple infection and variation from a cohort of highly susceptible, HIV(+) pregnant women using next-generation sequencing and an in-house pipeline for HPV full-length genome assembly. Seventy-two consensus sequences representing complete or near-complete (>97%) HPV genomes were assembled, spanning 28 different types. Genetic distance and phylogenetic analyses allowed us to propose the classification of novel HPV lineages and sublineages across nine HPV types, including two high-risk types. HPV diversity may be a hallmark of immunosuppressed patients upon HIV infection and AIDS progression.

The Asian-American variant of human papillomavirus type 16 exhibits higher activation of MAPK and PI3K/AKT signaling pathways, transformation, migration and invasion of primary human keratinocytes

Asian-American (AA) HPV-16 variants are associated with higher risk of cancer. Abnormal activation of intracellular signaling play a critical role in cancer development and progression. Our aim was to elucidate mechanisms underlying the higher oncogenic potential attributed to AA variant. We evaluated activation of MAPK and PI3K/AKT pathways in primary human keratinocytes (PHKs) transduced with E6/E7 of three HPV-16 variants: E-P, AA, E-350G. Phenotypes examined included migration, anchorage independent growth and invasion. AA PHKs presented the highest levels of active proteins involved in all cascades analyzed: MAPK-ERK, MAPK-p38 and PI3K-AKT. AA PHKs were more efficient in promoting anchorage independent growth, and in stimulating cell migration and invasion. MEK1 inhibition decreased migration. The mesenchymal phenotype marker vimentin was increased in AA PHKs. Our results suggest that MEK1, ERK2, AKT2 hyperactivation influence cellular behavior by means of GSK-3b inactivation and EMT induction prompting AA immortalized PHKs to more efficiently surpass carcinogenesis steps.

Genetic diversity of HPV16 and HPV18 in Brazilian patients with invasive cervical cancer

Cervical cancer is the fourth most common cancer among women, and ∼70-80% of these cancers are associated with two human papillomavirus types: HPV16 and HPV18. Several studies have reported that intra-type diversity is associated with the progression of infection to invasive cancer. Herein, we report the genetic diversity of HPV16 and HPV18 in a cohort of 594 Brazilian women with invasive cervical cancer and describe the prevalence of lineages and intra-type diversity prior to the implementation of the public immunization program in Brazil. HPV detection and genotyping were performed using PCR, PGMY/GP primers, and DNA extracted from fresh tumors. The HPV16 (378 women) and HPV18 (80 women) lineages were identified by PCR and sequencing of the LCR and E6 fragments, followed by SNV comparison and phylogenetic analysis. In our cohort, was found a higher frequency of the lineage A (in 217 women), followed by lineage D (in 97 women) and lineages B and C (in 10 women each) for HPV16; and a higher frequency of lineage A (in 56 women) followed by lineage B (in 15 women) in HPV18. The genetic diversity of HPV16 indicated a recent expansion of specific variants or a selective advantage that is associated with invasive cancer; this pattern was not observed for HPV18.

Distribution of human papillomavirus type 16 variants in Lithuanian women with cervical cancer

BACKGROUND AND OBJECTIVE

Cervical cancer usually is caused by HPV 16. However, HPV 16 varies within type; different genotypes are described as prototype or variants. Prevalence of different variants differ according the geographic regions and has an unequal impact for cervical cancer development. Our study aimed to identify which variant of HPV 16 was most prevalent in biological samples taken from Lithuanian women with cervical cancer.

MATERIALS AND METHODS

A total of 122 HPV 16 positive cervical samples (invasive cancer and cervical intraepithelial neoplasia) were investigated and sequenced to identify different variants. HPV 16 was detected using type specific PCR, exact sequence of the virus was obtained by viral DNA sequencing.

RESULTS

Adequate HPV sequence was detected in 106 cases from 122 (86.9% of all cases). After histological confirmation, 96 cases were included in the final analysis. In 33 cases (34.4%) HPV 16 prototype was detected; in 50 cases (52.1%), L83V variant; and in remaining 13 cases (13.5%), multivariant of HPV 16. The frequency of L83V variant in invasive cancer and carcinoma in situ samples was the same (66.7% and 62.0%, respectively; P=0.696). Of analyzed multivariants, 10 were attributed to the European phylogenetic line; 1, to the North American, and 1, to the Asian-American. One sample was not attributed to any of the known phylogenetic lines.

CONCLUSIONS

The European HPV 16 L83V variant is usually associated with high risk of cervical cancer among women. However, statistically significant difference was not achieved when comparing difference of L83V variants between investigated groups and in HPV 16 L83V variant and prototype distribution in CIN3/Ca in situ and cancer.

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.

Association Study between Cervical Lesions and Single or Multiple Vaccine-Target and Non-Vaccine Target Human Papillomavirus (HPV) Types in Women from Northeastern Brazil

We performed an association between high-grade squamous intraepithelial lesions (HSIL), low-grade squamous intraepithelial lesions (LSIL) and single or multiple vaccine-target as well as non-vaccine target Human papillomavirus (HPV) types. Using bead-based HPV genotyping, 594 gynecological samples were genotyped. An association between squamous intraepithelial lesion (SIL) and presence of HPV16, 18, 31, 58 and 56 types were calculated. The risk was estimated by using odds ratio (OR) and 95% of confidence intervals (CI). A total of 370 (62.3%) women were HPV positive. Among these, 157 (42.7%) presented a single HPV infection, and 212 (57.3%) were infected by more than one HPV type. HPV31 was the most prevalent genotype, regardless single and multiple HPV infections. Single infection with HPV31 was associated with LSIL (OR=2.32; 95%CI: 1.01 to 5.32; p=0.04); HPV31 was also associated with LSIL (OR=3.28; 95%CI: 1.74 to 6.19; p= 0.0002) and HSIL (OR=3.82; 95%CI: 2.10 to 6.97; p<0.001) in multiple HPV infections. Risk to harbor cervical lesions was observed in multiple HPV infections with regard to the HPV56 (OR=5.39; 95%CI: 2.44 to 11.90; p<0.001for LSIL; OR=5.37; 95%CI: 2.71 to 10.69; p<0.001) and HPV58 (OR=3.29; 95%CI: 1.34 to 8.09; p=0.0091 for LSIL; OR=3.55; 95%CI: 1.56 to 8.11; p=0.0026) genotypes. In addition, women coinfected with HPV16/31/56 types had 6 and 5-fold increased risk of HSIL (OR=6.46; 95%CI: 1.89 to 22.09; p=0.002) and LSIL (OR=5.22; 95%CI: 1.10 to 24.70; p=0.03), respectively. Multiple HPV infections without HPV16/18 has 2-fold increased risk of HSIL (OR=2.57; 95%CI: 1.41 to 4.70; p=0.002) and LSIL OR=2.03; 95%CI: 1.08 to 3.79; p=0.02). The results of this study suggest that single and multiple vaccine target as well as non-vaccine target HPV types are associated with LSIL and HSIL. These finding should be taken into consideration in the design of HPV vaccination strategies.

High-level of viral genomic diversity in cervical cancers: A Brazilian study on human papillomavirus type 16

Invasive cervical cancer (ICC) is the third most frequent cancer among women worldwide and is associated with persistent infection by carcinogenic human papillomaviruses (HPVs). The combination of large populations of viral progeny and decades of sustained infection may allow for the generation of intra-patient diversity, in spite of the assumedly low mutation rates of PVs. While the natural history of chronic HPVs infections has been comprehensively described, within-host viral diversity remains largely unexplored. In this study we have applied next generation sequencing to the analysis of intra-host genetic diversity in ten ICC and one condyloma cases associated to single HPV16 infection. We retrieved from all cases near full-length genomic sequences. All samples analyzed contained polymorphic sites, ranging from 3 to 125 polymorphic positions per genome, and the median probability of a viral genome picked at random to be identical to the consensus sequence in the lesion was only 40%. We have also identified two independent putative duplication events in two samples, spanning the L2 and the L1 gene, respectively. Finally, we have identified with good support a chimera of human and viral DNA. We propose that viral diversity generated during HPVs chronic infection may be fueled by innate and adaptive immune pressures. Further research will be needed to understand the dynamics of viral DNA variability, differentially in benign and malignant lesions, as well as in tissues with differential intensity of immune surveillance. Finally, the impact of intralesion viral diversity on the long-term oncogenic potential may deserve closer attention.

Sites of disruption within E1 and E2 genes of HPV16 and association with cervical dysplasia

Integration of HPV16 DNA into the host chromosome usually disrupts the E1 and/or E2 genes. The present study investigated the disruption of E1, E2 genes in a total of eighty four HPV16-positive precancerous and cervical cancer specimens derived from Greek women (seventeen paraffin-embedded cervical biopsies and sixty seven Thin Prep samples). Complete E2 and E1 genes were amplified using three and nine overlapping primer sets respectively, in order to define the sites of disruption. Extensive mapping analysis revealed that disruption/deletion events within E2 gene occurred in high grade and cervical cancer samples (x(2) test, P < 0.01), while no evidence of E2 gene disruption was documented among low grade cervical intraepithelial neoplasias. In addition, disruptions within the E1 gene occur both in high and low grade cervical intraepithelial neoplasia. This leads to the assumption that in low grade cervical intraepithelial neoplasias only E1 gene disruption was involved (Fisher's exact test, P < 0.05), while in high grade malignancies and cervical cancer cases deletions in both E1 and E2 genes occurred. Furthermore, the most prevalent site of disruption of E1 gene was located between nucleotides 1059 and 1323, while the most prevalent deleted region of the E2 gene was located between nucleotides 3172 and 3649 (E2 hinge region). Therefore, it is proposed that each population has its own profile of frequencies and sites of disruptions and extensive mapping analysis of E1 and E2 genes is mandatory in order to determine suitable markers for HPV16 DNA integration analysis in distinct populations.

Genetic variability in E6 and E7 oncogenes of human papillomavirus Type 16 from Congolese cervical cancer isolates

Background

The molecular epidemiological studies showed that some variants of HPV-16, distributed geographically, would present a higher risk of causing cervical cancer. This study aimed to analyze nucleotide changes of HPV-16 E6 and E7 genomic regions from infected Southwestern Congolese women.

Methods

DNA of twenty HPV-16 isolates was analyzed by amplifying the E6 and E7 genes using type-specific primers PCR and direct sequencing. The sequences obtained were aligned with the HPV-16 GenBank reference sequences.

Results

Thirteen (65.0%) out of 20 DNA-samples were successfully amplified. Genetic analysis revealed 18 and 4 nucleotide changes in E6 and E7 genomic regions respectively. The most frequently observed nucleotide variations were the missense C143G, G145T and C335T in E6 (100%), leading to the non-synonymous amino acid variation Q14D and H78Y. E7 genomic region was found to be highly conserved with two most common T789C and T795G (100%) silent variations. All HPV-16 variants identified belonged to the African lineage: 7 (53.8%) belonged to Af-1 lineage and 6 (46.1%) to Af-2 lineage. The missense mutation G622A (D21N) in the E7 region seems to be described for the first time in this study.

Conclusion

This study reported for the first time the distribution of HPV-16 E6 and E7 genetic variants in infected women from southwest Congo. The findings confirmed almost ascendancy of the African lineage in our study population.

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

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

Prevalence of human papillomavirus variants and genetic diversity in the L1 gene and long control region of HPV16, HPV31, and HPV58 found in North-East Brazil

This study showed the prevalence of human papillomavirus (HPV) variants as well as nucleotide changes within L1 gene and LCR of the HPV16, HPV31, and HPV58 found in cervical lesions of women from North-East Brazil.

The role of the globular heads of the C1q receptor in HPV-16 E2-induced human cervical squamous carcinoma cell apoptosis via a mitochondria-dependent pathway

Background

Human papillomavirus type-16 (HPV-16) E2 protein acts as a transcriptional modulator and plays a key role in regulating many biological responses. The purpose of this study was to investigate the relationship between HPV-16 E2, the receptor for the globular heads of human C1q (gC1qR) gene expression, mitochondrial dysfunction and apoptosis regulation in human cervical squamous carcinoma cells (C33a and SiHa).

Methods

HPV-16 E2 and gC1qR expression was examined using real-time PCR and western blot analysis. Apoptosis in C33a and SiHa cells was assessed by flow cytometry. Mitochondrial function was detected via ROS generation, the amount of cytosolic Ca²⁺, and changes in the mitochondrial membrane potential (Δψm).

Results

The expression of the HPV-16 E2 and gC1qR gene significantly decreased in human cervical squamous carcinoma samples relative to the non-cancerous cervix samples. C33a and SiHa cells that were transfected with a vector encoding HPV-16 E2 displayed significantly increased gC1qR gene expression and mitochondrial dysfunction as well as an up-regulation of cellular apoptosis, which was abrogated by the addition of gC1qR small-interfering RNA (siRNA).

Conclusions

These data support a mechanism whereby gC1qR plays an important role in HPV-16 E2-induced human cervical squamous carcinoma cell apoptosis via a mitochondria-dependent pathway.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-014-0286-y) contains supplementary material, which is available to authorized users.

Human papillomavirus induced transformation in cervical and head and neck cancers

Human papillomavirus (HPV) is one of the most widely publicized and researched pathogenic DNA viruses. For decades, HPV research has focused on transforming viral activities in cervical cancer. During the past 15 years, however, HPV has also emerged as a major etiological agent in cancers of the head and neck, in particular squamous cell carcinoma. Even with significant strides achieved towards the screening and treatment of cervical cancer, and preventive vaccines, cervical cancer remains the leading cause of cancer-associated deaths for women in developing countries. Furthermore, routine screens are not available for those at risk of head and neck cancer. The current expectation is that HPV vaccination will prevent not only cervical, but also head and neck cancers. In order to determine if previous cervical cancer models for HPV infection and transformation are directly applicable to head and neck cancer, clinical and molecular disease aspects must be carefully compared. In this review, we briefly discuss the cervical and head and neck cancer literature to highlight clinical and genomic commonalities. Differences in prognosis, staging and treatment, as well as comparisons of mutational profiles, viral integration patterns, and alterations in gene expression will be addressed.