Exploring the IRE1 interactome: From canonical signaling functions to unexpected roles.

The unfolded protein response is a mechanism aiming at restoring endoplasmic reticulum (ER) homeostasis and is likely involved in other adaptive pathways. The unfolded protein response is transduced by three proteins acting as sensors and triggering downstream signaling pathways. Among them, inositol-requiring enzyme 1 alpha (IRE1α) (referred to as IRE1 hereafter), an endoplasmic reticulum-resident type I transmembrane protein, exerts its function through both kinase and endoribonuclease activities, resulting in both X-box binding protein 1 mRNA splicing and RNA degradation (regulated ire1 dependent decay). An increasing number of studies have reported protein-protein interactions as regulators of these signaling mechanisms, and additionally, driving other noncanonical functions. In this review, we deliver evolutive and structural insights on IRE1 and further describe how this protein interaction network (interactome) regulates IRE1 signaling abilities or mediates other cellular processes through catalytic-independent mechanisms. Moreover, we focus on newly discovered targets of IRE1 kinase activity and discuss potentially novel IRE1 functions based on the nature of the interactome, thereby identifying new fields to explore regarding this protein's biological roles.

Cirscan: a shiny application to identify differentially active sponge mechanisms and visualize circRNA-miRNA-mRNA networks.

BACKGROUND: Non-coding RNAs represent a large part of the human transcriptome and have been shown to play an important role in disease such as cancer. However, their biological functions are still incompletely understood. Among non-coding RNAs, circular RNAs (circRNAs) have recently been identified for their microRNA (miRNA) sponge function which allows them to modulate the expression of miRNA target genes by taking on the role of competitive endogenous RNAs (ce-circRNAs). Today, most computational tools are not adapted to the search for ce-circRNAs or have not been developed for the search for ce-circRNAs from user's transcriptomic data. RESULTS: In this study, we present Cirscan (CIRcular RNA Sponge CANdidates), an interactive Shiny application that automatically infers circRNA-miRNA-mRNA networks from human multi-level transcript expression data from two biological conditions (e.g. tumor versus normal conditions in the case of cancer study) in order to identify on a large scale, potential sponge mechanisms active in a specific condition. Cirscan ranks each circRNA-miRNA-mRNA subnetwork according to a sponge score that integrates multiple criteria based on interaction reliability and expression level. Finally, the top ranked sponge mechanisms can be visualized as networks and an enrichment analysis is performed to help its biological interpretation. We showed on two real case studies that Cirscan is capable of retrieving sponge mechanisms previously described, as well as identifying potential novel circRNA sponge candidates. CONCLUSIONS: Cirscan can be considered as a companion tool for biologists, facilitating their ability to prioritize sponge mechanisms for experimental validations and identifying potential therapeutic targets. Cirscan is implemented in R, released under the license GPL-3 and accessible on GitLab ( https://gitlab.com/geobioinfo/cirscan_Rshiny ). The scripts used in this paper are also provided on Gitlab ( https://gitlab.com/geobioinfo/cirscan_paper ).

TGFß-induced long non-coding RNA LINC00313 activates Wnt signaling and promotes cholangiocarcinoma.

Cholangiocarcinoma is a devastating liver cancer characterized by high aggressiveness and therapy resistance, resulting in poor prognosis. Long non-coding RNAs and signals imposed by oncogenic pathways, such as transforming growth factor ß (TGFß), frequently contribute to cholangiocarcinogenesis. Here, we explore novel effectors of TGFß signalling in cholangiocarcinoma. LINC00313 is identified as a novel TGFß target gene. Gene expression and genome-wide chromatin accessibility profiling reveal that nuclear LINC00313 transcriptionally regulates genes involved in Wnt signalling, such as the transcriptional activator TCF7. LINC00313 gain-of-function enhances TCF/LEF-dependent transcription, promotes colony formation in vitro and accelerates tumour growth in vivo. Genes affected by LINC00313 over-expression in CCA tumours are associated with KRAS and TP53 mutations and reduce overall patient survival. Mechanistically, ACTL6A and BRG1, subunits of the SWI/SNF chromatin remodelling complex, interact with LINC00313 and affect TCF7 and SULF2 transcription. We propose a model whereby TGFß induces LINC00313 in order to regulate the expression of hallmark Wnt pathway genes, in co-operation with SWI/SNF. By modulating key genes of the Wnt pathway, LINC00313 fine-tunes Wnt/TCF/LEF-dependent transcriptional responses and promotes cholangiocarcinogenesis.

IRE1 endoribonuclease signaling promotes myeloid cell infiltration in glioblastoma.

BACKGROUND: Intrinsic or environmental stresses trigger the accumulation of improperly folded proteins in the endoplasmic reticulum (ER), leading to ER stress. To cope with this, cells have evolved an adaptive mechanism named the unfolded protein response (UPR) which is hijacked by tumor cells to develop malignant features. Glioblastoma (GB), the most aggressive and lethal primary brain tumor, relies on UPR to sustain growth. We recently showed that IRE1 alpha (referred to IRE1 hereafter), 1 of the UPR transducers, promotes GB invasion, angiogenesis, and infiltration by macrophage. Hence, high tumor IRE1 activity in tumor cells predicts a worse outcome. Herein, we characterized the IRE1-dependent signaling that shapes the immune microenvironment toward monocytes/macrophages and neutrophils. METHODS: We used human and mouse cellular models in which IRE1 was genetically or pharmacologically invalidated and which were tested in vivo. Publicly available datasets from GB patients were also analyzed to confirm our findings. RESULTS: We showed that IRE1 signaling, through both the transcription factor XBP1s and the regulated IRE1-dependent decay controls the expression of the ubiquitin-conjugating E2 enzyme UBE2D3. In turn, UBE2D3 activates the NFκB pathway, resulting in chemokine production and myeloid infiltration in tumors. CONCLUSIONS: Our work identifies a novel IRE1/UBE2D3 proinflammatory axis that plays an instrumental role in GB immune regulation.

Non-canonical miRNA-RNA base-pairing impedes tumor suppressor activity of miR-16.

Uveal melanoma (UM), the most common primary intraocular tumor in adults, has been extensively characterized by omics technologies during the last 5 yr. Despite the discovery of gene signatures, the molecular actors driving cancer aggressiveness are not fully understood, and UM is still associated with very poor overall survival (OS) at the metastatic stage. By defining the miR-16 interactome, we revealed that miR-16 mainly interacts via non-canonical base-pairing to a subset of RNAs, promoting their expression levels. Consequently, the canonical miR-16 activity, involved in the RNA decay of oncogenes, such as <i>cyclin D3</i>, is impaired. This non-canonical base-pairing can explain both the derepression of miR-16 targets and the promotion of oncogene expression observed in patients with poor OS in two cohorts. miR-16 activity, assessment using our RNA signature, discriminates the patient's OS as effectively as current methods. To the best of our knowledge, this is the first time that a predictive signature has been composed of genes belonging to the same mechanism (miR-16) in UM. Altogether, our results strongly suggest that UM is a miR-16 disease.

A Novel Pseudogene Methylation Signature to Predict Temozolomide Outcome in Non-G-CIMP Glioblastomas.

Objective: Alterations in the methylation state of pseudogenes may serve as clinically useful biomarkers of glioblastomas (GBMs) that do not have glioma-CpG island methylator phenotype (G-CIMP). Methods: Non-G-CIMP GBM datasets were included for evaluation, and a RISK-score signature was determined from the methylation state of pseudogene loci. Both bioinformatic and experimental analyses were performed for biological validation. Results: By integrating clinical information with DNA methylation microarray data, we screened a panel of eight CpGs from discovery cohorts of non-G-CIMP GBMs. Each CpG could accurately and independently predict the prognosis of patients under a treatment regime that combined radiotherapy (RT) and temozolomide (TMZ). The 8-CpG signature appeared to show opposite prognostic correlations between patients treated with RT/TMZ and those treated with RT monotherapy. The analyses further indicated that this signature had predictive value for TMZ efficacy because different survival benefits between RT/TMZ and RT therapies were observed in each risk subgroup. The incorporation of other risk factors, such as age and O-6-methylguanine-DNA methyltransferase (MGMT) promoter methylation status, with our pseudogene methylation signature could provide precise risk classification. In vitro experimental data revealed that two locus-specific pseudogenes (ZNF767P and CLEC4GP1) may modulate TMZ resistance via distinct mechanisms in GBM cells. Conclusion: The biologically and clinically relevant RISK-score signature, based on pseudogene methylation loci, may offer information for predicting TMZ responses of non-G-CIMP GBMs, that is independent from, but complementary to, MGMT-based approaches.

Integration of Raman spectra with transcriptome data in glioblastoma multiforme defines tumour subtypes and predicts patient outcome.

Raman spectroscopy is an imaging technique that has been applied to assess molecular compositions of living cells to characterize cell types and states. However, owing to the diverse molecular species in cells and challenges of assigning peaks to specific molecules, it has not been clear how to interpret cellular Raman spectra. Here, we provide firm evidence that cellular Raman spectra (RS) and transcriptomic profiles of glioblastoma can be computationally connected and thus interpreted. We find that the dimensions of high-dimensional RS and transcriptomes can be reduced and connected linearly through a shared low-dimensional subspace. Accordingly, we were able to predict global gene expression profiles by applying the calculated transformation matrix to Raman spectra and vice versa. From these analyses, we extract a minimal gene expression signature associated with specific RS profiles and predictive of disease outcome.

Integrative analysis identifies an immune-relevant epigenetic signature for prognostication of non-G-CIMP glioblastomas.

The clinical and molecular implications of DNA methylation alterations remain unclear among the majority of glioblastomas (GBMs) without glioma-CpGs island methylator phenotype (G-CIMP); integrative multi-level molecular profiling may provide useful information. Independent cohorts of non-G-CIMP GBMs or IDH wild type (wt) lower-grade gliomas (LGGs) from local and public databases with DNA methylation and gene expression microarray data were included for discovery and validation of a multimarker signature, combined using a RISK score model. Bioinformatic and in vitro functional analyses were employed for biological validation. Using a strict multistep selection approach, we identified eight CpGs, each of which was significantly correlated with overall survival (OS) of non-G-CIMP GBMs, independent of age, the O-6-methylguanine-DNA methyltransferase (MGMT) methylation status, treatments and other identified CpGs. An epigenetic RISK signature of the 8 CpGs was developed and validated to robustly and independently prognosticate prognosis in different cohorts of not only non-G-GIMP GBMs, but also IDHwt LGGs. It also showed good discriminating value in stratified cohorts by current clinical and molecular factors. Bioinformatic analysis revealed consistent correlation of the epigenetic signature to distinct immune-relevant transcriptional profiles of GBM bulks. Functional experiments showed that S100A2 appeared to be epigenetically regulated by one identified CpG and was associated with GBM cell proliferation, apoptosis, invasion, migration and immunosuppression. The prognostic 8-CpGs RISK score signature may be of promising value for refining current glioma risk classification, and its potential links to distinct immune phenotypes make it a promising biomarker candidate for predicting response to anti-glioma immunotherapy.

CRISPR screens identify tumor-promoting genes conferring melanoma cell plasticity and resistance.

Most genetic alterations that drive melanoma development and resistance to targeted therapy have been uncovered. In contrast, and despite their increasingly recognized contribution, little is known about the non-genetic mechanisms that drive these processes. Here, we performed in vivo gain-of-function CRISPR screens and identified SMAD3, BIRC3, and SLC9A5 as key actors of BRAFi resistance. We show that their expression levels increase during acquisition of BRAFi resistance and remain high in persister cells and during relapse. The upregulation of the SMAD3 transcriptional activity (SMAD3-signature) promotes a mesenchymal-like phenotype and BRAFi resistance by acting as an upstream transcriptional regulator of potent BRAFi-resistance genes such as EGFR and AXL. This SMAD3-signature predicts resistance to both current melanoma therapies in different cohorts. Critically, chemical inhibition of SMAD3 may constitute amenable target for melanoma since it efficiently abrogates persister cells survival. Interestingly, decrease of SMAD3 activity can also be reached by inhibiting the Aryl hydrocarbon Receptor (AhR), another druggable transcription factor governing SMAD3 expression level. Our work highlights novel drug vulnerabilities that can be exploited to develop long-lasting antimelanoma therapies.

EZH2 and KDM6B Expressions Are Associated with Specific Epigenetic Signatures during EMT in Non Small Cell Lung Carcinomas.

The role of Epigenetics in Epithelial Mesenchymal Transition (EMT) has recently emerged. Two epigenetic enzymes with paradoxical roles have previously been associated to EMT, EZH2 (Enhancer of Zeste 2 Polycomb Repressive Complex 2 (PRC2) Subunit), a lysine methyltranserase able to add the H3K27me3 mark, and the histone demethylase KDM6B (Lysine Demethylase 6B), which can remove the H3K27me3 mark. Nevertheless, it still remains unclear how these enzymes, with apparent opposite activities, could both promote EMT. In this study, we evaluated the function of these two enzymes using an EMT-inducible model, the lung cancer A549 cell line. ChIP-seq coupled with transcriptomic analysis showed that EZH2 and KDM6B were able to target and modulate the expression of different genes during EMT. Based on this analysis, we described INHBB, WTN5B, and ADAMTS6 as new EMT markers regulated by epigenetic modifications and directly implicated in EMT induction.

Local intracerebral inhibition of IRE1 by MKC8866 sensitizes glioblastoma to irradiation/chemotherapy in vivo.

Glioblastoma multiforme (GBM) is the most severe primary brain cancer. Despite an aggressive treatment comprising surgical resection and radio/chemotherapy, patient's survival post diagnosis remains short. A limitation for success in finding novel improved therapeutic options for such dismal disease partly lies in the lack of a relevant animal model that accurately recapitulates patient disease and standard of care. In the present study, we have developed an immunocompetent GBM model that includes tumor surgery and a radio/chemotherapy regimen resembling the Stupp protocol and we have used this model to test the impact of the pharmacological inhibition of the endoplasmic reticulum (ER) stress sensor IRE1, on treatment efficacy.

Transcriptomics in cancer revealed by Positron Emission Tomography radiomics.

Metabolic images from Positron Emission Tomography (PET) are used routinely for diagnosis, follow-up or treatment planning purposes of cancer patients. In this study we aimed at determining if radiomic features extracted from 18F-Fluoro Deoxy Glucose (FDG) PET images could mirror tumor transcriptomics. In this study we analyzed 45 patients with locally advanced head and neck cancer (H&N) that underwent FDG-PET scans at the time of diagnosis and transcriptome analysis using RNAs from both cancer and healthy tissues on microarrays. Association between PET radiomics and transcriptomics was carried out with the Genomica software and a functional annotation was used to associate PET radiomics, gene expression and altered biological pathways. We identified relationships between PET radiomics and genes involved in cell-cycle, disease, DNA repair, extracellular matrix organization, immune system, metabolism or signal transduction pathways, according to the Reactome classification. Our results suggest that these FDG PET radiomic features could be used to infer tissue gene expression and cellular pathway activity in H&N cancers. These observations strengthen the value of radiomics as a promising approach to personalize treatments through targeting tumor-specific molecular processes.

Cytosine methylation of mature microRNAs inhibits their functions and is associated with poor prognosis in glioblastoma multiforme.

BACKGROUND: Literature reports that mature microRNA (miRNA) can be methylated at adenosine, guanosine and cytosine. However, the molecular mechanisms involved in cytosine methylation of miRNAs have not yet been fully elucidated. Here we investigated the biological role and underlying mechanism of cytosine methylation in miRNAs in glioblastoma multiforme (GBM). METHODS: RNA immunoprecipitation with the anti-5methylcytosine (5mC) antibody followed by Array, ELISA, dot blot, incorporation of a radio-labelled methyl group in miRNA, and miRNA bisulfite sequencing were perfomred to detect the cytosine methylation in mature miRNA. Cross-Linking immunoprecipiation qPCR, transfection with methylation/unmethylated mimic miRNA, luciferase promoter reporter plasmid, Biotin-tagged 3'UTR/mRNA or miRNA experiments and in vivo assays were used to investigate the role of methylated miRNAs. Finally, the prognostic value of methylated miRNAs was analyzed in a cohorte of GBM pateints. RESULTS: Our study reveals that a significant fraction of miRNAs contains 5mC. Cellular experiments show that DNMT3A/AGO4 methylated miRNAs at cytosine residues inhibit the formation of miRNA/mRNA duplex and leading to the loss of their repressive function towards gene expression. In vivo experiments show that cytosine-methylation of miRNA abolishes the tumor suppressor function of miRNA-181a-5p miRNA for example. Our study also reveals that cytosine-methylation of miRNA-181a-5p results is associated a poor prognosis in GBM patients. CONCLUSION: Together, our results indicate that the DNMT3A/AGO4-mediated cytosine methylation of miRNA negatively.

Novel predictive epigenetic signature for temozolomide in non-G-CIMP glioblastomas.

OBJECTIVE: To identify novel epigenetic signatures that could provide predictive information that is complementary to promoter methylation status of the O-6-methylguanine-DNA methyltransferase (MGMT) gene for predicting temozolomide (TMZ) response, among glioblastomas (GBMs) without glioma-CpGs island methylator phenotype (G-CIMP) METHODS: Different cohorts of primary non-G-CIMP GBMs with genome-wide DNA methylation microarray data were included for discovery and validation of a multimarker signature, combined using a RISK score model. Different statistical analyses and functional experiments were performed for clinical and biological validation. RESULTS: By employing discovery cohorts with radiotherapy (RT) and TMZ versus RT alone and a strict multistep selection strategy, we identified seven CpGs, each of which was significantly correlated with overall survival (OS) of non-G-CIMP GBMs with RT/TMZ, independent of age, MGMT promoter methylation status, and other identified CpGs. A RISK score signature of the 7 CpGs was developed and validated to distinguish non-G-CIMP GBMs with differential survival outcomes to RT/TMZ, but not to RT alone. The interaction analyses also showed differential outcomes to RT/TMZ versus RT alone within the RISK score-based subgroups. The signature could also improve the risk classification by age and MGMT promoter methylation status. Functional experiments showed that HSBP2 appeared to be epigenetically regulated by one identified CpG and was associated with TMZ resistance, but it was not associated with cell proliferation or apoptosis in GBM cell lines. The predictive value of the single CpG methylation of HSBP2 by pyrosequencing was observed in a local cohort of isocitrate dehydrogenase 1 (IDH1) R132H wild-type GBMs. CONCLUSIONS: This novel epigenetic signature might be a promising predictive (but not a general prognostic) biomarker and be helpful for refining the MGMT-based guiding approach to TMZ usage in non-G-CIMP GBMs.

A five-CpG signature of microRNA methylation in non-G-CIMP glioblastoma.

AIMS: DNA methylation has been found to regulate microRNAs (miRNAs) expression, but the prognostic value of miRNA-related DNA methylation aberration remained largely elusive in cancers including glioblastomas (GBMs). This study aimed to investigate the clinical and biological feature of miRNA methylation in GBMs of non-glioma-CpG island methylator phenotype (non-G-CIMP). METHODS: Prognostic miRNA methylation loci were analyzed, with TCGA and Rennes cohort as training sets, and independent datasets of GBMs and low-grade gliomas (LGGs) were obtained as validation sets. Different statistical and bioinformatic analysis and experimental validations were performed to clinically and biologically characterize the signature. RESULTS: We identified and validated a risk score based on methylation status of five miRNA-associated CpGs which could predict survival of GBM patients in a series of training and validation sets. This signature was independent of age and O-6-methylguanine-DNA methyltransferase (MGMT) promoter methylation status. The risk subgroup was associated with angiogenesis and accordingly differential responses to bevacizumab-contained therapy. MiRNA target analysis and in vitro experiments further confirmed the accuracy of this signature. CONCLUSION: The five-CpG signature of miRNA methylation was biologically relevant and was of potential prognostic and predictive value for GBMs. It might be of help for improving individualized treatment.

EMT is associated with an epigenetic signature of ECM remodeling genes.

Type III epithelial-mesenchymal transition (EMT) has been previously associated with increased cell migration, invasion, metastasis, and therefore cancer aggressiveness. This reversible process is associated with an important gene expression reprogramming mainly due to epigenetic plasticity. Nevertheless, most of the studies describing the central role of epigenetic modifications during EMT were performed in a single-cell model and using only one mode of EMT induction. In our study, we studied the overall modulations of gene expression and epigenetic modifications in four different EMT-induced cell models issued from different tissues and using different inducers of EMT. Pangenomic analysis (transcriptome and ChIP-sequencing) validated our hypothesis that gene expression reprogramming during EMT is largely regulated by epigenetic modifications of a wide range of genes. Indeed, our results confirmed that each EMT model is unique and can be associated with a specific transcriptome profile and epigenetic program. However, we could select some genes or pathways that are similarly regulated in the different models and that could therefore be used as a common signature of all EMT models and become new biomarkers of the EMT phenotype. As an example, we can cite the regulation of gene-coding proteins involved in the degradation of the extracellular matrix (ECM), which are highly induced in all EMT models. Based on our investigations and results, we identified ADAM19 as a new biomarker of in vitro and in vivo EMT and we validated this biological new marker in a cohort of non-small lung carcinomas.

Rapid genetic and phenotypic changes in Pseudomonas aeruginosa clinical strains during ventilator-associated pneumonia.

Treatment with antibiotics leads to the selection of isolates with increased resistance. We investigated if evolution towards resistance was associated with virulence changes, in the context of P. aeruginosa ventilator-associated pneumonia (VAP). Four patients were selected because they had multiple VAP episodes during short periods (12 days to 5 weeks), with emergence of resistance. We performed whole-genome sequencing of 12 P. aeruginosa from bronchoalveolar lavages or blood culture (3 isolates per patient). Production of quorum sensing-dependent virulence factors, serum resistance, cytotoxicity against A549 cells, biofilm production, and twitching motility were studied. Each patient was infected with a unique strain. For all patients, resistance development was explained by genetic events in ampD, mexR or oprD. Additional variations were detected in virulence- and/or fitness-associated genes (algB, gacA, groEL, lasR, mpl, pilE, pilM, rhlR) depending on the strain. We noticed a convergence towards quorum sensing deficiency, correlated with a decrease of pyocyanin and protease production, survival in serum, twitching motility and cytotoxicity. In one patient, changes in pilM and pilE were related to enhanced twitching. We show that the emergence of resistance in P. aeruginosa is associated with virulence modification, even in acute infections. The consequences of this short-term pathoadaptation need to be explored.

The expression of EMX2 lead to cell cycle arrest in glioblastoma cell line.

BACKGROUND: Glioblastoma (GB) is a highly invasive primary brain tumor that nearly always systematically recurs at the site of resection despite aggressive radio-chemotherapy. Previously, we reported a gene expression signature related to tumor infiltration. Within this signature, the EMX2 gene encodes a homeodomain transcription factor that we found was down regulated in glioblastoma. As EMX2 is reported to play a role in carcinogenesis, we investigated the impact of EMX2 overexpression in glioma-related cell lines. METHODS: For that purpose, we constructed tetracycline-inducible EMX2 expression lines. Transfected cell phenotypes (proliferation, cell death and cell cycle) were assessed in time-course experiments. RESULTS: Restoration of EMX2 expression in U87 glioblastoma cells significantly inhibited cell proliferation. This inhibition was reversible after EMX2 removal from cells. EMX2-induced proliferative inhibition was very likely due to cell cycle arrest in G1/S transition and was not accompanied by signs of cell death. CONCLUSION: Our results suggest that EMX2 may constitute a putative therapeutic target for GB treatment. Further studies are required to decipher the gene networks and transduction signals involved in EMX2's effect on cell proliferation.

Integration of transcriptome and proteome profiles in glioblastoma: looking for the missing link.

BACKGROUND: Glioblastoma (GB) is the most common and aggressive tumor of the brain. Genotype-based approaches and independent analyses of the transcriptome or the proteome have led to progress in understanding the underlying biology of GB. Joint transcriptome and proteome profiling may reveal new biological insights, and identify pathogenic mechanisms or therapeutic targets for GB therapy. We present a comparison of transcriptome and proteome data from five GB biopsies (TZ) vs their corresponding peritumoral brain zone (PBZ). Omic analyses were performed using RNA microarray chips and the isotope-coded protein label method (ICPL). RESULTS: As described in other cancers, we found a poor correlation between transcriptome and proteome data in GB. We observed only two commonly deregulated mRNAs/proteins (neurofilament light polypeptide and synapsin 1) and 12 altered biological processes; they are related to cell communication, synaptic transmission and nervous system processes. This poor correlation may be a consequence of the techniques used to produce the omic profiles, the intrinsic properties of mRNA and proteins and/or of cancer- or GB-specific phenomena. Of interest, the analysis of the transcription factor binding sites present upstream from the open reading frames of all altered proteins identified by ICPL method shows a common binding site for the topoisomerase I and p53-binding protein TOPORS. Its expression was observed in 7/11 TZ samples and not in PBZ. Some findings suggest that TOPORS may function as a tumor suppressor; its implication in gliomagenesis should be examined in future studies. CONCLUSIONS: In this study, we showed a low correlation between transcriptome and proteome data for GB samples as described in other cancer tissues. We observed that NEFL, SYN1 and 12 biological processes were deregulated in both the transcriptome and proteome data. It will be important to analyze more specifically these processes and these two proteins to allow the identification of new theranostic markers or potential therapeutic targets for GB.

A novel prognostic six-CpG signature in glioblastomas.

AIMS: We aimed to identify a clinically useful biomarker using DNA methylation-based information to optimize individual treatment of patients with glioblastoma (GBM). METHODS: A six-CpG panel was identified by incorporating genome-wide DNA methylation data and clinical information of three distinct discovery sets and was combined using a risk-score model. Different validation sets of GBMs and lower-grade gliomas and different statistical methods were implemented for prognostic evaluation. An integrative analysis of multidimensional TCGA data was performed to molecularly characterize different risk tumors. RESULTS: The six-CpG risk-score signature robustly predicted overall survival (OS) in all discovery and validation cohorts and in a treatment-independent manner. It also predicted progression-free survival (PFS) in available patients. The multimarker epigenetic signature was demonstrated as an independent prognosticator and had better performance than known molecular indicators such as glioma-CpG island methylator phenotype (G-CIMP) and proneural subtype. The defined risk subgroups were molecularly distinct; high-risk tumors were biologically more aggressive with concordant activation of proangiogenic signaling at multimolecular levels. Accordingly, we observed better OS benefits of bevacizumab-contained therapy to high-risk patients in independent sets, supporting its implication in guiding usage of antiangiogenic therapy. Finally, the six-CpG signature refined the risk classification based on G-CIMP and MGMT methylation status. CONCLUSIONS: The novel six-CpG signature is a robust and independent prognostic indicator for GBMs and is of promising value to improve personalized management.