Unfurling the functional association between long intergenic noncoding RNAs (lincRNAs) and HPV16-related cervical cancer pathogenesis through weighted gene co-expression network analysis of differentially expressed lincRNAs and coding genes.

Long intergenic noncoding RNAs (lincRNAs) do not overlap annotated coding genes and are located in intergenic regions, as opposed to antisense and sense-intronic lncRNAs, located in genic regions. LincRNAs influence gene expression profiles and are thereby key to disease pathogenesis. In this study, we assessed the association between lincRNAs and HPV16-positive cervical cancer (CaCx) pathogenesis using weighted gene co-expression network analysis (WGCNA) with coding genes, comparing differentially expressed lincRNA and coding genes (DElincGs and DEcGs, respectively) in HPV16-positive patients with CaCx (n = 44) with those in HPV-negative healthy individuals (n = 34). Our analysis revealed five DElincG modules, co-expressing and correlating with DEcGs. We validated a substantial number of such module-specific correlations in the HPV16-positive cancer TCGA-CESC dataset. Four such modules, displayed significant correlations with patient traits, such as HPV16 physical status, lymph node involvement, and overall survival (OS), highlighting a collaborative effect of all genes within specific modules on traits. Using the DAVID bioinformatics knowledgebase, we identified the underlying biological processes associated with these modules as cancer development and progression-associated pathways. Next, we identified the top 10 DElincGs with the highest connectivity within each functional module. Focusing on the prognostic module hub genes, downregulated CTD-2619J13.13 expression was associated with poor patient OS. This lincRNA gene interacted with 25 coding genes of its module and was associated with such biological processes as keratinization loss and keratinocyte differentiation, reflecting severe disease phenotypes. This study has translational relevance in fighting various cancers with high mortality rates in underdeveloped countries.

MAL expression downregulation through suppressive H3K27me3 marks at the promoter in HPV16-related cervical cancers is prognostically relevant and manifested by the interplay of novel MAL antisense long noncoding RNA AC103563.8, E7 oncoprotein and EZH2.

BACKGROUND: MAL (T-lymphocyte maturation-associated protein) is highly downregulated in most cancers, including cervical cancer (CaCx), attributable to promoter hypermethylation. Long noncoding RNA genes (lncGs) play pivotal roles in CaCx pathogenesis, by interacting with human papillomavirus (HPV)-encoded oncoproteins, and epigenetically regulating coding gene expression. Hence, we attempted to decipher the impact and underlying mechanisms of MAL downregulation in HPV16-related CaCx pathogenesis, by interrogating the interactive roles of MAL antisense lncRNA AC103563.8, E7 oncoprotein and PRC2 complex protein, EZH2. RESULTS: Employing strand-specific RNA-sequencing, we confirmed the downregulated expression of MAL in association with poor overall survival of CaCx patients bearing HPV16, along with its antisense long noncoding RNA (lncRNA) AC103563.8. The strength of positive correlation between MAL and AC103563.8 was significantly high among patients compared to normal individuals. While downregulated expression of MAL was significantly associated with poor overall survival of CaCx patients bearing HPV16, AC103563.8 did not reveal any such association. We confirmed the enrichment of chromatin suppressive mark, H3K27me3 at MAL promoter, using ChIP-qPCR in HPV16-positive SiHa cells. Subsequent E7 knockdown in such cells significantly increased MAL expression, concomitant with decreased EZH2 expression and H3K27me3 marks at MAL promoter. In silico analysis revealed that both E7 and EZH2 bear the potential of interacting with AC103563.8, at the same binding domain. RNA immunoprecipitation with anti-EZH2 and anti-E7 antibodies, respectively, and subsequent quantitative PCR analysis in E7-silenced and unperturbed SiHa cells confirmed the interaction of AC103563.8 with EZH2 and E7, respectively. Apparently, AC103563.8 seems to preclude EZH2 and bind with E7, failing to block EZH2 function in patients. Thereby, enhanced EZH2 expression in the presence of E7 could potentially inactivate the MAL promoter through H3K27me3 marks, corroborating our previous results of MAL expression downregulation in patients. CONCLUSION: AC103563.8-E7-EZH2 axis, therefore, appears to crucially regulate the expression of MAL, through chromatin inactivation in HPV16-CaCx pathogenesis, warranting therapeutic strategy development.

Identification of HPV16 positive cervical cancer subsets characterized by divergent immune and oncogenic phenotypes with potential implications for immunotherapy.

BACKGROUND: Cervical cancers (CaCx), like many other cancer types, portray high molecular heterogeneity that affects response to therapy, including immunotherapy. In India and other developing countries, CaCx mortality rates are very high because women report to the clinics with advanced cancers in absence of organized screening programs. This calls for implementation of newer therapeutic regimens for CaCx, like immunotherapy, which is again not used commonly in such countries. OBJECTIVE: Therefore, we focused on dissecting tumour immune heterogeneity, if any, identify immune gene-based biomarkers of heterogeneity and subsets of such cancers with the potential for immunotherapy. We also attempted to characterize the cancer-associated phenotypes of such subsets, including viral load, to decipher the relationship of tumour immunogenicity with oncogenicity. METHODS: Employing RNA-seq analysis of 44 HPV16 positive CaCx patients, immune subtypes were identified by unsupervised hierarchical clustering of global immune-gene expression profiles. Proportions of tumor infiltrating immune cells in the tumor milieu were estimated, employing Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT), using gene expression data from RNA-seq. The oncogenic phenotypes of the immune subtypes of CaCx were deciphered through differential gene expression (DEGs) and pathway enrichment analysis. Viral load was estimated through TaqMan-based qRT-PCR analysis. RESULTS: Analysis revealed the presence of two immune subtypes of CaCx, A (26/44; 59.09%) and B (18/44; 40.90%). Compared to Subtype-A, Subtype-B portrayed overexpression of immune genes and high infiltration of immune cells, specifically CD8+ T cells (p < 0.0001). Besides, a significant correlation between PD-1 and PD-L1 co-expression among Subtype-B, as opposed to Subtype-A, confirmed the interactive roles of these immune checkpoint molecules in Subtype B. Stepwise discriminant analysis pin-pointed ten immune-genes that could classify 100% of the patients significantly (p < 0.0001) into the two immune subtypes and serve as potential biomarkers of CaCx immunity. Differential gene expression analysis between the subtypes unveiled that Subtype-B was more biologically aggressive than Subtype-A, reflecting loss of structural integrity and promotion of cancer progression. The viral load was significantly lower in Subtype-B (average viral load = 10.74/100 ng of genomic DNA) compared to Subtype-A (average viral load = 14.29/100 ng of genomic DNA). Thus viral load and the ten-gene panel underscore their association with immunogenicity and oncogenicity. CONCLUSION: Our study provides strong evidence that only a subset, about 41% of HPV16 positive CaCx patients in India, portray immune enrichment of the tumor milieu coupled with aggressive phenotypes. Such subtypes are therefore likely to benefit through checkpoint molecule-based or tumor infiltrating lymphocyte-based immunotherapy, which could be a leap forward in tackling aggressive forms of such CaCx in India and other developing countries.

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.

Prevalence and age-wise distribution of Human Papillomavirus type 16/18 infections among hospital screened women of a peri-urban area in West Bengal: Impact of socio-demographic factors.

BACKGROUND: We undertook the current study on cervical Human Papillomavirus (HPV) prevalence along with cytology in women visiting the Gynecology Out-patient Department of a hospital for common gynecological ailments, subsequent to our earlier population-based study on HPV prevalence from India. METHODS: We analyzed data on cervical-cytology (Pap smears) and PCR-based molecular detection of HPV infection along with socio-demographic variables (N = 696). RESULTS: We identified 36.84% HPV-positive women amongst whom, HPV16 and 18 together predominated (79.37%) over other HPV types (20.63%). Contrarily, only 6.4% women revealed abnormal cytological lesions, of which, 46.51% were HPV-positive and 95% of such women harbored HPV16/18, while 5% harbored other HPV types. Individuals with normal cytology portrayed 36.09% HPV infections, of which, 77.97% were HPV16/18-positive and 22.03% harbored other HPV types. Overall HPV prevalence decreased significantly (ptrend  = 0.047) with increase in age, but HPV16/18 infections were significantly over-represented compared to the other HPV types across all age-groups. Specifically, HPV16 prevalence increased (p trend < 0.01) with increase in severity of cervical lesions. HPV16 prevalence did not differ between the Hindus and Muslims but HPV18 was significantly higher among the cytologically normal Muslim women (24.14%, p = 0.02), compared to the Hindus (11.91%), specifically among those ≥ 30 years of age. There was a significant (p < 0.05) overrepresentation of HPV16 prevalence among women who were users of oral contraceptive-pills, irrespective of cytology. CONCLUSIONS: Our study highlights the need for HPV16/18-based screening of cervical cancers in India considering the immense socio-cultural and genetic diversity at the population level.

Clinical manifestations of nocardiosis: Study of risk factors and outcomes in a tertiary care hospital.

PURPOSE: The aim of this study is to evaluate the predisposing risk factors, clinical presentations, laboratory parameters, and treatments taken and outcomes in patients of nocardiosis in the span of 5 years in a tertiary care hospital. MATERIALS AND METHODS: The patients whose specimens showed Nocardia like organism in Gram-staining, Kinyoun staining and characteristic colonies in culture were included in the retrospective analysis study. Retrospective analysis of associated risk factors, clinical presentations, and radiological findings was performed. RESULTS: Of the thirteen patients, 11 (76.9%) had immunosuppressive pathologies including solid organ transplantation, autoimmune disease, use of steroids, and immunosuppressive drugs as important risk factors. Four types of clinical manifestations were observed, pulmonary (46.1%), cutaneous (23.07%), cerebral (15.3%), and bacteremia (15.3%). The most common presentation was pulmonary with steroid therapy as a significant risk factor. Consolidation and pleural effusion were the common radiological findings in these cases. In eight of the nine patients anti-nocrdial drugs were given. Cotrimoxazole as monotherapy was given in four cases (44.44%), cotrimoxazole in combination with meropenem in two cases (22.22%); minocycline and linezolid were given in one case each. The overall mortality was 36.36% and was seen in patients with pulmonary nocardiosis. CONCLUSIONS: The study indicates that Nocardial infections are re-emerging on account of an increase in numbers of immunocompromised patients due to increased organ transplants, autoimmune diseases, malignancies, and use of immunosuppressive drugs and steroids. The diagnosis is often missed/not suspected and delayed because of the clinical resemblance to many other infections. Nocardial infection should be suspected and assessed particularly in immunocompromised patients not responding to treatment/improving clinically.

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.

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.

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.

Polymerase chain reaction and deoxyribonucleic acid-sequencing based study on distribution of human papillomavirus 16/18 among histopathological types of cervical intra-epithelial neoplasia and primary invasive cervical carcinoma: A scenario in North Bengal, India.

INTRODUCTION: Human papillomavirus (HPV) types 16/18 are reportedly most common in cervical cancer (CaCx) with geographical variation of genotypes. HPV16 predominates both in squamous cell carcinoma (SCC) and adenocarcinoma in India, contrary to reported global predominance of HPV18 in the latter. Our study was aimed to determine the occurrence of HPV16/18 among histopathological types of cervical intra-epithelial neoplasia (CIN) and invasive CaCx from North Bengal, India and to identify any major deviation from the known Indian scenario of distribution of HPV16/18 genotypes in cases of SCC and adenocarcinoma. MATERIALS AND METHODS: This was a retrospective, cross-sectional, case-only type of study, in which 40 cases were histopathologically diagnosed as CIN/CaCx, on which polymerase chain reaction (PCR), deoxyribonucleic acid (DNA)-sequencing and bioinformatics by basic search local alignment tool were performed for HPV-genotyping. STATISTICAL ANALYSIS: The distribution of HPV genotypes among cases of SCC and adenocarcinoma was compared by Fisher's exact-test. RESULTS: HPV was detected in 97.5% (39/40) cases. HPV16-infected cases (32/39; 82.05%) predominated over HPV18-infected ones (7/39; 17.95%). However, HPV18-only infection was significantly (P = 0.0045, one-sided Fisher's exact test) more among adenocarcinoma (3/4; 75%) than SCC (2/26; 7.69%) contrary to HPV16-only infection (SCC = 24/26, 92.31%; adenocarcinoma = 1/4; 25%) whereas both CIN3 cases were HPV16-positive. CONCLUSION: Predominance of HPV18 over HPV16 in cases of adenocarcinoma in this region was contrasting to that of earlier Indian studies suggesting research on HPV18 related cervical carcinogenesis. PCR and DNA-sequencing could prove to be highly effective tools in HPV detection and genotyping. The study reported HPV16/18 infection in almost 98% of the cases, the knowledge about which might prove useful in future population based studies on HPV genotyping and designing of appropriate HPV-vaccines for this region.

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.

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.

Oncogenic HPV among HIV infected female population in West Bengal, India.

BACKGROUND: Prevalence of both cervical cancer and Human Immunodeficiency Virus (HIV) infection are very high in India. Natural history of Human Papilloma Virus (HPV) infection is known to be altered in HIV positive women and there is an increased possibility of persistence of HPV infections in this population. Therefore, this study was conducted to understand the epidemiology and circulating genotypes of oncogenic HPV among HIV positive and negative female population in West Bengal, India. METHODS: In this hospital-based cross-sectional study, 93 known HIV positive females attending a pre-ART registration clinic and 1106 HIV negative females attending a Reproductive and Child Health Care Clinic were subjected to study. Cervical cell samples collected from the study population were tested for the presence of HPV 16, 18 using specific primers. Roche PCR assay was used to detect other specific HPV genotypes in the cervical cells specimens of HIV positive cases only. RESULTS: Prevalence of HPV 16, 18 among HIV positive females (32.2%; n = 30) was higher than HIV negative females (9.1%; n = 101). About 53% (23/43) of cases with oncogenic HPV were infected with genotypes other than 16, 18 either as single/multiple infections. HPV 18 and HPV 16 were the predominant genotypes among HIV positive and HIV negative subjects respectively. Oncogenic HPV was not found to be associated with age and duration of sexual exposure. But the presence of HIV was found to a statistically significant predictor oncogenic HPV. CONCLUSION: The currently available HPV vaccines offer protection only against HPV 16 and 18 and some cross- protection to few associated genotypes. These vaccines are therefore less likely to offer protection against cervical cancer in HIV positive women a high percentage of who were infected with non-16 and non-18 oncogenic HPV genotypes. Additionally, there is a lack of sufficient evidence of immunogenicity in HIV infected individuals. Therefore, prevention of cervical cancer in HIV positive women must be focused towards early detection of oncogenic HPV & cervical cytological abnormality followed by an appropriate treatment.

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.

TP53 codon 72 polymorphism and cervical cancer: a pooled analysis of individual data from 49 studies.

BACKGROUND: Cervical cancer is caused primarily by human papillomaviruses (HPV). The polymorphism rs1042522 at codon 72 of the TP53 tumour-suppressor gene has been investigated as a genetic cofactor. More than 80 studies were done between 1998 and 2006, after it was initially reported that women who are homozygous for the arginine allele had a risk for cervical cancer seven times higher than women who were heterozygous for the allele. However, results have been inconsistent. Here we analyse pooled data from 49 studies to determine whether there is an association between TP53 codon 72 polymorphism and cervical cancer. METHODS: Individual data on 7946 cases and 7888 controls from 49 different studies worldwide were reanalysed. Odds ratios (OR) were estimated using logistic regression, stratifying by study and ethnic origin. Subgroup analyses were done for infection with HPV, ethnic origin, Hardy-Weinberg equilibrium, study quality, and the material used to determine TP53 genotype. FINDINGS: The pooled estimates (OR) for invasive cervical cancer were 1.22 (95% CI 1.08-1.39) for arginine homozygotes compared with heterozygotes, and 1.13 (0.94-1.35) for arginine homozygotes versus proline homozygotes. Subgroup analyses showed significant excess risks only in studies where controls were not in Hardy-Weinberg equilibrium (1.71 [1.21-2.42] for arginine homozygotes compared with heterozygotes), in non-epidemiological studies (1.35 [1.15-1.58] for arginine homozygotes compared with heterozygotes), and in studies where TP53 genotype was determined from tumour tissue (1.39 [1.13-1.73] for arginine homozygotes compared with heterozygotes). Null results were noted in studies with sound epidemiological design and conduct (1.06 [0.87-1.29] for arginine homozygotes compared with heterozygotes), and studies in which TP53 genotype was determined from white blood cells (1.06 [0.87-1.29] for arginine homozygotes compared with heterozygotes). INTERPRETATION: Subgroup analyses indicated that excess risks were most likely not due to clinical or biological factors, but to errors in study methods. No association was found between cervical cancer and TP53 codon 72 polymorphism when the analysis was restricted to methodologically sound studies. FUNDING: German Research Foundation (DFG).

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.