In-silico analysis of TCGA data showing multiple POLE-like favourable subgroups overlapping with TP53 mutated endometrial cancer: Implications for clinical practice in low and middle-income countries.

Introduction: The Cancer Genome Atlas cohort of endometrial carcinoma (TCGA-UCEC) includes almost 40% TP53-mutants encompassing missense and truncated variants. TCGA revealed 'POLE', characterized by POLE gene bearing exonuclease domain mutation (EDM), as the prognostically best molecular profile. The worst profile was characterized by TP53-mutated Type 2 cancer requiring adjuvant therapy having cost implications in low-resource settings. We aimed to find more 'POLE-like' favourable subgroups by searching TCGA cohort, especially within TP53 mutated risk group, that could eventually avoid adjuvant treatment in resource-poor settings. Method: Our study was an in-silico survival analysis performed on the TCGA-UCEC dataset using SPSS statistical package. TP53 and POLE mutations, microsatellite instability (MSI), time-to-event and clinicopathological parameters were compared among 512 endometrial cancer cases. Deleterious POLE-mutations were identified by Polyphen2. Progression free survival was studied using Kaplan-Meier plots keeping original 'POLE' as comparator. Result: In presence of wild type (WT)-TP53, other deleterious POLE-mutations behaved like POLE-EDM. Only truncated and not missense TP53 benefitted from POLE/MSI overlap. However, TP53 missense mutation, Y220C, was found to be as favourable as 'POLE'. Overlapping POLE, MSI and WT-TP53 also performed favourably. Truncated TP53 overlapped with POLE and/or MSI, TP53 Y220C alone and, WT-TP53 overlapped with POLE and MSI both, were named 'POLE-like' for prognostically behaving like the comparator 'POLE'. Conclusion: Obesity being a lesser frequent event in low and middle-income countries (LMICs), relative proportion of women with lower BMI and Type 2 endometrial cancers may be high. Identification of 'POLE-like' groups may facilitate therapeutic de-escalation in some TP53-mutated cases - a novel option. Instead of 5% (POLE-EDM), potential beneficiary would then comprise 10% (POLE-like) of TCGA-UCEC.

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.

Impact of genetic variations and transcriptional alterations of HLA class I genes on cervical cancer pathogenesis.

In a novel attempt to understand the variations in DNA sequences underlying HLA class I alleles associated with HPV16-related CaCx, we determined the alleles by reconstructing SNP-based haplotypes from resequencing of the most polymorphic exons 2 and 3 of HLA-A, HLA-B and HLA-C. We also determined the impact of SNPs and transcriptional alterations of the genes on CaCx. A high density of SNPs was identified from resequencing. HLA expression was determined by real-time PCR. We identified that even a single associated HLA allele had many underlying SNP-based haplotypes. Out of the most frequent (≥5%) HLA class I alleles, HLA-B*40:06 and HLA-B*15:02 respectively imparted significant risk towards and protection from CaCx as well as HPV16 infection. Employing median-joining networks to detect clusters of sequence-variations for specific HLA alleles, we found the protective SNP-based signature, GAATTTA, in all SNP-based haplotypes of HLA-B*15:02 allele. The signature was derived from seven SNPs within HLA-B which were newly associated with the disease. Contrarily, similarly derived risk-signature, TTGCGCC, mapped only to 52% of SNP-based haplotypes of HLA-B*40:06 allele. This indicated that all SNP-based haplotypes underlying a particular associated HLA allele might or might not have a single signature of risk/protection. HLA-A, HLA-B and HLA-C expressions were downregulated among CaCx cases compared to asymptomatic infections and HPV-negative controls. HLA-A and HLA-B were repressed in both cases harbouring episomal and integrated HPV16, whereas HLA-C in only the latter. Novel genetic variations and differential downregulation-patterns of HLA class I have a significant bearing on HPV16-related CaCx pathogenesis.

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.

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.

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.