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.

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.

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.

Expression signatures of HOX cluster genes in cervical cancer pathogenesis: Impact of human papillomavirus type 16 oncoprotein E7.

The Homeobox (HOX) genes encode important transcription factors showing deregulated expression in several cancers. However, their role in cervical cancer pathogenesis, remains largely unexplored. Herein, we studied their association with Human Papillomavirus type 16 (HPV16) mediated cervical cancers. Our previously published gene expression microarray data revealed a significant alteration of 12 out of 39 HOX cluster members among cervical cancer cases, in comparison to the histopathologically normal controls. Of these, we validated seven (HOXA10, HOXA13, HOXB13, HOXC8, HOXC9, HOXC11 and HOXD10) by quantitative real-time PCR. We identified decreased HOXA10 expression as opposed to the increased expression of the rest. Such decrease was independent of the integration status of HPV16 genome, but correlated negatively with E7 expression in clinical samples, that was confirmed in vitro. HOXA10 and HOXB13 revealed association with Epithelial-Mesenchymal Transition (EMT). While HOXA10 expression correlated positively with E-Cadherin and negatively with Vimentin expression, HOXB13 showed the reverse trend. Chromatin immunoprecipitation study in vitro revealed the ability of E7 to increase HOX gene expression by epigenetic regulation, affecting the H3K4me3 and H3K27me3 status of their promoters, resulting from a loss of PRC2-LSD1 complex activity. Thus, besides identifying the deregulated expression of HOX cluster members in HPV16 positive cervical cancer and their association with EMT, our study highlighted the mechanism of HPV16 E7-mediated epigenetic regulation of HOX genes in such cancers, potentially serving as bedrock for functional studies in the future.

Gestational trophoblastic tumor at a tertiary level cancer center: a retrospective study.

OBJECTIVE: To retrospectively analyze gestational trophoblastic disease at a tertiary level cancer center. STUDY DESIGN: Thirty-six cases of gestational trophoblastic tumor at a tertiary level cancer hospital during the period 1996-2004 were studied. Various epidemiologic features were studied, diagnosed an and patients were classified according to the International Federation of Gynecology and Obstetrics scoring system. Chemotherapy was the main mode of treatment. RESULTS: There was 100% success in low-risk cases and 90% success in high-risk cases. CONCLUSION: Gestational trophoblastic tumors have an excellent prognosis if diagnosed and treated in time, and the potential for childbearing can be maintained.

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.

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 [ORage-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.

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.