Intratype HPV16 sequence variation within LCR of isolates from asymptomatic carriers and cervical cancers

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.

Detailed Characteristics of Tonsillar Tumors with Extrachromosomal or Integrated Form of Human Papillomavirus

The human papillomavirus (HPV) integration, the critical step in viral carcinogenesis, most frequently occurs in the E2 gene, which results in its inactivation and in an increase of E6/E7 transcription. However, in a substantial number of tumors, the virus is present in an extrachromosomal form. For those tumors, the transformation mechanisms are not fully elucidated. Here we evaluated the possible mechanism of inactivating the E2 without interruption of the gene, methylation or mutation of the E2 binding sites (E2BSs) in HPV16-positive tonsillar tumors by next-generation and Sanger sequencing. Viral genome status was analyzed by the amplification of papillomavirus oncogene transcripts assay (APOT) and mRNA mapping, and expression of viral oncogenes was performed by qPCR. The methylation of E2BSs was significantly higher in tumors with an integrated, in comparison to extrachromosomal, form of the viral genome. No mutations were detected in the E2BSs. The viral oncogenes were equally expressed in samples with an integrated and extrachromosomal form of the virus. Only the nucleotide variants were identified in the E2 gene. No proposed mechanism of E2 inactivation was confirmed in tonsillar tumors with an extrachromosomal form of the HPV genome. The expression of E6/E7 genes seems to be sufficient to initiate and maintain the carcinogenic process.

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.

Unique LCR variations among lineages of HPV16, 18 and 45 isolates from women with normal cervical cytology in Ghana

Background

In addition to being useful for classification, sequence variations of human Papillomavirus (HPV) genotypes have been implicated in differential oncogenic potential and a differential association with the different histological forms of invasive cervical cancer. These associations have also been indicated for HPV genotype lineages and sub-lineages. In order to better understand the potential implications of lineage variation in the occurrence of cervical cancers in Ghana, we studied the lineages of the three most prevalent HPV genotypes among women with normal cytology as baseline to further studies.

Methods

Of previously collected self- and health personnel-collected cervical specimen, 54, which were positive for HPV16, 18 and 45, were selected and the long control region (LCR) of each HPV genotype was separately amplified by a nested PCR. DNA sequences of 41 isolates obtained with the forward and reverse primers by Sanger sequencing were analysed.

Results

Nucleotide sequence variations of the HPV16 genotypes were observed at 30 positions within the LCR (7460 – 7840). Of these, 19 were the known variations for the lineages B and C (African lineages), while the other 11 positions had variations unique to the HPV16 isolates of this study. For the HPV18 isolates, the variations were at 35 positions, 22 of which were known variations of Africa lineages and the other 13 were unique variations observed for the isolates obtained in this study (at positions 7799 and 7813). HPV45 isolates had variations at 35 positions and 2 (positions 7114 and 97) were unique to the isolates of this study.

Conclusion

This study provides the first data on the lineages of HPV 16, 18 and 45 isolates from Ghana. Although the study did not obtain full genome sequence data for a comprehensive comparison with known lineages, these genotypes were predominately of the Africa lineages and had some unique sequence variations at positions that suggest potential oncogenic implications. These data will be useful for comparison with lineages of these genotypes from women with cervical lesion and all the forms of invasive cervical cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-017-0755-z) contains supplementary material, which is available to authorized users.

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

•

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

•

98% of variants were of the EUR lineage.

•

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

•

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

Objective

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

Methods

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

Results

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

Conclusions

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

Function of YY1 in Long-Distance DNA Interactions

During B cell development, long-distance DNA interactions are needed for V(D)J somatic rearrangement of the immunoglobulin (Ig) loci to produce functional Ig genes, and for class switch recombination (CSR) needed for antibody maturation. The tissue-specificity and developmental timing of these mechanisms is a subject of active investigation. A small number of factors are implicated in controlling Ig locus long-distance interactions including Pax5, Yin Yang 1 (YY1), EZH2, IKAROS, CTCF, cohesin, and condensin proteins. Here we will focus on the role of YY1 in controlling these mechanisms. YY1 is a multifunctional transcription factor involved in transcriptional activation and repression, X chromosome inactivation, Polycomb Group (PcG) protein DNA recruitment, and recruitment of proteins required for epigenetic modifications (acetylation, deacetylation, methylation, ubiquitination, sumoylation, etc.). YY1 conditional knock-out indicated that YY1 is required for B cell development, at least in part, by controlling long-distance DNA interactions at the immunoglobulin heavy chain and Igκ loci. Our recent data show that YY1 is also required for CSR. The mechanisms implicated in YY1 control of long-distance DNA interactions include controlling non-coding antisense RNA transcripts, recruitment of PcG proteins to DNA, and interaction with complexes involved in long-distance DNA interactions including the cohesin and condensin complexes. Though common rearrangement mechanisms operate at all Ig loci, their distinct temporal activation along with the ubiquitous nature of YY1 poses challenges for determining the specific mechanisms of YY1 function in these processes, and their regulation at the tissue-specific and B cell stage-specific level. The large numbers of post-translational modifications that control YY1 functions are possible candidates for regulation.

Susceptibility to cervical cancer: an overview

Cervical cancer is the second most common cancer in females worldwide. It is well-established that Human Papillomavirus (HPV) infections play a critical role in the development of cervical cancer. However, a large number of women infected with oncogenic HPV types will never develop cervical cancer. Thus, there are several external environment and genetic factors involved in the progression of a precancerous lesion to invasive cancer. In this review article, we addressed possible susceptible phenotypes to cervical cancer, focusing on host genome and HPV DNA variability, multiple HPV infections, co-infection with other agents, circulating HPV DNA and lifestyle.

Human papillomavirus 16 variants from Zambian women with normal pap smears

Human papillomavirus (HPV) type 16 is the most prevalent high-risk viral genotype associated with cervical cancer. Six distinct phylogenetic clusters of HPVs have been identified and are distributed differently across five continents. HPV16 DNA was extracted from cervicolavage samples from women with normal pap smears. The LCR regions were amplified in triplicate, cloned, sequenced, and analyzed from a total of 11 recovered HPV16 positive samples [Ng'andwe et al. (2007): BMC Infect Dis 7:77] were analyzed for sequence variation. The HPV16 LCR variants were assessed for promoter activity by use of a luciferase reporter gene. Six novel HPV16 variants with nucleotide exchanges in the LCR region were identified. Five clones were classified as European group HPV16 variants and one as an African group variant. Two of these variants had relatively lower promoter activity, 30% of that of the wild-type strain. The decreased promoter activity of some HPV16 variants may decrease expression of viral oncogenes and may be linked with the development, phenotype and severity of the cervical lesions in women infected with these across HPV16 variants.

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

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

Cellular factor YY1 downregulates the human papillomavirus 16 E6/E7 promoter, P97, in vivo and in vitro from a negative element overlapping the transcription-initiation site

Cellular factors that bind to cis sequences in the human papillomavirus 16 (HPV-16) upstream regulatory region (URR) positively and negatively regulate the viral E6 and E7 oncogene promoter, P97. DNase I footprinting has revealed the binding of cellular proteins to two previously undetected cis elements overlapping and 3′ of the transcription-initiation site of the P97 promoter. Mutations within homologous motifs found in both of these cis elements abolished their negative function in vivo and the binding of the same cellular complex in vitro. This factor was identified as YY1 by complex mobility and binding specificity in comparison with vaccinia virus-expressed, purified recombinant YY1 protein and by antigenic reactivity with YY1 antisera. Cis mutations in the ‘initiator’ YY1 site activated the P97 promoter in vivo and in vitro. P97 was also activated threefold in vitro by depletion of endogenous YY1 with wild-type, but not mutant, YY1 oligonucleotides from the IgH kappa E3′ enhancer. Furthermore, increasing concentrations of exogenous, purified recombinant YY1 repressed wild-type P97 transcript levels by up to threefold, but did not influence the P97 promoter mutated in the ‘initiator’ YY1 site. Thus, the promoter-proximal YY1 site was not necessary for correct transcription initiation at the P97 promoter, but was found to be required for downregulation of P97 transcription in vivo and in vitro. In contrast to other viral and cellular promoters, where YY1 is thought to function as a positive transcription-‘initiator’ factor, HPV-16 P97 transcription is downregulated by YY1 from a critical motif overlapping the transcription start site.

Upstream regulatory region alterations found in human papillomavirus type 16 (HPV-16) isolates from cervical carcinomas increase transcription, ori function, and HPV immortalization capacity in culture

Human papillomavirus (HPV) DNAs isolated from cervical and head and neck carcinomas frequently contain nucleotide sequence alterations in the viral upstream regulatory region (URR). Our study has addressed the role such sequence changes may play in the efficiency of establishing HPV persistence and altered keratinocyte growth. Genomic mapping of integrated HPV type 16 (HPV-16) genomes from 32 cervical cancers revealed that the viral E6 and E7 oncogenes, as well as the L1 region/URR, were intact in all of them. The URR sequences from integrated and unintegrated viral DNA were found to harbor distinct sets of nucleotide substitutions. A subset of the altered URRs increased the potential of HPV-16 to establish persistent, cell growth-altering viral-genome replication in the cell. This aggressive phenotype in culture was not solely due to increased viral early gene transcription, but also to augmented initial amplification of the viral genome. As revealed in a novel ori-dependent HPV-16 plasmid amplification assay, the altered motifs that led to increased viral transcription from the intact genome also greatly augmented HPV-16 ori function. The nucleotide sequence changes correlate with those previously described in the distinct geographical North American type 1 and Asian-American variants that are associated with more aggressive disease in epidemiologic studies and encompass, but are not limited to, alterations in previously characterized sites for the negative regulatory protein YY1. Our results thus provide evidence that nucleotide alterations in HPV regulatory sequences could serve as potential prognostic markers of HPV-associated carcinogenesis.

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

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

Identification and functional analysis of sequence variants in the long control region and the E2 open reading frame of bovine papillomavirus type 1 isolated from equine sarcoids

BPV-1 DNA is the predominant viral type detected in equine sarcoids and represents the only reported natural cross species infection of papillomaviruses. In this study, nucleotide variations in the LCR and the E2 regions of equine sarcoid-associated BPV-1 were characterised by sequence analysis. Variants particular to sarcoid BPV-1 were identified in both the LCR and E2 sequence. The functionality of the most common LCR variant was examined in equine and bovine cells. These studies showed that the activity of the variant LCR was higher in equine cells than bovine cells; the activity of the variant LCR in the presence of the E2 variant was similar to the reference/wild-type sequences in equine cells, whereas in bovine cells the variant function was reduced by 50%. These data suggest the viral BPV variants commonly detected in sarcoids have an enhanced function in equine cells compared to their function in bovine cells.

New HPV-16 European and non-European variants in Central Brazil

HPV-16 is the most prevalent human papillomavirus genotype found in cervical intraepithelial neoplasias. The regulatory region of the HPV genome, LCR, has several binding sites for cellular and viral transcription factors, and nucleotide substitutions in this genomic region can interfere with the viral oncogenic expression. The present study aims to determine the LCR variability of European and non-European HPV-16 variants found in Brazil. Through automated sequencing, it was possible to characterize the LCR of ten non-European (eight Asian-American, one African 1, one African 2) and twelve European isolates. Among the 22 isolates analyzed, nine may be new variants of HPV-16, with different combinations of previously reported nucleotide substitutions, and three showed new substitutions not previously reported. Two new nucleotide substitutions, the insertion of T at position 7621 and the substitution of A to G at position 7836, were found in a single isolate, Bsb-14, a putative new African 1 variant. The characterization of the LCR of human papillomaviruses can be of pivotal importance to the understanding of the viral replication and pathogenicity.

HPV16 E2 gene disruption and polymorphisms of E2 and LCR: some significant associations with cervical cancer in Indian women

OBJECTIVES

We evaluated the status of the HPV16 E2 gene (disrupted or intact), nucleotide sequence alterations within intact E2 genes and LCR of HPV16 isolates in a group of CaCx cases (invasive squamous cell carcinomas, n = 81) and population controls (normal cervical scrapes, n = 27) from Indian women.

METHODS

E2 disruption was detected by amplifying the entire E2 gene with single set of primers, while overlapping primers were used to determine if any particular region got selectively disrupted. Nucleotide variations in E2 and LCR were analyzed by PCR amplification followed by bi-directional sequencing. The associations between the viral factors and CaCx were analyzed using Fisher's Exact or Chi-squared test and interpreted as OR (95% CI) and P values.

RESULTS

E2 disruption was significantly higher among the cases [3.38 (1.07-10.72); P = 0.02], which was maximum in the region between nucleotides 3650 and 3872 (DNA-binding region). The European (E) variant was found to be the prevalent subgroup (87.76% among cases and 96.30% among the controls), and the remaining samples were Asian-American variants. Among the E subgroup, variation at position 7450 (T > C) within the E2-binding site-IV was found to be significantly higher among the E2 undisrupted cases (21/37; 56.76%), compared to controls (5/18; 27.78%) [3.41 (1.01-11.55); P = 0.03].

CONCLUSIONS

Besides HPV16 E2 disruption, LCR 7450T > C variation within undisrupted E2 of E subgroup appears to be a major factor contributing to the risk of CaCx development in Indian women. Furthermore, polymorphisms in the E2 gene of HPV16 may not be significant for disease risk.

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

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