ABSTRACT
Human papillomavirus (HPV) drives cervical cancer (CaCx) pathogenesis and viral oncoproteins jeopardize global gene expression in such cancers. In this study, our aim was to identify differentially expressed coding (DEcGs) and long noncoding RNA genes (DElncGs) specifically sense intronic and Natural Antisense Transcripts as they are located in the genic regions and may have a direct influence on the expression pattern of their neighbouring coding genes. We compared HPV16-positive CaCx patients (N = 44) with HPV-negative normal individuals (N = 34) by employing strand-specific RNA-seq and determined the relationships between DEcGs and DElncGs and their clinical implications. By performing Gene set enrichment and protein-protein interaction (PPI) analyses of DEcGs, we identified enrichment of processes crucial for abortive virus life cycle and cancer progression. The DEcGs formed 16 gene clusters which we identified through Molecular Complex Detection (MCODE) plugin of Cytoscape. All the gene clusters portrayed cancer-related functions. We recorded significantly correlated expression levels of 79 DElncGs with DEcGs at proximal genomic loci based on Pearson's Correlation coefficients. Of these gene pairs, 24 pairs portrayed significantly altered correlation coefficients among patients, compared to normal individuals. Of these, 6 DEcGs of 6 such gene pairs, belonged to 5 of the identified gene clusters, one of which was survival-associated. Out of the 24 correlated DEcG: DElncG pairs, we identified 3 pairs, where expression of both members was significantly associated with patient overall survival. The findings justify the cooperative roles of these gene pairs, in patient prognostication, thereby bearing immense potential for translation. Thus, elucidation of correlative strengths between paired DElncGs and DEcGs in patient and normal samples, could serve as a foundation for identification of therapeutic and prognostic targets of HPV16-positive CaCx.
Subject(s)
Gene Expression Regulation, Neoplastic , Human papillomavirus 16 , Papillomavirus Infections , RNA, Long Noncoding , Uterine Cervical Neoplasms , Humans , Uterine Cervical Neoplasms/genetics , Uterine Cervical Neoplasms/virology , Uterine Cervical Neoplasms/pathology , Female , RNA, Long Noncoding/genetics , Human papillomavirus 16/genetics , Human papillomavirus 16/pathogenicity , Gene Expression Regulation, Neoplastic/genetics , Papillomavirus Infections/genetics , Papillomavirus Infections/virology , Middle Aged , Multigene Family/genetics , Adult , Clinical RelevanceABSTRACT
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.
Subject(s)
Gene Expression Regulation, Neoplastic , Gene Regulatory Networks , Human papillomavirus 16 , Papillomavirus Infections , RNA, Long Noncoding , Uterine Cervical Neoplasms , Humans , RNA, Long Noncoding/genetics , Female , Uterine Cervical Neoplasms/virology , Uterine Cervical Neoplasms/genetics , Uterine Cervical Neoplasms/pathology , Human papillomavirus 16/genetics , Human papillomavirus 16/pathogenicity , Papillomavirus Infections/virology , Papillomavirus Infections/genetics , Papillomavirus Infections/complications , Papillomavirus Infections/pathology , Gene Expression Profiling , Prognosis , Biomarkers, Tumor/geneticsABSTRACT
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.
Subject(s)
Uterine Cervical Neoplasms , Female , Humans , Uterine Cervical Neoplasms/genetics , Uterine Cervical Neoplasms/therapy , Human papillomavirus 16/genetics , Immunotherapy , Phenotype , CD8-Positive T-LymphocytesABSTRACT
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.
Subject(s)
Myelin and Lymphocyte-Associated Proteolipid Proteins , Oncogene Proteins, Viral , RNA, Long Noncoding , Uterine Cervical Neoplasms , Female , Humans , Chromatin/metabolism , DNA Methylation , Down-Regulation , Enhancer of Zeste Homolog 2 Protein/genetics , Enhancer of Zeste Homolog 2 Protein/metabolism , Histones/metabolism , Human papillomavirus 16/genetics , Oncogene Proteins, Viral/genetics , Oncogene Proteins, Viral/metabolism , Papillomavirus E7 Proteins/genetics , Papillomavirus E7 Proteins/metabolism , RNA, Long Noncoding/genetics , RNA, Long Noncoding/metabolism , Uterine Cervical Neoplasms/pathology , Myelin and Lymphocyte-Associated Proteolipid Proteins/genetics , Myelin and Lymphocyte-Associated Proteolipid Proteins/metabolismABSTRACT
BACKGROUND: Ovarian germ cell tumours constitute a heterogeneous group of neoplasm with malignant potential being seen in 5% of cases. There is limited data on treatment outcomes of patients with malignant ovarian germ cell tumours (MOGCT). Here, we present our hospital audit of patients with MOGCT. MATERIAL AND METHODS: This is a retrospective data review of patients with MOGCT treated between May 2011 and December 2019. Patients were treated with staging laparotomy and adjuvant chemotherapy, wherever applicable. Surveillance was allowed for those at low risk for recurrence. Clinicopathologic features and treatment details were recorded, and survival analysis was performed. RESULTS: Sixty-five patients with a median age of 25 years (range: 11-52 years) were treated during the study period. The most common histology was immature teratoma in 35.3% of cases. International Federation of Gynecology and Obstetrics stage IC was the most common stage of presentation (47%). Surveillance was advised for 12.3% of cases. Systemic therapy was given in 51 (78%) patients. At a median follow-up of 46 months (range: 1-109 months), the median progression-free survival (PFS) was not reached. Five-year PFS was 79.3% (95% CI: 65.8-88). The most common toxicity was febrile neutropenia (22%) among those who received systemic therapy. CONCLUSION: Immature teratoma was the most common histology in our series. The majority presented in the early stage. MOGCT is a highly curable disease with surgery and systemic therapy.
ABSTRACT
BACKGROUND: Women with response to primary treatment for advanced ovarian cancer are said to have progression if CA125 increases more than double the upper normal limit (70 IU/L) on follow-up. It was, however, noted that large section of women with CA125 > 35 IU/L had disease on imaging. OBJECTIVE: To compare values of CA125 rise at which radiological recurrence can be detected. METHODS: This is a retrospective observational study where women with advanced epithelial ovarian cancer who underwent interval debulking surgery and completed treatment at Tata Medical Center, Kolkata, India, from 2012 to 2016, and were followed up with Ca125. If CA125 doubled or exceeded 35 IU/L or increased to ≥ 70 IU/L, women were subjected to imaging. RESULTS: Among 142 women who underwent treatment, 64 women with response to primary treatment had recurrence. Recurrence was noted in two (3%) patients with doubling of Ca125 but ≤ 35 IU/, 18 (24%) patients with CA125 > 35 IU/L and 41 (64%) patients when CA125 was ≥ 70 IU/L. Three patients (5%) with normal CA125 had recurrence. Among the recurrence group, 45 women had R0 during surgery of which 27 (60%) had CA125 ≥ 70 IU/L and 14 (31%) had CA125 > 35 IU/L during recurrence. Sensitivity and specificity of value > 35 IU/L were 30.51% and 33.33%, respectively, with accuracy of 32.03%, while sensitivity and specificity at > 70 IU/L were 69.49% and 66.67%, respectively, with accuracy of 67.97%. CONCLUSION: CA125 value of ≥ 70 IU/L is a better predictor of recurrence; however, imaging done when value rises > 35 IU/L would be able to detect significant recurrences early thus allowing early treatment.