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

Molecular epidemiology of human papillomavirus variants in cervical cancer in India

Background & objectives Cervical cancer (CC) has been documented as the fourth most common cancer worldwide. Persistent infections with high-risk human papillomavirus (hr-HPV) have been suggested in the development of CC. Although prophylactic vaccines are available for the prevention of prevalent hr-HPV types, intra-type variations exist within a particular HPV type that has varying oncogenic potential as well as the mechanism of pathogenicity and varying neutralization by antibodies. Therefore, we carried out a systematic review to determine the distribution of HPV intra-typic variations in different geographical locations of India and their reported implications. Methods Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) guidelines were followed to retrieve relevant articles from the standard databases using appropriate keywords. Consequently, 17 articles were included in the current review after screening based on inclusion and exclusion criteria. Results The majority of articles included in this review reported variations within the HPV16 E6 gene, followed by the L1 and E7 genes. Analysis of available data indicated the differential regional distribution of some variations. These variations have also been reported to impact the biological functions of various viral proteins. Interpretation & conclusions The distribution of lineages varied with the different genomic regions sequenced. Additionally, there were certain unique and common variations in the HPV genome with respect to geographical regions. Hence, we suggest the identification of region-specific variations for the development of diagnostic and prognostic interventions.

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.

Predominance of genomically defined A lineage of HPV16 over D lineage in Indian patients from eastern India with squamous cell carcinoma of the cervix in association with distinct oncogenic phenotypes

•

The HPV16 lineage A viruses are predominantly associated with cervical squamous cell carcinoma (SCC) in the eastern region of India.

•

Among lineage A viruses, the E5(Y44L, I65V), E6(L83V) and LCR: C7577T variants are associated with SCC.

•

The risk alleles are distributed among 10 clades comprising of 64 HPV16 genomes.

•

Low integration, high episomal copy numbers and high E7 mRNA expression in SCC contrasts lineage A viruses from lineage D.

Human papillomavirus type-16 (HPV16) is classified into lineages, A, B, C and D and 10 sub-lineages portraying variable infectivity, persistence, and cytological outcomes, however, with geographical variations. Our objective was to delineate the distinctive features of lineages among cervical squamous cell carcinoma (SCC) in the eastern region of India. A total of 145 SCC cases and 24 non-malignant specimens, harboring episomal HPV16, were included. The presence of higher proportion of lineage A over D was observed among SCC cases (86.89% A1, 8.97% D1 and 4.14% D2), while only A1 sub-lineage viruses were found among control specimens. Among the A1 viruses, an association of variants in the E5 (Y44L, I65V), E6 (L83V) genes and LCR: C7577T with SCC, with combined Odd's ratio (95% CI) of 20.5(4.61–91.25) was observed. Network analyses revealed the presence of 10 clades of lineage A viruses comprising of 64 HPV16 genomes harboring the risk alleles. High episomal HPV16 DNA copy numbers (adjusted p-value= 0.0271) and E7 mRNA expression (p-value=0.000017) predominated in SCC with lineage A, over D. Our study highlights the distinctive modalities of oncogenicity among different HPV16 lineages.

Cervical cancer subtypes harbouring integrated and/or episomal HPV16 portray distinct molecular phenotypes based on transcriptome profiling of mRNAs and miRNAs

Heterogeneity in cervical cancers (CaCx) in terms of HPV16 physical status prompted us to investigate the mRNA and miRNA signatures among the different categories of CaCx samples. We performed microarray-based mRNA expression profiling and quantitative real-time PCR-based expression analysis of some prioritised miRNAs implicated in cancer-related pathways among various categories of cervical samples. Such samples included HPV16-positive CaCx cases that harboured either purely integrated HPV16 genomes (integrated) and those that harboured episomal viral genomes, either pure or concomitant with integrated viral genomes (episomal), which were compared with normal cervical samples that were either HPV negative or positive for HPV16. The mRNA expression profile differed characteristically between integrated and episomal CaCx cases for enriched biological pathways. miRNA expression profiles also differed among CaCx cases compared with controls (upregulation—miR-21, miR-16, miR-205, miR-323; downregulation—miR-143, miR-196b, miR-203, miR-34a; progressive upregulation—miR-21 and progressive downregulation—miR-143, miR-34a, miR-196b and miR-203) in the order of HPV-negative controls, HPV16-positive non-malignant samples and HPV16-positive CaCx cases. miR-200a was upregulated in HPV16-positive cervical tissues irrespective of histopathological status. Expression of majority of the predicted target genes was negatively correlated with their corresponding miRNAs, irrespective of the CaCx subtypes. E7 mRNA expression correlated positively with miR-323 expression among episomal cases and miR-203, among integrated cases. miR-181c expression was downregulated only among the episomal CaCx cases and negatively correlated with protein coding transcript of the proliferative target gene, CKS1B of the significantly enriched “G2/M DNA Damage Checkpoint Regulation” pathway among CaCx cases. Thus, the two CaCx subtypes are distinct entities at the molecular level, which could be differentially targeted for therapy. In fact, availability of a small molecule inhibitor of CKS1B, suggests that drugging CKS1B could be a potential avenue of treating the large majority of CaCx cases harbouring episomal HPV16.

Niche adaptation and viral transmission of human papillomaviruses from archaic hominins to modern humans

Recent discoveries on the origins of modern humans from multiple archaic hominin populations and the diversity of human papillomaviruses (HPVs) suggest a complex scenario of virus-host evolution. To evaluate the origin of HPV pathogenesis, we estimated the phylogeny, timing, and dispersal of HPV16 variants using a Bayesian Markov Chain Monte Carlo framework. To increase precision, we identified and characterized non-human primate papillomaviruses from New and Old World monkeys to set molecular clock models. We demonstrate specific host niche adaptation of primate papillomaviruses with subsequent coevolution with their primate hosts for at least 40 million years. Analyses of 212 HPV16 complete genomes and 3582 partial sequences estimated ancient divergence of HPV16 variants (between A and BCD lineages) from their most recent common ancestors around half a million years ago, roughly coinciding with the timing of the split between archaic Neanderthals and modern Homo sapiens, and nearly three times longer than divergence times of modern Homo sapiens. HPV16 A lineage variants were significantly underrepresented in present African populations, whereas the A sublineages were highly prevalent in European (A1-3) and Asian (A4) populations, indicative of viral sexual transmission from Neanderthals to modern non-African humans through multiple interbreeding events in the past 80 thousand years. Remarkably, the human leukocyte antigen B*07:02 and C*07:02 alleles associated with increased risk in cervix cancer represent introgressed regions from Neanderthals in present-day Eurasians. The archaic hominin-host-switch model was also supported by other HPV variants. Niche adaptation and virus-host codivergence appear to influence the pathogenesis of papillomaviruses.

The association of human papillomavirus type 16 E2 variations with cervical cancer in a Han Chinese population

Human papillomavirus type 16 (HPV16) is considered to be the primary pathogen related to cervical cancer. The HPV16 E2 protein plays an important role in tumourigenicity of cervical carcinoma. In the current study, we enrolled 121 HPV16-positive cervical cancer patients in the case group and 130 HPV16-positive asymptomatic individuals in the control group, and we investigated the association between HPV16 E2 gene variations and cervical cancer. The HPV16 E2 DNA was amplified and sequenced. We identified two HPV variants (EUR and As) in the control group; the As variant was predominant (68.5%), followed by the EUR variant (31.5%). In the case group, three HPV variants (EUR, As and AA) were observed; the As variant was predominant (72.7%), followed by the EUR variant (22.3%) and the AA variant (5.0%). Our results showed a significant difference in the distribution of the HPV16 variants between the case and control groups (P < 0.05). Moreover, in the HPV16 E2 gene variation analysis, the distribution of sixteen variations was significantly different between the case and control groups (P < 0.05), and all of these variations were present in the AA variant. In the subgroup analysis, the frequency of the T3575G (S274A) variation in the EUR variant was significantly different between the case and control groups (P = 0.029); however, there was no significant difference in the frequency of the variations in the As variant between the case and control groups. Our findings in the current study could provide a better understanding of the relationship between HPV16 variants, E2 gene variations and cervical cancer.

Impact of viral and host DNA methylations on HPV16-related cervical cancer pathogenesis

Epigenetic alterations within human papillomavirus (HPV) and host cellular genomes are known to occur during cervical carcinogenesis. Our objective was to analyse the influence of (1) methylation within two immunostimulatory CpG motifs within HPV16 E6 and E7 genes around the viral late promoter and their correlation, if any, with expression deregulation of host receptor (TLR9) and DNA methyltransferases (DNMT1, DNMT3A and DNMT3B) and (2) global DNA methylation levels within CpGs of the repetitive Alu sequences, on cervical cancer (CaCx) pathogenesis. Significantly higher proportions of CaCx samples portrayed methylation in immunostimulatory CpG motifs, compared to HPV16-positive non-malignant samples, with cases harbouring episomal HPV16 showing decreased methylation compared to those with viral integration. A significant linear trend of TLR9 upregulation was recorded in the order of HPV-negative controls < HPV16-positive non-malignant samples < HPV16-positive CaCx cases. TLR9 upregulation in cases with episomal HPV16 was again higher among those with non-methylated immunostimulatory CpG motifs. Comparison of cases with HPV-negative controls revealed that DNMT3A was significantly downregulated only among integrated cases, DNMT3B was significantly overexpressed among both categories of cases, although at variable levels, while DNMT1 failed to show any deregulated expression among the cases. Global host DNA hypomethylation, also showed a significant linear increasing trend through the progressive CaCx development stages mentioned above and was most prominently higher among cases with episomal HPV16 as opposed to viral integration. Thus, HPV16 and host methylations appear to influence CaCx pathogenesis, with differential molecular signatures among CaCx cases with episomal and integrated HPV16.

Disparities in Cervical Cancer Incidence and Mortality: Can Epigenetics Contribute to Eliminating Disparities?

Screening for uterine cervical intraepithelial neoplasia (CIN) followed by aggressive treatment has reduced invasive cervical cancer (ICC) incidence and mortality. However, ICC cases and carcinoma in situ (CIS) continue to be diagnosed annually in the United States, with minorities bearing the brunt of this burden. Because ICC peak incidence and mortality are 10-15 years earlier than other solid cancers, the number of potential years of life lost to this cancer is substantial. Screening for early signs of CIN is still the mainstay of many cervical cancer control programs. However, the accuracy of existing screening tests remains suboptimal. Changes in epigenetic patterns that occur as a result of human papillomavirus infection contribute to CIN progression to cancer, and can be harnessed to improve existing screening tests. However, this requires a concerted effort to identify the epigenomic landscape that is reliably altered by HPV infection specific to ICC, distinct from transient changes.

Identification of genetic variation in the lncRNA HOTAIR associated with HPV16-related cervical cancer pathogenesis

PURPOSE

Previously, over-expression of the long noncoding RNA (lncRNA) HOTAIR has been found to be associated with the invasive and metastatic capacities of several epithelial cancers, including cervical cancer (CaCx). Here, we aimed at identifying functionally relevant genetic variants that may be employed to differentiate between clinically distinct CaCx subtypes, i.e., those exhibiting high HOTAIR levels and molecular signatures of metastasis and those lacking such signatures in the presence of low HOTAIR expression levels.

METHODS

Genomic DNA isolated from various cervical tissue samples (characterized by histopathology and HPV status) was used for HOTAIR amplicon sequencing, followed by validation of the findings by Sanger sequencing. The impact of the genetic variants found on the secondary structure of HOTAIR and the concomitant alterations in miRNA binding sites were determined through in silico analysis, followed by miRNA expression analysis by quantitative real-time PCR and confirmation of miRNA binding using a luciferase reporter assay.

RESULTS

We found that rs2366152C was over-represented [OR_age-adjusted = 2.58 (1.23-5.57); p = 0.014] in low HOTAIR expressing HPV positive CaCx cases compared to HPV negative controls. This genetic variant showed the propensity of a secondary structure alteration and gain of a miR-22 binding site in HOTAIR, which was found to be concordant with miR-22 over-expression in low HOTAIR CaCx cases compared to controls. We found that miR-22 expression negatively correlated with HOTAIR and E7 expression in HPV16 positive cases and in an E7 transfected HPV negative CaCx-derived cell line (C33A), but was not altered in high HOTAIR cases compared to controls. Reduced luciferase activity of a HOTAIR rs2366152C expression plasmid in C33A cells through miR-22 co-transfection confirmed the ability of miR-22 to specifically target rs2366152C-harbouring HOTAIR lncRNA in CaCx cells, ultimately leading to its down-regulation.

CONCLUSIONS

Our data indicate that rs2366152C not only has the potential to serve as a marker for singling out CaCx cases lacking metastatic molecular signatures, but also to explain the functional inactivation of HOTAIR in these cases, including the mechanism of its down-regulation.

Complete Genome Sequences of Eight Human Papillomavirus Type 16 Asian American and European Variant Isolates from Cervical Biopsies and Lesions in Indian Women

Human papillomavirus type 16 (HPV16), a member of the Papillomaviridae family, is the primary etiological agent of cervical cancer. Here, we report the complete genome sequences of four HPV16 Asian American variants and four European variants, isolated from cervical biopsies and scrapings in India.

Bridging Links between Long Noncoding RNA HOTAIR and HPV Oncoprotein E7 in Cervical Cancer Pathogenesis

Human Papillomavirus (HPV) type 16 oncoprotein E7 plays a major role in cervical carcinogenesis by interacting with and functionally inactivating various host regulatory molecules. Long noncoding RNA (lncRNA) HOTAIR is one such regulator that recruits chromatin remodelling complex PRC2, creating gene silencing H3K27 me3 marks. Hence, we hypothesized that HOTAIR could be a potential target of E7, in HPV16 related cervical cancers (CaCx). We identified significant linear trend of progressive HOTAIR down-regulation through HPV negative controls, HPV16 positive non-malignants and CaCx samples. Majority of CaCx cases portrayed HOTAIR down-regulation in comparison to HPV negative controls, with corresponding up-regulation of HOTAIR target, HOXD10, and enrichment of cancer related pathways. However, a small subset had significantly higher HOTAIR expression, concomitant with high E7 expression and enrichment of metastatic pathways. Expression of HOTAIR and PRC2-complex members (EZH2 and SUZ12), showed significant positive correlation with E7 expression in CaCx cases and E7 transfected C33A cell line, suggestive of interplay between E7 and HOTAIR. Functional inactivation of HOTAIR by direct interaction with E7 could also be predicted by in silico analysis and confirmed by RNA-Immunoprecipitation. Our study depicts one of the causal mechanisms of cervical carcinogenesis by HPV16 E7, through modulation of HOTAIR expression and function.

Situation of HPV16 E2 gene status during radiotherapy treatment of cervical carcinoma

BACKGROUND

Human papillomavirus (HPV) integration within the E2 gene has been proposed as a critical event in cervical carcinogenesis. This study concerned whether HPV16 status and E2 gene intactness are predictive of radiation response in patients with cervical cancer.

MATERIALS AND METHODS

Biopsies of 44 patients with cervical cancer were collected before or after radiotherapy. The presence of HPV16 was assessed by polymerase chain reaction (PCR) using specific primers for the L1 region. E2 disruption was detected by amplifying the entire E2 gene.

RESULTS

HPV16 DNA was found in 54.5% of the clinical samples. Overall, 62.5% of the HPV16 positive tumors had integrated viral genome and 37.5% had episomal genome. There was a tendency of increase of HPV16 E2 negative tumors compared with HPV16 L1 ones in advanced stages (75% versus 20% in stage III respectively). Detection of E2 gene appeared influenced by the radiotherapy treatment, as the percentage of samples containing an intact HPV16 E2 was more frequent in pretreated patients compared to radiotherapy treated patients (66.6% versus 20%). The radiation therapy caused an eight-fold [OR= 8; CI=1.22-52.25; p=0.03] increase in the risk of HPV16 genome disruption. The integration status is influenced by the irradiation modalities, interestingly E2 disruption being found widely after radiotherapy treatment (75%) with a total fractioned dose of 50 Gy.

CONCLUSIONS

This study reveals that the status of the viral DNA may be used as a marker to optimize the radiation treatment.

Viral load, integration and methylation of E2BS3 and 4 in human papilloma virus (HPV) 16-positive vaginal and vulvar carcinomas

Objective

To investigate if viral load, integration and methylation of E2BS3 and 4 represent different ways of tumor transformation in vaginal and vulvar carcinoma and to elucidate its clinical impact.

Methods

Fifty-seven samples, positive for HPV16, were selected for the study. Detection of viral load was made with realtime-PCR using copy numbers of E6 and integration was calculated from comparing E2 to E6-copies. Methylation of E2BS3 and 4 was analysed using bisulphite treatment of tumor DNA, followed by PCR and pyrosequencing.

Results

Vaginal tumors were found to have a higher viral load (p = 0.024) compared to vulvar tumors but a high copy number (> median value, 15 000) as well as high methylation (>50%) was significantly (p = 0.010 and p = 0.045) associated with a worse cancer-specific survival rate in vulvar carcinoma, but not in vaginal carcinoma. Four groups could be defined for the complete series using a Cluster Two step analysis; (1) tumors holding episomal viral DNA, viral load below 150 000 copies not highly methylated (n = 25, 46.3%); (2) tumors harboring episomal viral DNA and being highly methylated (>50%; n = 6, 11.1%); (3) tumors with viral DNA fully integrated (n = 11, 20.4%), and (4) tumors harboring episomal viral DNA and being medium- or unmethylated (<50%) and having a high viral load (> total mean value 150 000; n = 12, 22.2%). The completely integrated tumors were found to be distinct group, whilst some overlap between the groups with high methylation and high viral load was observed.

Conclusion

HPV16- related integration, methylation in E2BS3 and 4 and viral load may represent different viral characteristics driving vaginal and vulvar carcinogenesis. HPV16- related parameters were found to be of clinical importance in the vulvar series only.

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.

HPV-16 E2 physical status and molecular evolution in vivo in cervical carcinomas

A key event in the development of cervical carcinoma is the deregulated expression of high-risk human papillomavirus (HR-HPV) oncogenes, most commonly due to HPV integration into host DNA. Here we explored whether HPV-16 E2 gene integrity is a biomarker of progressive disease with oncogenes expression. HPV-16 genome disruption was assessed by amplification of the entire E2 gene, while mRNA expression patterns of the E1, E2, E6, and E7 genes were evaluated by reverse transcription PCR (RT-PCR). As expected, E2 disruption was significantly higher among patients with cervical cancers than subjects with benign lesions (p=0.02). The status of the E2 gene correlated with tumorogenesis, and seemed also to correlate with the stage of the carcinomas, since integrated HPV-16 DNA was frequently detected in patients with advanced cancer stages (75% of stage III vs 60% stages I and II). In bivariate analysis, the lesions&rsquo; grade was most significantly associated with HPV-16 DNA disruption (p&lt;0.05). In cervical carcinoma the deletion pattern involved more frequently the E2 gene rather than the E1 gene (62.5% vs 45.8%). The prevalence of the E6/E7 HPV-16 transcripts in cervical carcinoma specimens and in benign cervical lesions were detected with frequencies of, respectively, 91.6% and 45.4%. The mRNA levels of the HPV-16 E6/E7 genes were expressed at approximately the same levels in each physical state. We consistently observed that E6/E7 were absent or weakly detectable in the presence of E2. However, in the absence of E2 the levels of E6/E7 markedly increased (p&lt;0.05). This study underscores the significance of investigating alternative mechanisms of E2 expression and oncogenes E6/E7 transcripts in vivo as biomarkers for disease severity in cervical carcinomas.

Differential expression of HPV16 L2 gene in cervical cancers harboring episomal HPV16 genomes: influence of synonymous and non-coding region variations

We tested the hypothesis that (i) synonymous variations within the coding regions, and (ii) variations within the non-coding regions of HPV, influence cervical cancer (CaCx) pathogenesis under the impact of intact HPV16 genomes. Whole genome sequence analysis of HPV16 isolates within 70 CaCx cases and 25 non-malignant samples revealed that synonymous variations were significantly higher within the E6 (p = 0.014), E5 (p = 0.001) and L2 (p = 0.0002) genes of HPV16 isolates within cases, compared to isolates within non-malignant samples. All of the 25 (100%) humanized codons identified within L2 ORF of the samples analyzed, were harbored by CaCx cases, while 8 out of 25 (32%) were harbored by HPV16 positive non-malignant samples (p = 3.87105E-07). L2 (mRNA and protein) expression was evident only among cases with episomal viral genomes and L2 mRNA expression correlated significantly with E2 gene copy numbers suggesting expression from all episomal genomes. Among such cases, Asian American (AA) isolates portrayed all of the humanized codons (100%; 4–6/sample) recorded within L2, which was significantly higher (p = 2.02E-7) compared to the European (E) isolates (22.8%; none or 1–2/sample). Additionally, majority of E variant isolates within cases (54/57; 94.7%) portrayed a variation (T4228C) within the short non-coding region (NCR2) between E5 and L2 genes, which portrays a weak promoter activity specific for L2 mRNA expression. This resulted in loss of 9 out of 14 miRNA binding sites (hsa-miR-548 family), despite the significant overexpression of miR548a-5p and miR548d-5p among such cases (28.64 and 36.25 folds, respectively), in comparison to HPV negative control samples. The findings exemplify the biological relevance of sequence variations in HPV16 genomes and highlight that episomal HPV16 in CaCx cases employ multiple mechanisms to sustain L2 expression, thereby justifying the potential role of L2 in such cancers, as opposed to those harboring viral integration.

Differential methylation of E2 binding sites in episomal and integrated HPV 16 genomes in preinvasive and invasive cervical lesions

Enhanced expression of the HPV 16 E6-E7 oncogenes may trigger neoplastic transformation of the squamous epithelial cells at the uterine cervix. The HPV E2 protein is a key transcriptional regulator of the E6-E7 genes. It binds to four E2 binding sites (E2BSs 1-4) in the viral upstream regulatory region (URR). Modification of E2 functions, for example, by methylation of E2BSs is hypothesized to trigger enhanced expression of the viral E6-E7 oncogenes. In the majority of HPV-transformed premalignant lesions and about half of cervical carcinomas HPV genomes persist in an extra-chromosomal, episomal state, whereas they are integrated into host cells chromosomes in the remaining lesions. Here we compared the methylation profile of E2BSs 1-4 of the HPV 16 URR in a series of 18 HPV16-positive premalignant lesions and 33 invasive cervical cancers. CpGs within the E2BSs 1, 3, and 4 were higher methylated in all lesions with only episomal HPV16 genomes compared with lesions displaying single integrated copies. Samples with multiple HPV16 integrated copies displayed high methylation levels for all CpGs suggesting that the majority of multiple copies were silenced by extensive methylation. These data support the hypothesis that differential methylation of the E2BSs 1, 3 and 4 is related to the activation of viral oncogene expression in cervical lesions as long as the viral genome remains in the episomal state. Once the virus becomes integrated into host cell chromosomes these methylation patterns may be substantially altered due to complex epigenetic changes of integrated HPV genomes.

Whole genome analysis of human papillomavirus type 16 multiple infection in cervical cancer patients

The characterization of the whole genome of human papillomavirus type 16 (HPV16) from cervical cancer specimens with multiple infections in comparison with single infection samples as the oncogenic potential of the virus may differ. Cervical carcinoma specimens positive for HPV16 by PCR and INNO-LiPA were randomly selected for whole genome characterization. Two HPV16 single infection and six HPV16 multiple infection specimens were subjected to whole genome analysis by using conserved primers and subsequent sequencing. All HPV16 whole genomes from single infection samples clustered in the European (E) lineage while all multiple infection specimens belonged to the non-European lineage. The variations in nucleotide sequences in E6, E7, E2, L1 and Long control region (LCR) were evaluated. In the E6 region, amino acid changes at L83V were related to increased cancer progression. An amino acid variation N29S within the E7 oncoprotein significantly associated with severity of lesion was also discovered. In all three domains of the E2 gene non synonymous mutations were found. The L1 region showed various mutations which may be related to conformation changes of viral epitopes. Some transcription factor binding sites in the LCR region correlated to virulence were shown on GRE/1, TEF- 1, YY14 and Oct-1. HPV16 European variant prone to single infection may harbor a major variation at L83V which significantly increases the risk for developing cervical carcinoma. HPV16 non-European variants prone to multiple infections may require many polymorphisms to enhance the risk of cervical cancer development.

Some novel insights on HPV16 related cervical cancer pathogenesis based on analyses of LCR methylation, viral load, E7 and E2/E4 expressions

This study was undertaken to decipher the interdependent roles of (i) methylation within E2 binding site I and II (E2BS-I/II) and replication origin (nt 7862) in the long control region (LCR), (ii) expression of viral oncogene E7, (iii) expression of the transcript (E7-E1/E4) that encodes E2 repressor protein and (iv) viral load, in human papillomavirus 16 (HPV16) related cervical cancer (CaCx) pathogenesis. The results revealed over-representation (p<0.001) of methylation at nucleotide 58 of E2BS-I among E2-intact CaCx cases compared to E2-disrupted cases. Bisulphite sequencing of LCR revealed overrepresentation of methylation at nucleotide 58 or other CpGs in E2BS-I/II, among E2-intact cases than E2-disrupted cases and lack of methylation at replication origin in case of both. The viral transcript (E7-E1/E4) that produces the repressor E2 was analyzed by APOT (amplification of papillomavirus oncogenic transcript)-coupled-quantitative-RT-PCR (of E7 and E4 genes) to distinguish episomal (pure or concomitant with integrated) from purely integrated viral genomes based on the ratio, E7 C(T)/E4 C(T). Relative quantification based on comparative C(T) (threshold cycle) method revealed 75.087 folds higher E7 mRNA expression in episomal cases over purely integrated cases. Viral load and E2 gene copy numbers were negatively correlated with E7 C(T) (p = 0.007) and E2 C(T) (p<0.0001), respectively, each normalized with ACTB C(T), among episomal cases only. The k-means clustering analysis considering E7 C(T) from APOT-coupled-quantitative-RT-PCR assay, in conjunction with viral load, revealed immense heterogeneity among the HPV16 positive CaCx cases portraying integrated viral genomes. The findings provide novel insights into HPV16 related CaCx pathogenesis and highlight that CaCx cases that harbour episomal HPV16 genomes with intact E2 are likely to be distinct biologically, from the purely integrated viral genomes in terms of host genes and/or pathways involved in cervical carcinogenesis.

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

BACKGROUND

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

METHODOLOGY/PRINCIPAL FINDINGS

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

CONCLUSIONS/SIGNIFICANCE

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

Correlation between ebv co-infection and HPV16 genome integrity in Tunisian cervical cancer patients

Infection with high risk Human papillomavirus (HR-HPV) is necessary but not sufficient to cause cervical carcinoma. This study explored whether multiple HR-HPV or coinfection with Epstein-Barr virus (EBV) influence the integration status of HPV16 genome. The presence and typing of HPV in a series of 125 cervical specimens were assessed by polymerase chain reaction (PCR) using the specific primers for the HPV L1 region. As for EBV infection, the viral EBNA1 gene was used for its detection through PCR amplification. Disruption of the HPV E2 gene was assessed by amplification of the entire E2 gene with single set of primers, while E2 transcripts were evaluated by a reverse transcription PCR method (RT-PCR). The overall prevalence of HPVDNA was of 81.8% in cervical cancers versus 26.9% in benign lesions. In HPV positive cases, HPV16 and HPV18 were the most prevalent types, followed by HPV types 33, 31. EBV EBNA1 prevalence was statistically more frequent in cervical carcinomas than in benign lesions (29.5%, vs 9.6%; P=0.01). No viral infection was detected in healthy control women. The uninterrupted E2 gene was correlated with the presence of E2 transcripts originating from the HPV episomal forms. It was observed that integration was more common in HPV18 and EBV coinfection. The presence of EBV caused a five-fold [OR= 5; CI= 1.15-21.8; P = 0.04] increase in the risk of HPV16 genome integration in the host genome. This study indicates that EBV infection is acting as a cofactor for induction of cervical cancer by favoring HPVDNA integration.

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.

Association of viral load with HPV16 positive cervical cancer pathogenesis: causal relevance in isolates harboring intact viral E2 gene

We tested the hypothesis that cervical cancers (CaCx) harbor high HPV16 viral load compared to controls and this is influenced by E2 status and age of subjects. Viral load (natural log transformed values) per 100ng genomic DNA was estimated (152 cases and 87 controls) by Taqman assay. Median viral load was significantly higher (Mann-Whitney U test) among cases (17.21) compared to controls (9.86), irrespective of E2 status or upon considering E2 status as a covariate in logistic regression model (p<0.001). Viral load of E2 intact cases (17.80) was significantly higher (p<0.001) compared to those with disrupted E2 (9.78). At equivalent probability of being a case, viral load was higher among individuals (i) of lower age, irrespective of E2 status, and (ii) with intact E2 but of similar age as those with disrupted E2. Thus viral load in association with E2 status and/or age might be of causal relevance in CaCx pathogenesis.

Human papilloma virus genotyping, variants and viral load in tumors, squamous intraepithelial lesions, and controls in a north Indian population subset

A study of human papilloma virus (HPV) types and variants is important for developing preventive protocols and appropriate intervention targets. The presence of HPV types, their variants, and viral load in a population subset from North India was studied. Polymerase chain reaction (PCR) and line blots were used for HPV genotyping; HPV 16 and 18 viral loads were measured using real-time PCR. Variant analysis was done by sequencing of the PCR-amplified E6/E7 regions of HPV 16 and the long control region and E6/E7 regions of HPV 18. The 93.6%, 78.6%, and 10% of tumors, squamous intraepithelial lesions (SILs), and controls were HPV-positive, respectively. The most commonly observed type was HPV 16. Human papilloma virus 73 which is uncommonly observed was seen in 2 tumors. Multiple infections were more common in controls and SILs than tumors. The majority (86.4%) of the HPV 16-positive and all of the HPV 18-positive samples belonged to the European variant class. Five novel nonsynonymous changes were seen in the HPV 16-positive and 2 in HPV 18-positive samples. There was a significant increase in viral loads from controls through SILs to tumors, but no significant differences in viral loads were observed between different stages of cancer. In tumors, a significant increase in HPV 16 viral loads was seen with increasing age. The study shows a similar HPV type and variant distribution to European studies, with some differences in type distribution. Viral load does not appear to be good marker for stage wise progression and intralesional variability may affect its use as a differentiating parameter between high-grade squamous intraepithelial lesion and low-grade squamous intraepithelial lesions.

Addition of a single E2 binding site to the human papillomavirus (HPV) type 16 long control region enhances killing of HPV positive cells via HPV E2 protein-regulated herpes simplex virus type 1 thymidine kinase-mediated suicide gene therapy

Human papillomavirus type 16 (HPV16) is associated with the development of anogenital cancers and their precursor lesions, intraepithelial neoplasia. Treatment strategies against HPV-induced intraepithelial neoplasia are not HPV specific and mostly consist of physical removal or ablation of lesions. We had previously designed an HPV-specific approach to kill HPV-infected cells by the herpes simplex virus type 1 thymidine kinase (TK) gene driven by HPV E2 binding to E2-binding sites (E2BS) in the native HPV16 long control region. E2-induced TK expression renders the cells sensitive to the prodrug ganciclovir. To optimize this therapeutic approach, we modified the native long control region by adding variable numbers of E2BS adjacent to E2BS4, resulting in greatly increased cell death in HPV-positive cell lines with variable levels of E2 protein expression and no reduction in HPV specificity. Our results showed maximum increase in TK expression and cell killing when one additional E2BS was added adjacent to E2BS. As HPV-infected patients also exhibit variable E2 expression across lesions and within a lesion, this approach may potentiate the clinical utility of the herpes simplex virus type 1 TK/ganciclovir therapeutic approach.

Cell-type specific transcriptional activities among different papillomavirus long control regions and their regulation by E2

This study systematically examined the viral long control region (LCR) activities and their responses to E2 for human papillomavirus (HPV) types 11, 16, and 18 as well as bovine papillomavirus 1 (BPV1) in a number of different cell types, including human cervical cancer cell lines, human oral keratinocytes, BJ fibroblasts, as well as CV1 cells. The study revealed cell- and virus-type specific differences among the individual LCRs and their regulation by E2. In addition, the integration of the LCR into the host genome was identified as a critical determinant for LCR activity and its response to E2. Collectively, these data indicate a more complex level of transcriptional regulation of the LCR by cellular and viral factors than previously appreciated, including a comparatively low LCR activity and poor E2 responsiveness for HPV16 in most human cells. This study should provide a valuable framework for future transcriptional studies in the papillomavirus field.

Characterization of sequence variations within HPV16 isolates among Indian women: prediction of causal role of rare non-synonymous variations within intact isolates in cervical cancer pathogenesis

We re-sequenced HPV16 genome (~6 kb) implicated in cervical carcinogenesis (LCR, E2, E5, E6, E7, L1, L2) to prioritize sequence variants for functional validation as biomarkers, using CaCx cases (n=74) and asymptomatic controls (n=24). Of the nucleotide variations recorded (n=271), non-synonymous changes in L2 region were significantly higher (p=0.005) among cases (2.67%) compared to controls (1.27%). Using SIFT database, 29 non-synonymous changes (frequency=0.01-0.03) predicted as deleterious to protein functions were identified. Haplotype analysis considering 110 polymorphic variations (frequency> or =0.05) within intact viral isolates (53 CaCx cases and 21 controls) using NETWORK software, confirmed Asian-American (AA, 14.86%) and European (E, 85.14%) variants, differing at 78 positions. The E-variants portrayed thirty-six haplotypes, of which, E-12 was most prevalent within cases (38.1%; 16/42) and controls (28.57%; 6/21) harboring polymorphic variations at 10 positions, in contrast to HPV16R. Cases of the E-12 haplotype harbored 7 deleterious mutations distributed within L1 (n=1), E2 (n=1), E5 (n=1), and L2 (n=4), while none within similar controls. Thus rare deleterious variations within genes implicated in productive infection over the E-12 haplotype background of intact HPV16 isolates might be of causal relevance for CaCx development.

The natural history of cervical HPV infection: unresolved issues

The identification of high-risk human papillomavirus (HPV) types as a necessary cause of cervical cancer offers the prospect of effective primary prevention and the possibility of improving the efficiency of cervical screening programmes. However, for these opportunities to be realized, a more complete understanding of the natural history of HPV infection, and its relationship to the development of epithelial abnormalities of the cervix, is required. We discuss areas of uncertainty, and their possible effect on disease prevention strategies.

CpG methylation of HPV 16 LCR at E2 binding site proximal to P97 is associated with cervical cancer in presence of intact E2

Human papillomavirus type 16 (HPV-16) E2 protein negatively regulates transcription of the E6 and E7 genes. This study was done to test the hypothesis that methylation of the HPV 16 long control region (LCR) is overrepresented among cervical cancer (CaCx) cases compared to cytologically normal controls harboring intact E2 gene. Methylation of the E2 binding site (E2BS-I), proximal to the P97 promoter, was assessed by HpaII/ MspI restriction digestion while McrBC digestion was used to assess LCR-E6 (7289-540) for 57 CaCx samples and 15 normal controls. E2BS-I methylation was found to be significantly higher (56.14%) in cases compared to (20%) controls [OR(age-adjusted) (95% CI): 4.53 (1.05-19.43) p=0.042]. The difference between cases (54.39%) and controls (40%) with respect to LCR-E6 methylation status [OR(age-adjusted) (95% CI): 1.77(0.5-6.3); p=0.38] was not significant. Sequencing of a randomly selected set of 13 methylated malignant samples revealed absence or rare presence, of methylation at CpGs 7579, 7535, 7683 and 7862 respectively. Methylation was found to be more at CpGs within E2 binding sites proximal to the P97 promoter. These results indicate the involvement of E2 binding site methylation in presence of intact E2, leading to loss of E2 repressor activity in CaCx.