Uracil as a biomarker for spatial pyrimidine metabolism in the development of gingivobuccal oral squamous cell carcinoma.

No biomarker has yet been identified that allows accurate diagnosis and prognosis of oral cancers. In this study, we investigated the presence of key metabolites in oral cancer using proton nuclear magnetic resonance (NMR) spectroscopy to identify metabolic biomarkers of gingivobuccal oral squamous cell carcinoma (GB-OSCC). NMR spectroscopy revealed that uracil was expressed in 83.09% of tumor tissues and pyrimidine metabolism was active in GB-OSCC; these results correlated well with immunohistochemistry (IHC) and RNA sequencing data. Based on further gene and protein analyses, we proposed a pathway for the production of uracil in GB-OSCC tissues. Uridinetriphosphate (UTP) is hydrolyzed to uridine diphosphate (UDP) by CD39 in the tumor microenvironment (TME). We hypothesized that UDP enters the cell with the help of the UDP-specific P2Y6 receptor for further processing by ENTPD4/5 to produce uracil. As the ATP reserves diminish, the weakened immune cells in the TME utilize pyrimidine metabolism as fuel for antitumor activity, and the same mechanism is hijacked by the tumor cells to promote their survival. Correspondingly, the differential expression of ENTPD4 and ENTPD5 in immune and tumor cells, respectively, indicatedtheir involvement in disease progression. Furthermore, higher uracil levels were detected in patients with lymph node metastasis, indicating that metastatic potential is increased in the presence of uracil. The presence of uracil and/or expression patterns of intermediate molecules in purine and pyrimidine pathways, such asCD39, CD73, and P2Y6 receptors together with ENTPD4 and ENTPD5, hold promise as biomarker(s) for oral cancer diagnosis and prognosis.

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.

Integrative analysis of genomic and transcriptomic data of normal, tumour, and co-occurring leukoplakia tissue triads drawn from patients with gingivobuccal oral cancer identifies signatures of tumour initiation and progression.

A thickened, white patch - leukoplakia - in the oral cavity is usually benign, but sometimes (in ~9% of individuals) it progresses to malignant tumour. Because the genomic basis of this progression is poorly understood, we undertook this study and collected samples of four tissues - leukoplakia, tumour, adjacent normal, and blood - from each of 28 patients suffering from gingivobuccal oral cancer. We performed multiomics analysis of the 112 collected tissues (four tissues per patient from 28 patients) and integrated information on progressive changes in the mutational and transcriptional profiles of each patient to create this genomic narrative. Additionally, we generated and analysed whole-exome sequence data from leukoplakia tissues collected from 11 individuals not suffering from oral cancer. Nonsynonymous somatic mutations in the CASP8 gene were identified as the likely events to initiate malignant transformation, since these were frequently shared between tumour and co-occurring leukoplakia. CASP8 alterations were also shown to enhance expressions of genes that favour lateral spread of mutant cells. During malignant transformation, additional pathogenic mutations are acquired in key genes (TP53, NOTCH1, HRAS) (41% of patients); chromosomal-instability (arm-level deletions of 19p and q, focal-deletion of DNA-repair pathway genes and NOTCH1, amplification of EGFR) (77%), and increased APOBEC-activity (23%) are also observed. These additional alterations were present singly (18% of patients) or in combination (68%). Some of these alterations likely impact immune-dynamics of the evolving transformed tissue; progression to malignancy is associated with immune suppression through infiltration of regulatory T-cells (56%), depletion of cytotoxic T-cells (68%), and antigen-presenting dendritic cells (72%), with a concomitant increase in inflammation (92%). Patients can be grouped into three clusters by the estimated time to development of cancer from precancer by acquiring additional mutations (range: 4-10 years). Our findings provide deep molecular insights into the evolutionary processes and trajectories of oral cancer initiation and progression. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

dbGENVOC: database of GENomic Variants of Oral Cancer, with special reference to India.

Oral cancer is highly prevalent in India and is the most frequent cancer type among Indian males. It is also very common in southeast Asia. India has participated in the International Cancer Genome Consortium (ICGC) and some national initiatives to generate large-scale genomic data on oral cancer patients and analyze to identify associations and systematically catalog the associated variants. We have now created an open, web-accessible database of these variants found significantly associated with Indian oral cancer patients, with a user-friendly interface to enable easy mining. We have value added to this database by including relevant data collated from various sources on other global populations, thereby providing opportunities of comparative geographical and/or ethnic analyses. Currently, no other database of similar nature is available on oral cancer. We have developed Database of GENomic Variants of Oral Cancer, a browsable online database framework for storage, retrieval and analysis of large-scale data on genomic variants and make it freely accessible to the scientific community. Presently, the web-accessible database allows potential users to mine data on ∼24 million clinically relevant somatic and germline variants derived from exomes (n = 100) and whole genomes (n = 5) of Indian oral cancer patients; all generated by us. Variant data from The Cancer Genome Atlas and data manually curated from peer-reviewed publications were also incorporated into the database for comparative analyses. It allows users to query the database by a single gene, multiple genes, multiple variant sites, genomic region, patient ID and pathway identities. Database URL: http://research.nibmg.ac.in/dbcares/dbgenvoc/.

Application of Random Forest and data integration identifies three dysregulated genes and enrichment of Central Carbon Metabolism pathway in Oral Cancer.

BACKGROUND: Studies of epigenomic alterations associated with diseases primarily focus on methylation profiles of promoter regions of genes, but not of other genomic regions. In our past work (Das et al. 2019) on patients suffering from gingivo-buccal oral cancer - the most prevalent form of cancer among males in India - we have also focused on promoter methylation changes and resultant impact on transcription profiles. Here, we have investigated alterations in non-promoter (gene-body) methylation profiles and have carried out an integrative analysis of gene-body methylation and transcriptomic data of oral cancer patients. METHODS: Tumor and adjacent normal tissue samples were collected from 40 patients. Data on methylation in the non-promoter (gene-body) regions of genes and transcriptome profiles were generated and analyzed. Because of high dimensionality and highly correlated nature of these data, we have used Random Forest (RF) and other data-analytical methods. RESULTS: Integrative analysis of non-promoter methylation and transcriptome data revealed significant methylation-driven alterations in some genes that also significantly impact on their transcription levels. These changes result in enrichment of the Central Carbon Metabolism (CCM) pathway, primarily by dysregulation of (a) NTRK3, which plays a dual role as an oncogene and a tumor suppressor; (b) SLC7A5 (LAT1) which is a transporter dedicated to essential amino acids, and is overexpressed in cancer cells to meet the increased demand for nutrients that include glucose and essential amino acids; and, (c) EGFR which has been earlier implicated in progression, recurrence, and stemness of oral cancer, but we provide evidence of epigenetic impact on overexpression of this gene for the first time. CONCLUSIONS: In rapidly dividing cancer cells, metabolic reprogramming from normal cells takes place to enable enhanced proliferation. Here, we have identified that among oral cancer patients, genes in the CCM pathway - that plays a fundamental role in metabolic reprogramming - are significantly dysregulated because of perturbation of methylation in non-promoter regions of the genome. This result compliments our previous result that perturbation of promoter methylation results in significant changes in key genes that regulate the feedback process of DNA methylation for the maintenance of normal cell division.

Study of Caspase 8 mutation in oral cancer and adjacent precancer tissues and implication in progression.

It is hypothesized that same driver gene mutations should be present in both oral leukoplakia and cancer tissues. So, we attempted to find out mutations at one of the driver genes, CASP8, in cancer and adjacent leukoplakia tissues. Patients (n = 27), affected by both of cancer and adjacent leukoplakia, were recruited for the study. Blood and tissue DNA samples were used to identify somatic mutations at CASP8 by next generation sequencing method. In total, 56% (15 out of 27) cancer and 30% (8 out of 27) leukoplakia tissues had CASP8 somatic mutations. In 8 patients, both cancer and adjacent leukoplakia tissues, located within 2-5 cm of tumor sites, had identical somatic mutations. But, in 7 patients, cancer samples had somatic mutations but none of the leukoplakia tissues, located beyond 5cm of tumor sites, had somatic mutations. Mutated allele frequencies at CASP8 were found to be more in cancer compared to adjacent leukoplakia tissues. This study provides mutational evidence that oral cancer might have progressed from previously grown leukoplakia lesion. Leukoplakia tissues, located beyond 5cm of cancer sites, were free from mutation. The study implies that CASP8 mutation could be one of the signatures for some of the leukoplakia to progress to oral cancer.

A novel hotspot and rare somatic mutation p.A138V, at TP53 is associated with poor survival of pancreatic ductal and periampullary adenocarcinoma patients.

BACKGROUND: Pancreatic Ductal Adenocarcinoma (PDAC) is a cancer of the exocrine pancreas and 5-year survival rates remain constant at 7%. Along with PDAC, Periampullary Adenocarcinoma (PAC) accounts for 0.5-2% of all gastrointestinal malignancies. Genomic observations were well concluded for PDAC and PACs in western countries but no reports are available from India till now. METHODS: Targeted Next Generation Sequencing were performed in 8 (5 PDAC and 3 PAC) tumour normal pairs, using a panel of 412 cancer related genes. Primary findings were replicated in 85 tumour samples (31 PDAC and 54 PAC) using the Sanger sequencing. Mutations were also validated by ASPCR, RFLP, and Ion Torrent sequencing. IHC along with molecular dynamics and docking studies were performed for the p.A138V mutant of TP53. Key polymorphisms at TP53 and its associated genes were genotyped by PCR-RFLP method and association with somatic mutations were evaluated. All survival analysis was done using the Kaplan-Meier survival method which revealed that the survival rates varied significantly depending on the somatic mutations the patients harboured. RESULTS: Among the total 114 detected somatic mutations, TP53 was the most frequently mutated (41%) gene, followed by KRAS, SMAD4, CTNNB1, and ERBB3. We identified a novel hotspot TP53 mutation (p.A138V, in 17% of all patients). Low frequency of KRAS mutation (33%) was detected in these samples compared to patients from Western counties. Molecular Dynamics (MD) simulation and DNA-protein docking analysis predicted p.A138V to have oncogenic characteristics. Patients with p.A138V mutation showed poorer overall survival (p = 0.01). So, our finding highlights elevated prevalence of the p53p.A138V somatic mutation in PDAC and pancreatobiliary PAC patients. CONCLUSION: Detection of p.A138V somatic variant in TP53 might serve as a prognostic marker to classify patients. It might also have a role in determining treatment regimes. In addition, low frequency of KRAS hotspot mutation mostly in Indian PDAC patient cohort indicates presence of other early drivers in malignant transformation.

Epigenomic dysregulation-mediated alterations of key biological pathways and tumor immune evasion are hallmarks of gingivo-buccal oral cancer.

BACKGROUND: Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. Although genomic drivers of OSCC-GB have been identified by whole exome and whole genome sequencing, no epigenome-wide study has been conducted in OSCC-GB; our study has filled this gap, and has discovered and validated epigenomic hallmarks of gingivobuccal oral cancer. METHODS: We have carried out integrative analysis of epigenomic (n = 87) and transcriptomic (n = 72) profiles of paired tumor-normal tissues collected from OSCC-GB patients from India. Genome-wide DNA methylation assays and RNA-sequencing were performed on high-throughput platforms (Illumina) using a half-sample of randomly selected patients to discover significantly differentially methylated probes (DMPs), which were validated on the remaining half-sample of patients. RESULTS: About 200 genes showed significant inverse correlation between promoter methylation and expression, of which the most significant genes included genes that act as transcription factors and genes associated with other cancer types. Novel findings of this study include identification of (a) potential immunosuppressive effect in OSCC-GB due to significant promoter hypomethylation driven upregulation of CD274 and CD80, (b) significant dysregulation by epigenetic modification of DNMT3B (upregulation) and TET1 (downregulation); and (c) known drugs that can reverse the direction of dysregulation of gene expression caused by promoter methylation. CONCLUSIONS: In OSCC-GB patients, there are significant alterations in expression of key genes that (a) regulate normal cell division by maintenance of balanced DNA methylation and transcription process, (b) maintain normal physiological signaling (PPAR, B cell receptor) and metabolism (arachidonic acid) pathways, and (c) provide immune protection against antigens, including tumor cells. These findings indicate novel therapeutic targets, including immunotherapeutic, for treatment of OSCC-GB.

Lymph node metastasis in oral cancer is strongly associated with chromosomal instability and DNA repair defects.

Oral squamous cell carcinoma (OSCC) is highly prevalent in south and southeast Asia. Many (30-50%) OSCC patients develop lymph node metastasis (LNM), which is the most important prognostic factor in OSCC. To identify genomic correlates of LNM, we compared exome sequences and copy number variation data of blood and tumor DNA from highly contrasting subgroups of patients to reduce false inferences-(i) patients with LNM and (ii) patients with late stage disease but without LNM. We found that LNM is associated with (i) specific hotspot somatic mutations in TP53 and CASP8; (ii) rare nonsilent germline mutations in BRCA2 and FAT1; (iii) mutations in mito-G2/M and nonhomologous end joining (NHEJ) pathways; (iv) recurrent deletion of genes for DNA repair by homologous recombination; and (v) chromosomal instability. LN+ patients with NHEJ pathway mutations have longer disease-free survival. Five genomic features have a high predictive value of LNM.

Profiling of genomic alterations of mitochondrial DNA in gingivobuccal oral squamous cell carcinoma: Implications for disease progress.

We have identified 164 somatic mutations in mitochondrial DNA in gingivobuccal oral cancer by deep sequencing the mitochondrial genome from paired tumor and blood DNA samples from 89 patients. We have found evidence of positive selection of somatic nonsynonymous mutations. Non-synonymous mutations in mitochondrial respiratory genes were found to increase the risk of lymph node metastasis (P = 0.0028). We have observed a significant reduction in mitochondrial DNA copy number in tumor DNA of these patients compared to the DNA from adjacent normal tissue samples (P < 1 × 10-6). Analysis of transcriptome data of tumor and adjacent normal tissue revealed patients harboring mutations in mitochondrial protein-coding genes exhibited reduced expression of mitochondrial transcripts.

Analysis of the whole transcriptome from gingivo-buccal squamous cell carcinoma reveals deregulated immune landscape and suggests targets for immunotherapy.

BACKGROUND: Gingivo-buccal squamous cell carcinoma (GBSCC) is one of the most common oral cavity cancers in India with less than 50% patients surviving past 5 years. Here, we report a whole transcriptome profile on a batch of GBSCC tumours with diverse tobacco usage habits. The study provides an entire landscape of altered expression with an emphasis on searching for targets with therapeutic potential. METHODS: Whole transcriptomes of 12 GBSCC tumours and adjacent normal tissues were sequenced and analysed to explore differential expression of genes. Expression changes were further compared with those in TCGA head and neck cohort (n = 263) data base and validated in an independent set of 10GBSCC samples. RESULTS: Differentially expressed genes (n = 2176) were used to cluster the patients based on their tobacco habits, resulting in 3 subgroups. Immune response was observed to be significantly aberrant, along with cell adhesion and lipid metabolism processes. Different modes of immune evasion were seen across 12 tumours with up-regulation or consistent expression of CD47, unlike other immune evasion genes such as PDL1, FUT4, CTLA4 and BTLA which were downregulated in a few samples. Variation in infiltrating immune cell signatures across tumours also indicates heterogeneity in immune evasion strategies. A few actionable genes such as ITGA4, TGFB1 and PTGS1/COX1 were over expressed in most samples. CONCLUSION: This study found expression deregulation of key immune evasion genes, such as CD47 and PDL1, and reasserts their potential as effective immunotherapeutic targets for GBSCC, which requires further clinical studies. Present findings reiterate the idea of using transcriptome profiling to guide precision therapeutic strategies.

Hedgehog signaling pathway is active in GBM with GLI1 mRNA expression showing a single continuous distribution rather than discrete high/low clusters.

Hedgehog (Hh) signaling pathway is a valid therapeutic target in a wide range of malignancies. We focus here on glioblastoma multiforme (GBM), a lethal malignancy of the central nervous system (CNS). By analyzing RNA-sequencing based transcriptomics data on 149 clinical cases of TCGA-GBM database we show here a strong correlation (r = 0.7) between GLI1 and PTCH1 mRNA expression--as a hallmark of the canonical Hh-pathway activity in this malignancy. GLI1 mRNA expression varied in 3 orders of magnitude among the GBM patients of the same cohort showing a single continuous distribution-unlike the discrete high/low-GLI1 mRNA expressing clusters of medulloblastoma (MB). When compared with MB as a reference, the median GLI1 mRNA expression in GBM appeared 14.8 fold lower than that of the "high-Hh" cluster of MB but 5.6 fold higher than that of the "low-Hh" cluster of MB. Next, we demonstrated statistically significant up- and down-regulation of GLI1 mRNA expressions in GBM patient-derived low-passage neurospheres in vitro by sonic hedgehog ligand-enriched conditioned media (shh-CM) and by Hh-inhibitor drug vismodegib respectively. We also showed clinically achievable dose (50 µM) of vismodegib alone to be sufficient to induce apoptosis and cell cycle arrest in these low-passage GBM neurospheres in vitro. Vismodegib showed an effect on the neurospheres, both by down-regulating GLI1 mRNA expression and by inducing apoptosis/cell cycle arrest, irrespective of their relative endogenous levels of GLI1 mRNA expression. We conclude from our study that this single continuous distribution pattern of GLI1 mRNA expression technically puts almost all GBM patients in a single group rather than discrete high- or low-clusters in terms of Hh-pathway activity. That is suggestive of therapies with Hh-pathway inhibitor drugs in this malignancy without a need for further stratification of patients on the basis of relative levels of Hh-pathway activity among them.

Variant allele frequency enrichment analysis in vitro reveals sonic hedgehog pathway to impede sustained temozolomide response in GBM.

Neoplastic cells of Glioblastoma multiforme (GBM) may or may not show sustained response to temozolomide (TMZ) chemotherapy. We hypothesize that TMZ chemotherapy response in GBM is predetermined in its neoplastic clones via a specific set of mutations that alter relevant pathways. We describe exome-wide enrichment of variant allele frequencies (VAFs) in neurospheres displaying contrasting phenotypes of sustained versus reversible TMZ-responses in vitro. Enrichment of VAFs was found on genes ST5, RP6KA1 and PRKDC in cells showing sustained TMZ-effect whereas on genes FREM2, AASDH and STK36, in cells showing reversible TMZ-effect. Ingenuity pathway analysis (IPA) revealed that these genes alter cell-cycle, G2/M-checkpoint-regulation and NHEJ pathways in sustained TMZ-effect cells whereas the lysine-II&V/phenylalanine degradation and sonic hedgehog (Hh) pathways in reversible TMZ-effect cells. Next, we validated the likely involvement of the Hh-pathway in TMZ-response on additional GBM neurospheres as well as on GBM patients, by extracting RNA-sequencing-based gene expression data from the TCGA-GBM database. Finally, we demonstrated TMZ-sensitization of a TMZ non-responder neurosphere in vitro by treating them with the FDA-approved pharmacological Hh-pathway inhibitor vismodegib. Altogether, our results indicate that the Hh-pathway impedes sustained TMZ-response in GBM and could be a potential therapeutic target to enhance TMZ-response in this malignancy.

Somatic mutations in arachidonic acid metabolism pathway genes enhance oral cancer post-treatment disease-free survival.

The arachidonic acid metabolism (AAM) pathway promotes tumour progression. Chemical inhibitors of AAM pathway prolong post-treatment survival of cancer patients. Here we test whether non-synonymous somatic mutations in genes of this pathway, acting as natural inhibitors, increase post-treatment survival. We identify loss-of-function somatic mutations in 15 (18%) of 84 treatment-naïve oral cancer patients by whole-exome sequencing, which we map to genes of AAM pathway. Patients (n = 53) who survived ≥ 12 months after surgery without recurrence have significantly (P = 0.007) higher proportion (26% versus 3%) of mutations than those who did not (n = 31). Patients with mutations have a significantly (P = 0.003) longer median disease-free survival (24 months) than those without (13 months). Compared with the presence of a mutation, absence of any mutation increases the hazard ratio for death (11.3) significantly (P = 0.018). The inferences are strengthened when we pool our data with The Cancer Genome Atlas (TCGA) data. In patients with AAM pathway mutations, some downstream pathways, such as the PI3K-Akt pathway, are downregulated.

Exome sequencing reveals the likely involvement of SOX10 in uveal melanoma.

PURPOSE: To identify the spectrum of somatic mutations in an Asian Indian patient with uveal melanoma (UM) without metastasis using exome sequencing. CASE REPORT: A 49-year-old man from India was diagnosed as having cilio-choroidal (uveal) melanoma (UM), without metastasis, in his right eye with the help of magnetic resonance imaging. This was later confirmed by histopathological evaluation. Two individuals from India with non-neoplastic blind eyes were recruited as controls. The affected eyes from the UM patient and the two control individuals were enucleated, and uveal tissues were collected. DNA was extracted from uveal tissue, and the matched blood sample from each of the three individuals was followed by exome sequencing. Statistical and bioinformatic analyses were done to identify somatic mutations and their putative associations with UM. Thirty-one somatic mutations (25 amino acid altering) in protein-coding (exonic) regions were detected in the UM patient. Of the amino acid-altering somatic mutations, 16 mutations were predicted to be candidate mutations relevant to UM. Somatic mutations, putatively causal for UM, were identified in GNAQ, SF3B1, and SOX10. CONCLUSIONS: Somatic mutations in GNAQ and SF3B1 genes were probable drivers of UM in the Indian patient; these were also reported earlier in some White patients. In addition, a frameshift deletion of 20 base pairs has been identified in SOX10 in the UM patient. Somatic mutations in SOX10, a transcription factor, which acts upstream of microphthalmia-associated transcription factor and synergizes with microphthalmia-associated transcription factor, was identified in some melanoma cell lines. The transcription factor SOX10 was found to have an essential role in melanocyte development and pigmentation. Our finding of the frameshift deletion (p.H387fs) in exon 4 of SOX10 in UM provides an important insight and complements earlier findings of mutations in GNAQ and SF3B1 on the genomic basis of UM.

Increased risk of oral cancer in relation to common Indian mitochondrial polymorphisms and Autosomal GSTP1 locus.

BACKGROUND: Polymorphisms at mitochondrial (mt) loci could modulate the risk of diseases including cancers. Here the mtDNA polymorphisms at 12,308 nucleotide pairs (np), 11,467 np, 10,400 np, and 10,398 np were studied to examine the association with the risk of oral cancer and leukoplakia, alone and in combination with polymorphisms at the GST loci. METHODS: Polymorphisms at mt loci were screened in 310 cancer, 224 leukoplakia, and 389 control individuals by polymerase chain reaction (PCR) restriction length polymorphism (RFLP) and most of the GST genotype data were taken from previously published reports. Data were analyzed to determine the risk of the diseases. RESULTS: The major allele, A, at 12,308 np on tRNA(Leu) (CUN), increased the risk of cancer (odd ratio [OR] of 1.7; 95% confidence interval [95% CI], 1.1-2.6) but not that of leukoplakia. The same allele also appeared to increase the risk of cancer in smokers (OR of 4.0; 95% CI, 1.1-14.4), who are mostly males (OR of 1.8; 95% CI, 1.1-3-2), but not in smokeless tobacco users, who are mostly females. The major allele A at 11467 np demonstrated identical results as the major allele, A, at 12,308 np. The major alleles G at 10,398 np and T at 10,400 np (ie, M-haplogroup) increased the risk of cancer significantly in smokers (OR of 2.6; 95% CI, 1.2-5.7 and OR of 2.4; 95% CI, 1.1-5.1, respectively). The risk-risk genotype-allele combination at GSTP1 and mt12308 np loci increased the risk of cancer (OR of 2.6; 95% CI, 1.4-4.9) when compared with the nonrisk-nonrisk combination in leukoplakia patients. CONCLUSIONS: Polymorphisms at the mt loci alone and in combination with the risk genotype at GSTP1 increased the risk of oral cancer. Thus, risk genotypes from 2 different organelles may work in combination to increase the risk of oral cancer.